Next Article in Journal
Sustainable and Nutritional Recommendations for the Development of Menus by School Food Services in Spain
Next Article in Special Issue
Microwave Spectroscopy Investigation of Carasau Bread Doughs: Effects of Composition up to 8.5 GHz
Previous Article in Journal
Effects of Far-Infrared Radiation Drying on Starch Digestibility and the Content of Bioactive Compounds in Differently Pigmented Rice Varieties
Previous Article in Special Issue
Optimization of Cellulosic Fiber Extraction from Parsley Stalks and Utilization as Filler in Composite Biobased Films
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling Tool for Studying the Influence of Operating Conditions on the Enzymatic Hydrolysis of Milk Proteins

by
Pedro Valencia
1,
Karen Espinoza
1,
Carolina Astudillo-Castro
2 and
Fernando Salazar
2,*
1
Department of Chemical and Environmental Engineering, Universidad Técnica Federico Santa María, Avenida España 1680, Valparaíso 2390123, Chile
2
School of Food Engineering, Pontificia Universidad Católica de Valparaíso, Waddington 716, Valparaíso 2360100, Chile
*
Author to whom correspondence should be addressed.
Foods 2022, 11(24), 4080; https://doi.org/10.3390/foods11244080
Submission received: 15 November 2022 / Revised: 8 December 2022 / Accepted: 12 December 2022 / Published: 16 December 2022
(This article belongs to the Special Issue Modeling of Food Systems and Design of Experiments)

Abstract

:
Systematic modeling of the enzymatic hydrolysis of milk proteins is needed to assist the study and production of partially hydrolyzed milk. The enzymatic hydrolysis of milk proteins was characterized and evaluated as a function of the temperature and protease concentration using Alcalase, Neutrase and Protamex. Modeling was based on the combination of two empirical models formed by a logarithmic and a polynomial equation to correlate the kinetic constants and the operating conditions. The logarithmic equation fitted with high accuracy to the experimental hydrolysis curves with the three proteases (R2 > 0.99). The kinetic constants were correlated with the operating conditions (R2 > 0.97) using polynomial equations. The temperature and protease concentration significantly affected the initial rate of hydrolysis, i.e., the kinetic constant a, while the kinetic constant b was not significantly affected. The values for the kinetic constant a were predicted according to the operating conditions and they were strongly correlated with the experimental data (R2 = 0.95). The model allowed for a high-quality prediction of the hydrolysis curves of milk proteins. This modeling tool can be used in future research to test the correlation between the degree of hydrolysis and the functional properties of milk hydrolysates.

Graphical Abstract

1. Introduction

The use of hydrolyzed milk proteins is a very attractive practice directed to children’s food not only for the hypoallergenicity of this product, but also for the immunomodulatory action of enzyme-generated peptides [1]. The hydrolyzed proteins are obtained by the proteolytic action of proteases, which reduces the molecular weight of proteins hydrolyzing the peptide bonds, thus decreasing the allergenicity caused by epitopes [2,3]. An increase in the degree of hydrolysis and the levels of low molecular weight peptides are generally correlated with the reduced antigenicity of milk proteins [4]. Wroblewska, et al. [5] found that an increase in the degree of hydrolysis (DH) from 15.3% to 15.9% after a two-step process with Alcalase and papain showed a more effective decrease in immunoreactivity despite the nonsignificant increase in the DH. This is a key control parameter during the production of partially hydrolyzed formulas, which could be preferred over extensively hydrolyzed formulas because of their better taste and nutritional value [2]. Partially hydrolyzed formulas are more orally tolerated and less expensive than extensively hydrolyzed formulas [6]. The DH is a useful parameter to characterize the progression of the proteolytic reaction and it is influenced by the operating conditions, traditionally pH, temperature, protease dose and protease-to-substrate ratio [7]. As mentioned above, control of the DH is a major concern during the formulation of hypoallergenic milk. The correlation of DH with hypoallergenicity is a common strategy to determine the adequate DH to be achieved. In addition, the DH achieved depends on the operating conditions and the protease used. The characterization of the hydrolysis process regarding the operating conditions must be evaluated before allergenicity tests. Theoretical modeling by the Michaelis-Menten equation has been proposed by various authors with different protein substrates and proteases [8,9,10,11,12,13,14,15]. The logarithmic equation obtained by Márquez-Moreno and Fernández-Cuadrado [16] has been used by different authors [17,18,19,20,21,22,23,24] due to its simplicity and successful fitting to the experimental data. Empirical modeling based on response surface methodology (RSM) has been applied to protein hydrolysis [25,26,27,28,29]. In a recent publication, a two-level structured model was used as a new methodology to characterize the enzymatic hydrolysis of proteins [30]. The advantages of this methodology are (i) a good description of the hydrolysis process, (ii) an operating definition of kinetic constants, (iii) excellent fitting to experimental data, (iv) the simultaneous exploration of multiple variables (operating conditions) and (v) very good predictability.
The objective of this work was to characterize the enzymatic hydrolysis of milk proteins regarding the operating conditions with different proteases using this recently published modeling methodology. The proposed methodology can be used to correlate the DH with immunoreactivity reduction during the hydrolysis of milk proteins. The information obtained will be useful to produce hypoallergenic milk by designing and setting the operating conditions, thus achieving the proper DH during the process.

2. Materials and Methods

2.1. Materials

Milk powder from Hormel Foods USA was used. The milk contained 33.7% (w/w) protein and 5.3% (w/w) moisture. Commercial proteases Alcalase 2.4 L (serine endoprotease from Bacillus licheniformis), Neutrase (metallo-endoprotease from Bacillus amyloliquefaciens) and Protamex (a protease with endoprotease and exopeptidase activity from Bacillus licheniformis and Bacillus amyloliquefaciens) supplied by Novozymes (Bagsvaerd, Denmark) were used. Analytical grade quality reagents were used in all experiments.

2.2. Hydrolysis Curves

The milk was prepared by mixing the low-heat milk powder with deionized water at a 1:10 proportion, resulting in a 3.1% (w/w) protein solution. The mixture was subjected to magnetic agitation and the homogenized mixture is hereafter delineated milk. The milk was loaded in a glass vessel and thermoregulated by a water bath at a constant temperature. After thermal equilibrium, the pH was fixed at 6.5. The protease previously diluted with deionized water was added to start the hydrolysis reaction. Different protease amounts were added to 40 g of milk. The release of free α-amino groups was followed by the pH-stat technique (G20 Compact Titrator, Mettler-Toledo) to obtain the hydrolysis progress. The pH was constantly monitored and 0.5 N NaOH was added to maintain the pH of 6.5 during the 60 min of reaction. Calculations were performed according to Equation (1):
α -NH   mM = V · N α · V T
In the above equation, V is the added NaOH volume, N is the NaOH concentration, VT is the total reaction volume and α is the average degree of dissociation of the α-NH groups (an average value of 0.2 was used for all conditions). The total α-amino groups in milk proteins were determined by quantification with the o-phthaldialdehyde method (OPA) reported by Nielsen [31] after total hydrolysis in 6 N HCl at 110 °C for 24 h. The DH in each experiment was calculated according to Equation (2):
D H = V · N α · h T · M P
where hT is the content of α-amino groups per protein mass in milk (10.7 meq/gprot) and MP is the mass of protein in the reaction mixture (1.24 g). The release of α-amino groups against time was plotted to build the hydrolysis curves.

2.3. Experimental Design and Statistics

A central composite circumscribed (CCC) design was used to set the number and conditions of each hydrolysis experiment. Each experiment corresponded to a hydrolysis curve with fixed reaction conditions of temperature (T) and protease concentration (E). The T and E values resulting from the design were between 48–62 °C and 19–231 mAU, respectively (Table 1). A total of 10 experiments resulted from the CCC design shown in Table 2. Each experiment was randomly run as one replicate. The central point was replicated in Experiments 5 and 10 (Table 2). Each experiment consisted of a hydrolysis curve where the product concentration (α-NH) was registered against the reaction time. The autotitrator registered one data point per second, resulting in 3600 points of product concentration (α-NH) as a function of time for a 1 h experiment.
The first modeling level corresponds to the description of each hydrolysis curve performed by the estimation of the kinetics constants a and b from Equation (3). The logarithmic equation proposed by Márquez-Moreno and Fernández-Cuadrado [16] in terms of DH can be rewritten in terms of product concentration:
P = 1 b ln ( a b t + 1 )
where P is the concentration of released α-amino groups during milk protein hydrolysis, t is the reaction time and a and b are the kinetic constants of the model. These constants were estimated by the nonlinear least-squares method by fitting Equation (3) to the experimental data. The determination coefficient (R2) and the standard error (se) were used to characterize the fitting quality according to Equation (4):
s e = σ ^ 2 · C j j ,
where σ ^ 2 is the estimator of σ2 calculated from the mean sum of squares and Cjj is the diagonal element of the dispersion matrix (XTX)−1. The second modeling level correlated each of the kinetic constants a and b with the reaction conditions. Equation (5) resulted from considering the variables T (x1) and E (x2) (operating conditions) and the combination of linear, square and interaction components.
y = β 0 + β 1 x 1 + β 2 x 2 + β 11 x 1 2 + β 22 x 2 2 + β 12 x 1 x 2
A multivariable regression was performed to estimate the regression coefficients (βi). A t test was used to calculate the statistical significance of each individual regression coefficient. This information was used to eliminate the nonsignificant regression coefficients and reduce the polynomial model. An analysis of variance (ANOVA) was used to test the significance of the regression (Table 3, Table 4 and Table 5). Additional hydrolysis experiments were performed to test the predictability of the model (Table 6). The prediction range was calculated from the confidence intervals according to Equation (6):
y ^ x 0 ± t α / 2 , n p σ ^ 2 · x 0 T ( X T X ) 1 x 0
where x0 corresponds to the vector of the experimental points considered in the validation experiment and x 0 T ( X T X ) 1 x 0 is the prediction variance used to calculate the prediction error of the kinetic constants a and b.

3. Results and Discussion

The description of the enzymatic hydrolysis of milk proteins was assessed through a two-level mathematical model using a previously published methodology [30]. The first-level model corresponds to the logarithmic equation used to model a wide spectrum of hydrolysis curves with different protein sources and proteases. The effect of E and T on the reaction performance can be observed in Figure 1. The hydrolysis curves generated similar results for Alcalase and Protamex. On the other hand, Neutrase exhibited high efficiency at the initial phase of the hydrolysis curve (high a values) and a sharp decrease in the hydrolysis rate after 5 or 10 min of reaction. Consequently, lower degrees of hydrolysis were obtained with Neutrase compared to Alcalase and Protamex.
These results agree with previous observations concerning Neutrase thermal stability. Neutrase activity has an optimal temperature range of 45 °C to 55 °C and its activity sharply decreases at temperatures over 55 °C [32]. The half-life was estimated to be approximately 10 min at 60 °C and a pH of 7. This important decrease in Neutrase activity explains the observed shape of the hydrolysis curve. The logarithmic equation was fit to determine the kinetic constants a and b in each experiment. These results are shown in Table 2 for Alcalase, Neutrase and Protamex. A high fitting quality was obtained for all the hydrolysis curves (R2 > 0.99). This result agrees with previous fittings of hydrolysis curves [22,30]. As stated previously, the kinetic constant a corresponds to the initial slope of the hydrolysis curve, i.e., d P / d t at t = 0 [30,33]. Thus, it can be predicted that the kinetic constant a increases when T and E increase. It can be observed in Figure 1 that faster hydrolysis rates and larger reaction extensions were obtained at higher temperatures and protease concentrations. The highest values for the kinetic constant a and the lowest values for the kinetic constant b were obtained with Neutrase. This is reflected in a high initial rate and an abrupt decrease after 5 or 10 min of reaction. Despite this high initial rate, Neutrase produced the lowest DHs. On the other hand, Alcalase and Protamex produced lower initial rates but higher DHs than Neutrase. This result can be explained by the higher thermal stability of Alcalase and Protamex. Valencia et al. [15] reported an Alcalase activity loss of 20% after 3 h of reaction during the hydrolysis of 1.7% (w/w) salmon muscle proteins.
According to the obtained results, the kinetic constants a and b can characterize the different catalytic efficiencies that proteases exhibit during milk protein hydrolysis. The second modeling level was based on the correlation between the kinetic constants (a, b) and the reaction conditions T and E. The polynomial equation contained T and E as the independent variables. The multivariable regression results are shown in Table 3, Table 4 and Table 5 for Alcalase, Neutrase and Protamex, respectively. The nonsignificant coefficients (βi) were eliminated from the polynomial equation to obtain a reduced model for each protease (Table 3, Table 4 and Table 5). The reaction conditions significantly affected the kinetic constant a for all proteases. A lower correlation between the kinetic constant a and operating conditions was obtained with Neutrase (R2 = 0.8746) compared to Alcalase and Protamex. This can be explained by the lower thermal stability of Neutrase. A better correlation can be observed at lower temperatures for Neutrase (40–50 °C). The kinetic constant b was poorly correlated with the operating conditions when Alcalase and Protamex were used (R2 = 0.7972 and R2 = 0.8693, respectively). For this reason, the kinetic constant b was fixed at 0.0410 and 0.0423 for Alcalase and Protamex, respectively, corresponding to the central value. The reaction conditions impacted the kinetic constant b to a lesser extent than the kinetic constant a. Average values of b for different operating conditions have been used in previous publications [16,19,21,22,24]. This result is consistent with previous findings, where the kinetic constant b depended exclusively on the substrate concentration [30]. In the present study, the protein concentration was constant because the same milk formulation was used in all experiments. This explains the poor correlation of the kinetic constant b with the T and E. Regression analysis and the analysis of variance, shown in Table 3, Table 4 and Table 5. The response surfaces from the reduced models for the kinetic constant a were plotted against T and E in Figure 2.
The ANOVA results shown in Table 3, Table 4 and Table 5 indicated that these models were significant in describing the kinetic constant a as a function of operating conditions, with the exception of the kinetic constant a for Neutrase (p > 0.01). This result was explained by the thermal inactivation of Neutrase. A better fitting can be achieved at lower temperatures for Neutrase. Considering the results obtained with Alcalase and Protamex, the polynomial models can be used to calculate the values of the kinetic constant a from a set of operating conditions (T, E) and plot the hydrolysis curves considering the fixed values for the kinetic constant b. The predictability of these models was tested with additional experiments, which are presented in Table 6. The models for kinetic constant a were used to calculate their predicted values, while the experimental values were calculated from the fitting of the hydrolysis curves from validation experiments. The predicted vales of the kinetic constant a were correlated against their experimental values in Figure 3. Good agreement was obtained for both Alcalase (R2 = 0.9830) and Protamex (R2 = 0.9295) independently and with both proteases together (R2 = 0.9528). For Neutrase, high prediction errors were obtained, thus, this protease was excluded from validation experiments. This error is explained by significant thermal inactivation of Neutrase. The predicted hydrolysis curves are plotted in Figure 4 and Figure 5 for Alcalase and Protamex, respectively.
The predicted hydrolysis curves were in good agreement with the experimental curves, despite the errors observed in Table 6 for the kinetic constant a, reaching approximately 20% in some cases. Indeed, the predicted DH after 60 min of reaction resulted in less than 5% error for most of the cases (Table 7). This result is an important source of validation for the present model because it can not only be used to predict the DH, but also to study the reactor performance and design to achieve a certain DH value. As noted in a previous article [30], the prediction quality of this model can still be improved by using a different experimental design and/or expanding the studied range of the operating conditions. Furthermore, this methodology allowed us to compare the performance of three different commercial proteases. It was observed that in the range of temperatures studied, Neutrase was rapidly inactivated, achieving a lower DH after 60 min of reaction compared to Alcalase and Protamex. Protamex showed an initial rate twice that observed for Alcalase. This performance is reflected in the values of the kinetic constant a (Table 2). Nevertheless, Alcalase reached almost the same values of DH as Protamex after 60 min of reaction. This performance is reflected in the lower values of the kinetic constant b obtained with Alcalase (Table 2). The performance of Protamex can be explained by a possible higher thermal inactivation and/or product inhibition than Alcalase. In resolution, this methodology can be used to evaluate the performance of any protease under different operating conditions. Different milk sources or formats can be evaluated using this methodology. Considerations about the nitrogen content must be prevented for the proper calculation of the DH. In the case that autotitration is not available, a different analytical technique can be used for the quantification of released α-NH, such as the TNBS [34] or OPA [31] methods. Sample withdrawal is required in these methods in order for the sample to be mixed with the proper reagent. If parameter conditions other than protease dose needs to be tested, such as milk concentration, pH or E/S, a new factor can be added to the polynomial equation after including this parameter in the experimental design. Moreover, different experimental designs can be tried if better statistical properties are achieved in the experimental configuration. Thus, the actual methodology offers not only reliability, but also versatility.
The logarithmic equation has largely been demonstrated to be a reliable model because of its good description, explanation, and prediction. The proposed modeling tool will improve the study of the correlation between DH and allergenicity during the enzymatic hydrolysis of milk proteins. An intermediate correlation between DH and allergenicity is further needed. In this way, DH can be used as an objective and quality control parameter during the hydrolysis process, thus allowing the reaction to be stopped at the desired allergenicity. This improvement will allow us to establish the operating conditions to produce a specific DH in the milk protein hydrolysate and, consequently, a specific functional property, for example, a specific allergenicity reduction.

4. Conclusions

The performance of the enzymatic hydrolysis of milk proteins was characterized and mathematically modeled against the temperature and protease concentration. The enzymatic hydrolysis performance was readily and accurately predicted using this modeling methodology to compare the protease efficiency and temperature effect. The temperature and protease concentration significantly affected the initial rate of hydrolysis, while the reaction extension was not affected. The model allows the possibility of studying the effect of additional variables, such as operating conditions, by simply adding variables to the polynomial equation and using an adequate experimental design. This methodology can be used in future research to test the effect of hydrolysis conditions on the reduction in antigenicity of hydrolysates. In this way, antigenicity can be correlated with operating conditions and the process can be designed to obtain hypoallergenic milk products.

Author Contributions

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

Funding

This research was funded by ANID FONDECYT/Iniciación 11110249 and 11110402, by Euro-Chile, French Embassy in Chile and Institut Français Project 2012–2013 and by INNOVA-CORFO I + D 11IDL1-10462.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data is contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gauthier, S.F.; Pouliot, Y.; Saint-Sauveur, D. Immunomodulatory peptides obtained by the enzymatic hydrolysis of whey proteins. Int. Dairy J. 2006, 16, 1315–1323. [Google Scholar] [CrossRef]
  2. Bu, G.H.; Luo, Y.K.; Chen, F.S.; Liu, K.L.; Zhu, T.W. Milk processing as a tool to reduce cow’s milk allergenicity: A mini-review. Dairy Sci. Technol. 2013, 93, 211–223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Dąbrowska, A.; Bajzert, J.; Babij, K.; Szołtysik, M.; Stefaniak, T.; Willak-Janc, E.; Chrzanowska, J. Reduced IgE and IgG antigenic response to milk proteins hydrolysates obtained with the use of noncommercial serine protease from Yarrowia lipolytica. Food Chem. 2020, 302, 125350. [Google Scholar] [CrossRef] [PubMed]
  4. Liang, X.; Yang, H.; Sun, J.; Cheng, J.; Luo, X.; Wang, Z.; Yang, M.; Bing Tao, D.; Yue, X.; Zheng, Y. Effects of enzymatic treatments on the hydrolysis and antigenicity reduction of natural cow milk. Food Sci. Nutr. 2021, 9, 985–993. [Google Scholar] [CrossRef] [PubMed]
  5. Wroblewska, B.; Jedrychowski, L.; Szabo, E.; Hajos, G. The reduction of cow milk proteins immunoreactivity by two-step enzymatic hydrolysis. Acta Aliment. 2005, 34, 307–315. [Google Scholar] [CrossRef]
  6. Abd El-Salam, M.H.; El-Shibiny, S. Preparation, properties, and uses of enzymatic milk protein hydrolysates. Crit. Rev. Food Sci. Nutr. 2017, 57, 1119–1132. [Google Scholar] [CrossRef]
  7. Adler-Nissen, J. Enzymic Hydrolysis of Food Proteins; Elsevier Applied Science Publishers: Barking, UK, 1986. [Google Scholar]
  8. Apar, D.K.; Özbek, B. A Kinetic Study on Corn Gluten Hydrolysis. Chem. Eng. Technol. 2009, 32, 673–675. [Google Scholar] [CrossRef]
  9. Barros, R.M.; Xavier Malcata, F. A kinetic model for hydrolysis of whey proteins by cardosin A extracted from Cynara cardunculus. Food Chem. 2004, 88, 351–359. [Google Scholar] [CrossRef]
  10. Demirhan, E.; Apar, D.K.; Ozbek, B. A kinetic study on sesame cake protein hydrolysis by Alcalase. J. Food Sci. 2011, 76, C64–C67. [Google Scholar] [CrossRef]
  11. Sousa, R.; Lopes, G.P.; Tardioli, P.W.; Giordano, R.L.C.; Almeida, P.I.F.; Giordano, R.C. Kinetic model for whey protein hydrolysis by alcalase multipoint-immobilized on agarose gel particles. Braz. J. Chem. Eng. 2004, 21, 147–153. [Google Scholar] [CrossRef]
  12. Tardioli, P.W.; Sousa, R.; Giordano, R.C.; Giordano, R.L.C. Kinetic model of the hydrolysis of polypeptides catalyzed by Alcalase® immobilized on 10% glyoxyl-agarose. Enzym. Microb. Technol. 2005, 36, 555–564. [Google Scholar] [CrossRef]
  13. Trusek-Holownia, A. Production of protein hydrolysates in an enzymatic membrane reactor. Biochem. Eng. J. 2008, 39, 221–229. [Google Scholar] [CrossRef]
  14. O’Meara, G.M.; Munro, P.A. Kinetics of the hydrolysis of lean meat protein by alcalase: Derivation of two alternative rate equations and their fit to experimental data. Biotechnol. Bioeng. 1985, 27, 861–869. [Google Scholar] [CrossRef] [PubMed]
  15. Valencia, P.; Pinto, M.; Almonacid, S. Identification of the key mechanisms involved in the hydrolysis of fish protein by Alcalase. Process Biochem. 2014, 49, 258–264. [Google Scholar] [CrossRef]
  16. Márquez-Moreno, M.C.; Fernández-Cuadrado, V. Enzymic hydrolysis of vegetable proteins: Mechanism and kinetics. Process Biochem. 1993, 28, 481–490. [Google Scholar] [CrossRef]
  17. Apar, D.K.; Ozbek, B. Hydrolysis and solubilization of corn gluten by Neutrase. J. Chem. Technol. Biotechnol. 2007, 82, 1107–1114. [Google Scholar] [CrossRef]
  18. Apar, D.K.; Ozbek, B. Corn gluten hydrolysis by Alcalase: Effects of process parameters on hydrolysis solubilization and enzyme inactivation. Chem. Biochem. Eng. Q. 2008, 22, 203–212. [Google Scholar]
  19. Demirhan, E.; Apar, D.K.; Özbek, B. Sesame cake protein hydrolysis by alcalase: Effects of process parameters on hydrolysis, solubilization, and enzyme inactivation. Korean J. Chem. Eng. 2011, 28, 195–202. [Google Scholar] [CrossRef]
  20. González-Tello, P.; Camacho, F.; Jurado, E.; Páez, M.P.; Guadix, E.M. Enzymatic hydrolysis of whey proteins: I. Kinetic models. Biotechnol. Bioeng. 1994, 44, 523–528. [Google Scholar] [CrossRef]
  21. Márquez, M.C.; Vázquez, M.A. Modeling of enzymatic protein hydrolysis. Process Biochem. 1999, 35, 111–117. [Google Scholar] [CrossRef]
  22. Qi, W.; He, Z. Enzymatic hydrolysis of protein: Mechanism and kinetic model. Front. Chem. China 2006, 1, 308–314. [Google Scholar] [CrossRef]
  23. Martinez-Araiza, G.; Castano-Tostado, E.; Amaya-Llano, S.L.; Regalado-Gonzalez, C.; Martinez-Vera, C.; Ozimek, L. Modeling of Enzymatic Hydrolysis of Whey Proteins. Food Bioprocess Technol. 2012, 5, 2596–2601. [Google Scholar] [CrossRef]
  24. Ruan, C.Q.; Chi, Y.J.; Zhang, R.D. Kinetics of Hydrolysis of Egg White Protein by Pepsin. Czech J. Food Sci. 2010, 28, 355–363. [Google Scholar] [CrossRef] [Green Version]
  25. Bhaskar, N.; Benila, T.; Radha, C.; Lalitha, R.G. Optimization of enzymatic hydrolysis of visceral waste proteins of Catla (Catla catla) for preparing protein hydrolysate using a commercial protease. Bioresour. Technol. 2008, 99, 335–343. [Google Scholar] [CrossRef] [PubMed]
  26. Bhaskar, N.; Mahendrakar, N.S. Protein hydrolysate from visceral waste proteins of Catla (Catla catla): Optimization of hydrolysis conditions for a commercial neutral protease. Bioresour. Technol. 2008, 99, 4105–4111. [Google Scholar] [CrossRef]
  27. Chabeaud, A.; Dutournie, P.; Guerard, F.; Vandanjon, L.; Bourseau, P. Application of response surface methodology to optimize the antioxidant activity of a saithe (Pollachius virens) hydrolysate. Mar. Biotechnol. 2009, 11, 445–455. [Google Scholar] [CrossRef]
  28. Liaset, B.; Nortvedt, R.; Lied, E.; Espe, M. Studies on the nitrogen recovery in enzymic hydrolysis of Atlantic salmon (Salmo salar, L.) frames by Protamex™ protease. Process Biochem. 2002, 37, 1263–1269. [Google Scholar] [CrossRef]
  29. Ovissipour, M.; Kenari, A.A.; Motamedzadegan, A.; Nazari, R.M. Optimization of Enzymatic Hydrolysis of Visceral Waste Proteins of Yellowfin Tuna (Thunnus albacares). Food Bioprocess Technol. 2012, 5, 696–705. [Google Scholar] [CrossRef]
  30. Valencia, P.; Espinoza, K.; Ceballos, A.; Pinto, M.; Almonacid, S. Novel modeling methodology for the characterization of enzymatic hydrolysis of proteins. Process Biochem. 2015, 50, 589–597. [Google Scholar] [CrossRef]
  31. Nielsen, P.M.; Petersen, D.; Dambmann, C. Improved Method for Determining Food Protein Degree of Hydrolysis. J. Food Sci. 2001, 66, 642–646. [Google Scholar] [CrossRef]
  32. Benjakul, S.; Morrissey, M.T. Protein Hydrolysates from Pacific Whiting Solid Wastes. J. Agric. Food Chem. 1997, 45, 3423–3430. [Google Scholar] [CrossRef]
  33. Valencia, P.L.; Solis, T.; Rojas, P.; Ibañez, F.; Astudillo-Castro, C.; Pinto, M.; Almonacid, S. Proteolytic susceptibility of food byproduct proteins: An evaluation by means of a quantitative index. Process Biochem. 2019, 77, 63–69. [Google Scholar] [CrossRef]
  34. Adler-Nissen, J. Determination of the degree of hydrolysis of food protein hydrolysates by trinitrobenzenesulfonic acid. J. Agric. Food Chem. 1979, 27, 1256–1262. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Experimental curves of milk protein hydrolysis with (a) Alcalase, (b) Neutrase and (c) Protamex at different temperatures (125 mAU) and protease concentrations (55 °C).
Figure 1. Experimental curves of milk protein hydrolysis with (a) Alcalase, (b) Neutrase and (c) Protamex at different temperatures (125 mAU) and protease concentrations (55 °C).
Foods 11 04080 g001
Figure 2. Response surfaces for kinetic constant a against the operating conditions T and E for the hydrolysis of milk proteins with (a) Alcalase and (b) Protamex.
Figure 2. Response surfaces for kinetic constant a against the operating conditions T and E for the hydrolysis of milk proteins with (a) Alcalase and (b) Protamex.
Foods 11 04080 g002
Figure 3. Correlation between predicted and experimental values of the kinetic constant a for the hydrolysis curves obtained in validation experiments with Alcalase (R2 = 0.9830) and Protamex (R2 = 0.9295). The correlation for both data is R2 = 0.9528.
Figure 3. Correlation between predicted and experimental values of the kinetic constant a for the hydrolysis curves obtained in validation experiments with Alcalase (R2 = 0.9830) and Protamex (R2 = 0.9295). The correlation for both data is R2 = 0.9528.
Foods 11 04080 g003
Figure 4. Predicted and experimental hydrolysis curves for the hydrolysis of milk proteins with Alcalase. Experimental data (circles), predicted curves (continuous line) and confidence interval of predicted curves at 95% probability (dashed lines) for experiments presented in Table 6.
Figure 4. Predicted and experimental hydrolysis curves for the hydrolysis of milk proteins with Alcalase. Experimental data (circles), predicted curves (continuous line) and confidence interval of predicted curves at 95% probability (dashed lines) for experiments presented in Table 6.
Foods 11 04080 g004
Figure 5. Predicted and experimental hydrolysis curves for the hydrolysis of milk proteins with Protamex. Experimental data (circles), predicted curve (continuous line) and confidence interval of predicted curves at 95% probability (dashed lines) for experiments presented in Table 6.
Figure 5. Predicted and experimental hydrolysis curves for the hydrolysis of milk proteins with Protamex. Experimental data (circles), predicted curve (continuous line) and confidence interval of predicted curves at 95% probability (dashed lines) for experiments presented in Table 6.
Foods 11 04080 g005
Table 1. Codified values for operating conditions of temperature and protease concentration.
Table 1. Codified values for operating conditions of temperature and protease concentration.
VariableLevels
−1.4−1.00+1.0+1.4
E (mUA)1950125200231
T (°C)4850556062
Table 2. Experimental CCC design and kinetic constant values obtained from fitting to hydrolysis curves and %DH after 60 min of reaction.
Table 2. Experimental CCC design and kinetic constant values obtained from fitting to hydrolysis curves and %DH after 60 min of reaction.
ExpVariablesAlcalaseNeutraseProtamex
x1 (S)x2 (E)x3 (T)a ± se
(mM/min)
b ± se
(mM−1)
R2a ± se
(mM/min)
b ± se
(mM−1)
R2a ± se
(mM/min)
b ± se
(mM−1)
R2
1−1.40.0−1.44.06 ± 0.0150.0494 ± 1.36 × 10−40.994025.6 ± 0.1370.135 ± 2.02 × 10−40.99679.15 ± 0.0300.0631 ± 9.98 × 10−50.9966
2−1.01.0−1.07.56 ± 0.0270.0432 ± 9.01 × 10−50.995530.2 ± 0.1960.104 ± 1.93 × 10−40.992316.5± 0.0420.0509 ± 5.38 × 10−50.9983
3−1.0−1.0−1.02.89 ± 0.0090.0498 ± 1.34 × 10−40.99558.76 ± 0.0110.105 ± 5.28 × 10−50.99964.36 ± 0.0100.0542 ± 8.37 × 10−50.9980
40.01.40.015.2 ± 0.0380.0356 ± 4.39 × 10−50.998042.7 ± 0.1220.080 ± 6.38 × 10−50.998526.6 ± 0.0500.0409 ± 2.99 × 10−50.9991
50.00.00.09.46 ± 0.0420.0359 ± 9.36 × 10−50.992826.3 ± 0.0480.082 ± 4.68 × 10−50.999316.1 ± 0.0300.0407 ± 3.44 × 10−50.9990
60.0−1.40.02.15 ± 0.0030.0508 ± 8.14 × 10−50.99874.96 ± 0.0070.082 ± 6.66 × 10−50.99932.93 ± 0.0060.0498 ± 8.65 × 10−50.9981
71.01.01.021.7 ± 0.0420.0361 ± 2.96 × 10−50.999086.5 ± 0.5930.078 ± 1.28 × 10−40.993340.7 ± 0.1260.0386 ± 4.12 × 10−50.9979
81.0−1.01.05.76 ± 0.0150.0353 ± 6.53 × 10−50.997422.0 ± 0.1240.085 ± 1.58 × 10−40.99349.44 ± 0.0140.0392 ± 3.19 × 10−50.9993
91.40.01.417.8 ± 0.0450.0343 ± 4.01 × 10−50.9981105.8 ± 1.080.090 ± 2.04 × 10−40.986630.2 ± 0.1130.0411 ± 5.69 × 10−50.9967
100.00.00.010.4 ± 0.0350.0401 ± 7.31 × 10−50.996237.7 ± 0.0910.084 ± 5.79 × 10−50.998916.5 ± 0.0310.0399 ± 3.38 × 10−50.9990
Table 3. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Alcalase.
Table 3. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Alcalase.
Regression analysis for kinetic constant a = 9.69 + 4.56 x 1 + 4.88 x 2 + 2.82 x 1 x 2
PredictorCoefficientsetcalcP
β09.6940.35027.681.47 × 10−7
β14.5630.39211.652.41 × 10−5
β24.8800.39212.461.63 × 10−5
β122.8200.5545.0930.00223
R2 = 0.9814R2 aj = 0.9721
Analysis of variance for kinetic constant a
Source of
variation
Degrees of
freedom
Sum of
squares
Mean squareFP
Regression4388.861697.215466.040.000162
Error57.35991.47198
Total9396.2215
Regression analysis for kinetic constant b = 0.0410 0.00537 x 1 0.00344 x 2
PredictorCoefficientsetcalcP
β00.041030.0010937.752.38 × 10−9
β1−0.005370.001224.420.00310
β2−0.003440.001222.830.0253
R2 = 0.7972R2 aj = 0.7393
Analysis of variance kinetic constant b
Source of
variation
Degrees of
freedom
Sum of
squares
Mean squareFP
Regression30.00032510.0001087.8640.01679
Error60.00008270.000014
Total90.0004078
Table 4. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Neutrase.
Table 4. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Neutrase.
Regression analysis for kinetic constant a = 25.07 + 22.86 x 1 + 17.41 x 2 + 17.49 x 1 2
PredictorCoefficientsetcalcP
β025.076.533.840.00855
β122.864.984.590.00375
β217.414.983.490.01291
β1117.495.962.940.02610
R2 = 0.8746R2 aj = 0.8119
Analysis of variance for kinetic constant a
Source of
variation
Degrees of freedomSum of
squares
Mean squareFP
Regression48322.22080.558.720.0177
Error51193.2238.63
Total99515.4
Regression analysis for kinetic constant b = 0.0809 0.0135 x 1 + 0.0147 x 1 2
PredictorCoefficientsetcalcP
β00.080910.0018842.989.64 × 10−10
β1−0.013520.001449.403.21 × 10−5
β110.014700.001728.565.92 × 10−5
R2 = 0.9585R2 aj = 0.9466
Analysis of variance kinetic constant b
Source of
variation
Degrees of freedomSum of
squares
Mean squareFP
Regression30.0026730.00089146.170.000154
Error60.0001160.000019
Total90.002788
Table 5. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Protamex.
Table 5. Regression and variance analysis for RSM models for kinetic constants a and b for the enzymatic hydrolysis of milk proteins with Protamex.
Regression analysis for kinetic constant a = 17.25 + 7.40 x 1 + 9.60 x 2 + 4.79 x 1 x 2
PredictorCoefficientsetcalcP
β017.250.81821.097.41 × 10−7
β17.400.9148.090.000191
β29.600.91410.504.38 × 10−5
β124.791.2933.710.0100
R2 = 0.9693R2 aj = 0.9540
Analysis of variance kinetic constant a
Source of
variation
Degrees of freedomSum of
squares
Mean squareFP
Regression41267.3316.8339.470.000565
Error540.18.03
Total91307.4
Regression analysis for kinetic constant b = 0.0423 0.00731 x 1 + 0.00436 x 1 2
PredictorCoefficientsetcalcP
β00.04230.0015727.012.44 × 10−8
β1−0.007310.001206.100.000489
β110.004360.001433.050.0186
R2 = 0.8693R2 aj = 0.8320
Analysis of variance kinetic constant b
Source of
variation
Degrees of freedomSum of
squares
Mean squareFP
Regression30.0005340.00017813.300.004637
Error60.0000800.000013
Total90.000614
Table 6. Predicted and experimental values for kinetic constants a and b obtained from validation experiments with Alcalase and Protamex.
Table 6. Predicted and experimental values for kinetic constants a and b obtained from validation experiments with Alcalase and Protamex.
ExpVariablesAlcalaseProtamex
x1 (T)x2 (E)a ± sea ± se
PredictedExperimentalError (%)PredictedExperimentalError (%)
11−0.507.41 ± 0.449.10 ± 0.0218.613.6 ± 2.5117.0 ± 0.0120.4
120.5012.0 ± 0.4410.2 ± 0.0116.821.0 ± 2.5116.9 ± 0.0324.2
130−14.81 ± 0.583.92 ± 0.0122.87.65 ± 3.297.44 ± 0.012.9
140114.6 ± 0.5813.8 ± 0.025.726.8 ± 3.2925.9 ± 0.033.8
Table 7. Predicted and experimental values for DH after 60 min of milk protein hydrolysis from validation experiments with Alcalase and Protamex.
Table 7. Predicted and experimental values for DH after 60 min of milk protein hydrolysis from validation experiments with Alcalase and Protamex.
ExpVariablesDH
x1
(T)
x2
(E)
AlcalaseProtamex
PredictedExperimentalError (%)PredictedExperimentalError (%)
11−0.5021.825.715.325.529.614.0
120.5025.224.81.528.527.53.8
130−118.819.63.921.623.47.9
140126.627.22.230.330.10.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Valencia, P.; Espinoza, K.; Astudillo-Castro, C.; Salazar, F. Modeling Tool for Studying the Influence of Operating Conditions on the Enzymatic Hydrolysis of Milk Proteins. Foods 2022, 11, 4080. https://doi.org/10.3390/foods11244080

AMA Style

Valencia P, Espinoza K, Astudillo-Castro C, Salazar F. Modeling Tool for Studying the Influence of Operating Conditions on the Enzymatic Hydrolysis of Milk Proteins. Foods. 2022; 11(24):4080. https://doi.org/10.3390/foods11244080

Chicago/Turabian Style

Valencia, Pedro, Karen Espinoza, Carolina Astudillo-Castro, and Fernando Salazar. 2022. "Modeling Tool for Studying the Influence of Operating Conditions on the Enzymatic Hydrolysis of Milk Proteins" Foods 11, no. 24: 4080. https://doi.org/10.3390/foods11244080

APA Style

Valencia, P., Espinoza, K., Astudillo-Castro, C., & Salazar, F. (2022). Modeling Tool for Studying the Influence of Operating Conditions on the Enzymatic Hydrolysis of Milk Proteins. Foods, 11(24), 4080. https://doi.org/10.3390/foods11244080

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