Cancer immunotherapy is a rapidly developing field, with numerous drugs and 15 therapy combinations waiting to be tested in pre-clinical and clinical settings. However, the costly and time-consuming trial-and-error approach to development of new treatment paradigms creates a research bottleneck, motivating the development of complementary approaches. Computational modelling is a compelling candidate for this task, however, difficulties associated with the validation of such models limits their use in pre-clinical 20 and clinical settings. Here we propose a bottom-up deterministic computational model to simulate tumour response to treatment with anti-programmed-death-1 antibodies (anti-PD-1). The model was built with validation in mind, and so contains minimum number of parameters, and only 4 free parameters. Moreover, all model parameters can be measured experimentally. Free parameters were tuned by fitting the model to 25 experimental data from the literature, using B16-F10 murine melanoma implanted into wild type (C57BL/6), as well as into immunodeficient (NSG) mice strains, and treated with anti-PD-1 antibodies. The model's predictive ability was verified on two independent datasets from literature with different but well-known inputs. To identify possible biomarkers of response to anti-PD-1 immunotherapy, sensitivity study of key model 30 parameters was performed. Good agreement between the simulated tumour growth curves and the experimental data was achieved, with mean relative deviations in the range of 13