Abstract

Kinetic modeling of biomass and lactic acid production by Enterococcus faecalis SLT13 have been developed during batch culture in M17 and Hydrolyzed Cheese Whey (HCW) in 2 L and 20 L bioreactors. The specific growth rate μmax was higher in 20 L bioreactor (1.09 h−1); however, the maximum specific lactic acid production rate qpmax and maximum specific sugar utilization rate qsmax were higher in 2 L bioreactor. Biomass and sugar utilization were affected by lactic acid inhibition in HCW. No effects of substrate inhibition have been observed. Substrate limitation of biomass has been observed on HCW in 20 L bioreactor; the substrate limitation constant for biomass Ksx was 4.229 g/L. Substrate limitation of sugar consumption has been observed on M17 in 2 L bioreactor; the substrate limitation constant for sugar consumption Kss was 2.73 g/L. Compared to experimental data, the model provided good predictions for biomass, sugar consumption, and lactic acid production.

1. Introduction

Actually, the lactic acid and its derivative polylactic acid (PLA) market appears among the most important markets of the chemical industry. This market is in clear progress and it is predictable to attain USD 3.82 billion and 5.16 billion on 2020 for the lactic acid and the PLA, respectively. This growth is due to several reasons: the increasing request for the eco-friendly products, the increasing prices of products derived from the petrochemistry, and the limited reserves of fossil energy [1].

The lactic acid presents several applications, essentially the food industries (preservative and acidulant), textile, and pharmaceutical industries. Moreover, the expansion of PLA market is the main reason for lactic acid market rising. Indeed, the lactic acid is used as monomer for PLA production [2]. For its bioavailability and biodegradability, the PLA is used not only in the food industries for packaging of sensitive foods (i.e., eggs) and tablets in the pharmaceutical industry, but also in electronics and textile, which can be substituted for synthetic polymers derived from petrochemistry.

Two processes are used to produce industrially lactic acid: chemical synthesis and microbial fermentation. However, microbial fermentation presents numerous advantages compared to chemical synthesis; these advantages are essentially environmental; fermentation is a clean process; also the lactic acid obtained by fermentation is optically pure (L (+)-lactic acid) compared to the racemic mixture (50% L-lactic acid and 50% D-lactic acid) obtained chemically [3]. Indeed, PLA with L (+)-lactic acid content higher than 90% is a crystalline polymer and is preferred for commercial uses than the amorphous polymer with high content of the D-isomer [4].

For advantageous commercial process, choosing an appropriate carbon source for lactic acid production is crucial. Indeed, a promising substrate must be of low cost, does not contain contaminants, and allows fast growth of the lactic acid producer microorganism as well as high recovery in lactic acid [5]. Thus, several substrates have been used for the production of lactic acid through microbial fermentation. Among them, cheese whey, a major biowaste of the dairy industry, is a good candidate for this purpose. This byproduct contains 5–8% (w/w) of dry matter, in which lactose is approximately 60–80%, proteins 10–20%, minerals, vitamins, fat, lactic acid, and trace elements [6]. Cheese whey has been widely used for various productions including organic acids, single-cell proteins, enzymes, and ethanol [7].

Many microorganisms are involved in lactic acid production from varied feedstock including Lactic Acid Bacteria (LAB), yeasts, and fungi. LAB species used for lactic acid production include Lactobacillus, Streptococcus, and Enterococcus producing L(+)-lactic acid, while Leuconostoc and Lactobacillus bulgaricus produce D(−)-lactic acid [8]. Genetically modified Saccharomyces cerevisiae was used for the production of lactic acid in continuous fermentation [9]. Rhizopus is the major fungi used for lactic acid production; the species R. oryzae and R. microsporus have been previously used and high amounts of lactic acid were formed [10]. Others bacterial strains, Bacillus and Escherichia coli, were also used for the production of both lactate isomers [11, 12].

Among lactic acid bacteria involved in lactic acid production, Enterococcus faecalis has been described in many studies as lactic acid producer. Enterococcus faecalis RKY1 is able to produce lactic acid from agricultural feedstock, such as wheat, barley, corn, hydrol, soybean curd residues, malt, and from single and mixed sugars [1315]. The same strain was used by Nandsana and Kumar [16] to produce lactic acid from molasses. Enterococcus faecium No. 78 has been used for lactic acid production from liquefied sago starch [17]. Enterococcus faecalis QU11 is able to produce lactic acid from glycerol fermentation [18].

With increasing interest in the industrial application of lactic acid and PLA, different kinetic models have been developed in order to study the growth and lactic acid production by different species of LAB. Indeed, kinetic modeling is an essential step to optimize a fermentation process, as though models help in process control, reducing process costs and time, and improving product quality [19].

Earlier studies have described lactic acid production by Lactobacillus sp. using Luedeking–Piret model. This kind of model relates cell growth and lactic acid production by including two coefficients: growth-associated coefficient (α) and a non-growth-associated coefficient (β). Unstructured kinetic models such as the Gompertz equation have been also used [20]. Further developments have proposed models which included Monod equation to describe the relation between a limiting substrate and biomass growth and then product formation [21, 22]. Monod equation relates the specific growth rate (µ) to the concentration of a single growth-limiting substrate (S) via two parameters, the maximum specific growth rate (µmax) and the substrate affinity constant (KS). Other authors have proposed the inhibitory effects of high initial lactose concentrations, lactose limitation, and lactate inhibition [23, 24], while others developed models showing the inhibition effects of both associated and dissociated forms of lactic acid [25, 26].

This study is conducted to find the most suitable model that describes biomass and lactic acid production by Enterococcus faecalis SLT13 in batch cultures in a medium prepared with hydrolyzed cheese whey and in a synthetic medium (M17). The effects of culture media and scale-up from 2 to 20 L bioreactors on the kinetic parameters were studied.

2. Materials and Methods

2.1. Bacterial Strain and Culture Conditions

The strain Enterococcus faecalis SLT13 used in this work was isolated from Tunisian traditional fermented milk “Leben” [27]. Stock cultures of this strain were stored in M17 broth containing 20% glycerol (Merck) at −80°C.

Enterococcus faecalis cells were grown in M17 broth (Merck), containing glucose as main carbon source or in papaïn Hydrolyzed Cheese Whey (HCW). Cheese whey powder (BHA, Belgium SA) was rehydrated in distilled water (6% w/v). All fermentations were supplemented with (/L): 2 g yeast extract (YE); 0.5 g NH4Cl; 2 g K2HPO4; and 25 mg MnSO4. Papaïn was provided by the Walloon Center of Industrial Biology (CWBI, Belgium) and had an activity of 22 000 IU. The enzyme : substrate ratio was 0.5 :100 (w/w). The whey hydrolysis was carried out at pH 5.0 and 60°C. The enzyme was inactivated by thermal treatment (90°C for 5 min) and cooled in an ice bath. Hydrolysis time was 30 min. After hydrolysis, the pH was adjusted to 7.0.

Batch cultures using M17 or HCW media were carried out in 2 L bioreactor Biostat B (B. Braun Biotech International, Melsungen, Germany). A working volume of 1.5 L was used, which comprised 1.4 L growth medium and 100 mL of inoculum. Inoculum was prepared by adding 100 μL of the preserved strain in 250 mL Erlenmeyer flask containing 100 mL of M17 broth (approx. 1 × 107 CFU/mL) and then incubated at 30°C. The media used were sterilized in the bioreactor at 121°C for 15 min. The set point of pH was 7.0 ± 0.2, temperature 30°C, and stirrer speed 100 rpm. The pH was adjusted with 1 N HCl or 1 N NaOH. The fermentation was run for 10 h in M17 and 24 h in HCW. Samples were withdrawn every 2–4 h for determination of glucose, lactose, and lactic acid concentrations and viable cells count (CFU/mL).

Batch culture in 20 L bioreactor was conducted as follows: 100 μL of the preserved strain was inoculated into 20 mL vial containing 10 mL M17 broth and then incubated at 30°C for 12 h and transferred again in 250 mL Erlenmeyer flask containing 200 mL M17 broth (preculture P1) that was incubated at 30°C for 16 h. The preculture P1 from the exponential growth phase was inoculated into 1 L Erlenmeyer flask containing 800 mL M17 broth (preculture P2). Batch culture was performed in 20 L stirred bioreactor (Biolafitte, France) in HCW. The media used were sterilized in the bioreactor at 121°C for 15 min. The preculture P2 in the exponential growth phase prepared in 1 L Erlenmeyer flask was added to the culture medium. The temperature of the bioreactor was maintained at 30°C and the pH was controlled at 7 ± 0.2 by the addition of NaOH (3N). Samples were withdrawn every 2–4 h for determination of lactose and lactic acid concentrations and viable cells count (CFU/mL). All experiments were conducted in triplicate.

2.2. Biomass Determination

For biomass determination, a calibration curve relating dry cell weight with the cell concentration (estimated by viable cells counting method) was used. For cells counting, 0.1 mL of diluted sample was spread over M17-agar surface. The plates were then incubated for one night at 30°C. The number of colonies was counted and expressed in CFU/mL (Colony-Forming Units/mL).

2.3. Analytical Analysis

Glucose, lactose, and lactic acid concentrations in the fermentation broth were measured by HPLC (Agilent Technologies) using an ion-exclusion column (SupelCo Gel C-610H) operated at 35°C. Components were eluted with 0.1% phosphoric acid at a flow rate of 0.5 mL·min−1. Detection was accomplished by a refractive index detector.

2.4. Kinetic Model
2.4.1. Model Development

The model used in this work is the one developed by Nandasana and Kumar [16] for lactic acid production by Enterococcus faecalis RKY1 and is as follows:

Specific growth rate:

The biomass production rate:

The substrate consumption rate:

The product production rate:

For (3), the usual assumption proposed by Nandasana and Kumar [16] was not used, due to the nature of the process; a minus sign was included. According to the same authors, these assumptions are taken in consideration for the simplification of the modeling step:(i)Kss = Ksp, low lactose concentration affects in the same way lactose uptake and lactic acid production;(ii)Kis = Kip, high lactose concentration inhibits in the same way lactose uptake and lactic acid production;(iii)Kps = Kpp, lactate concentration affects in the same way lactose uptake and lactic acid production.

The method used to solve model differential equations for a given set of kinetic parameters was ODE113s solver of MATLAB [28], in order to obtain X, S, and values. The initial conditions were the experimental values of X, S, and at time zero (t = t0).

2.4.2. Model Integration

In order to estimate kinetic parameters, a search of these parameters values, for which predicted values of X, S, and are closer to the experimental one, needs to be applied. For this process, an optimization method was used. The nonlinear optimization method chosen for this problem was the interior-point approach to constrained minimization using the “fmincon” function in MATLAB [29, 30]. This function is a local search algorithm that starts with an initial guess and finds a constrained minimum.

The objective function (OF) used in this minimization was in the form of

Thus, it was considered in the estimation of parameters that all variables are equally important and given equal weight to the approximation error for each one of them. The procedure to find the parameters was repeated at least 30 times, with random initial guesses, in order to assure a good solution.

3. Results and Discussion

3.1. Effect of Culture Media and Scale-Up on Lactic Acid and Biomass Productivity

This study is conducted to evaluate the effect of scale-up on growth and lactic acid production by Enterococcus faecalis SLT13. The fermentation was first carried out at laboratory scale, in 2 L bioreactor, than at semipilot scale, in 20 L bioreactor. All the studied fermentations were conducted in Stirred Tank Reactor (STR). Moreover, two distinct culture media were analyzed including the synthetic medium M17 and the papaïn Hydrolyzed Cheese Whey (HCW) in order to characterize the growth kinetics of Enterococcus faecalis SLT13. 2 L bioreactor cultures were conducted in M17 and HCW under controlled conditions, while only HCW was used for 20 L bioreactor; the same parameters (pH, temperature, and stirrer) were maintained.

The key kinetic parameters during the culture of E. faecalis in the conditions mentioned below are estimated experimentally and summarized in Table 1.

Whereas the final biomass is higher in synthetic media (M17), lactic acid concentration is more important when E. faecalis was grown on HCW. Lactic acid yield on lactose (Yps) increases significantly in HCW (from 0.314 to 0.98 g/g). Lactic acid yield obtained for the strain E. faecalis RKY1 based on sugar consumption varies from 0.90 to 0.99 for different initial sugar concentrations [2].

The ability of Enterococcus faecalis SLT13 to conduct lactose-to-lactic acid production is enhanced in hydrolyzed whey. This suggests that the carbon/nitrogen ratio in this kind of culture media is more favorable to lactic acid production than biomass growth. This fact is confirmed by the decrease of biomass yield on lactose (Yxs), which suggests that the carbon source available in the medium is rather directed towards lactic acid production. Values of Yxs are in the same range of magnitude than those reported for others species: Enterococcus faecalis CBRDO1 cultivated in LA5 medium supplemented with various concentration of glucose (5–20 g·L−1) presented Yxs from 0.03 to 0.06 g/g [31]. Enterococcus faecalis RKY1 grown on cane-sugar molasses had 0.07–0.19 g/g [2]. Ziadi et al. [21] presented values varying between 0.18 and 0.20 g/g for two strains of Lactococcus lactis SLT6 and SLT10 cultured on hydrolyzed cheese whey.

The volumetric lactic acid productivities are higher in HCW and reach 1.33 g/L·h in 2 L bioreactor. The ability of E. faecalis to produce L(+)- lactic acid from glucose or lactose by homolactic fermentation was earlier demonstrated in the literature. Lactate is generated through pyruvate reduction and two enzymes are involved in this reaction; the cytosolic NAD + -dependent L-(+)-lactate dehydrogenases: ldh-1 and ldh-2 [32].

Enterococcus faecalis is able to produce lactic acid from different substrates; Oh et al. [15] reported productivities higher than 0.8 g/L·h for E. faecalis RKY1 using barley and wheat as nutrient source.

While biomass is enhanced by scale-up, lactic acid production is negatively affected although the scale-up ratio 1 : 10 is respected. This decrease can be explained by many reasons. During industrial fermentation, scale-up causes changes in the geometry and consequently hydrodynamic of the bioreactor, especially for the STR. Two important changes occurred and can affect microbial survival conducting to a decrease of productivities: the mass transfer and the shear stress. The changed geometric conditions affect mixing behavior and consequently mass transfer which lead to a less availability of nutriments for the cells. The increasing shear stress during scale-up affects intrinsic resistance of the lactic-acid-producer microorganism [33].

Indeed, in mechanical agitated bioreactors and independently from microorganism type, the stirrer is the main dispersing tool allowing contact between both phases, biotic and abiotic, of the system. For an optimal process with optimal kinetic parameters, mass transfer between cells and culture media is crucial. Thus, choosing the appropriate stirrer speed and design is essential [34]. Several stirrer designs are available; the most used on microbial fermentation is the flat blade turbine type Rushton (containing 4, 6, or 8 blades) assuring a radial flow. This type of turbine is adapted to the agitation of Newtonian fluids such as culture media. However, shear stress caused by this type of turbine is important, particularly in the case of fragile microorganisms. Knowing that the impeller used in both 2 L and 20 L bioreactors consists of 2 and 3 turbine type Rushton, respectively, the decrease in lactic acid productivity during scale-up can be attributed to the high shear stress caused by this kind of turbine. The effect of shear stress on LAB was earlier demonstrated when studying the cell metabolism of Lactobacillus delbrueckii subsp. bulgaricus; cells cultivated at 72 Pa were affected by shearing forces; however, when cultivated at 36 Pa, metabolism was improved [35].

Thus, a good comprehension of kinetic behavior during bioreactor scale-up of lactic acid production is of crucial importance.

3.2. Estimation of Parameters

The experimental data obtained from batch fermentations of Enterococcus faecalis SLT13 on M17 and Hydrolyzed Cheese Whey in 2 L and 20 L bioreactors were used for developing the kinetic models and estimation of kinetic parameters. Results are presented in Table 2.

Differences in the maximum specific growth rate (μmax) values obtained from the parameters estimation are observed. Higher growth rate is achieved when the strain is cultivated in HCW in 20 L bioreactor (1.09 h−1). Similar values have been reported for Enterococcus faecalis RKY1 grown on cane-sugar molasses 1.6 h−1 [16]. The low growth rate in M17 was found previously for Enterococcus faecalis EF37 grown on the same medium [20]. Enterococcus faecalis CBRDO1 grown on LA5 medium supplemented with various concentration of glucose (5–20 g·L−1) showed μmax values of 0.59–0.64 h−1 [31].

The substrate limitation constants Ksx, Kss and Ksp (Monod constants) for biomass production, sugar utilization and lactic acid production respectively are estimated. Although no effects of substrate limitation are observed for fermentation in HCW in 2 L bioreactor, substrate limitation of biomass is observed when Enterococcus faecalis SLT13 was cultured on HCW in 20 L bioreactor (the substrate limitation constant for growth of biomass Ksx is 4.229 g/L). Substrate limitation of sugar consumption is observed when the strain is cultured on M17 in 2 L bioreactor (the substrate limitation constant for sugar consumption Kss was 2.73 g/L).

In this study, substrate inhibition is not observed; the key constants for substrate inhibition (Kix, Kis and Kip) are ≥100 g/L (Table 2). Meanwhile, initial substrate concentration does not exceed 45 g/L.

The product inhibition key constants (Kpx, Kps and Kpp) are also estimated and listed in Table 2. Although no substrate inhibition is observed, lactate inhibition of biomass growth occurs in HCW since the values of Kpx are 3.77 and 5.00 g/L for 2 L and 20 L bioreactors respectively. The values of Kpx, Kpp and Kps reported for Enterococcus faecalis RKY1, cultured in molasses, were 17.074 and 29.1664 g/L respectively [16]. The inhibition by lactic acid can be caused by both undissociated and dissociated forms. Seen that lactic acid is completely dissociated at pH ≥ 6 and the fermentations occurs at pH 7, product inhibition is mainly due to completely dissociated form.

For sugar consumption, the higher maximum sugar uptake rate (qsmax) is estimated to be 4.99 g/(g·h) for HCW in 2 L bioreactor, in spite of the values are close for the three fermentations. The maximum specific lactic acid production rate (qpmax) is higher for HCW and reach 2.04 g/(g·h) in 2 L bioreactor.

The growth-associated term (α) is very low for all the fermentations. In spite of that lactic acid production is both growth-associated and non-growth-associated and seen the value of α, the production of lactic acid is essentially occurred in the stationary phase. Higher growth-associated parameter (α) was reported by Nandasana and Kumar [16] for Enterococcus faecalis molasses and was 0.26 g/g. Lactococcus lactis NZ133 grown on lactose showed a highest value 0.36 g/g [36].

Evaluation of the cell death coefficient or endogenous decay constant (Kd) of Enterococcus faecalis is conducted for the three fermentations. The estimated decay constants are very low and equal to 0.0001 h−1 on HCW. Nevertheless, decay constant is relatively important for the fermentation in M17 and is equal to 0.013 h−1. Death coefficient is strain dependent and different Kd values were reported in the literature. Lactococcus lactis SLT6 and SLT10 grown in the same media used in this study (HCW) presented values of 0.093–0.26 h−1respectively [21]. Enterococcus faecalis RKY1 grown on cane-sugar molasses had Kd value of 0.00318 h−1 [16]. The lower Kd value for HCW can be ascribed to the protective effect of cheese whey proteins. Indeed, whey proteins are commonly used as protective media during the spray-drying and storage of lactic acid bacteria [37].

3.3. Comparison of Predicted and Experimental Parameters

The model used in this study and developed by Nandasana and Kumar [16] reflects closely the kinetic behavior of Enterococcus faecalis SLT13. The comparison of experimental data and predicted one is shown in Figures 13.

The value of correlation coefficient (R2) is presented in Table 2 and is used to determine if the model fits well the experimental data. A good linearity between model and experimental data is observed, especially for the fermentations occurring in 2 L M17 and 20 L HCW.

For all the fermentations, no lag phase is observed proving the ability of enterococci to adapt to different growth conditions. Exponential growth phase lasts 6 h when the strain was cultured in M17 followed by a short stationary phase. When cultured in HCW, exponential growth phase lasts only 5 h, while the stationary phase is relatively long. Glucose is totally consumed after 6 h of fermentation which confirmed that in M17 medium, the main factor affecting sugar utilization is substrate limitation (the substrate limitation constant for sugar consumption Kss was 2.73 g/L). In HCW, while few quantities of lactose remain at the end of fermentation in 2 L bioreactor, it is almost totally consumed in 20 L bioreactor after about 10 h of fermentation (1.4 g/L) which confirms that the growth of biomass is limited by lactose concentration since Ksx is 4.229 g/L. Furthermore, the initial concentration of lactose is higher in 2 L bioreactor. Indeed, in cheese whey powder, concentration of lactose can vary from batch to batch.

In all fermentations, lactate production began earlier in the exponential growth phase and continues during the stationary phase. Indeed, lactic acid production is growth-associated as well as non-growth-associated. Lactic acid production remains constant after 7 h and 10 h of fermentation in M17 and in HCW in 20 L bioreactor, respectively. However, on HCW in 2 L bioreactor, lactate concentration continues to increase with the availability of the substrate.

4. Conclusions

The overall results obtained in our study show that cheese whey is suitable for lactic acid production; volumetric productivities are higher than the one obtained in M17. The mathematical modeling is useful and allows characterizing growth and lactic acid production of Enterococcus faecalis SLT13 in two different media and two different volume bioreactors. The scale-up affects lactic acid productivity but enhances biomass growth. Substrate inhibition did not impact significantly the maximum lactic acid productivities attainable in the both 2 L and 20 L bioreactors. However, product inhibition seems to play a major role in the two studied media.

Abbreviations

Kd:Death coefficient (h−1)
Kip:Substrate inhibition constant for lactic acid production (g/L)
Kis:Substrate inhibition constant for sugar consumption (g/L)
Kix:Substrate inhibition constant for growth of biomass (g/L)
Kpp:Product inhibition constant for lactic acid production (g/L)
Kps:Product inhibition constant for sugar consumption (g/L)
Kpx:Product inhibition constant for growth of biomass (g/L)
Ksp:Substrate limitation constant for lactic acid production (g/L)
Kss:Substrate limitation constant for sugar consumption (g/L)
Ksx:Substrate limitation constant for growth of biomass (g/L)
P:Lactic acid concentration (g/L)
qpmax:Maximum specific lactic acid production rate (g/(g h))
qsmax:Maximum specific sugar utilization rate (g/(g h))
S:Sugar concentration (g/L)
X:Biomass concentration (g/L)
Yp/x:Lactic acid yield on growth of biomass (g/g)
Yp/s:Lactic acid yield on sugar consumption (g/g)
Yx/s:Biomass yield on sugar consumption (g/g).
Greek Symbols
α:Growth-associated constant in Luedeking–Piret model (g/g)
μmax:Maximum specific growth rate (h−1) of lactic acid
μ:Specific growth rate (h−1).

Data Availability

All the numerical data used to support the findings of this study are available from the corresponding author upon reasonable request.

Additional Points

Kinetic modeling of lactic acid and biomass production by Enterococcus faecalis SLT13 was realized on synthetic media and Hydrolyzed Cheese Whey. Hydrolyzed Cheese Whey enhances lactic acid production. Scale-up from 2 L bioreactor to 20 L bioreactor increases biomass but not lactic acid production. Biomass growth and sugar utilization are affected by lactic acid inhibition in Hydrolyzed Cheese Whey.

Conflicts of Interest

The authors declare no financial or commercial conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the Ministry of Higher Education and Scientific Research (Tunisia) and the Walloon Region (Belgium) during the project: Technological Platform for the Development of Bioindustries.

Supplementary Materials

Graphical abstract representing lactic acid production by Enterococcus faecalis SLT13 during batch culture on Hydrolyzed Cheese Whey in 2 L and 20 L bioreactors. (Supplementary Materials)