- Open Access
Evaluation of toxic effects of several carboxylic acids on bacterial growth by toxicodynamic modelling
Microbial Cell Factories volume 10, Article number: 100 (2011)
Effects of organic acids on microbial fermentation are commonly tested in investigations about metabolic behaviour of bacteria. However, they typically provide only descriptive information without modelling the influence of acid concentrations on bacterial kinetics.
We developed and applied a mathematical model (secondary model) to capture the toxicological effects of those chemicals on kinetic parameters that define the growth of bacteria in batch cultures. Thus, dose-response kinetics were performed with different bacteria (Leuconostoc mesenteroides, Carnobacterium pisicola, Escherichia coli, Bacillus subtilis and Listonella anguillarum) exposed at increasing concentrations of individual carboxylic acids (formic, acetic, propionic, butyric and lactic). In all bioassays the acids affected the maximum bacterial load (X m ) and the maximum growth rate (v m ) but only in specific cases the lag phase (λ) was modified. Significance of the parameters was always high and in all fermentations the toxicodynamic equation was statistically consistent and had good predictability. The differences between D and L-lactic acid effects were significant for the growth of E. coli, L. mesenteroides and C. piscicola. In addition, a global parameter (EC50,τ) was used to compare toxic effects and provided a realistic characterization of antimicrobial agents using a single value.
The effect of several organic acids on the growth of different bacteria was accurately studied and perfectly characterized by a bivariate equation which combines the basis of dose-response theory with microbial growth kinetics (secondary model). The toxicity of carboxylic acids was lower with the increase of the molecular weight of these chemicals.
The study of the inhibitory capacity of antibacterial chemicals on microbial growth is generally based on point estimates of the effect, even though microbial exposure to them may be associated with complex kinetic profiles . Evaluation of the whole time course of observed effects using a toxicodynamic analysis would be a more interesting approach in characterizing the corresponding bacterial responses . Among antibacterial agents, carboxylic acids are commonly applied as preservative of foods, disinfectant of materials and surfaces, agent for the control of fermentations, extracting solvent of biological compounds and substrates for biopolymer manufacturing [3–5]. Different features as susceptibility, adaptability, tolerance, resistance and survival of several bacteria to those weak acids have been extensively reported [6–9]. Thereby, Escherichia coli is one of the most commonly studied bacteria because it is a well-known food-poisoning pathogen [10, 11]. Listonella anguillarum (also know as Vibrio anguillarum) is another bacteria that has attracted scientific interest due to its association with high mortalities in aquaculture [12–14] in which probiotics and organic acid treatments are recommended [15–17]. Consequently, the accurate modelling of the acid effects on the microbial growth is necessary and indispensable to describe and compare the efficiencies of different treatments, optimize their inhibitory properties and dosing strategies as well as develop standardized protocols of application.
The toxicodynamic analysis using dose-response (DR) bioassay is a powerful tool widely used in different experimental contexts [18, 19]. Based on this perspective, a bivariate model formed by a logistic equation predicting growth profiles and another sigmoid equation simulating DR tendencies was previously investigated to evaluate the simultaneous effects of detergent , alkyl esters  and heavy metals  on microbial growth or mortality. In this type of models the most representative kinetic parameters (maximum bacterial load, lag phase and maximum growth rate) are non-linearly affected by the concentration or dose of a chemical. Furthermore, simultaneous fits of all experimental data from control and toxic-dosed cultures must be accomplished for correct modelling . This proposal is experimentally more realistic than the conventional toxicological assessment that is focused on the estimation of the specific growth rate from biomass quantified at two data in the exponential phase .
The aim of the present work is to evaluate the effect of carboxylic acids on different bacteria using a toxicodynamic model. This mathematical model was formulated by the combination of the Weibull equation for dose-response description and the logistic equation for bacterial growth in a bivariate model. Data from growth curve studies of Leuconostoc mesenteroides (Ln), Carnobacterium pisicola (Cb), E. coli (Ec), Bacillus subtilis (Bs) and L. anguillarum (La) exposed to a wide range of concentrations of five carboxylic acids (formic, acetic, propionic, butyric and lactic) is used for model validation. The accuracy and suitability of the proposed model is extensively demonstrated in all bacterial kinetics and the numerical estimates allow us a complete toxic characterization of carboxylic acid effects on growth parameters.
Selected representation of the tested combinations is depicted in Figure 1 where experimental data from each culture were simultaneously fitted to equation (1). In all cases, surfaces predicted by the model showed good agreement with data obtained in the bioassays. Graphical responses indicated the difference of bacterial sensitivity to the acids exposed. For instance, formic acid was more toxic for La growth than for the production of biomass by Bs. Table 1 summarizes the parametric estimates and statistical analysis of the effects generated by formic and acetic acids on the growth of the five bacteria tested. In most of cultures affected by those acids the effects on lag phase were statistically not significant. Only in the cultures of Cb and Ln dosed with formic, significant effects for λ-parameter were observed.
The numerical results for propionic and butyric acids are shown in Table 2. In all bacterial cultures, these weak acids involved significant changes in two parameters from the logistic equation (maximum bacterial load Xm and maximum growth rate vm) but not in lag phase (λ). However, L+lactic affected the three parameters of growth on Cb and Ec fermentations (Table 3). The lack of significance in the lag phase was clear in the rest of bacteria. In addition, the comparative effect between the two isomeric forms of lactic acid (D- and L+) was also assayed for these last bacteria. Two of these cases, Ec and Ln, are depicted in Figure 2.
The fitting of results was always satisfactory both graphically and statistically. The mathematical equations were robust and consistent (p-values < 0.001 from Fisher's F test), the residuals were randomly distributed and autocorrelations were not observed by Durbin-Watson test (data not shown). Furthermore, all the adjusted coefficients of multiple determination between predicted and observed values were always R2adj > 0.95, with a wide majority of the fittings superior at 0.99. Bias and accuracy factors (Bf and Af) also indicated the lack of bias and high accuracy of equation (1) to describe experimental effects of acids on bacterial growth.
For the global description of the carboxylic acid toxicity, a single index -EC 50,τ - was used (Table 4). This parameter is a summary of all the effects on the biomass produced at a given time -τ- or time required to achieve the semi-maximum biomass . This value could be the main parameter with practical interest in terms of environmental assessment and operational applications of chemicals. Thereby, comparison of toxicity among different agents and bacterial sensitivity can be also evaluated using this parameter. Our results revealed that the toxicity of carboxylic acids almost always decreased (higher value of EC 50,τ ) with increasing of their molecular weights but not depending on pKa value. Thus, formic and acetic acids led to greater toxic values than butyric and lactic acids. In all treatments, La was the microorganism most susceptible to the acids effects and Ln and Bs the most resistant. In addition, L+lactic was more toxic than D-lactic on Ec growth whereas the effect of D-lactic on Ln was higher than its optical isomer.
The validation and generalization of equation (1) was studied using relative viable cell count as dependent variable of bacterial load. Four cases were selected and mathematical description was elaborated (Figure 3). The parameter estimations were in agreement with previously reported using optical density at 700 nm (Table 5). Lag phase of Cb and Ln was not significantly affected by acetic, formic and propionic acids. However, both parameters, X m and v m , were significantly modified by those acids. The statistical consistency, goodness-of-fit, lack of autocorrelations for residuals and high predictability of the model was clearly showed.
The combination of sigmoid equations for modelling DR relationships of agents affecting the most important parameters from microbial growth has been scarcely reported. However, excellent results were obtained when it was applied to the study of detergent  and heavy metals  on algae and bacterial kinetics, respectively. From the results of acid-dosed cultures performed in this study, it can be seen that the bacterial system showed different profiles of growth and responses to the toxic effect of carboxylic acids. All 3D-surfaces of experimental data were perfectly modelled by the proposed equation (1) and good predictability of this bivariate model was also demonstrated. Statistical analysis confirmed these results. In all cases, distribution of residuals was randomly scattered around zero and grouped data and autocorrelations were not observed. Adjusted coefficients of multiple determination indicated the goodness of fit and p-values from Fisher's F test showed the robustness and consistence of mathematical model. Bf and Af factors also revealed the good agreement between experimental and simulated data. In addition, we have shown that the parameter -EC 50,τ - used in the global characterization of toxic effects by carboxylic acids produces an interesting resource to compare toxicities between chemicals and sensibility of microorganisms to them.
Although absorbance measurements have been less used within microbiological modelling, it is a simple and inexpensive technique, in comparison with laborious and expensive plate count methodology, for generation of several growth data (i.e., bivariate experiments) that allow the formulation of accuracy mathematical models as equation (1). This equation was also validated using relative population size  as growth response. With present data, the use of relative viable cell count as dependent variable avoided the problem reported by Dalgaard and Koutsoumanis , in where the differences of lag phase and rate parameters obtained by two kinds of bacterial load quantification (count plate and absorbance) limited the ability of comparison among different experimental conditions. Indeed, the units of parameters and the numerical value of them are different (Tables 6 and 5) but behaviour among two approaches was easily comparable.
In the cell measurement, the parameter of rate from equation (1) was defined as maximum specific growth rate (units of h-1) meanwhile, in optical quantification, it was denominated as maximum growth rate (units of AU h-1). Both parameters are related in the description of the growth rate of bacteria but those estimates can not be directly compared. However, the tendencies predicted by proposed equation in both quantifications of growth were similar and statistically significant.
Previous works have reported the importance of the concentration of carboxylic acid exposures in suppressing the growth of bacteria [7, 26–28]. Conventional studies of bacterial resistance , adaptation  and membrane permeabilization  have demonstrated the dependence of metabolic responses with the microorganism and acid concentration evaluated. Nonetheless, these results have not been used optimally to characterize mechanisms and modes of action underlying the effects, to guide the selection of dosing ranges and to compare the effectiveness of acid treatments. These descriptions are very important when weak acids must be applied to reduce the levels of food-borne pathogens on food processing, control industrial fermentations or avoid contaminations on surfaces and devices. Consequently, our modelling approach provides a simple and consistent tool for such descriptions.
In all cases, we have found that carboxylic acids significantly affected the kinetic parameters, maximum biomass and maximum growth rate, but only in specific cultures the lag phase was altered by these chemicals. These heterogeneous responses in the lag phase duration have been also observed using heavy metals [22, 31], lactic acid  and salts and bacteriocins  as effectors. Cherrington et al.  suggested that direct comparisons of different organic acids in relation to their antimicrobial activities are difficult due to the variation in physical characteristics. However, we have demonstrated that the effects and toxicity of carboxylic acids can be successfully compared and significant kinetic parameters can be obtained. Thus, the toxicity was associated to the molecular weight and the sensitivity of bacteria was very varied being La the most clearly sensitive to acids (Table 4). Regarding with stereoisomers of lactic acid, antimicrobial capacity of D-lactic on Cb and Ec fermentations was lower than L+lactic. This finding is in agreement with previously reported by Leitch and Stewart  in which the susceptibility of several strains of E. coli was greater to L+lactic than that observed with D-lactic, but not depending on transmembrane pH gradient.
In summary, a general bivariate model has been proposed for the characterization of the effects from five carboxylic acids (formic, acetic propionic, butyric and lactic) against five bacteria (C. piscicola, E. coli, B. subtilis, L. mesenteroides and L. anguillarum). This equation could be easily applied to other chemicals and microorganisms and may support a resource for exhaustive comparison and formal description of agent effects on bacterial growth. Toxicity of acids decreased with the increase of molecular weight and L. anguillarum was the least resistant bacterium. In addition, L+lactic was more toxic than D-lactic for the growth of C. piscicola and E. coli.
Microbiological methods and chemicals
The studied bacteria are shown in Table 7. These microorganisms were chosen by their different features about cell wall structure (Gram-positive and negative), behaviour (free, opportunistic parasite, probiotic), habitats (marine and terrestrial) and metabolic properties (homo and heterofermentative). L. anguillarum and L. mesenteroides were kindly provided by Dr. Harry Birkbeck (University of Glasgow, UK) and Dr. B. Ray (University of Wyoming, Laramie, USA), respectively.
Stock cultures of bacteria were kept at -80°C in commercial MRS (L. mesenteroides HD-IIM_1 and C. piscicola CECT 4020), marine (L. anguillarum 90-11-287) and tryptone soy broths (B. subtilis CECT 35 and E. coli CECT 731) with 25% glycerol [36, 37]. Marine medium were provided by Difco (Becton, Dickinson and Company, MD, USA), MRS (de Man, Rogosa and Sharpe) medium by Pronadisa (Hispanlab S.A., Spain) and tryptone soy broth (TSB) by Panreac (Panreac Química SA, Barcelona, Spain). Culture media were prepared as indicating on commercial formulation and sterilized at 121°C for 15 min. Carboxylic acids were in all cases purchased from Sigma-Aldrich (St. Louis, MO, USA). Concentrated solutions of these acids were separately prepared and sterilized with steam flow at 101°C for 1 h.
Bioassay and culture conditions
Fermentations were carried out in triplicate using methods which were described in detail in previous reports [38, 39]. To prepare the microbial suspensions, 12-h cultures in commercial media were centrifuged at 4,000 g for 15 min and the sediments washed with phosphate buffer 0.005 M. Subsequently, the sediments were again centrifuged and resuspended in the corresponding buffed medium for each bacterium: 0.05 M pH = 6.0 biphtalate-NaOH buffer in MRS (Ln and Cb), 0.05 M pH = 7.4 Tris-HCl buffer in marine medium (La) and 0.05 M pH = 7.4 Tris-HCl buffer in TSB (Bs and Ec) and adjusted to an absorbance (700 nm, A700) of 0.200. For bioassay, equal volumes of individual concentrations of each acid (Table 8) and microbial suspension were mixed in 15 ml tubes. Thus, under these conditions the inocula concentrations for the tested bacteria were established among 106 and 107 cfu ml-1. All incubations were performed with orbital shaking at 200 rpm and 22°C (La), 30°C (Ln, Cb) and 37°C (Bs, Ec).
At pre-established times, samples were centrifuged at 4,000 g for 15 min. Sediments (bacterial biomass) were washed and resuspended in distilled water to the appropriate dilution for measuring the bacterial growth by A700. For comparative purposes, some sediments were used for quantify viable cells by means of plate count technique on commercial media with agar. Serial, 10-fold dilutions were prepared in peptone-buffered solutions and 0.1 ml samples were plated in triplicate, incubated at previously indicated temperatures and manually counted (after 3 days of incubation). Results were expressed as relative viable cell count [24, 37]: ln (N/N0); where N is the colony-forming units per ml (cfu ml-1) and N0 is the initial cfu ml-1.
Recently, we have proposed a bivariate equation (secondary model) that modelled satisfactorily the effect of three heavy metals on the growth parameters of several bacteria . This equation is based on the combination of Weibull function as DR model [40, 41] affecting the most important parameters of the reparametrized logistic equation used for growth description :
This equation can be easily and accurately used when the dependent variable or response (growth and population size) is described by means of optical density (absorbance at 700 nm) or relative viable cell count (ln (N/N0)). In the present work, X is the growth measured as optical density or relative population size, v m is the maximum growth rate, X m is the maximum bacterial load (growth), λ is the lag phase and C is the acid concentration. The meaning of the rest of symbolic notations and corresponding units are summarized in Table 6.
In addition, a global parameter (EC 50,τ ) was also defined for the overall description of chemical effects on kinetic studies. This parameter was defined as the dose of agent (in mM) that reduces the biomass by 50% compared to that produced by the control at time (τ) that reduces the biomass by 50% .
Fitting procedures and parametric estimates from the experimental results were performed by minimisation of the sum of quadratic differences between observed and model-predicted values, using the nonlinear least-squares (quasi-Newton) method provided by the macro 'Solver' of Microsoft Excel 2003 spreadsheet. Firstly, growth curves without acid were fitted and subsequently effect of acid on parameters was calculated with fixed values of X m , v m and λ. For relative viable cell count, an assumption was established: when (N/N0) data were < 1 then values of ln (N/N0) were forced to zero. Subsequently, confidence intervals from the parametric estimates (Student's t test), consistence of mathematical models (Fisher's F test) and residual analysis (Durbin-Watson test) were determined with 'SolverAid' macro, which is freely available from de Levie's Excellaneous website: http://www.bowdoin.edu/~rdelevie/excellaneous/. In addition, bias and accuracy factors of the equation (1) were calculated to evaluate the fitting of that model to experimental data :
where log (predicted/observed) is the logarithmic relation between the predicted and the experimental values, and n is the number of data. The nearer the values of Bf (bias factor) and Af (accuracy factor) to 1 indicate the better the fitting of the models to experimental data (a value of 1 indicates that there is perfect agreement among predicted and observed data).
Tam VH, Schilling AN, Nikolaou M: Modelling time-kill studies to discern the pharmacodynamics of meropenem. J Antimicrob Chemother. 2005, 55: 699-706.
Murado MA, Vázquez JA: Biphasic toxicodynamic features of some antimicrobial agents on microbial growth: a dynamic mathematical model and its implications on hormesis. BMC Microbiol. 2010, 10: 220-
Beek AT, Brul S: To kill or not to kill Bacilli: Opportunities for food biotechnology. Cur Opin Biotechnol. 2010, 21: 168-174.
Hirshfield IN, Terzulli S, O'Byrne C: Weak organic acids: a panoply of effects on bacteria. Sci Progr. 2003, 86: 245-269.
Lu D, Zhang X, Zhou T, Ren Z, Wang S, Lei Z: Biodegradable poly (lactic acid) copolymers. Progr Chem. 2008, 20: 339-350.
Choi YM, Kim OY, Kim KH, Kim BC, Rhee MS: Combined effect of organic acids and supercritical carbon dioxide treatments against nonpathogenic Escherichia coli, Listeria monocytogenes, Salmonella typhimurium and E. coli O157:H7 in fresh pork. Lett Appl Microbiol. 2009, 49: 510-515.
Hsiao CP, Siebert KJ: Modeling the inhibitory effects of organic acids on bacteria. Int J Food Microbiol. 1999, 47: 189-201.
Hua Q, Joyce AR, Palsson BØ, Fong SS: Metabolic characterization of Escherichia coli strains adapted to growth on lactate. Appl Environ Microbiol. 2007, 73: 4639-4647.
Leyer GJ, Wang LL, Johnson EA: Acid adaptation of Escherichia coli O157:H7 increases survival in acidic foods. Appl Environ Microbiol. 1995, 61: 3752-3755.
Bjornsdottir K, Breidt F, McFeeters RF: Protective effects of organic acids on survival of Escherichia coli O157:H7 in acidic environments. Appl Environ Microbiol. 2006, 72: 660-664.
Presser KA, Ross T, Ratkowsky DA: Function of temperature, pH, lactic acid concentration and water activity. Appl Environ Microbiol. 1998, 64: 1773-1779.
Engelsen AR, Sandlund N, Fiksdal IU, Bergh Ø: Immunohistochemistry of Atlantic cod larvae Gadus morhua experimentally challenged with Vibrio anguillarum. Dis Aquat Org. 2008, 80: 13-20.
Grisez L, Chair M, Sorgeloos P, Ollevier F: Mode of infection and spread of Vibrio anguillarum in turbot Scophthalmus maximus larvae after oral challenge through live feed. Dis Aquat Org. 1996, 26: 181-187.
Planas M, Pérez-Lorenzo M, Vázquez JA, Pintado J: A model for experimental infections with Vibrio (Listonella) anguillarum in first feeding turbot (Scophthalmus maximus L.) larvae under hatchery conditions. Aquaculture. 2005, 250: 232-243.
Balcázar JL, Blas ID, Ruiz-Zarzuela I, Cunningham D, Vendrell D, Múzquiz JL: The role of probiotics in aquaculture. Veterin Microbiol. 2006, 114: 173-186.
Gatesoupe FJ: Updating the importance of lactic acid bacteria in fish farming: Natural occurrence and probiotic treatments. J Molec Microbiol Biotechnol. 2007, 14: 107-114.
Vázquez JA, González MP, Murado MA: Effects of acid lactic bacteria cultures on pathogenic microbiota from fish. Aquaculture. 2005, 245: 149-161.
Liu P, Rand KH, Obermann B, Derendorf H: Pharmacokinetic-pharmacodynamic modelling of antibacterial activity of cefpodoxime and cefixime in in vitro kinetic models. Int J Antimicrob Agents. 2005, 25: 120-129.
Riobó P, Paz B, Franco JM, Vázquez JA, Murado MA: 2008. Proposal for a simple and sensitive haemolytic assay for palytoxin. Toxicological dynamics, kinetics, ouabain inhibition and thermal stability. Harmful Algae. 2008, 7: 415-429.
Kooijman SALM, Hanstveit AO, Oldersma H: Parametric analyses of population growth in bio-assays. Water Res. 1983, 17: 527-538.
Bowmer CT, Hooftman RN, Hanstveit AO, Venderbosch PWM, van der Hoeven N: The ecotoxicity and the biodegradability of lactic acid, alkyl lactate esters and lactate salts. Chemosphere. 1998, 37: 1317-1333.
Rial D, Vázquez JA, Murado MA: Effects of three heavy metals on the bacteria growth kinetics: a bivariate model for toxicological assessment. Appl Microbiol Biotechnol. 2011, 90: 1095-1109.
Strotmann UJ, Pagga U: A growth inhibition test with sewage bacteria -results of an international ring test 1995. Chemosphere. 1996, 32: 921-933.
Zwietering MH, Jongenburger I, Rombouts FM, Van't Riet K: Modeling of the bacterial growth curve. Appl Environ Microbiol. 1990, 56: 1875-1881.
Dalgaard P, Koutsoumanis K: Comparison of maximum specific growth rates and lag times estimated from absorbance and viable count data by different mathematical models. J Microbiol Methods. 2001, 43: 183-196.
Conner DE, Kotrola JS: Growth and survival of Escherichia coli O157:H7 under acidic conditions. Appl Environ Microbiol. 1995, 61: 382-385.
Leitch ECM, Stewart CS: Escherichia coli O157 and non-O157 isolates are more susceptible to L-lactate than to D-lactate. Appl Environ Microbiol. 2002, 68: 4676-4678.
Ray B, Sandine WE: Acetic, propionic, and lactic acids of starter culture bacteria as biopreservatives. Food preservatives of microbial origin. Edited by: Ray B, Daeschel M. 1992, 103-136. CRC Press, Boca Raton, Fla,
Bruno-Bárcena JM, Azcárate-Peril A, Hassan HH: Role of antioxidant enzymes in bacterial resistance to organic acids. Appl Environ Microbiol. 2010, 76: 2747-2753.
Alakomi HL, Skytta E, Saarela M, Mattila-Sandholm T, Latva-Kala K, Helander IM: Lactic acid permeabilizes gram-negative bacteria by disrupting the outer membrane. Appl Environ Microbiol. 2000, 66: 2001-2005.
Giotta L, Agostiano A, Italiano F, Milano F, Trotta M: Heavy metal ion influence on the photosynthetic growth of Rhodobacter sphaeroides. Chemosphere. 2006, 62: 1490-1499.
Cárdenas FC, Giannuzzi L, Zaritzky NE: Mathematical modelling of microbial growth in ground beef from Argentina. Effect of lactic acid addition, temperature and packaging film. Meat Sci. 2008, 79: 509-520.
Boziaris IS, Skandamis PN, Anastasiadi M, Nychas GJE: Effect of NaCl and KCl on fate and growth/no growth interfaces of Listeria monocytogenes Scott A at different pH and nisin concentrations. J Appl Microbiol. 2007, 102: 796-805.
Cherrington CA, Hinton M, Mead GS, Chopra I: Organic acids: chemistry, antibacterial activity and practical applications. Adv Microbiol Physiol. 1991, 32: 87-108.
Skov MN, Pedersen K, Larsen JL: Comparison of pulsed-field gel electrophoresis, ribotyping, and plasmid profiling for typing of Vibrio anguillarum serovar O1. Appl Environ Microbiol. 1995, 61: 1540-1545.
Cabo ML, Murado MA, González MP, Vázquez JA, Pastoriza L: An empirical model for describing the effects of nitrogen sources on nisin production. Lett Appl Microbiol. 2001, 33: 425-429.
Vázquez JA, González MP, Murado MA: A new marine medium. Use of the different fish peptones and comparative study of the growth of selected species of marine bacteria. Enz Microb Technol. 2004, 35: 385-392.
Cabo ML, Murado MA, González MP, Pastoriza L: A method for bacteriocin quantification. J Appl Microbiol. 1999, 87: 907-914.
Murado MA, González MP, Vázquez JA: Dose-reponse relationships. An overview a generative model and its application to the verification of descriptive models. Enz Microb Technol. 2002, 31: 439-455.
Murado MA, Vázquez JA, Rial D, Beiras R: Dose-response modeling with two agents: Application to the bioassay of oil and shoreline cleaning agents. J Hazard Mater. 2011, 185: 807-817.
Riobó P, Paz B, Franco JM, Vázquez JA, Murado MA, Cacho E: Mouse bioassay for palytoxin. Specific symptoms and dose-response against dose-death time relationships. Food Chem Toxicol. 2008, 46: 2639-2647.
Vázquez JA, Murado MA: Mathematical tools for objective comparison of microbial cultures. Application to evaluation of 15 peptones for lactic acid bacteria productions. Biochem Eng J. 2008, 39: 276-287.
Ross T: Indices for performance evaluation of predictive models in food microbiology. J Appl Bacteriol. 1996, 81: 501-508.
We thank Margarita Nogueira and Araceli Menduiña for their excellent technical assistance. Diego Rial Conde was awarded with Isabel Barreto contract and we are greateful for this financial support (Dirección Xeral de Investigación, Desenvolvemento e Innovación; Xunta de Galicia). Miguel Angel Prieto had a predoctoral contract (JAE-predoctoral CSIC-grant co-financed by FEDER funds from European Union). This study was partially funded by Xunta de Galicia (Programa de Consolidación para estructuración das unidades de investigación do sistema galego I+D+I 2008-2010, IN845B-2010/004) and Ministerio de Ciencia e Innovación (CTM2010-18411, financed with FEDER funds by European Union).
The authors declare that they have no competing interests.
JAV performed the experiments, worked on mathematical modelling of experimental data and wrote the manuscript. AD, IRA, MAP and DR helped in the analytical determinations. MAM developed the mathematical models and designed the present study. All authors read and approved the manuscript.