Mathematical Problems in Engineering

Volume 2014 (2014), Article ID 392054, 7 pages

http://dx.doi.org/10.1155/2014/392054

Research Article

Comparison of Two Mechanistic Microbial Growth Models to Estimate Shelf Life of Perishable Food Package under Dynamic Temperature Conditions

Dong Sun Lee

Department of Food Science and Biotechnology, Kyungnam University, 7 Kyungnamdaehak-ro, Masanhappo-gu, Changwon 631-701, Republic of Korea

Received 14 July 2014; Accepted 17 September 2014; Published 1 October 2014

Academic Editor: Kit Keith L. Yam

Copyright © 2014 Dong Sun Lee. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Two mechanistic microbial growth models (Huang’s model and model of Baranyi and Roberts) given in differential and integrated equation forms were compared in predicting the microbial growth and shelf life under dynamic temperature storage and distribution conditions. Literatures consistently reporting the microbial growth data under constant and changing temperature conditions were selected to obtain the primary model parameters, set up the secondary models, and apply them to predict the microbial growth and shelf life under fluctuating temperatures. When evaluated by general estimation behavior, bias factor, accuracy factor, and root-mean-square error, Huang’s model was comparable to Baranyi and Roberts’ model in the capability to estimate microbial growth under dynamic temperature conditions. Its simple form of single differential equation incorporating directly the growth rate and lag time may work as an advantage to be used in online shelf life estimation by using the electronic device.

1. Introduction

Microbial growth or spoilage is often the primary quality factor determining the shelf life of many perishable foods [1–3]. Sometimes it also causes health hazard with pathogenic organisms. Microbial growth model is very useful to decide the shelf life of the packaged perishable foods for any condition of food storage and distribution [2]. Even though the shelf life determination is commonly conducted based on the most probable storage or distribution temperature, temperature abuse or fluctuation occurring in the food supply chain makes it difficult to apply a single shelf life to practical situation [4, 5]. Dynamic or real-time estimation of the microbial food quality is suggested as a way to ensure the quality level delivered to the consumers and can be realized by devices such as time temperature integrator (TTI) and temperature-monitoring sensor attached on the food packages [6–9]. Recent advances in sensor and communication engineering can handle dynamic storage data for effective shelf life determination of the packaged foods [10–12].

Any dynamic shelf life controls can be possible with appropriate prediction model to estimate the microbial growth under temperature-variant conditions. Several attempts or methods have been reported for this purpose [13–15]. Among microbial growth models proposed so far by several researchers, models given in differential equation(s) forms  are easier and more powerful to deal with fluctuating temperature conditions [16, 17]. Model of Baranyi and Roberts [18] seems to be the most widely used one for this purpose after its introduction. Recently Huang [19] developed a new mechanistic model to describe the microbial growth in a differential equation form. He also proposed a new functional relationship to express the temperature dependence of growth rate instead of square root model. However, his model has not been tried for predicting microbial growth under dynamic temperature conditions and thus needs to be examined or tested for its capability.

This study therefore aims to apply Huang’s mechanistic microbial growth model to dynamic temperature conditions and compare it with model of Baranyi and Roberts in the ability to predict the microbial growth of the packaged perishable food in the food supply chain.

2. Materials and Methods

2.1. Huang’s Model versus Model of Baranyi and Roberts

Huang’s model to describe increase rate of cell number against time is written in a differential equation given by [19]

where , , and are maximum specific growth rate, lag time, and maximum cell density, respectively.

Equation (1) can be written in an integrated format of primary model for a constant temperature condition:

where is initial microbial load and .

Huang’s model needs four parameters of , , , and to characterize a microbial growth curve under constant temperature storage condition.

Model of Baranyi and Roberts [18] is written in two differential equations:

where is a hypothetical physiological state of the cell population representing normalized concentration of unknown substance critically needed for cell growth.

Equation (4) in combination with (3) can be formulated into an integrated format of primary model for a constant environmental condition:

where is defined by and is initial state of at time zero.

Model of Baranyi and Roberts (will be named as Baranyi model later for simplicity) also needs four parameters of , , , and to characterize a microbial growth curve at constant temperature condition. In their model and can be related to obtain lag time as

While square root model (see (7) and (8)) as secondary model has mostly been used to describe the temperature dependence of growth rate and lag time in many predictive models including that of Baranyi and Roberts [13, 20], Huang [19] proposed using a different form of (9) for growth rate. Consider the following:

where and are constant and minimum growth temperature for each equation, respectively.

2.2. Literature Data for Testing

Usually in dynamic shelf life modeling, it is assumed that kinetic parameters are obtained from growth experiments at constant temperatures and their temperature dependence is established to be applied to temperature-variant conditions. All those functional relationships should be verified in temperature-changing storage and distribution of food packages. Thus literatures reporting microbial growth data sufficiently under constant and fluctuating temperature conditions in consistent way were looked for to test and compare the two models of microbial growth. Growth data of pseudomonads on seabream from Koutsoumanis [21], those of Escherichia coli O157:H7 and Salmonella spp. on lettuce from Koseki and Isobe [22], and of lactic acid bacterial growth data on cured meat from Mataragas et al. [16] were collected and used in this study. Even though pathogenic bacteria such as E. coli and Salmonella are not usually used for shelf life determination, the data from Koseki and Isobe [22] were included for the analysis in this study due to their consistency and relevance in relation to product supply chain to the consumers. Primary model parameters under constant temperature conditions were obtained by nonlinear regression of the data in (2) and (5) except that those given by Koseki and Isobe [22] were directly adopted for the Baranyi model. Secondary model parameters for (7), (8), and (9) were also determined by linear or nonlinear regression.

3. Results and Discussion

3.1. Adequacy Evaluation of the Simulation Models Using Literature Data

Table 1 shows parameters of primary model, for both Huang and Baranyi models, which were determined or taken for the literature sources. From the primary model parameters, secondary models to express the functional relationship of temperature effect on microbial growth have been presented in the literature or can be formulated following (7), (8), and/or (9). Summarized results are given in Table 2, and those relationships could be connected to the solution of differential equation (1) for Huang model and to that of (3) and (4) for Baranyi model. Maximum cell density is assumed to be indifferent to temperature or be in a simple linear function of temperature depending on the microbial strains (Table 2). For the Baranyi model using (3), at initial state is required, and average values from constant temperature conditions have been used as in many other studies [13, 23] and thus are given in Table 2.

Table 1: Parameters of the primary models for the microbial growth data under constant temperature conditions.

Microorganism
Temperature (°C)
Huang
Baranyi
Source of data

(d)
(d⁻¹)

(d)
(d⁻¹)

Pseudomonads
0
1.38
1.21
20.3
1.62
1.25
20.2
Koutsoumanis [21]

2
1.07
1.50
18.5
1.12
1.55
18.5

5
0.28
2.27
19.0
0.30
2.31
18.9

8
0.35
3.53
18.6
0.34
3.40
18.6

10
0.45
4.17
18.5
0.46
4.38
18.4

15
0.21
6.50
18.6
0.21
6.68
18.6

E. coli O157:H7
10
1.42
2.08
13.6
1.36
1.66
11.6
Koseki and Isobe [22]

15
0.34
8.02
14.5
0.31
6.63
13.9

20
0.23
19.61
14.6
0.22
14.37
14.6

25
0.14
25.00
15.6
0.19
24.32
15.7

Salmonella spp.
10
1.65
1.88
11.6
1.46
1.11
11.6
Koseki and Isobe [22]

15
0.28
8.63
12.6
0.27
7.19
12.4

20
0.23
16.84
14.0
0.22
14.37
13.8

25
0.15
32.89
14.6
0.12
22.66
15.0

Lactic acid bacteria
0
2.83
0.32
19.5
1.89
0.32
19.5
Mataragas et al. [16]

4
1.40
0.58
19.8
0.76
0.57
19.8

8
1.41
1.33
18.4
1.52
1.36
18.4

12
0.30
1.86
20.0
0.17
1.73
20.0

Table 2: Secondary models used for the microbial growth estimation.

Microorganism
Parameter
Huang
Baranyi
Source of data

Pseudomonads
(d) or

Koutsoumanis [21]

(d⁻¹)

E. coli O157:H7
(d) or

Koseki and Isobe [22]

(d⁻¹)

Salmonella spp.
(d) or

Koseki and Isobe [22]

(d⁻¹)

Lactic acid bacteria
(d) or

Mataragas et al. [16]

(d⁻¹)

Figures 1–4 show the results of microbial growth estimation for pseudomonads, E. coli O157:H7, Salmonella, and lactic acid bacteria under dynamic temperature conditions. All the estimated microbial evolutions through time generally agreed with those worked by the respective researchers [16, 21, 22]. Huang model gave predictions pretty close to the experimental data of pseudomonads on fish and almost identical to those given by Baranyi model (Figure 1). Only slight discrepancy from the experimental data happened after passing lag time before reaching the stationary phase in Figure 1(b) in the same way with Baranyi model. In predicting E. coli growth on the lettuce under dynamic temperature conditions, Huang model seems to be better than Baranyi model in the estimation except for the last duration of storage (Figure 2). Prediction of higher count than experimental data may be due to the slight decrease of the microbial load in the extended storage of the lettuce, which is not common in usual food storage and is rarely accounted for in most microbial growth modeling. The outcome of the higher prediction is thought to be a kind of fail-safe phenomenon. The similar trend was observed for the estimation of Salmonella growth (Figure 3). Both models were poor with prediction of higher Salmonella growth than the real data under dynamic temperature conditions of Figure 3(b). Almost same degree of estimations was observed between two models for lactic acid bacteria growth under two nonisothermal conditions (Figure 4). A slight but clear difference between two models is that Huang model shows more distinct ending of lag time for all the conditions.

(a)

(b)

(c)

(d)

(e)

Figure 1: Estimation of pseudomonads growth under several fluctuating temperature conditions. Experimental data are from Koutsoumanis [21]. Thin solid lines are the temperature, and thick solid and dotted lines are microbial counts estimated by Huang and Barayi models, respectively. ▲: experimental microbial count.

(a)

(b)

Figure 2: Estimation of E. coli O157:H7 growth under two fluctuating temperature conditions. Experimental data are from Koseki and Isobe [22]. Thin solid lines are the temperature, and thick solid and dotted lines are microbial counts estimated by Huang and Barayi models, respectively. ▲: experimental microbial count.

(a)

(b)

Figure 3: Estimation of Salmonella growth under two fluctuating temperature conditions. Experimental data are from Koseki and Isobe [22]. Thin solid lines are the temperature, and thick solid and dotted lines are microbial counts estimated by Huang and Barayi models, respectively. ▲: experimental microbial count.

(a)

(b)

Figure 4: Estimation of lactic acid bacteria growth under two fluctuating temperature conditions. Experimental data are from Mataragas et al. [16]. Thin solid lines are the temperature, and thick solid and dotted lines are microbial counts estimated by Huang and Barayi models, respectively. ▲: experimental microbial count.

Table 3 summarizes the comparison between two models in terms of the prediction parameters of bias factor , accuracy factor , and root-mean-square error , which were defined as follows:

where is the predicted microbial count (CFU/g), is the experimental count (CFU/g), and is number of data.

Table 3: Evaluation of the prediction of microbial growth under dynamic temperature conditions.

Conditions and microorganism
BF
AF
RMSE ( CFU/g)

Huang
Baranyi
Huang
Baranyi
Huang
Baranyi

Pseudomonads in Figure 1(a)
1.037
1.039
1.056
1.058
0.511
0.515

Pseudomonads in Figure 1(b)
1.081
1.083
1.089
1.087
0.633
0.607

Pseudomonads in Figure 1(c)
1.013
1.038
1.044
1.062
0.324
0.416

Pseudomonads in Figure 1(d)
1.009
1.027
1.045
1.056
0.325
0.387

Pseudomonads in Figure 1(e)
0.918
0.982
1.102
1.062
0.480
0.373

E. coli O157:H7 in Figure 2(a)
0.984
0.967
1.078
1.042
0.356
0.182

E. coli O157:H7 in Figure 2(b)
1.041
1.110
1.041
1.110
0.236
0.577

Salmonella spp. in Figure 3(a)
1.044
1.059
1.068
1.066
0.351
0.299

Salmonella spp. in Figure 3(b)
0.998
1.099
1.181
1.099
0.858
0.496

Lactic acid bacteria in Figure 4(a)
1.073
1.097
1.082
1.098
0.506
0.585

Lactic acid bacteria in Figure 4(b)
1.018
1.033
1.075
1.068
0.496
0.480

Overall in terms of the prediction parameters, Huang and Baranyi models showed a similar degree of prediction capacity in predicting the microbial growth under fluctuating temperature conditions (Table 3). For the used literature data sets, for Haung model ranged from 0.92 to 1.08 and that for Baranyi model was from 0.93 to 1.11. was from 1.04 to 1.18 for Haung model and from 1.04 to 1.11 for Baranyi model. was from 0.24 to 0.86 for Haung model and from 0.17 to 0.61 for Baranyi model. These ranges of and are considered adequate or acceptable [16, 23]. and values from both models in pseudomonads were closer to 1.0 than those obtained originally by Koutsoumanis [21] using four-parameter logistic primary model and square root secondary model. This is a verification showing the advantage of these two mechanistic models over the latter model. and values of Baranyi analysis for lactic acid bacteria on cured meat in Table 3 are close to those obtained by Mataragas et al. [16].

As a way to validate the microbial growth models, time to reach an acceptable limit of microbial load (10⁷ cells/g for pseudomonads and 108.4 cells/g for lactic acid bacteria) was calculated as shelf life estimate. The shelf life estimates based on the criteria set in the literatures were almost the same between Huang and Baranyi models and also similar to the predicted values given in the sources (Table 4).

Table 4: Shelf life prediction by Huang and Baranyi models.

Spoilage organism and temperature condition
Shelf life predicted (days)
Shelf life predicted in the source (days)
Criterion of spoilage level at end of shelf life

Huang
Baranyi

Pseudomonads in Figure 1(a)
2.7
2.7
3.0
10⁷/g

Pseudomonads in Figure 1(b)
3.6
3.7
3.8
10⁷/g

Pseudomonads in Figure 1(c)
4.3
4.1
4.2
10⁷/g

Pseudomonads in Figure 1(d)
4.0
3.8
3.9
10⁷/g

Pseudomonads in Figure 1(e)
7.2
6.7
6.9
10⁷/g

Lactic acid bacteria in Figure 4(a)
28.3
28.2
28
108.4/g

Lactic acid bacteria in Figure 4(b)
25.7
25.6
24
108.4/g

Huang model has advantage of directly determining the growth rate and lag time from the fitting of its primary model to the microbial growth data at constant temperature, which can also be used for building secondary model. Relationships of growth rate and lag time as function of temperature may affect the model’s capability to predict the microbial growth under nonisothermal conditions. Huang [19] proposed (9) for growth rate and (8) for the lag time. Because it is difficult to predict lag phase duration and express its temperature dependence in simple square root model [24], there may still remain a possible limitation in Huang model’s capability to predict the microbial growth under varying temperature conditions. Baranyi model has been set up based on theoretical hypothesis of microbial growth and usually needs only (7) for secondary model defining the temperature dependence of the growth rate [13, 16, 22]. However, selection of proper value or initial physiological cell state is known to be critical for successful estimation of microbial growth under dynamic temperature conditions by using two simultaneous differential equations (3) and (4) [13, 22]. Even though values from different temperature are averaged for use in fluctuating temperature storage, too much different values from several temperature conditions make it difficult to determine the proper value [7, 19]. From the finding that Huang model is comparable to Baranyi model in estimating the microbial growth under dynamic temperature conditions, its simplicity in using single differential equation (1) with direct incorporation of the growth rate and lag time may work as an advantage for application of plugging into quality estimation devices such as RFID or other electronic sensors, which are being introduced recently for the shelf life control [10, 11].

Even with some limitations, Huang and Baranyi models are found to work effectively to predict the microbial growth dynamics under changing temperature conditions being useful for shelf life estimation and control in dynamic food supply chain.

4. Conclusions

Huang model was comparable to Baranyi model in the capability to estimate microbial growth under dynamic temperature conditions. Its simple form of single differential equation incorporating directly the growth rate and lag time may work as an advantage to be used in online shelf life estimation by using the electronic device.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by Kyungnam University Foundation Grant, 2013. The author thanks Yun Hee Jo, Min-Ji Kwon, and Ji Hye Lee for reading the data from the literatures.

References

T. A. McMeekin and T. Ross, “Shelf life prediction: status and future possibilities,” International Journal of Food Microbiology, vol. 33, no. 1, pp. 65–83, 1996.

K. Koutsoumanis and G.-J. E. Nychas, “Application of a systematic experimental procedure to develop a microbial model for rapid fish shelf life predictions,” International Journal of Food Microbiology, vol. 60, no. 2-3, pp. 171–184, 2000.

D. S. Lee, “Packaging and the microbial shelf life of food,” in Food Packaging and Shelf Life, G. L. Robertson, Ed., pp. 55–79, CRC Press, Boca Raton, Fla, USA, 2009.

T. Tsironi, E. Gogou, E. Velliou, and P. S. Taoukis, “Application and validation of the TTI based chill chain management system SMAS (Safety Monitoring and Assurance System) on shelf life optimization of vacuum packed chilled tuna,” International Journal of Food Microbiology, vol. 128, no. 1, pp. 108–115, 2008.

S. Bruckner, B. Petersen, and J. Kreyenschmidt, “Influence of cold chain interruptions on the shelf life of fresh pork and poultry,” International Journal of Food Science and Technology, vol. 47, no. 8, pp. 1639–1646, 2012.

P. S. Taoukis, K. Koutsoumanis, and G. J. E. Nychas, “Use of time-temperature integrators and predictive modelling for shelf life control of chilled fish under dynamic storage conditions,” International Journal of Food Microbiology, vol. 53, no. 1, pp. 21–31, 1999.

D. S. Lee, K.-J. Hwang, D. S. An, J. P. Park, and H. J. Lee, “Model on the microbial quality change of seasoned soybean sprouts for on-line shelf life prediction,” International Journal of Food Microbiology, vol. 118, no. 3, pp. 285–293, 2007.

I. Seo, D. S. An, and D. S. Lee, “Development of convenient software for online shelf-life decisions for Korean prepared side dishes based on microbial spoilage,” Food Science and Biotechnology, vol. 18, no. 5, pp. 1243–1252, 2009.

J. K. Heising, M. Dekker, P. V. Bartels, and M. A. J. S. (Tiny) Van Boekel, “Monitoring the quality of perishable foods: opportunities for intelligent packaging,” Critical Reviews in Food Science and Nutrition, vol. 54, no. 5, pp. 645–654, 2014.

L. Qi, M. Xu, Z. Fu, T. Mira, and X. Zhang, “C²SLDS: a WSN-based perishable food shelf-life prediction and LSFO strategy decision support system in cold chain logistics,” Food Control, vol. 38, no. 1, pp. 19–29, 2014.

M. Kärkkäinen, “Increasing efficiency in the supply chain for short shelf life goods using RFID tagging,” International Journal of Retail & Distribution Management, vol. 31, no. 10, pp. 529–536, 2003.

K. L. Yam, P. T. Takhistov, and J. Miltz, “Intelligent packaging: concepts and applications,” Journal of Food Science, vol. 70, no. 1, pp. R1–R10, 2005.

J. Baranyi, T. P. Robinson, A. Kaloti, and B. M. Mackey, “Predicting growth of Brochothrix thermosphacta at changing temperature,” International Journal of Food Microbiology, vol. 27, no. 1, pp. 61–75, 1995.

L. Huang, “Estimation of growth of Clostridium perfringens in cooked beef under fluctuating temperature conditions,” Food Microbiology, vol. 20, no. 5, pp. 549–559, 2003.

K. Koutsoumanis, P. S. Taoukis, and G. J. E. Nychas, “Development of a safety monitoring and assurance system for chilled food products,” International Journal of Food Microbiology, vol. 100, no. 1–3, pp. 253–260, 2005.

M. Mataragas, E. H. Drosinos, A. Vaidanis, and I. Metaxopoulos, “Development of a predictive model for spoilage of cooked cured meat products and its validation under constant and dynamic temperature storage conditions,” Journal of Food Science, vol. 71, no. 6, pp. M157–M167, 2006.

J. F. van Impe, F. Poschet, A. H. Geeraerd, and K. M. Vereecken, “Towards a novel class of predictive microbial growth models,” International Journal of Food Microbiology, vol. 100, no. 1–3, pp. 97–105, 2005.

J. Baranyi, “Mathematics of predictive food microbiology,” International Journal of Food Microbiology, vol. 26, no. 2, pp. 199–218, 1995.

L. Huang, “A new mechanistic growth model for simultaneous determination of lag phase duration and exponential growth rate and a new Bělehdrádek-type model for evaluating the effect of temperature on growth rate,” Food Microbiology, vol. 28, no. 4, pp. 770–776, 2011.

K. Bernaerts, E. Dens, K. Vereecken et al., “Modeling microbial dynamics under time-varying conditions,” in Modeling Microbial Responses in Food, R. C. McKellar and X. Lu, Eds., pp. 243–261, CRC Press, Boca Raton, Fla, USA, 2004.

K. Koutsoumanis, “Predictive modeling of the shelf life of fish under nonisothermal conditions,” Applied and Environmental Microbiology, vol. 67, no. 4, pp. 1821–1829, 2001.

S. Koseki and S. Isobe, “Prediction of pathogen growth on iceberg lettuce under real temperature history during distribution from farm to table,” International Journal of Food Microbiology, vol. 104, no. 3, pp. 239–248, 2005.

D. A. Longhi, F. Dalcanton, G. M. de Aragão, B. A. Carciofi, and J. B. Laurindo, “Assessing the prediction ability of different mathematical models for the growth of Lactobacillus plantarum under non-isothermal conditions,” Journal of Theoretical Biology, vol. 335, no. 1, pp. 88–96, 2013.

I. A. M. Swinnen, K. Bernaerts, E. J. J. Dens, A. H. Geeraerd, and J. F. van Impe, “Predictive modelling of the microbial lag phase: a review,” International Journal of Food Microbiology, vol. 94, no. 2, pp. 137–159, 2004.