The main aim of the current work was to build up a new mathematical modelling approach in predictive food microbiology field for the prediction of growth kinetics of microorganisms. For this purpose, the bacterial growth data of Pseudomonas spp. in whole fish (gilt-head seabream) subjected to isothermal and non-isothermal storage temperatures were collected from previously published growth curves. Maximum specific growth rate (1/h) and lag phase duration (h) were described as a function of storage temperature using the direct two-step, direct one-step and inverse dynamic modelling approaches based on various meta-heuristic optimization algorithms. The fitting capability of the modelling approaches and employed optimization algorithms was separately compared, and the one-step modelling approach for the direct methods and the Bayesian optimization method for the used algorithms provided the best goodness of fit results. These two were then further processed in validation step. The inverse dynamic modelling approach based on the Bayesian optimization algorithm yielded satisfactorily statistical indexes (1.02 > Bias factor > 1.09 and 1.07 > Accuracy factor > 1.13), which indicates it can be reliably used as an alternative way of describing the growth behaviour of Pseudomonas spp. in fish in a fast and efficient manner with minimum labour effort.