Next Article in Journal
Impact of Packaging Film and Beef Trimmings on Ground Beef Shelf Life
Next Article in Special Issue
Dynamic Thermal Properties Estimation Using Sensitivity Coefficients for Rapid Heating Process
Previous Article in Journal
Food Waste of Italian Families: Proportion in Quantity and Monetary Value of Food Purchases
Previous Article in Special Issue
Dynamic Modeling of the Impact of Temperature Changes on CO2 Production during Milk Fermentation in Batch Bioreactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Modeling of Carnobacterium maltaromaticum CNCM I-3298 Growth and Metabolite Production and Model-Based Process Optimization

by
Cristian Puentes
1,2,
Amélie Girardeau
1,
Stephanie Passot
1,
Fernanda Fonseca
1 and
Ioan-Cristian Trelea
1,*
1
INRAE, AgroParisTech, UMR SayFood, Université Paris-Saclay, F-78850 Thiverval-Grignon, France
2
CentraleSupélec, LGPM, Université Paris-Saclay, F-91192 Gif-sur-Yvette, France
*
Author to whom correspondence should be addressed.
Foods 2021, 10(8), 1922; https://doi.org/10.3390/foods10081922
Submission received: 16 July 2021 / Revised: 10 August 2021 / Accepted: 17 August 2021 / Published: 19 August 2021
(This article belongs to the Special Issue Dynamic Modelling and Simulation of Food Systems)

Abstract

:
Carnobacterium maltaromaticum is a species of lactic acid bacteria found in dairy, meat, and fish, with technological properties useful in food biopreservation and flavor development. In more recent years, it has also proven to be a key element of biological time–temperature integrators for tracking temperature variations experienced by perishable foods along the cold-chain. A dynamic model for the growth of C. maltaromaticum CNCM I-3298 and production of four metabolites (formic acid, acetic acid, lactic acid, and ethanol) from trehalose in batch culture was developed using the reaction scheme formalism. The dependence of the specific growth and production rates as well as the product inhibition parameters on the operating conditions were described by the response surface method. The parameters of the model were calibrated from eight experiments, covering a broad spectrum of culture conditions (temperatures between 20 and 37 °C; pH between 6.0 and 9.5). The model was validated against another set of eight independent experiments performed under different conditions selected in the same range. The model correctly predicted the growth kinetics of C. maltaromaticum CNCM I-3298 as well as the dynamics of the carbon source conversion, with a mean relative error of 10% for biomass and 14% for trehalose and the metabolites. The paper illustrates that the proposed model is a valuable tool for optimizing the culture of C. maltaromaticum CNCM I-3298 by determining operating conditions that favor the production of biomass or selected metabolites. Model-based optimization may thus reduce the number of experiments and substantially speed up the process development, with potential applications in food technology for producing starters and improving the yield and productivity of the fermentation of sugars into metabolites of industrial interest.

Graphical Abstract

1. Introduction

Carnobacterium maltaromaticum is a psychotropic species of lactic acid bacteria widely found in food such as dairy products, fish, and meat. It is a Gram-positive, facultative anaerobic bacterium, able to grow at alkaline pH (up to 9.6) [1,2].
In the food industry, C. maltaromaticum has potential applications related to health protection and organoleptic properties. These include the biopreservation of food, by inhibiting the growth of foodborne pathogens such as Listeria sp. in cold conditions, and the development of flavor in ripened cheese varieties [2,3,4].
This lactic acid bacterium may also be used as a biological indicator in time–temperature integrators (TTI): ‘smart-labels’ that monitor the time–temperature history of chilled products throughout the cold-chain [5,6]. Concentrates of the strain CNCM I-3298 have been selected as inoculum for TopCryo® labels, the only biological TTI that has been taken to market to date. A pH decline of the label medium, associated with bacterial growth and acid production, produces an irreversible color change from green to red as an indication to the consumer about the spoilage of the food to which the TTI is attached [7].
In these applications, C. maltaromaticum concentrates produced by fermentation may be used alone or in association with other microorganisms. Some experimental studies on C. maltaromaticum fermentation under different culture conditions have been reported in the literature [3,4,6,7,8,9,10,11,12,13,14]. The effect of temperature and pH on the acidifying activity was evaluated and modelled by Girardeau et al. [7]. However, there is a lack of knowledge on the characterization and optimization of Carnobacteria growth and production of various metabolites such as acids or flavor compounds in a bioreactor.
Carnobacteria are considered to be homofermentative lactic acid bacteria that produce lactic acid from glucose, with pyruvate as a central metabolic intermediate (via the Embden–Meyerhof pathway) [15,16,17]. However, pyruvate may be alternatively converted to acetate, ethanol, formate, and CO2 [16,18] under anaerobic conditions and substrate limitation [19], arising for example at the end of fermentation [20]. The production of organic acids by Carnobacteria is also strain-dependent [8,16,21]. A recent study reported that lactic, formic, and acetic acids are key organic acids produced by C. maltaromaticum in a meat juice medium [22], indicating that this microorganism has the enzymatic machinery to perform mixed-acid fermentation (Figure 1).
For optimization purposes, modeling has proven to be a powerful tool, enabling the exploration of a wider range of operating conditions while minimizing cost, compared with the experimental approach [24,25,26,27,28,29]. To our knowledge, the only dynamic model dealing with C. maltaromaticum strains has been published by Ellouze et al. [6]. That research was oriented towards a biological TTI setting associated with a sausage-like packaging instead of a bioreactor and taking into account lactic acid as the single metabolite.
The aim of this study was thus to develop and validate a dynamic model predicting the impact of fermentation conditions (temperature and pH) on the growth and bioconversion fermentation dynamics of C. maltaromaticum CNCM I-3298 using trehalose as a carbon source and considering the four main identified metabolites: formic acid, acetic acid, lactic acid, and ethanol. This study was conducted as part of a research project on the production and conservation of C. maltaromaticum concentrates. In that context, the growth of C. maltaromaticum was tested in different sugars: glucose, maltose, mannitol, and trehalose, with similar growth rates. Trehalose was chosen in this study because this molecule is known for its ability to protect cells during bacterial stabilization processes (freeze-drying in particular). Therefore, the residual trehalose (not consumed during fermentation) could be used as cryoprotectant after production of bacterial concentrates.
The model development involved four major steps, presented in Section 3: derivation of the main governing equations based on the known mixed-acid fermentation pathway, mass balances, and kinetic rate expressions (Section 3.1); parameter identification for each fermentation experiment (Section 3.2); construction of response surfaces of the calibrated parameters as a function of temperature and pH (Section 3.3); and final validation of the complete model. The resulting model is shown to be a useful tool in determining the optimal conditions for producing bacterium concentrates in bioreactors and for assessing the productivity of the bioconversion fermentation of sugars into metabolites of potential industrial interest (Section 4.4).

2. Materials and Methods

Data used to calibrate and validate the model were partially reported in a previous study, in which a modified central composite experimental design was carried out to study the effect of operating conditions on the technological properties of C. maltaromaticum CNCM I-3298 [7]. Sixteen lab-scale fermentations (hereafter named F01 to F16) were performed using a wide range of regulated operating conditions (Figure 2): temperature between 20 and 37 °C and pH between 6.0 and 9.5.
Fermentation durations varied between 20 h and 45 h, and the initial conditions were: for biomass (X0) 0.077 molC·L−1, trehalose (S0) between 0.091 mol·L−1 and 0.107 mol·L−1, and medium volume (V0) 3.5 L.
The main fermentation settings and the kinetic measurements are reported below.

2.1. Fermentation

2.1.1. Culture Medium and Bacterial Strain

The fermentation medium was composed of the following ingredients for 1 kg of final solution: 40 g of trehalose (Treha™; Tokyo Japan); 10 g of proteose peptone (Oxoid; Waltham, MA, USA); 5 g of yeast extract (Humeau; La-Chapelle-sur-Erdre, France); 5 g of Tween 80 (VWR; Leuven, Belgium); 0.41 g of MgSO4 (Merck; Darmstadt, Germany); 0.056 g MnSO4 (Merck; Darmstadt, Germany); and water to reach a total of 1 kg of solution. All medium components were sterilized together at 121 °C for 20 min. Fermentations were carried out on C. maltaromaticum CNCM I-3298 pre-cultures. Pre-cultures were prepared by inoculating 10 mL of sterilized fermentation medium with 100 µL of C. maltaromaticum CNCM I-3298 stock culture and were incubated for 13 to 16 h at 30 °C. An amount of 1 mL of the resulting culture was transferred into 50 mL of fresh medium and then incubated again for 11 h under the same conditions. The resulting culture was then used to inoculate the bioreactor. Inoculation was performed at an initial concentration of approximately 107 CFU mL−1.

2.1.2. Bioreactor and Parameter Control

The bioreactor (Minifors, Infors HT, Bottmingen, Switzerland) had a total volume of 5 L and was equipped with a heat mantle and a cryostat for temperature control. It contained 3.5 L of fermentation medium, inoculated with an initial cell concentration of approximately 107 CFU·mL−1. Initial pH was adjusted to the desired value with 5 M NaOH or 0.01 M H2SO4 solutions. During fermentation, pH was controlled to the desired setpoint for each investigated condition (Figure 2) by automatic addition of 5 M NaOH. Culture homogenization was performed with an agitation device set at 150 rpm. Temperature was set according to the investigated operating conditions mentioned above (Figure 2).

2.2. Kinetic Measurements

2.2.1. Cell Growth

Cell growth was monitored using an infrared probe (Excell210, CellD, Roquemaure, France) continuously measuring absorbance at 880 nm and storing data every minute. The absorbance data were calibrated in dry weight. Dry cell weight was determined by filtering 10 mL of bacterial suspension (straight out of the bioreactor) through a 0.20 µm polyethersulfone membrane (Supor®, PALL Biotech, Saint-Germain-en-Laye, France). The filter was then dried for 24 h at 80 °C. Measurements were obtained in triplicate. Mass concentrations were finally converted to molC L−1 (carbon-mol of biomass per liter), assuming the simplified unit-carbon biomass formula CH1.8O0.5 [30].

2.2.2. Total Acid Production

Total acid production was determined according to the volume of NaOH solution injected into the bioreactor to maintain a constant pH. The pH was regulated/controlled to set values using the IRIS NT V5 software (Infors, AG, Bottmingen, Switzerland).

2.2.3. Substrate Consumption and Metabolite Production

Trehalose consumption and metabolite production were determined using high-performance liquid chromatography (HPLC, Waters Associates, Millipore; Molsheim, France). HPLC was performed on culture media samples of a few mL, aseptically retrieved from the bioreactor at different times during fermentation and filtered through 0.22 µm pores (Sartorius stedim, Biotech; Göttigen, Germany). Analyses were made using a cation exchange column (Aminex Ion Exclusion HPX-87 300 × 7.8 mm, Bio-Rad, Richmond, VA, USA) at 35 °C. Mobile phase was 0.0005 M H2SO4, and flow rate was set at 0.6 mL·min−1 (LC-6A pump, Shimadzu, Courtaboeuf, France).
HPLC analysis showed that C. maltaromaticum CNCM I-3298 produced not only lactic acid but also formic acid, acetic acid, and ethanol in variable proportions according to the fermentation conditions.

3. Dynamic Model

The mathematical model was a set of ordinary differential equations implemented in MATLAB R2018b (the MathWorks Inc. Natick, MA, USA). Model parameters and response surface coefficients were identified by nonlinear regression analysis using the Statistic and Machine Learning Toolbox of MATLAB.

3.1. Model Formulation

The dynamic model developed in this study combined biochemical knowledge about the metabolism of the selected bacterium and mass balances of the main compounds: substrate, biomass, and identified metabolites. Expressions of specific growth and metabolite production rates included substrate limitation, product inhibition phenomena, and time lags due to microbial metabolism adaptation [31]. The surface response method was used to express the empiric dependence of some model parameters on operating conditions. The model assumed the bioreactor was perfectly stirred and there were no differences between individual cells. It was thus unsegregated and zero-dimensional, predicting average spatial concentrations [32].
Seven state variables were considered: six volume concentrations (biomass [X], trehalose [S], formic acid [F], acetic acid [A], lactic acid [L], and ethanol [E], Figure 1) and the culture medium volume (V). This latter variable varied continuously with the addition of base (NaOH) for pH control but also changed in a discrete way due to periodic sampling for biological and chemical analysis.
Mass balances for the considered metabolites resulted in the following set of differential equations:
d X dt = μ X X Q V X
d F dt = π F X Q V F
d A dt = π A X Q V A
d L dt = π L X Q V L
d E dt = π E X Q V E
d A T dt = d F dt + d A dt + d L dt
d S dt = μ X Y X / S   + π F Y F / S   + π A Y A / S   + π L Y L / S   + π E Y E / S   X Q V S
dV dt = Q
Here, μX is the specific growth rate (h−1); πF, πA, πL, and πE are the specific production rates of four metabolites (h−1); and YX/S, YF/S, YA/S, YL/S, and YE/S are the yield of biomass and metabolites with respect to the substrate (mol.mol−1). Q is the experimentally measured rate of NaOH solution (L.h−1) added for pH control throughout fermentation.
In Equation (6), [AT] is the total acid concentration, defined as the sum of formic, acetic, and lactic acid concentrations. These compounds are assumed to be mainly responsible for the pH change of the liquid medium.
Specific growth and production rates were defined using the Monod law to account for substrate limitation, modified with product inhibition and enzymatic activation factors [33,34,35]:
μ X = μ max , X I X E A S K SX + S
π m = π max , m I m E A S K Sm + S         m = F ,   A ,   L ,   E
In these equations, IX and Im are inhibition factors that depend on the inhibitor concentration. They vary between 1 and 0. Inhibition increases with the inhibitor concentration, and its effect on the specific rate is maximal when the corresponding factor is 0. In this model, progressive inhibition factors of the following form were used [36,37]:
I X = 1 1 + C I K IX n
I m = 1 1 + C I K Im p         m = F ,   A ,   L ,   E
KIX and KIm represent characteristic concentrations of the inhibitors (mol L−1) such that the corresponding rates (μX and πm) are reduced by a factor of 2 compared with the absence of inhibitor, n and p are shape factors, and CI is the concentration of the inhibitor. Since all the metabolites were produced in similar proportions and no biochemical knowledge about their relative inhibiting nature was available, CI was simply defined as the sum of the four metabolite concentrations:
C I = F + A + L + E
To illustrate the role of the shape factor n, Figure 3a depicts the evolution of IX with CI for different n values and a lag-time of 5 h. A more or less sharp change in the inhibition factor occurs around the characteristic inhibitor concentration, CI = KIX. The significance of the shape factor p is similar.
The enzymatic adaptation factor EA is an empirical representation of the lag time, a period of adaptation to the culture environment where the microorganism produces new enzymatic machinery [38,39,40]. Based on the shape of experimental data, the following equation was proposed:
E A = 1 1 + exp r t t lag
where tlag (h) is the lag time experimentally observed. Figure 3b shows that EA is an increasing function of time, tending to 1 when t ≫ tlag. In analogy with n, r is a shape factor that describes the gradual transition from the lag phase to the active phase of growth. A higher value of r implies a steeper change of EA around t = tlag.
To illustrate the features of the proposed model, a representation of the dimensionless specific growth and production rates (μ/μmax and π/πmax) over time is depicted in Figure 4. The dynamic behavior of both variables is similar given the similarity of Equations (1)–(5). The specific rates achieve a maximum value in the active growth phase, and they are zero when t ≪ tlag and when the substrate is depleted. The shape of the curve is defined by three factors: in the increasing region (0 to 10 h in Figure 4), the dominant effect is enzyme activation EA (Equation (14)); in the slowly decreasing region (10 to 20 h), the rate is controlled by inhibition (Equation (11) or (12)), whereas in the sharply decreasing region (20 to 22 h) it is controlled by substrate limitation, corresponding to the Monod-like factor in Equation (9) or (10).

3.2. Model Parameter Identification

The system of kinetic equations for a single fermentation experiment included 24 parameters: five yield coefficients, five inhibition parameters, five growth/production rates, five Monod-like saturation constants, three shape factors, and one lag time. Due to a limited number of experimental data and to facilitate the identification procedure, a single value was adopted for the inhibition parameter (KIm) and the Monod saturation constant (KSm) of the four identified metabolites. Moreover, 10 parameters were fixed for all experiments: the shape factors, the yield coefficients, and the Monod saturation constants (KSX and KSm). For each fermentation, lag time was determined by graphical readout. This simplification of fixing parameters independent of operating conditions is supported by two assumptions often used in the literature: (1) metabolite production yields are constant and therefore independent of culture conditions [41] and (2) the saturation constant of the Monod model depends only on the nature of the substrate [33,38], which was the same in all experiments of this study.
The remaining group of seven parameters (μmax,X, πmax,F, πmax,A, πmax,L, πmax,E, KIX, KIm) were identified for each fermentation of the experimental design by nonlinear regression. Here, the Levenberg–Marquardt algorithm [42] was used to minimize the sum of squares of the errors between experimental and predicted concentrations. However, since the ranges and the number of measurements were slightly different among the metabolites, the values compared in the least squares function were normalized by dividing by their maximum value and were weighted by the relevant number of experimental measurements.
The quality of the model representation was quantified with two error indicators, defined as follows:
Root mean square error:
RMSE = 1 N i = 1 N C model , i   -   C exp , i 2 1 / 2
Relative mean error (as a percentage):
RME = 1 N i = 1 N C model , i   -   C exp , i C exp , max   -   C exp , min   · 100 %
where N is the number of available measurements, Cmodel and Cexp are respectively the values of the concentration variables calculated with the model and measured experimentally.

3.3. Response Surface Model for Parameter Dependence on Fermentation Conditions

Nonlinear regression was performed to model the relationship between the seven parameters of the dynamic model specific to each experiment and the fermentation operating conditions—namely, temperature (T) and pH. The regression model had a similar form for all parameters, the logarithm of the parameter being expressed as a second-order polynomial with interaction:
log 10 Par i = β 0 i + β 1 i T + β 2 i pH + β 3 i T 2 + β 4 i pH 2 + β 5 i TpH
The regression coefficients (β) for all parameters depending on operating conditions (μmax,X, πmax,F, πmax,A, πmax,L, πmax,E, KIX, KIm) were simultaneously computed by least-squares optimization based on all available concentration measurements. In this way, the accuracy and standard errors of the coefficients were statistically acceptable, due to a large number of degrees of freedom: several hundreds of concentration data were used to estimate 42 coefficients. Initial guesses for these coefficients were obtained using Equation (17), and parameter values were determined separately for each experiment.
In this procedure, two sets of data from the experimental design were defined as indicated in Figure 2: eight calibration experiments, located in extreme positions of the experimental domain, used simultaneously for coefficients (β) estimation, and eight validation experiments, only used a posteriori to verify the accuracy of the complete dynamic model.

4. Results and Discussion

4.1. Model Parameter Identification

The values of the parameters that are independent of operating conditions, summarized in Table 1, were determined from the experimental data of experiment F10. This run was placed in a central position in the composite experimental design (T = 30 °C, pH = 8) (Figure 2). Monod saturation constants are usually difficult to determine from batch experiments because the number of measurements is typically very low in the substrate limitation zone. Saturation constants were thus fixed to a common value with a typical order of magnitude [43]. As for yields, they were found to differ from the theoretical ones defined through standard stoichiometric reactions of anabolism and catabolism. These differences can be due to other reactions involving the carbon substrate, whose products were not analytically measured and were not considered in the model.
After fixing the parameters in Table 1 for the whole set of experiments, the group of seven adjustable parameters of the model (μmax,X, πmax,F, πmax,A, πmax,L, πmax,E, KIX, KIm) were identified for each run by nonlinear regression.
The parameters obtained by this procedure are summarized in Table 2. Standard errors were computed from the variance–covariance matrix of the nonlinear optimization algorithm. These errors represented between 5% and 13% of the value of the identified parameters, a reasonable uncertainty level for a biological model.
For the whole set of experiments, the prediction errors are reported in Appendix A Table A1. Except for some runs for variables S, F, and A, all RME were lower than 15%. Additionally, the average RMSE and RME values for each concentration were of the same order of magnitude as the experimental variability, here defined as the biological repeatability for run F01, for which three independent replicates were performed. These results validate the formulation and accuracy of the proposed model under the operating conditions included in the experimental design.
In the specific case of reference run F10, a comparison between the model simulation (using the corresponding parameters from Table 2) and experimental data is illustrated in Figure 5.
Three growth phases are apparent in Figure 5: a lag phase (phase 1, between 0 and 10 h); a phase of active growth, substrate consumption, and metabolite production (phase 2, between 10 h and 21 h); and a final phase where concentrations do not change over time, owing to the depletion of the carbon source or growth inhibition by metabolites (phase 3, after 21 h). Regarding culture volume evolution, as already mentioned, the discrete variations at regular intervals were due to sampling for analysis of the culture medium and the gradual increase was due to NaOH addition for pH control. One can also observe that the four metabolites were produced simultaneously, with no gap for the growth dynamics. The metabolites were thus primary end products generated during a single trophophase [44]. This justifies the choice of a global inhibitor concentration (Equation (13)), which included four correlated concentrations.
In consideration of the visual fit from Figure 5, the model representation is reasonably satisfactory. The most pronounced discrepancy between the model and experimental data appears for lactic acid, for which the model predicted a lower concentration before substrate depletion. This is related to a slightly underestimated yield factor YL/S.

4.2. Response Surface Model for Parameter Dependence on Fermentation Conditions

Model parameters were expressed as a function of temperature and pH, according to the surface model (Equation (17)). The values of the β regression coefficients were adjusted globally using the whole set of calibration data.
The resulting response surfaces for the seven model parameters are plotted in Figure 6. For the five kinetic parameters, (i.e., the maximum specific growth and production rates), the response surfaces have the same convex shape, with a well-defined maximum value at intermediate T and pH conditions. These maxima likely indicate the optimal temperatures and pH for cellular growth, as well as the enzymatic activity catalyzing each of the reactions, leading to the production of the different metabolites (Figure 1).
Concerning the inhibition concentrations, the response surface for KIm has a concave shape with a local minimum, whereas that of KIX resembles a saddle surface. For this latter case, the surface shape indicates that for every pH there is a T where KIX is minimal, and for every T there is a pH where KIX is maximal. Both KIm and KIX represent the combined effect of several inhibiting metabolites (Equations (11)–(13)) with potentially different inhibition mechanisms.
For completeness, the final values of the regression coefficients of Equation (17) for the seven adjustable parameters of the dynamic model are reported in Appendix A Table A2. All coefficients in Equation (17) for each model parameter were significantly different from zero at a 0.05 level. A comparison between the parameter values determined for each experiment (Section 4.1) and the parameter values computed with Equation (17) (from globally adjusted β coefficients) is depicted in Appendix A Figure A1. The goodness of the fit was assessed through the coefficient of determination, R2. This coefficient is higher than 0.89 for six out of seven model parameters, which is a high threshold for biological data. In the case of KIm, only 66% of the variance of this parameter was explained by variables T and pH. The remaining 34% could be associated with inherent experimental variability and factors not included in the model, for instance transient variability of the inhibition and kinetics parameters and actual dependence of the fixed parameters (Table 1) with T and pH [45]. From a more general point of view, differences from experimental data could be due to features that were not represented by the mathematical model, such as population segregation, internal pH variability, and concentration gradients in the culture medium [46,47].

4.3. Model Validation

The ability of the dynamic model including the parameters calculated from operating conditions (Equation (17)) to predict data of independent experiments was assessed with a set of validation experiments.
A comparison between the average RMSE values obtained in Section 4.1 (determined for each experiment) and Section 4.2 (calculated from operating conditions) for calibration and validation sets is depicted in Figure 7. In most cases, RMSE values were higher than the corresponding experimental variabilities, indicating that more complex models could capture additional phenomena not included in the present model, such as dependence of yields, saturation constants, or lag time (Table 1) on operating conditions. As one might expect, RMSE was generally lower for the calibration experiments than for the validation experiments, not used for parameter determination. However, the relative difference remained small (less than 30%), indicating a satisfactory ability of the developed model to predict time evolution of the considered biomass, substrate, and metabolites under new conditions within the explored experimental range.
It also appears in Figure 7 that average RMSE values with parameters given by the response surface model (Table A3) are about 50% higher than with parameters determined separately for each experiment (Table A1), for both calibration and validation sets. This result could be expected since in the global calibration step, data from eight independent experiments were combined as a whole for the least squares estimation, with a detrimental effect on the individual representation of each experiment. However, results with the parameters calculated from operating conditions are the most useful in engineering purposes since they enable a quick prediction of growth and metabolites production dynamics, based on the selected combination of temperature and pH.
In light of this quantitative analysis, the prediction accuracy of the empirical dynamic model coupled to the regression model may be considered satisfactory within the operating domain covered in this study.

4.4. Model-Based Optimization of Fermentation Operating Conditions for Industrial Use

Optimal conditions for growth and metabolite production of C. maltaromaticum calculated using the developed model are summarized in Table 3. Two optimization criteria were considered: final concentrations and final productivities calculated for a 99.9% substrate consumption.
For a detailed representation of the evolution of final concentrations and productivities for biomass and metabolites with temperature and pH, the reader is referred to Appendix A Figure A2 and Figure A3. As a general trend, the highest productivities were obtained around 35 °C and pH 7.5, although the exact optimal conditions depended on the considered metabolite (Table 3). No general trend was readily apparent for the maximization of the final concentrations.
These data can be useful in optimizing industrial processes involving the growth of C. maltaromaticum cells in a trehalose-based substrate. A first application consists of producing C. maltaromaticum concentrates, regardless of metabolite production. In this case two conditions of cultivation appear advisable: 20 °C and pH 7.8 to maximize concentration (227 mmolC·L−1) or 33.5 °C and pH 7.5 in order to maximize productivity (7.49 mmolC·L−1·h−1) and thus the biomass production per unit of time, at the expense of a 12% reduction of the final biomass concentration (199 mmolC·L−1).
A second application deals with the development and parametrization of time–temperature integrators (TTI), labels in which a pH decline, associated with acids synthesis, entails an irreversible color change from green to red. Modulating the acidifying activity of C. maltaromaticum thus allows a reliable shelf-life estimation of different food products. Long shelf-lives can be tracked using TTI composed of concentrates exhibiting low acidifying activities (minimal production of total acids), while short shelf-lives can be tracked using concentrates exhibiting high acidifying activities. In the scenario of maximizing acidifying activity, the production of total acids must be favored, and thus fermentation should be carried out under two possible conditions: 37.0 °C and pH 6.0 to maximize their final concentration (433 mmol·L−1) or 35.5 °C and pH 7.7 to maximize their productivity (14.90 mmol·L−1·h−1). Under these conditions, the biomass production decreases respectively by 20% and 4% with respect to its optimal values. If the objective is, on the contrary, to minimize acidifying activity, two conditions can be envisaged to favor the lowest production of total acids: 37.0 °C and pH 9.5 for a final concentration of 352 mmol L−1 or 25.0 °C and pH 9.5 for a final productivity of 4.53 mmol·L−1·h−1. Under these conditions, the mean biomass production would decrease respectively by 48% and 77% with respect to the maximal values.
Data from Table 3 show that the conditions to minimize the total acids concentration (37.0 °C and pH 9.5) coincide with those to maximize the ethanol concentration (the non-acidifying metabolite, 139 mmol L−1) and lead to a lactic acid concentration close to its minimal value (128 mmol L−1 versus the minimum around 120 mmol L−1). Conversely, when the production of total acids is maximized, the lactic acid concentration is also maximal (296 mmol L−1) and that of ethanol is close to its minimum (72 mmol L−1 versus 68 mmol L−1).
Furthermore, it should be noted that the condition 27 °C and pH 7.6 leads both to a good biomass productivity (6.93 mmol·L−1·h−1 versus the maximum 7.49 mmol·L−1·h−1) and a low total acids concentration (372 mmol L−1 versus the minimum 352 mmol L−1). Cultivation under this condition turns out be advantageous to ally a high biomass production and a relatively low total acidification.

5. Conclusions

The dynamic model developed in this study is able to predict with satisfactory accuracy the growth of C. maltaromaticum CNCM I-3298 (average error of 10%) as well as the conversion of trehalose into four primary metabolites (average error of 14%) under a wide range of conditions of temperature and pH. The interpolation capability of the model was verified with a set of eight independent validation experiments, for which the average relative error was 13%.
This model constitutes a useful tool for optimizing C. maltaromaticum cultures. Based on two easily controllable parameters, pH and temperature, it could be implemented in industrial applications of food technology to define optimal growth and metabolite production conditions with various objectives, such as the maximization of biomass for production of bacterium concentrates or the maximization or minimization of the acidifying activity. A typical operating condition for this bacterium could be, for instance, 30.0 °C and pH 7.0. If the goal is to produce bacterium concentrates, to maximize final biomass concentration, our results suggest that a quite different condition should be selected (20.0 °C and pH 7.8), while for maximum biomass productivity, 33.5 °C and pH 7.5 is most appropriate. Such results are quite difficult to anticipate from the qualitative knowledge of the bacterium alone, and a large number of time-consuming experiments would be required to locate these optimal conditions experimentally without constructing a dynamic model of the process.
The effort of developing the model is especially cost effective when a variety of scenarios are explored. If the goal is to develop time–temperature integrators (TTI) to track the cold-chain of food products, a set of labels with specific shelf-lives has to be designed for various target products. The range of desired shelf-lives can be as large as 1 to 30 days, requiring very different TTI designs. In a traditional approach, for each desired shelf-life duration, a range of factors such as the initial bacterium concentration and the buffer properties of the medium have to be explored in a series of relatively time-consuming experiments. In such an environment, temperature varies in an arbitrary but known way, and pH depends on the produced acids. The presented dynamic model can be extended to predict the moment when a specific amount of acids is produced, corresponding to the pH-induced color change of the TTI label and hence to the desired shelf-life. Model-based design of the TTI labels is expected to be faster and more accurate than a trial and error procedure.
On a more fundamental level, further work is required to incorporate the effect of other culture parameters, such as aeration, nutrient concentrations, or the use of a different carbon source, which may modify growth kinetics and metabolite production. Additionally, it would be relevant to deepen the understanding of inhibition mechanisms of the metabolites to give more biological significance to the associated parameters in the model.

Author Contributions

Conceptualization, C.P., A.G., F.F. and I.-C.T.; methodology, C.P., A.G., F.F., and I.-C.T.; software, C.P. and I.-C.T.; validation, C.P., A.G., S.P., F.F. and I.-C.T.; formal analysis, C.P. and I.-C.T.; investigation, C.P. and A.G.; resources, F.F. and I.-C.T.; data curation, C.P. and A.G.; writing—original draft preparation, C.P. and A.G.; writing—review and editing, C.P., A.G., S.P., F.F. and I.-C.T.; visualization, C.P., A.G. and I.-C.T.; supervision, S.P., F.F. and I.-C.T.; project administration, S.P., F.F. and I.-C.T.; funding acquisition, S.P. and F.F. All authors have read and agreed to the published version of the manuscript.

Funding

This study received funding from the European Union’s Horizon 2020 Marie Skłodowska-Curie research and innovation program under grant agreement no. 777657.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to an ongoing research project.

Acknowledgments

The authors thank Marie-Nöelle Lecercq-Perlat for HPLC analysis.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Nomenclature

AAcetic acid
EEthanol
FFormic acid
LLactic acid
SCarbon substrate (trehalose)
XBiomass
EAEnzymatic activation factor
C (mol·L−1)Concentration (in the calculation of errors and the definition of the inhibition factors)
[i] (mol·L−1)Concentration of species i (substrate, metabolite, biomass) in the culture medium (in the system of differential equations)
ImProduction inhibition factor of metabolite m
IXBiomass growth inhibition factor
KIm (mol·L−1)Concentration for 50% production rate inhibition of metabolite m
KIX (molC·L−1)Concentration for 50% growth rate inhibition of biomass
KSm (mol·L−1)Concentration of production rate saturation of metabolite m
KSX (molC·L−1)Concentration of biomass growth rate saturation
molCCarbon-mol of biomass
nShape factor of the growth inhibition function
pShape factor of the production inhibition function
pHPotential of hydrogen
Q (L·h−1)Rate of base addition for pH control
RShape factor of the enzymatic activation function
T (K)Temperature
TTITime–temperature indicator
RMSERoot-mean square error
RMERelative mean error
SEStandard error
t (h)Time
tlag (h)Lag time
V (L)Culture medium volume
Yi/S (mol·mol−1)Yield of product i on substrate S
YX/S (molC·mol−1)Biomass yield on substrate S
μX (h−1)Specific growth rate
μmax,X (h−1)Maximal specific growth rate
πm (h−1)Specific production rate of metabolite m
πmax,m (h−1)Maximum specific production rate of metabolite m

Appendix A

Appendix A.1. Model Parameter Identification

Table A1. Residual modelling error with model parameters determined for each experiment and summarized in Table 2.
Table A1. Residual modelling error with model parameters determined for each experiment and summarized in Table 2.
RMSE (Mol L−1)RME (%)
XSFALEA TotalXSFALEA Total
F010.0080.0060.0100.0070.0220.0070.033447136107
F020.0050.0070.0060.0030.0040.0030.00541688873
F030.0050.0060.0030.0010.0030.0010.0064932321
F040.0060.0060.0030.0010.0050.0020.00571355663
F050.0080.0050.0140.0070.0270.0070.034449111098
F060.0060.0040.0120.0050.0240.0060.0125277883
F070.0060.0070.0110.0060.0090.0060.0074514334112
F080.0060.0070.0110.0050.0300.0040.003451723871
F090.0040.0050.0020.0010.0020.0020.00371767562
F100.0050.0090.0060.0040.0280.0040.0263747956
F110.0030.0180.0090.0050.0110.0060.0323279109116
F120.0050.0130.0020.0010.0030.0010.004136134822
F130.0020.0060.0120.0050.0240.0050.0072476652
F140.0030.0030.0090.0050.0120.0060.0154316205124
F150.0090.0050.0090.0060.0230.0050.01154713882
F160.0030.0010.0010.0010.0080.0010.0142124623
Mean0.0050.0070.0080.0040.0150.0040.014511811774
Response surface model for parameter dependence on fermentation conditions.
Table A2. Response surface coefficients fitted to experimental data by multiple regression.
Table A2. Response surface coefficients fitted to experimental data by multiple regression.
VariableCoefficientµmaxX (h−1)πmaxF (h−1)πmaxA (h−1)πmaxL (h−1)πmaxE (h−1)KIX (Mol L−1)KIm (Mol L−1)
ConstantValue (β0)−1.38 × 101−1.64 × 101−1.74 × 101−1.21 × 101−1.12 × 101−8.61 × 10−11.19 × 101
Standard error0.05 × 1010.07 × 1010.09 × 1010.09 × 1010.07 × 1010.05 × 10−10.03 × 101
TValue (β1)2.63 × 10−12.66 × 10−12.71 × 10−11.82 × 10−11.27 × 10−1−1.87 × 10−1−5.80 × 10−1
Standard error0.13 × 10−10.26 × 10−10.28 × 10−10.23 × 10−10.14 × 10−10.05 × 10−10.17 × 10−1
pH Value (β2)2.42 × 1002.89 × 1003.01 × 1002.27 × 1002.02 × 1006.65 × 10−1−1.27 × 100
Standard error0.11 × 1000.14 × 1000.19 × 1000.16 × 1000.13 × 1000.45 × 10−10.16 × 100
Value (β3)−4.26 × 10−3−3.45 × 10−3−3.26 × 10−3−2.05 × 10−3−1.92 × 10−33.97 × 10−37.95 × 10−3
Standard error0.13 × 10−30.32 × 10−30.33 × 10−30.81 × 10−30.22 × 10−30.09 × 10−30.32 × 10−3
pH²Value (β4)−1.46 × 10−1−1.63 × 10−1−1.66 × 10−1−1.34 × 10−1−1.22 × 10−1−4.03 × 10−24.35 × 10−2
Standard error0.06 × 10−10.10 × 10−10.12 × 10−10.10 × 10−10.09 × 10−10.31 × 10−20.55 × 10−2
T·pHValue (β5)−3.53 × 10−3−7.58 × 10−3−9.17 × 10−3−7.22 × 10−3−4.68 × 10−3−3.54 × 10−31.86 × 10−2
Standard error0.11 × 10−30.97 × 10−31.17 × 10−30.91 × 10−30.05 × 10−30.06 × 10−30.01 × 10−2
Figure A1. Comparison between model parameters determined for each experiment (Table 2) and parameters computed with the response surface models (Equation (17) using coefficients in Table A2).
Figure A1. Comparison between model parameters determined for each experiment (Table 2) and parameters computed with the response surface models (Equation (17) using coefficients in Table A2).
Foods 10 01922 g0a1

Appendix A.2. Model Validation

Error indicators RMSE and RME for calibration and validation sets are summarized in Appendix A Table A3. The RME values varied from 1% to 53%, with an average of 15% for calibration set and 13% for validation set. Likewise, the average RME was lower than 15% for most of experiments, except for runs F12 (17%), F15 (18%), F09 (20%), and F16 (34%). Considering both calibration and validation sets, the average RME was 10% for biomass and 14% for substrate and metabolites.
Table A3. Quality of fit of the model with parameters computed with the response surface models.
Table A3. Quality of fit of the model with parameters computed with the response surface models.
FermentationRMSE (Mol L−1)RME (%)
XSFALEAtXSFALEAt
CalibrationF010.0100.0100.0210.0100.0410.0140.03356142010187
F030.0210.0090.0080.0050.0080.0030.012151479754
F040.0050.0100.0090.0060.0070.0050.01552618239179
F070.0120.0090.0100.0060.0150.0060.0099813315103
F120.0030.0110.0080.0050.0030.0040.00985313168116
F130.0100.0070.0060.0030.0360.0030.0117544842
F140.0150.0020.0130.0070.0040.0090.01220223282163
F160.0110.0230.0330.0140.0640.0180.0859335045483818
Mean0.0110.0100.0140.0070.0220.0080.0231018182212157
ValidationF020.0050.0030.0170.0090.0080.0090.00758202215214
F050.0110.0080.0070.0060.0450.0050.01677481384
F060.0080.0090.0130.0050.0440.0070.03265761196
F080.0090.0100.0120.0060.0260.0080.0037819258141
F090.0110.0060.0050.0040.0180.0040.0091922152435157
F100.0210.0140.0250.0110.0630.0140.04411101718202010
F110.0180.0130.0070.0030.0170.0040.0501620761289
F150.0280.0230.0260.0100.0470.0140.0361520222117217
Mean0.0140.0110.0140.0070.0330.0080.0251113141616156
Model-based optimization of fermentation conditions.
Figure A2. Evolution of final concentrations (left) and batch-average productivities (right) with temperature and pH for biomass, formic acid, and acetic acid.
Figure A2. Evolution of final concentrations (left) and batch-average productivities (right) with temperature and pH for biomass, formic acid, and acetic acid.
Foods 10 01922 g0a2
Figure A3. Evolution of final concentrations (left) and batch-average productivities (right) with temperature and pH for lactic acid, ethanol, and total acids.
Figure A3. Evolution of final concentrations (left) and batch-average productivities (right) with temperature and pH for lactic acid, ethanol, and total acids.
Foods 10 01922 g0a3

References

  1. Edima, H.; Cailliez-Grimal, C.; Revol-Junelles, A.-M.; Tonti, L.; Linder, M.; Millière, J.B. A selective enumeration medium for Carnobacterium maltaromaticum. J. Microbiol. Methods 2007, 68, 516–521. [Google Scholar] [CrossRef]
  2. Afzal, M.; Jacquet, T.; Delaunay, S.; Borges, F.; Millière, J.B.; Revol-Junelles, A.M.; Cailliez-Grimal, C. Carnobacterium maltaromaticum: Identification, isolation tools, ecology and technological aspects in dairy products. Food Microbiol. 2010, 27, 573–579. [Google Scholar] [CrossRef]
  3. Cailliez-Grimal, C.; Edima, H.; Revol-Junelles, A.-M.; Millière, J.-B. Carnobacterium maltaromaticum: The Only Carnobacterium Species in French Ripened Soft Cheeses as Revealed by Polymerase Chain Reaction Detection. J. Dairy Sci. 2007, 90, 1133–1138. [Google Scholar] [CrossRef]
  4. Edima, H.; Cailliez-Grimal, C.; Revol-Junelles, A.M.; Rondags, E.; Millière, J.B. Impact of pH and Temperature on the Acidifying Activity of Carnobacterium Maltaromaticum. J. Dairy Sci. 2008, 91, 3806–3813. [Google Scholar] [CrossRef] [Green Version]
  5. Taoukis, P.; Labuza, T. Applicability of Time-Temperature Indicators as Shelf Life Monitors of Food Products. J. Food Sci. 1989, 54, 783–788. [Google Scholar] [CrossRef]
  6. Ellouze, M.; Pichaud, M.; Bonaiti, C.; Coroller, L.; Couvert, O.; Thuault, D.; Vaillant, R. Modelling pH evolution and lactic acid production in the growth medium of a lactic acid bacterium: Application to set a biological TTI. Int. J. Food Microbiol. 2008, 128, 101–107. [Google Scholar] [CrossRef]
  7. Girardeau, A.; Puentes, C.; Keravec, S.; Peteuil, P.; Trelea, I.C.; Fonseca, F. Influence of culture conditions on the technological properties of Carnobacterium maltaromaticum CNCM I-3298 starters. J. Appl. Microbiol. 2019, 126, 1468–1479. [Google Scholar] [PubMed]
  8. Borch, E.; Molin, G. The aerobic growth and product formation of Lactobacillus, Leuconostoc, Brochothrix, and Carnobacterium in batch cultures. Appl. Microbiol. Biotechnol. 1989, 30, 81–88. [Google Scholar] [CrossRef]
  9. dos Reis, F.B.; de Souza, V.M.; Thomaz, M.R.S.; Fernandes, L.P.; de Oliveira, W.P.; de Martinis, E.C.P. Use of Carnobacterium maltaromaticum cultures and hydroalcoholic extract of Lippia sidoides Cham. Against Listeria monocytogenes in fish model systems. Int. J. Food Microbiol. 2011, 146, 228–234. [Google Scholar] [CrossRef] [PubMed]
  10. Afzal, M.I.; Delaunay, S.; Paris, C.; Borges, F.; Revol-Junelles, A.M.; Cailliez-Grimal, C. Identification of metabolic pathways involved in the biosynthesis of flavor compound 3-methylbutanal from leucine catabolism by Carnobacterium maltaromaticum LMA 28. Int. J. Food Microbiol. 2012, 157, 332–339. [Google Scholar] [CrossRef]
  11. Afzal, M.I.; Boulahya, K.A.; Paris, C.; Delaunay, S.; Cailliez-Grimal, C. Effect of oxygen on the biosynthesis of flavor compound 3-methylbutanal from leucine catabolism during batch culture in Carnobacterium maltaromaticum LMA 28. J. Dairy Sci. 2013, 96, 352–359. [Google Scholar] [CrossRef]
  12. Rahman, A.; Gleinser, M.; Lanhers, M.-C.; Riedel, C.U.; Foligne, B.; Hanse, M.; Yen, F.T.; Klouj, A.; Afzal, M.I.; Back, A.; et al. Adaptation of the lactic acid bacterium Carnobacterium maltaromaticum LMA 28 to the mammalian gastrointestinal tract: From survival in mice to interactions with human cells. Int. Dairy J. 2014, 34, 93–99. [Google Scholar]
  13. Rahman, A.; El Kheir, S.M.; Back, A.; Mangavel, C.; Revoljunelles, A.-M.; Borges, F. Repeat-based bequence typing of Carnobacterium maltaromaticum. Int. J. Food Microbiol. 2016, 226, 1–4. [Google Scholar] [CrossRef] [PubMed]
  14. Zhang, P.; Badoni, M.; Gänzle, M.; Yang, X. Growth of Carnobacterium spp. isolated from chilled vacuum-packaged meat under relevant acidic conditions. Int. J. Food Microbiol. 2018, 286, 120–127. [Google Scholar]
  15. De Bruyn, I.N.; Louw, A.I.; Visser, L.; Holzapfel, W.H. Lactobacillus divergens is a homofermentative organism. Syst. Appl. Microbiol. 1987, 9, 173–175. [Google Scholar] [CrossRef]
  16. De Bruyn, I.N.; Holzapfel, W.H.; Visser, L.; Louw, A.I. Glucose metabolism by Lactobacillus divergens. J. Gen. Microbiol. 1988, 134, 2103–2109. [Google Scholar] [CrossRef] [Green Version]
  17. Leisner, J.J.; Laursen, B.G.; Prévost, H.; Drider, D.; Dalgaard, P. Carnobacterium: Positive and negative effects in the environment and in foods. FEMS Microbiol. Rev. 2007, 31, 592–613. [Google Scholar] [CrossRef] [Green Version]
  18. Schillinger, U.; Holzapfel, W.H. The genus Carnobacterium. In The Genera of Lactic Acid Bacteria. The Lactic Acid Bacteria; Wood, B.J.B., Holzapfel, W.H., Eds.; Springer: Boston, MA, USA, 1995; Volume 2, pp. 307–326. [Google Scholar]
  19. Gänzle, M.G. Lactic metabolism revisited: Metabolism of lactic acid bacteria in food fermentations and food spoilage. Curr. Opin. Food Sci. 2015, 2, 106–117. [Google Scholar] [CrossRef]
  20. Ward, B. Bacterial energy metabolism. In Molecular Medical Microbiology, 2nd ed.; Tang, Y.W., Liu, D., Schwartzman, J., Sussman, M., Poxton, I., Eds.; Elsevier: Waltham, MA, USA, 2015; pp. 201–233. [Google Scholar]
  21. Laursen, B.G.; Leisner, J.J.; Dalgaard, P. Carnobacterium Species: Effect of Metabolic Activity and Interaction with Brochothrix thermosphacta on Sensory Characteristics of Modified Atmosphere Packed Shrimp. J. Agric. Food Chem. 2006, 54, 3604–3611. [Google Scholar] [CrossRef]
  22. Zhang, P.; Gänzle, M.; Yang, X. Complementary Antibacterial Effects of Bacteriocins and Organic Acids as Revealed by Comparative Analysis of Carnobacterium spp. from Meat. Appl. Environ. Microbiol. 2019, 85, e01227-19. [Google Scholar] [CrossRef] [PubMed]
  23. Andersson, U.; Rådström, P. Beta-Glucose 1-phosphate-interconverting enzymes in maltose- and trehalose-fermenting lactic acid bacteria. Environ. Microbiol. 2002, 4, 81–88. [Google Scholar] [CrossRef]
  24. Levisauskas, D.; Galvanauskas, V.; Henrich, S.; Wilhelm, K.; Volk, N.; Lübbert, A. Model-based optimization of viral capsid protein production in fed-batch culture of recombinant Escherichia coli. Bioprocess. Biosyst. Eng. 2003, 25, 255–262. [Google Scholar] [CrossRef]
  25. Antelo, L.; Passot, S.; Fonseca, F.; Trelea, I.C.; Alonso, A.A. Towards optimal operation conditions of freeze-drying processes via a multi-level approach. Dry. Technol. 2012, 30, 1432–1448. [Google Scholar] [CrossRef] [Green Version]
  26. Kiparissides, A.; Pistikopoulos, E.N.; Mantalaris, A. On the model-based optimization of secreting mammalian cell (GS-NS0) cultures. Biotechnol. Bioeng. 2015, 112, 536–548. [Google Scholar] [CrossRef] [PubMed]
  27. Balsa-Canto, E.; Alonso, A.A.; Arias-Méndez, A.; Garcia, M.; López-Núñez, A.; Mosquera-Fernández, M.; Vázquez, C.; Vilas, C. Modeling and optimization techniques with applications in food processes, bio-processes and bio-systems. In Numerical Simulation in Physics and Engineering; Higueras, I., Roldán, T., Torrens, J., Eds.; Springer: Cham, Switzerland, 2016; Volume 9, pp. 187–216. [Google Scholar]
  28. Sanchez-Castaneda, A.; Athès, V.; Moussa, M.; Lopez-Miranda, J.; Paez-Lerma, J.; Soto-Cruz, N.; Trelea, I.-C. Modeling of isoamyl acetate production by fermentation with Pichia fermentans in an aerated system coupled to in situ extraction. Process. Biochem. 2018, 65, 11–20. [Google Scholar] [CrossRef]
  29. Sinner, P.; Kager, J.; Daume, S.; Herwig, C. Model-based analysis and optimisation of a continuous Corynebacterium glutamicum Bioprocess Utilizing Lignocellulosic Waste. IFAC-Pap. 2019, 52, 181–186. [Google Scholar] [CrossRef]
  30. Battley, E.H. Growth-Reaction Equations for Saccharomyces cerevisiae. Physiol. Plant. 1960, 13, 192–203. [Google Scholar] [CrossRef]
  31. Bastin, G.; Dochain, D. On-Line Estimation and Adaptive Control of Bioreactors; Elsevier: Amsterdam, The Netherlands, 1990; pp. 1–82. [Google Scholar]
  32. Bailey, J.E. Mathematical modeling and analysis in biochemical engineering: Past accomplishments and future opportunities. Biotechnol. Prog. 1998, 14, 8–20. [Google Scholar] [CrossRef]
  33. Monod, J. Recherches sur la Croissance des Cultures Bactériennes; Hermann & Cie: Paris, France, 1942. [Google Scholar]
  34. Levenspiel, O. The monod equation: A revisit and a generalization to product inhibition situations. Biotechnol. Bioeng. 1980, 22, 1671–1687. [Google Scholar] [CrossRef]
  35. Baranyi, J.; Roberts, T.A. A dynamic approach to predicting bacterial growth in food. Int. J. Food Microbiol. 1994, 23, 277–294. [Google Scholar] [CrossRef]
  36. Aiba, S.; Shoda, M. Reassessment of the product inhibition in alcohol fermentation. J. Ferment. Technol. 1969, 47, 790–803. [Google Scholar]
  37. Claret, C.; Bories, A.; Soucaille, P. Inhibitory effect of dihydroxyacetone on Gluconobacter oxydans: Kinetic aspects and expression by mathematical equations. J. Ind. Microbiol. 1993, 11, 105–112. [Google Scholar] [CrossRef]
  38. Monod, J. The growth of bacterial cultures. Annu. Rev. Microbiol. 1949, 3, 371–394. [Google Scholar] [CrossRef] [Green Version]
  39. Baranyi, J. Comparison of Stochastic and Deterministic Concepts of Bacterial Lag. J. Theor. Biol. 1998, 192, 403–408. [Google Scholar] [CrossRef] [PubMed]
  40. Rolfe, M.; Rice, C.; Lucchini, S.; Pin, C.; Thompson, A.; Cameron, A.; Alston, M.; Michael, F.; Stringer, M.; Betts, R.; et al. Lag Phase Is a Distinct Growth Phase that Prepares Bacteria for Exponential Growth and Involves Transient Metal Accumulation. J. Bacteriol. 2012, 194, 686–701. [Google Scholar] [CrossRef] [Green Version]
  41. Bouville, M. Fermentation kinetics including product and substrate inhibitions plus biomass death: A mathematical analysis. Biotechnol. Lett. 2007, 29, 737–741. [Google Scholar] [CrossRef] [Green Version]
  42. Seber, G.; Wild, C. Nonlinear Regression, 2nd ed.; Wiley: Hoboken, NJ, USA, 2003. [Google Scholar]
  43. Owens, J.D.; Legan, J.D. Determination of the Monod substrate saturation constant for microbial growth. FEMS Microbiol. Rev. 1987, 3, 419–432. [Google Scholar] [CrossRef]
  44. Brock, T.D.; Madigan, M.T. Biology of Microorganisms, 6th ed.; Prentice Hall: London, UK, 1991; p. 353. [Google Scholar]
  45. Rodríguez, J.; Lema, J.M.; Kleerebezem, R. Energy-based models for environmental biotechnology. Trends Biotechnol. 2008, 26, 366–374. [Google Scholar] [CrossRef]
  46. Bellgardt, K.H. Bioprocess Models. In Bioreaction Engineering: Modeling and Control; Schügerl, K., Bellgardt, K.H., Eds.; Springer: Berlin/Heidelberg, Germany, 2000; pp. 44–105. [Google Scholar]
  47. González-Cabaleiro, R.; Lema, J.M.; Rodríguez, J. Metabolic Energy-Based Modelling Explains Product Yielding in Anaerobic Mixed Culture Fermentations. PLoS ONE 2015, 10, e0126739. [Google Scholar] [CrossRef]
Figure 1. Mixed-acid fermentation pathway likely used by C. maltaromaticum to ferment trehalose. End products are shown in blue. ACK, acetate kinase; ADH, acetaldehyde dehydrogenase; LDH, lactate dehydrogenase; PFL, pyruvate formate lyase; PTA, phosphate acetyltransferase; PYK, pyruvate kinase; TreH, neutral trehalose. Adapted from [19,20,21,22,23].
Figure 1. Mixed-acid fermentation pathway likely used by C. maltaromaticum to ferment trehalose. End products are shown in blue. ACK, acetate kinase; ADH, acetaldehyde dehydrogenase; LDH, lactate dehydrogenase; PFL, pyruvate formate lyase; PTA, phosphate acetyltransferase; PYK, pyruvate kinase; TreH, neutral trehalose. Adapted from [19,20,21,22,23].
Foods 10 01922 g001
Figure 2. Operating conditions of C. maltaromaticum CNCM I-3298 fermentation experiments.
Figure 2. Operating conditions of C. maltaromaticum CNCM I-3298 fermentation experiments.
Foods 10 01922 g002
Figure 3. (a) Example of inhibition factor IX as a function of CI for different n values and KIX = 0.2 mol L−1. (b) Example of enzymatic activation factor EA as a function of t for different r values, and tlag = 5 h.
Figure 3. (a) Example of inhibition factor IX as a function of CI for different n values and KIX = 0.2 mol L−1. (b) Example of enzymatic activation factor EA as a function of t for different r values, and tlag = 5 h.
Foods 10 01922 g003
Figure 4. Typical evolution of the relative production rate over time.
Figure 4. Typical evolution of the relative production rate over time.
Foods 10 01922 g004
Figure 5. Evolution of concentrations over time for experiment F10 (T = 30 °C, pH = 8). Comparison between model (continuous line, using parameters from Table 1 and Table 2) and experimental data (symbols). The error bars for data represent the biological standard deviation, calculated from three independent repetitions of the run F01.
Figure 5. Evolution of concentrations over time for experiment F10 (T = 30 °C, pH = 8). Comparison between model (continuous line, using parameters from Table 1 and Table 2) and experimental data (symbols). The error bars for data represent the biological standard deviation, calculated from three independent repetitions of the run F01.
Foods 10 01922 g005
Figure 6. Response surfaces for model parameters, calculated with globally adjusted β coefficients in Equation (17).
Figure 6. Response surfaces for model parameters, calculated with globally adjusted β coefficients in Equation (17).
Foods 10 01922 g006
Figure 7. Comparison between experimental variability and average RMSE values for concentrations computed using parameters determined for each experiment (Table 2) and the response surface models (Table A2 and Equation (17)).
Figure 7. Comparison between experimental variability and average RMSE values for concentrations computed using parameters determined for each experiment (Table 2) and the response surface models (Table A2 and Equation (17)).
Foods 10 01922 g007
Table 1. Model parameters independent of operating conditions, determined from the experimental data of experiment F10 (T = 30 °C, pH = 8) with tlag = 10 h.
Table 1. Model parameters independent of operating conditions, determined from the experimental data of experiment F10 (T = 30 °C, pH = 8) with tlag = 10 h.
ParameterConstant Value
YX/S (molC.mol−1)6.9
YF/S (mol.mol−1)5.6
YA/S (mol.mol−1)3.8
YL/S (mol.mol−1)7.0
YE/S (mol.mol−1)4.7
KSX (mol L−1)0.001
KSm (mol L−1)0.001
n3
p1
r (h−1)0.8
Table 2. Model parameters determined for each experiment by nonlinear regression.
Table 2. Model parameters determined for each experiment by nonlinear regression.
FermentationµmaxX (h−1)πmaxF (h−1)πmaxA (h−1)πmaxL (h−1)πmaxE (h−1)KIX (Mol.L−1)KIm (Mol.L−1)
F01Value0.2240.1520.0640.2150.0780.1170.069
Standard error0.0030.0030.0030.0050.0010.0040.003
F02Value0.0960.1220.0640.0530.0640.1430.142
Standard error0.0210.0120.0070.0060.0060.0830.036
F03Value0.1640.1370.0640.1310.0720.0990.102
Standard error0.0120.0090.0040.0090.0050.0060.013
F04Value0.0780.0630.0310.0720.0340.1440.092
Standard error0.0120.0040.0020.0040.0020.0470.016
F05Value0.1300.1300.0600.1960.0700.1600.100
Standard error0.0050.0020.0020.0040.0010.0040.004
F06Value0.0940.1290.0630.1590.0630.1630.092
Standard error0.0010.0020.0010.0030.0010.0030.002
F07Value0.0740.0280.0100.0890.0220.1220.139
Standard error0.0020.0010.0010.0040.0020.0050.004
F08Value0.0600.0550.0240.1260.0480.1490.060
Standard error0.0020.0020.0000.0030.0020.0010.002
F09Value0.0710.0970.0480.1140.0710.1080.044
Standard error0.0090.0070.0040.0100.0060.0210.007
F10Value0.2200.2300.1000.3000.1200.1930.091
Standard error0.0510.0130.0070.0220.0070.0440.015
F11Value0.1210.1270.0600.1160.0660.1790.133
Standard error0.0040.0070.0040.0050.0030.0060.003
F12Value0.1320.1550.0820.0760.0860.0430.059
Standard error0.0100.0260.0030.0020.0060.0020.007
F13Value0.1470.1620.0770.2190.0920.1640.106
Standard error0.0060.0060.0020.0080.0040.0090.003
F14Value0.0470.0450.0200.1120.0320.1400.130
Standard error0.0050.0040.0010.0060.0030.0160.007
F15Value0.1600.1800.0800.3300.1100.2600.170
Standard error0.0560.0130.0080.0300.0100.0220.010
F16Value0.1100.1100.0500.2000.0700.3200.280
Standard error0.0030.0040.0020.0070.0030.0510.104
Table 3. Optimal conditions for growth and production of metabolites according to the developed model. In bold: targeted metabolite for each set of operating conditions. Final concentrations and productivities calculated with initial conditions: [X]0 = 80 mmol L−1, [S]0 = 100 mmol L−1, [F, A, L, E]0 = 0.
Table 3. Optimal conditions for growth and production of metabolites according to the developed model. In bold: targeted metabolite for each set of operating conditions. Final concentrations and productivities calculated with initial conditions: [X]0 = 80 mmol L−1, [S]0 = 100 mmol L−1, [F, A, L, E]0 = 0.
A. TargetT (°C)pHFinal Concentrations (mmol L−1)Final Productivities (mmol L−1.h−1)
XFALEATXFALEAT
BiomassBconc.20.07.822712956177733636.563.731.615.112.0910.46
Formic acidFconc.28.09.513317689123943871.682.201.121.541.184.86
Acetic acidAconc.28.09.513317689123943871.682.201.121.541.184.86
Lactic acidLconc.37.06.01809542296724331.682.201.121.541.184.86
EthanolEconc.37.09.5118147771281393522.352.901.542.542.766.98
EthanolEconc.27.07.621714362166683726.934.571.985.312.1611.86
Total acidsATconc.37.06.01809542296724333.331.760.775.481.338.01
Total acidsATconc.37.09.5118147771281393522.352.901.542.542.766.98
BiomassBprod.33.57.519913964178733827.495.252.416.712.7714.38
Formic acidFprod.34.58.017814871163823816.775.612.706.183.1214.48
Acetic acidAprod.35.08.117214872161853816.495.582.716.063.2114.35
Lactic acidLprod.37.07.118812358215783956.604.312.047.542.7313.89
EthanolEprod.37.08.315814473162973785.575.082.565.703.4313.34
EthanolEprod.28.06.017510841278784272.141.320.503.380.955.20
Total acidsATprod.35.57.718413967177813847.165.412.596.903.1314.90
Total acidsATprod.25.09.514617388126873861.712.031.031.481.024.53
maximization, minimization, conc.: final concentration, prod.: batch-average productivity.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Puentes, C.; Girardeau, A.; Passot, S.; Fonseca, F.; Trelea, I.-C. Dynamic Modeling of Carnobacterium maltaromaticum CNCM I-3298 Growth and Metabolite Production and Model-Based Process Optimization. Foods 2021, 10, 1922. https://doi.org/10.3390/foods10081922

AMA Style

Puentes C, Girardeau A, Passot S, Fonseca F, Trelea I-C. Dynamic Modeling of Carnobacterium maltaromaticum CNCM I-3298 Growth and Metabolite Production and Model-Based Process Optimization. Foods. 2021; 10(8):1922. https://doi.org/10.3390/foods10081922

Chicago/Turabian Style

Puentes, Cristian, Amélie Girardeau, Stephanie Passot, Fernanda Fonseca, and Ioan-Cristian Trelea. 2021. "Dynamic Modeling of Carnobacterium maltaromaticum CNCM I-3298 Growth and Metabolite Production and Model-Based Process Optimization" Foods 10, no. 8: 1922. https://doi.org/10.3390/foods10081922

APA Style

Puentes, C., Girardeau, A., Passot, S., Fonseca, F., & Trelea, I. -C. (2021). Dynamic Modeling of Carnobacterium maltaromaticum CNCM I-3298 Growth and Metabolite Production and Model-Based Process Optimization. Foods, 10(8), 1922. https://doi.org/10.3390/foods10081922

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop