In this study, we developed kinetic models to predict the growth of pathogenic Escherichia coli on cheeses during storage at constant and changing temperatures. A five-strain mixture of pathogenic E. coli was inoculated onto natural cheeses (Brie and Camembert) and processed cheeses (sliced Mozzarella and sliced Cheddar) at 3 to 4 log CFU/g. The inoculated cheeses were stored at 4, 10, 15, 25, and $30^{circ}C$ for 1 to 320 h, with a different storage time being used for each temperature. Total bacteria and E. coli cells were enumerated on tryptic soy agar and MacConkey sorbitol agar, respectively. E. coli growth data were fitted to the Baranyi model to calculate the maximum specific growth rate (${mu}_{max}$; log CFU/g/h), lag phase duration (LPD; h), lower asymptote (log CFU/g), and upper asymptote (log CFU/g). The kinetic parameters were then analyzed as a function of storage temperature, using the square root model, polynomial equation, and linear equation. A dynamic model was also developed for varying temperature. The model performance was evaluated against observed data, and the root mean square error (RMSE) was calculated. At $4^{circ}C$, E. coli cell growth was not observed on any cheese. However, E. coli growth was observed at $10{circ}C$ to $30^{circ}C$C with a ${mu}_{max}$ of 0.01 to 1.03 log CFU/g/h, depending on the cheese. The ${mu}_{max}$ values increased as temperature increased, while LPD values decreased, and ${mu}_{max}$ and LPD values were different among the four types of cheese. The developed models showed adequate performance (RMSE = 0.176-0.337), indicating that these models should be useful for describing the growth kinetics of E. coli on various cheeses.