Open access peer-reviewed chapter

Model-Based Evolutionary Operation Design for Batch and Fed- Batch Antibiotic Production Bioprocesses

Written By

Samuel Conceição de Oliveira

Submitted: October 26th, 2016 Reviewed: April 21st, 2017 Published: December 20th, 2017

DOI: 10.5772/intechopen.69395

Chapter metrics overview

1,233 Chapter Downloads

View Full Metrics


The control policy determination for batch and fed-batch antibiotic production bioprocesses is an important practical issue due to the high added value of these bioproducts. Since it is highly desirable to optimize the antibiotic production, several methods have been proposed aimed at this objective. Once having a mathematical model for the bioprocess, the optimization problem can be formulated within the framework of Pontryagin's maximum principle and of the optimal control theory to determinate the best control trajectory for certain key manipulated variables, such as temperature, pH, and substrate feed rate. In this chapter, applications of these model-based techniques to optimize and control antibiotics production bioprocesses are reviewed and new aspects are emphasized. The cases analyzed included the optimization of the substrate feed rate in a fed-batch reactor and of the temperature in a batch reactor during penicillin fermentations. The main contributions of this study were: (i) the proposition of a different procedure for calculating the second switching time of substrate feed rate, (ii) the application of simpler numerical methods to solve the two-point boundary-value problem associated with the temperature profile optimization, and (iii) the demonstration that the non-isothermal operation is more productive in antibiotic than the operation under constant temperature.


  • modeling
  • optimization
  • evolutionary operation
  • bioprocesses
  • antibiotic fermentation
  • batch bioreactors
  • fed-batch bioreactors

1. Introduction

Improvement in the productivity of many submerged fermentation processes is carried out by manipulating nutritional and physical parameters such as medium composition, agitation speed, aeration rate, pH, and temperature [1, 2]. Although the attainment of optimal conditions for a multivariable fermentation process is often tedious, it is possible to undertake a rational procedure by using statistical experimental designs [2].

Experimental designs can be divided in two distinct groups [3, 4, 5, 6]: (i) model-based experimental designs and (ii) statistical experimental designs. In model-based experimental designs, predictions of a mathematical model are used to determine how an experiment or process should be performed, whereas, with statistical experimental designs, these model predictions are not explicitly required.

The optimization and operation of fermentation processes play a key role in the biotechnology industry due to heavy competition among companies. Secondary metabolites, such as antibiotics and other pharmaceutical products, represent an important added value; therefore, improvements in the production of these bioproducts are of great interest to industries. To achieve high-performance operations, the optimization of manipulated variables that affect the fermentation process becomes a significant task.

In general, optimization problems can be classified in two categories: set-point and profile optimizations [7]. Set-point optimization problems involve finding the best set of values of manipulated variables that lead to the maximization of performance indexes [7]. Profile optimization consists of determining temporal or spatial functions (profiles), rather than a point in n-dimensional space, which lend an optimal value to the performance index [7].

In antibiotic fermentation, it is well known that the temperature and pH for the maximum rate of antibiotic production are different from those for the maximum rate of cell growth [8]. In this sense, the implementation of temperature and pH profiles plays an important role in significant improvements in antibiotic production bioprocesses [8].

Since the primary goal of a fermentation process is the cost-effective production of bioproducts, it is important to select the more appropriate operating mode that allows the production of the desired product at a high concentration with a high productivity and yield [9]. Fed-batch bioprocesses have been widely employed for the production of various bioproducts, including primary and secondary metabolites [9]. In the particular case of secondary metabolites, such as antibiotics, the interaction between growth metabolism and product biosynthesis is critically affected by growth-limiting nutrient concentrations. Since both the underfeeding and the overfeeding of nutrients are detrimental to cell growth and product formation, due to the occurrence of phenomena such as cell starvation and catabolite repression, establishing a suitable feeding strategy is crucial in fed-batch bioprocesses [9, 10].

A particular time sequence of control variables may be required in order to conduct the bioprocess over time in a trajectory that provides the greatest productivity. This can lead to complex optimal time profiles for the control variables, which are sometimes impossible to be determined purely experimentally. Thus, appropriate mathematical and numerical methods can be applied for the determination of these profiles in order to reduce the experimental effort and the required time for optimization.

The search for the optimal pH, temperature, and substrate feed-rate profiles in batch and fed-batch antibiotic fermentation is a typical problem of optimization and evolutionary operations for which the use of kinetic models and powerful mathematical techniques is essential for their solution [7, 10, 11]. According to Rani and Rao [12], several approaches for the determination of optimal time profiles for control variables have been reported in the literature [13, 14, 15]. In these reports, the optimization problem is generally formulated on the basis of Pontryagin’s maximum principle, taking as a starting point a phenomenological mathematical model of the bioprocess. For simple mathematical models, the problem can be solved analytically, from the Hamiltonian of the system, by applying an iterative scheme on the control variable to determine the optimal control profile [8, 16, 17, 18, 19].

In this chapter, two studies on the optimization and the evolutionary operation of antibiotic production bioprocesses are revisited, and new results are obtained and highlighted. Such studies report mathematical models of bioprocesses, in conjunction with Pontryagin’s maximum principle, to optimize the substrate feed-rate profile for a fed-batch bioreactor and the temperature profile for a batch fermentation in order to maximize the production of antibiotic. The fundamentals of Pontryagin’s maximum principle, when applied to the cases analyzed, are also presented.

The aim is to provide a theoretical basis for the application of a model-based methodology that can be used for the optimization and control of bioprocesses from other antibiotics and secondary metabolites with a broad structural diversity and therapeutic activity, including antibacterial, antifungal, antiviral, antitumor, immunosuppressive, antihypertensive, and antihypercholesterolemic compounds.


2. Case studies: batch and fed-batch antibiotic production bioprocesses

During batch and fed-batch bioprocesses, the state variables (cell, substrate, oxygen and product concentrations, temperature, and pH) change significantly, from initial to final values. This dynamic behavior motivates the development of optimization methods to find the optimal time trajectories for the control variables in order to improve the performance of these bioprocesses.

Two case studies on the optimization of control variables in batch and fed-batch antibiotic production bioprocesses are revisited, and additional results are obtained and presented. The cases studied are those reported by Costa [20] and Constantinides and Mostouffi [17], concerning the optimization of the substrate feed rate in a fed-batch reactor and the temperature in a batch reactor, respectively. These cases are presented and detailed in the following sections.

2.1. Case study #1: determination of the optimal substrate feed-rate profile in fed-batch bioreactor for penicillin production

In this case study, the bioprocess of penicillin production by Penicillium chrysogenumis described by the mathematical model presented by Costa [20], which is based on the classical model proposed by Bajpai-Reuss for penicillin fermentation. In this model, the specific growth rate (μ) takes into account diffusional limitations that occur in the filamentous fungal biomass, as described by the Contois model. The specific rate of product formation (π) considers that penicillin production is repressed by high substrate concentrations (catabolic repression), being modeled by the Andrews equation. Penicillin degradation by hydrolysis is also considered, assuming first-order kinetics for this reaction. The specific rate of substrate consumption (σ) is represented by the Herbert-Pirt generalized model, whereby the substrate is consumed for cell growth and maintenance and for product formation. The equations of the full mathematical model are as follows (note that if the dilution rate is used as the control variable, the total mass-balance equation is not required for the optimization problem formulation, thus reducing the equation system dimension):

dXdt=μXXu; μ=μmSBX+SE1-2
dSdt=σX+(SfS)u; σ=μYX/S+πYP/S+mE3-4
dPdt=(πXkhP)Pu; π=πmSkm+S+S2kiE5-6


  • X, S, and Pdenote the concentrations of cell, substrate, and product, respectively;

  • μ, σ, and πare the specific rates of cell growth, substrate consumption, and product formation;

  • μm, B, YX/S, YP/S, m, πm, km, and kiare the parameters of the mathematical model, including kinetic and yield parameters;

  • u=D=F/V:u=control variable;D=dilution rate;F=feed rate;V=culture volumeE7

In matrix notation:



X¯=[XSP]; f¯(X¯)=[μXσXπXkhP]; g¯(X¯)=[X(SfS)P]E9-11

The constraints imposed on the control variable uare as follows:


where umin = 0 and umax = Fmax / V.

Another constraint concerns the maximum volume of culture (final volume), i.e., V(tf) = Vmax = Vf, where tfis the final processing time.

The initial conditions are given by

X(0)=X0; S(0)=S0; P(0)=P0; V(0)=V0E13

The objective of the optimization/control problem is to determine the optimal time profile for the control variable that maximizes the antibiotic concentration at the end of the bioprocess.

According to the Pontryagin’s maximum principle, the optimal profile must maximize the Hamiltonian, given by

H=λ¯T [f¯(X¯)+g¯(X¯)u];  λ¯T=[λ1  λ2  λ3]E14


H=H0(t)+φ(t)u  {H0(t)=λ¯Tf¯(X¯)φ(t)=λ¯Tg¯(X¯)E15

Since f¯(X¯)=[f1f2f3]and g¯(X¯)=[g1g2g3], the following equations are obtained:


For optimal control, it is established that

  • If φ(t)>0, then u=umax.

  • If φ(t)<0, then u=umin.

  • If φ(t)=0,then u=usin.

In the singular interval: H(t)=0; ϕ(t)=0; H0(t)=0E18-20

Since λ¯=λ¯(t), the following equations can be developed:



Substituting Eqs. (21) and (24) into Eq. (22) gives


By developing the matrices indicated in the previous equation, one obtains

dφdt=[λ1  λ2  λ3]*[g1X1   g1X2   g1X3g2X1   g2X2   g2X3g3X1   g3X2   g3X3] *[f1f2f3][f1X1   f1X2   f1X3f2X1   f2X2   f2X3f3X1   f3X2   f3X3] *[g1g2g3]E27


g1=X; g2=(SfS); g3=P; X1=X; X2=S; X3=PE28-33
g1X1=1; g1X2=0; g1X3=0E34-36
g2X1=0; g2X2=1; g2X3=0E37-39
g3X1=0; g3X2=0; g3X3=1E40-42
f1=μX; f2=σX; f3=πXkhP; μ=μ(S,X); σ=σ(S,X); π=π(S)E43-48
f1X1=μ+XμX; f1X2=XμS; f1X3=0E49-51
f2X1=(σ+Xσ'X); f2X2=XσS; f2X3=0E52-54
f3X1=π; f3X2=XπS; f3X3=khE55-57


dϕdt=[λ1  λ2  λ3]*[100010001] *[f1f2f3][μ+XμXXμS0(σ+XσX)XσS0πXπSkh]*[g1g2g3]E58
dϕdt=[λ1  λ2  λ3]*[f1f2f3][g1(μ+XμX)+g2XμSg1(σ+XσX)g2XσS g1π+g2XπSg3kh] E59

In the singular interval:


From Eq. (64):


Substituting Eq. (66) into Eq. (65) results in the following equation:


By introducing the expression of λ2 into the expression of λ1, one obtains


Substituting λ1=αλ3and λ2=βλ3into the expression of dϕdtin the singular interval provides


The equation Q(X¯)=0is the expression of the singular arc, which is independent of the adjoint variables λ1, λ2, λ3. In the expression of Q(X¯), the indicated derivatives are given by

μ=μmSBX+SμX=μmBS(BX+S)2; μS=(BX+S)μmμmS(BX+S)2E72-73
π=πmSkm+S+S2ki  πS=(km+S+S2ki)πmπmS(1+2kiS)(km+S+S2ki)2E74
 σ=μYX/S+πYP/S+mσX=μXYX/S; σS=μSYX/S+πSYP/SE75-76

For the determination of the singular dilution rate (usin), one starts from the first derivative of ϕas follows:






Substituting the expression of 3/dtinto the expression of the second derivative of φgives


When u = usin, the second derivative of ϕis zero, i.e., d2ϕdt2=0. Thus,


The expression of usin determines how the dilution rate must be manipulated (varied) during the singular interval, being a function only of the state variables. When developed, the final expression obtained for usin is quite complex and extensive, which was implemented in computational programming language to evaluate the value of this variable over the singular interval.

The condition for stopping the integration of mass-balance equations during the period following the singular interval, which is conducted in batch mode (u= 0), is determined from the fact that for free final time (tf) problems, how is the case, H(tf) = 0. Thus


As the final conditions of the adjoint variables are λ1 = 0, λ2 = 0, and λ3 = 1, the stopping condition at t = tfis dP/dt= 0.

Thus, the problem-solving algorithm consisted of the following steps:

  1. Integrate the mass-balance equations with u = 0(batch operation) until the values of the state variables satisfy Q(X¯)=0, being the instant that this occurs, the first switching time t1;

  2. From instant t1, make u = usin until the reactor volume reaches its maximum value, which corresponds to the second switching time t2;

  3. Starting from time t2, return the operation of the reactor with u = 0until the stop condition is reached in tf.

From the data reported by Costa [20] and summarized in Table 1, the mass-balance equations were numerically integrated to determinate t1, i.e., the instant at which Q(X¯)=0. According to the data presented in Table 2, this time instant can be determined as t1 = 28.74 h, since Q(X¯)changes its signal around this time value.

Using the computational program developed for the calculation of usin, it was verified that this control variable was practically constant during the singular interval, being 0.00213 h-1 a value quite representative.

Kinetic parameters/operating variablesValues
μm(h-1); B(g-S/g-X)1.1×10-1; 6×10-3
πm(g g-1 h-1); km(g/L); ki(g/L); kh(h-1)4.0×10-3; 1.0×10-4; 1.0×10-1; 1.0×10-2
YX/S(g/g); YP/S(g/g); m(g g-1 h-1)0.47; 1.2; 2.9×10-2
X0(g/L); S0(g/L); P0(g/L); Sf(g/L)1.3; 69.0; 0.0; 500.0
V0(L); Vf(L)8.121; 10.0

Table 1.

Values of kinetic parameters and operating variables used in the case study of the optimization of penicillin production in fed-batch reactor (source: [20]).

With respect to the second switching time (t2), this was determined from usin and the initial and final (maximum) volumes as follows:


Thus, the singular interval duration is of 97.71 h and the second switching time (t2) is t2 = t1 + ∆t= 28.74 + 97.71 = 126.45 h. The mass-balance equations were then integrated until this time, and it was verified that, at the end of the integration, the stop condition (dP/dt=0) had been satisfied, dispensing with the complementary period of batch operation and reaching a final product concentration (P) of 6.35 g/L.

2.1.1. Simulation of the penicillin production bioprocess in fed-batch reactor under optimized conditions

The simulation of the bioprocess under optimized conditions was performed using a computer program in FORTRAN language. For the numerical integration of the ordinary differential equations corresponding to the mass balances, the variable-step fourth-order Runge-Kutta-Gill method was used [17]. The full profiles of the state variables during the bioprocess are shown in Figures 1 and 2.

Figure 1.

Temporal profiles of theS(dotted line),X(dashed line), andP(solid line) state variables during a fed-batch penicillin production bioprocess.

Figure 2.

Temporal profile of the state variableV(dotted line) and of the control variableu(solid line) during a fed-batch penicillin production bioprocess.

Due to the decoupling between biomass growth and product synthesis, this type of fermentation behaves as a biphasic process. Therefore, characteristic profiles of penicillin fermentation were obtained for the state variables as shown in Figure 1, i.e., a first phase of accumulation of the cell is observed in which the substrate is almost entirely consumed for this purpose, without associated product formation (trophophase). After this growth phase, the fed substrate is practically used for penicillin production since there is no further catabolic repression of the antibiotic synthesis due to the low levels of substrate concentration established in the reactor in this second phase (idiophase). In addition, the kinetic pattern observed is in agreement with that expected for a secondary metabolite, i.e., the production occurs mostly after cell growth.

Regarding the fermentation medium volume in the bioreactor, the behavior of this variable shown in Figure 2 was already expected since, during the batch operation, this volume is constant because there is no addition or removal of fermentation medium to or from the bioreactor. In the fed-batch operation with continuous feed of unfermented medium to the bioreactor at a constant flow rate, the volume increases linearly over time, as shown in Figure 2. The temporal profile exhibited by the control variable u(dilution rate) in a step format is due to the change in the bioreactor operation, from batch (u= 0) to fed-batch (u≠ 0) mode.

2.2. Case study #2: determination of the optimal temperature profile in a batch bioreactor for penicillin production

In this case study, the fungal growth is described by the logistic law, a substrate-independent model for microorganism population dynamics. In addition, the production of penicillin is also modeled considering that the formation of antibiotic is not associated with cell growth, and that the product is degraded by hydrolysis according to a first-order kinetics. The mathematical model, comprising two ordinary differential equations corresponding to the mass balances of cell and product in a batch bioreactor, and containing four parameters, is represented by (more information about this model can be found in Ref. [21]):


where tis the time, Xis the cell concentration, Pis the antibiotic concentration, rXis the cell growth rate, rPis the antibiotic production rate, rhis the product’s hydrolysis rate, and μm, Xm, β, and khare the parameters of the model, according to the following meanings: μmis the maximum specific growth rate, Xmis the maximum possible cell concentration to be achieved, βis the constant of product formation not associated with growth, and khis the rate constant for the antibiotic hydrolysis reaction.

For the application of the Pontryagin’s maximum principle, the model variables were dimensionless and expressions describing the kinetic parameters (bi) as a function of temperature (θ) were incorporated in order to extend the validity range of the model to non-isothermal conditions. These functions have shapes typical of those found in microbial or enzyme-catalyzed reactions (concave down parabolas). The hydrolysis of the antibiotic is neglected in this dimensionless version of the model, with this version given by the following equations [17]:

dy1dτ=b1y1b1b2y12, y1(0)=0.03E87
dy2dτ=b3y1, y2(0)=0.0E88


  • y1= dimensionless concentration of cell (-); y2= dimensionless concentration of product (-); τ= dimensionless time, 0τ1(-)

    b1=w1[1.0w2(θw3)21.0w2(25w3)2]; b2=w4[1.0w2(θw3)21.0w2(25w3)2]; b3=w5[1.0w2(θw6)21.0w2(25w6)2]E89-91

  • w1=13.1;w2=0.005;w3=30°C

  • w4=0.94;w5=1.71;w6=20°C

As in this case, g¯(X¯)=0and u= 0, due to the reactor being operated in a batch mode, the mass-balances equations are simplified to



X¯=[y1y2]; f¯(X¯)=[f1f2]=[b1y1b1b2y12b3y1]E93-94

As previously established in the first case study, the Hamiltonian is given by

H=λ¯T(f¯(X¯)+g¯(X¯) u0)H=λ¯Tf¯(X¯);λ¯T=[λ1λ2]H=[λ1λ2][f1f2]=λ1f1+λ2f2=λ1(b1y1b1b2y12)+λ2(b3y1)E95

The temporal variation rates of the adjoint variables λ1 and λ2 are formulated as

dλ¯dτ=HX¯ddτ[λ1λ2]=[λ1b1+2b1b2y1λ2b3                0]E96

From the previous equation, the following equations can be derived


The necessary condition for the optimization of the bioprocess is


From the expressions bi= bi(θ), the following derivatives can be obtained

b1θ=[2w1w2(θw3)1.0w2(25w3)2]; (b1/b2)θ=0; b3θ=[2w5w2(θw6)1.0w2(25w6)2]E100-102

By inserting the derivatives of the parameters with respect to the temperature into the expression of ∂H/∂θ= 0, the expression of the optimal temperature profile (θopt) is obtained as follows:

θopt=[2λ1y1w1w2w31.0w2(25w3)2+ 2y1w5w2w61.0w2(25w6)2]/[2λ1y1w1w21.0w2(25w3)2+ 2y1w5w21.0w2(25w6)2] E103

As previously demonstrated, when the objective is to maximize the antibiotic concentration at the end of the bioprocess, it is necessary that λ1(1) = 0 and λ2(1) = 1. Since 2/dt= 0, the second condition requires that λ2 be constant and equal to 1.0 in the entire time domain, i.e., λ2 = 1.0 for 0 ≤ τ≤ 1.

Several numerical methods have been developed to solve this two-point boundary-value problem arising from the application of the maximum principle of Pontryagin to a batch penicillin production bioprocess. Constantinides and Mostoufi [17] used the orthogonal collocation method to solve this problem, justifying that this method is more accurate than the finite difference method. The problem was solved here using a much simpler numerical method than that of the orthogonal collocation to integrate the differential equations, which is the variable-step fourth-order Runge-Kutta-Gill method [17]. Thus, the algorithm for solving the problem consisted of the following steps:

  1. Assignment of an initial value for λ1(0);

  2. Integration of the system of ODEs from τ= 0 to τ= 1 and verification if λ1(1) = 0. If not, assign a new value to λ1(0) until the final condition is satisfied.

In order to make the computational algorithm autonomous for the determination of λ1(0), the Newton-Raphson method [17] was coupled to the numerical integration method, solving the following non-linear algebraic equation:


2.2.1. Simulation of the penicillin production bioprocess in batch reactor under optimized non-isothermal conditions

The proposed algorithm was implemented in FORTRAN programming language, and the profiles of the state variables (y1, y2, and θ) are presented in Figures 35. These profiles are in strict agreement with those reported by Constantinides and Mostoufi [17] when they used the orthogonal collocation method to solve this problem. The value determined for λ1 at τ= 0 was λ1(0) = 3.61210035.

Figure 3.

Cell dimensionless concentration profile during a non-isothermal penicillin fermentation.

Figure 4.

Product dimensionless concentration profile during a non-isothermal penicillin fermentation.

Figure 5.

Exact and approximate optimal temperature profiles for a non-isothermal penicillin fermentation.

The cell concentration profile shown in Figure 3 depicts the main phases involved in a typical microbial growth curve, i.e., the exponential, stationary, and decline phases. The decline phase is attributed to the negative effects of low temperatures on cell growth. In Figure 4, concerning the penicillin production dynamics, an initial short lag phase can be observed, followed by a transition phase in which penicillin production is initiated, until a final linear production phase is achieved.

According to the presented formulation (bioprocess model and Pontryagin’s maximum principle), the optimum temperature profile varies between 20 and 30°C following the curve (a) shown in Figure 5. This profile suggests a variable operating temperature during the growth and penicillin production phases, contradicting the standard industrial operating procedure of maintaining a constant temperature throughout the bioprocess. Particularly during the penicillin production phase, the profile prescribes a decrease in the operating temperature so that a high concentration of antibiotic is reached at the end of the bioprocess. The temperature profile (a) shown in Figure 5 may bring some practical difficulty to its programming/execution. In this context, an approximate profile derived from the exact profile, such as that represented by the curve (b) in Figure 5, may make the temperature programming strategy more feasible. This proposal is in agreement with that reported by Bailey and Ollis [22], i.e., the temperature schedule predicted by these calculations can be closely approximated in industrial practice with little added cost. The approximate profile was built from the following equations, which were based on the analysis of the temperature data generated by the exact profile:

  • 0.0τ0.4θ(°C)=θ(τ=0.0)+(θ(τ=0.4)θ(τ=0.0)0.40.0)(τ0.0)θ=29.6511.625τ

  • 0.4<τ<0.9θ(°C)=25.0

  • 0.9τ1.0θ(°C)=θ(τ=0.9)+(θ(τ=1.0)θ(τ=0.9)1.00.9)(τ0.9)θ=25.050(τ0.9)

2.2.2. Simulation of the penicillin production bioprocess in batch reactor under non-optimized isothermal conditions

A pertinent simulation to be performed is one under isothermal conditions to verify whether this thermal operation mode is, in fact, less productive in penicillin than non-isothermal mode following an optimal temperature profile. For this purpose, simulations were performed for constant temperatures of 20, 25, and 30°C and the results were compared with those obtained with the optimized temperature profile (Figure 6). Figure 6(a) illustrates the well-known fact that high temperatures (30°C) favor the growth of the fungus, while low temperatures (20°C) favor the synthesis of the antibiotic since relative to the amount of penicillin produced, more and less biomass was accumulated at these respective temperature levels [22]. It is observed in Figure 6(b) that the isothermal operation at an intermediate temperature to those investigated (θ= 25°C) presents a performance very similar to the operation with optimized temperature profile, becoming a viable alternative if the variable temperature strategy cannot be implemented, although with lower productivity, according to the data presented in Table 3. Despite the fact that operation at a fixed temperature between 24 and 25°C predominates in current industrial practice, the benefits of temperature programming during batch antibiotic fermentations are clear.

28.74330.093.06 ×10-31.40×10-2−86336.84

Table 2.

Data of the numerical integration of the mass-balance equations used to determine the first switching time (t1) during the penicillin production in fed-batch reactor.

θ(°C)y1 (τ=1)y2 (τ=1)(y2/y1)τ=1

Table 3.

Final values of the dimensionless concentration of cells (y1) and product (y2) in isothermal and non-isothermal penicillin fermentations.

Figure 6.

Dimensionless concentration profiles of cell and product during isothermal penicillin fermentations at different temperatures.


3. Conclusions

In this chapter, the usefulness of the Pontryagin’s maximum principle has been demonstrated for the optimization and operation of complex antibiotic production bioprocesses such as those conducted in batch and fed-batch reactors under isothermal/non-isothermal conditions. By applying this principle, it was possible to determine the optimal profile of temperature in batch reactors and substrate feed rate in fed-batch reactors that maximize the antibiotic concentration at the end of the bioprocess. Although having a rather complex mathematical formulation, the Pontryagin’s maximum principle can be classified as a powerful and suitable tool for the optimization, control, and model-driven operation of bioprocesses aiming at maximum productivity of bioproducts. However, for the application of this principle, it is necessary to dispose a mathematical model, preferably phenomenological and representative of the bioprocess, in order to evaluate whether or not the solution found for a given problem is feasible. In the present study, two classical phenomenological models of penicillin production bioprocesses were used, together with the Pontryagin’s maximum principle, aiming to determine the optimal operating conditions for the production of antibiotic, and the solutions found are considered feasible and can be implemented in real cases. However, a more complete mathematical model, incorporating the medium oxygenation state, could provide better bioprocess control, since the productivity in penicillin fermentations is highly dependent upon dissolved oxygen concentration, with its critical level being around 30% of saturation. In the models used here, the dissolved oxygen concentration was implicitly assumed to be non-limiting of the bioprocess, making this a rather restrictive hypothesis.



The author wishes to thank CNPq for their financial support (protocol number: 455487/2014-6).


  1. 1. Sircar A, Sridhar P, Das PK. Optimization of solid state medium for the production of clavulanic acid byStreptomyces clavuligerus. Process Biochemistry. 1998;33(3):283-289
  2. 2. Chen WC. Medium improvement for β-fructofuranosidase production byAspergillus japonicus. Process Biochemistry. 1998;33(3):267-271
  3. 3. Berkholz R, Guthke R. Model based sequential experimental design for bioprocess optimisation—An overview. In: Hofman M, Thonart P, editors. Engineering and Manufacturing for Biotechnology. Heidelberg, the Netherlands: Springer; 2002. pp. 129-141. DOI: 10.1007/0-306-46889-1
  4. 4. Mandenius CF, Brundin A. Bioprocess optimization using design-of-experiments methodology. Biotechnology Progress. 2008;24:1191-1203
  5. 5. Núñez EGF, Véliz RV, Costa BLV, Rezende AG, Tonso A. Using statistical tools for improving bioprocesses. Asian Journal of Biotechnology. 2013;5:1-20
  6. 6. Kennedy M, Krouse D. Strategies for improving fermentation medium performance: A review. Journal of Industrial Microbiology & Biotechnology. 1999;23:456-475
  7. 7. Lim HC, Lee KS. Process control and optimization. In: Pons MN, editor. Bioprocess Monitoring and Control. Munique: Hanser Publishers; 1992. pp. 159-222
  8. 8. Chu WBZ, Constantinides A. Modeling, optimization, and computer control of the cephalosporin C fermentation process. Biotechnology and Bioengineering. 1988;32(3):277-288
  9. 9. Lee J, Lee SY, Park S, Middelberg APJ. Control of fed-batch fermentations. Biotechnology Advances. 1999;17:29-48
  10. 10. Chaudhuri B, Modak JM. Optimization of fed-batch bioreactor using neural network model. Bioprocess Engineering. 1998;1:71-79
  11. 11. Zhang H, Zhang Z, Lan LH. Evolutionary optimization of a fed-batch penicillin fermentation process. In: 2010 International Symposium on Computer, Communication, Control and Automation; 5-7 May 2010; Taina. IEEE; 2010. pp. 403-406
  12. 12. Rani KY, Rao VSR. Control of fermenters—A review. Bioprocess Engineering. 1999;21:77-88
  13. 13. Ashoori A, Moshiri B, Khaki-Sedigh A, Bakhtiari MR. Optimal control of a nonlinear fed-batch fermentation process using model predictive approach. Journal of Process Control. 2009;19:1162-1173
  14. 14. Rocha M, Mendes R, Rocha O, Rocha I, Ferreira EC. Optimization of fed-batch fermentation processes with bio-inspired algorithms. Expert Systems with Applications. 2014;41:2186-2195
  15. 15. Roeva O, Tzonkov S. A genetic algorithm for feeding trajectory optimisation of fed-batch fermentation processes. Bioautomation. 2009;12:1-12
  16. 16. Guthke R, Knorre WA. Optimal substrate profile for antibiotic fermentations. Biotechnology and Bioengineering. 1981;23:2771-2777
  17. 17. Constantinides A, Mostoufi N. Numerical Methods for Chemical Engineers with MATLAB Applications. Upper Saddle River: Prentice Hall PTR; 1999. p. 560
  18. 18. Van Impe JF, Nicolai BM, Vanrolleghem PA, Spriet JA, De Moor B, Vandewalle J. Optimal control of the penicillin G fed-batch fermentation: An analysis of the model of Heijnenet al. Optimal Control Applications & Methods. 1994;15:13-34
  19. 19. Skolpap W, Scharer JM, Douglas PL, Moo-Young M. Optimal feed rate profiles for fed-batch culture in penicillin production. Songklanakarin Journal of Science and Technology. 2005;27(5):1057-1064
  20. 20. Costa AC. Singular control in bioreactors [thesis]. Rio de Janeiro: Federal University of Rio de Janeiro (UFRJ); 1996
  21. 21. Oliveira SC. Mathematical modeling of batch antibiotic production: Comparative analysis between phenomenological and empirical approaches. Journal of Bioprocessing & Biotechniques. 2017;7(1):298. DOI: 10.4172/2155-9821.1000298
  22. 22. Bailley JE, Ollis DF. Biochemical Engineering Fundamentals. 2nd ed. New York, NY: McGraw-Hill; 1986. p. 984

Written By

Samuel Conceição de Oliveira

Submitted: October 26th, 2016 Reviewed: April 21st, 2017 Published: December 20th, 2017