On the use of in-silico simulations to support experimental design: A case study in microbial inactivation of foods

The mathematical models used in predictive microbiology contain parameters that must be estimated based on experimental data. Due to experimental uncertainty and variability, they cannot be known exactly and must be reported with a measure of uncertainty (usually a standard deviation). In order to increase precision (i.e. reduce the standard deviation), it is usual to add extra sampling points. However, recent studies have shown that precision can also be increased without adding extra sampling points by using Optimal Experiment Design, which applies optimization and information theory to identify the most informative experiment under a set of constraints. Nevertheless, to date, there has been scarce contributions to know a priori whether an experimental design is likely to provide the desired precision in the parameter estimates. In this article, two complementary methodologies to predict the parameter precision for a given experimental design are proposed. Both approaches are based on in silico simulations, so they can be performed before any experimental work. The first one applies Monte Carlo simulations to estimate the standard deviation of the model parameters, whereas the second one applies the properties of the Fisher Information Matrix to estimate the volume of the confidence ellipsoids. The application of these methods to a case study of dynamic microbial inactivation, showing how they can be used to compare experimental designs and assess their precision, is illustrated. The results show that, as expected, the optimal experimental design is more accurate than the uniform design with the same number of data points. Furthermore, it is demonstrated that, for some heating profiles, the uniform design does not ensure that a higher number of sampling points increases precision. Therefore, optimal experimental designs are highly recommended in predictive microbiology.

[1]  Heidi den Besten,et al.  Quantifying variability on thermal resistance of Listeria monocytogenes. , 2015, International journal of food microbiology.

[2]  G. Milliken Nonlinear Regression Analysis and Its Applications , 1990 .

[3]  K J Versyck,et al.  Introducing optimal experimental design in predictive modeling: a motivating example. , 1999, International journal of food microbiology.

[4]  Sandro Macchietto,et al.  Model-based design of experiments for parameter precision: State of the art , 2008 .

[5]  Kirk D. Dolan,et al.  Parameter estimation in food science. , 2013, Annual review of food science and technology.

[6]  J. F. Van Impe,et al.  Sensitivity analysis of microbial growth parameter distributions with respect to data quality and quantity by using Monte Carlo analysis , 2004, Math. Comput. Simul..

[7]  Ryuei Nishii,et al.  Optimality of experimental designs , 1993, Discret. Math..

[8]  B. Carciofi,et al.  Optimal experimental design for improving the estimation of growth parameters of Lactobacillus viridescens from data under non-isothermal conditions. , 2017, International journal of food microbiology.

[9]  Joseph G. Pigeon,et al.  Statistics for Experimenters: Design, Innovation and Discovery , 2006, Technometrics.

[10]  Jose Lucas Peñalver-Soto,et al.  Optimal characterization of thermal microbial inactivation simulating non-isothermal processes. , 2018, Food research international.

[11]  Alberto Garre,et al.  Bioinactivation: Software for modelling dynamic microbial inactivation. , 2017, Food research international.

[12]  Maarten Nauta,et al.  Separation of uncertainty and variability in quantitative microbial risk assessment models. , 2000 .

[13]  M. S. Khots,et al.  D-optimal designs , 1995 .

[14]  Xiangzhong Xie,et al.  The Impact of Global Sensitivities and Design Measures in Model-Based Optimal Experimental Design , 2018 .

[15]  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.

[16]  J. R. Esty,et al.  The Thermal Death Point in Relation to Time of Typical Thermophilic Organisms , 1920 .

[17]  M H Zwietering,et al.  Quantifying strain variability in modeling growth of Listeria monocytogenes. , 2015, International journal of food microbiology.

[18]  J. Claes,et al.  Practical Identification of Unstructured Growth Kinetics by Application of Optimal Experimental Design , 1997 .

[19]  J. V. Van Impe,et al.  Parameter estimations in predictive microbiology: Statistically sound modelling of the microbial growth rate. , 2017, Food research international.

[20]  E. L. Lehmann,et al.  Theory of point estimation , 1950 .

[21]  Alberto Garre,et al.  Bioinactivation FE: A free web application for modelling isothermal and dynamic microbial inactivation. , 2018, Food research international.

[22]  G. Aragão,et al.  Optimal experimental design to model spoilage bacteria growth in vacuum-packaged ham , 2018 .

[23]  P. Periago,et al.  When nanoemulsified, d-limonene reduces Listeria monocytogenes heat resistance about one hundred times , 2016 .

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

[25]  I Stamati,et al.  Optimal experimental design for discriminating between microbial growth models as function of suboptimal temperature: From in silico to in vivo. , 2016, Food research international.

[26]  Eva Balsa-Canto,et al.  Parameter estimation and optimal experimental design. , 2008, Essays in biochemistry.

[27]  Eva Balsa-Canto,et al.  Optimal design of dynamic experiments for improved estimation of kinetic parameters of thermal degradation , 2007 .

[28]  M. Ros-Chumillas,et al.  Nanoemulsified D-Limonene Reduces the Heat Resistance of Salmonella Senftenberg over 50 Times , 2017, Nanomaterials.

[29]  P. Fernández,et al.  Mathematical quantification of the induced stress resistance of microbial populations during non-isothermal stresses. , 2018, International journal of food microbiology.

[30]  Míriam R. García,et al.  Quality and shelf-life prediction for retail fresh hake (Merluccius merluccius). , 2015, International journal of food microbiology.

[31]  Ashim K. Datta,et al.  Error estimates for approximate kinetic parameters used in food literature , 1993 .

[32]  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.

[33]  Marcio Schwaab,et al.  Optimum reference temperature for reparameterization of the Arrhenius equation. Part 1 : Problems involving one kinetic constant , 2007 .

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

[35]  Philippe Nimmegeers,et al.  Comparing design of experiments and optimal experimental design techniques for modelling the microbial growth rate under static environmental conditions. , 2018, Food microbiology.

[36]  K. Koutsoumanis,et al.  Predictive Modeling of the Shelf Life of Fish under Nonisothermal Conditions , 2001, Applied and Environmental Microbiology.

[37]  Eva Balsa-Canto,et al.  An iterative identification procedure for dynamic modeling of biochemical networks , 2010, BMC Systems Biology.

[38]  Peter A. Vanrolleghem,et al.  Estimating the parameters of the Baranyi model for bacterial growth , 1999 .

[39]  P. Periago,et al.  Estimation of the non-isothermal inactivation patterns of Bacillus sporothermodurans IC4 spores in soups from their isothermal survival data. , 2004, International journal of food microbiology.

[40]  P. Fernández,et al.  Effect of storage conditions in the response of Listeria monocytogenes in a fresh purple vegetable smoothie compared with an acidified TSB medium. , 2018, Food microbiology.

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

[42]  J. Banga,et al.  Computational procedures for optimal experimental design in biological systems. , 2008, IET systems biology.

[43]  Tilman Barz,et al.  Monte Carlo Simulations for the Analysis of Non-linear Parameter Confidence Intervals in Optimal Experimental Design , 2019, Front. Bioeng. Biotechnol..