Bioinactivation: Software for modelling dynamic microbial inactivation.

This contribution presents the bioinactivation software, which implements functions for the modelling of isothermal and non-isothermal microbial inactivation. This software offers features such as user-friendliness, modelling of dynamic conditions, possibility to choose the fitting algorithm and generation of prediction intervals. The software is offered in two different formats: Bioinactivation core and Bioinactivation SE. Bioinactivation core is a package for the R programming language, which includes features for the generation of predictions and for the fitting of models to inactivation experiments using non-linear regression or a Markov Chain Monte Carlo algorithm (MCMC). The calculations are based on inactivation models common in academia and industry (Bigelow, Peleg, Mafart and Geeraerd). Bioinactivation SE supplies a user-friendly interface to selected functions of Bioinactivation core, namely the model fitting of non-isothermal experiments and the generation of prediction intervals. The capabilities of bioinactivation are presented in this paper through a case study, modelling the non-isothermal inactivation of Bacillus sporothermodurans. This study has provided a full characterization of the response of the bacteria to dynamic temperature conditions, including confidence intervals for the model parameters and a prediction interval of the survivor curve. We conclude that the MCMC algorithm produces a better characterization of the biological uncertainty and variability than non-linear regression. The bioinactivation software can be relevant to the food and pharmaceutical industry, as well as to regulatory agencies, as part of a (quantitative) microbial risk assessment.

[1]  Karline Soetaert,et al.  Inverse Modelling, Sensitivity and Monte Carlo Analysis in R Using Package FME , 2010 .

[2]  Heikki Haario,et al.  DRAM: Efficient adaptive MCMC , 2006, Stat. Comput..

[3]  Kirk D. Dolan,et al.  Parameter estimation for dynamic microbial inactivation: which model, which precision? , 2013 .

[4]  I. Leguerinel,et al.  On calculating sterility in thermal preservation methods: application of the Weibull frequency distribution model. , 2001, International journal of food microbiology.

[5]  J. Rose,et al.  Quantitative Microbial Risk Assessment , 1999 .

[6]  Kirk D. Dolan,et al.  Estimation of Kinetic Parameters for Nonisothermal Food Processes , 2003 .

[7]  R. C. Whiting,et al.  Microbial modeling in foods. , 1995, Critical reviews in food science and nutrition.

[8]  Zafiro Aspridou,et al.  Individual cell heterogeneity as variability source in population dynamics of microbial inactivation. , 2015, Food microbiology.

[9]  C. R. Stumbo,et al.  Thermobacteriology in Food Processing , 1973 .

[10]  Antonio Valero,et al.  Predictive Microbiology in Foods , 2012 .

[11]  A H Geeraerd,et al.  Structural model requirements to describe microbial inactivation during a mild heat treatment. , 2000, International journal of food microbiology.

[12]  Charles P. Gerba,et al.  Quantitative Microbial Risk Assessment: Second Edition , 2014 .

[13]  M Peleg,et al.  Reinterpretation of microbial survival curves. , 1998, Critical reviews in food science and nutrition.

[14]  A. Esnoz,et al.  Nonisothermal heat resistance determinations with the thermoresistometer Mastia , 2009, Journal of applied microbiology.

[15]  P. Fernández,et al.  Thermal inactivation kinetics of Bacillus stearothermophilus spores using a linear temperature program. , 1999, Journal of food protection.

[16]  A H Geeraerd,et al.  Inactivation model equations and their associated parameter values obtained under static acid stress conditions cannot be used directly for predicting inactivation under dynamic conditions. , 2008, International journal of food microbiology.

[17]  József Baranyi,et al.  Predictive modelling of Salmonella: From cell cycle measurements to e-models , 2012 .

[18]  W. D. Bigelow,et al.  The logarithmic nature of thermal death time curves , 1921 .

[19]  Karline Soetaert,et al.  Solving Differential Equations in R: Package deSolve , 2010 .

[20]  C. A. Ferreira,et al.  One-step global parameter estimation of kinetic inactivation parameters for Bacillus sporothermodurans spores under static and dynamic thermal processes. , 2016, Food research international.

[21]  Jorge J. Moré,et al.  The Levenberg-Marquardt algo-rithm: Implementation and theory , 1977 .

[22]  S. Condón,et al.  Predicting microbial heat inactivation under nonisothermal treatments. , 2007, Journal of food protection.

[23]  Arturo Esnoz,et al.  High Heating Rates Affect Greatly the Inactivation Rate of Escherichia coli , 2016, Front. Microbiol..

[24]  Mariem Ellouze,et al.  Software for predictive microbiology and risk assessment: a description and comparison of tools presented at the ICPMF8 Software Fair. , 2015, Food microbiology.

[25]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[26]  Robert Tibshirani,et al.  An Introduction to the Bootstrap , 1994 .

[27]  Alejandro Amézquita,et al.  Characterization of Bacillus sporothermodurans IC4 spores; putative indicator microorganism for optimisation of thermal processes in food sterilisation , 2010 .