Prostate cancer (PCa) is a major world-wide health concern. Current diagnostic methods involve Prostate-Specific Antigen (PSA) blood tests, biopsies, and Magnetic Resonance Imaging (MRI) to assess cancer aggressiveness and guide treatment decisions. MRI aligns with in silico medicine, as patient-specific image biomarkers can be obtained, contributing towards the development of digital twins for clinical practice. This work presents a novel framework to create a personalized PCa model by integrating clinical MRI data, such as the prostate and tumour geometry, the initial distribution of cells and the vasculature, so a full representation of the whole prostate is obtained. On top of the personalized model construction, our approach simulates and predicts temporal tumour growth in the prostate through the Finite Element Method, coupling the dynamics of tumour growth and the transport of oxygen, and incorporating cellular processes such as proliferation, differentiation, and apoptosis. In addition, our approach includes the simulation of the PSA dynamics, which allows to evaluate tumour growth through the PSA patient’s levels. To obtain the model parameters, a multi-objective optimization process is performed to adjust the best parameters for two patients simultaneously. This framework is validated by means of data from four patients with several MRI follow-ups. The diagnosis MRI allows the model creation and initialization, while subsequent MRI-based data provide additional information to validate computational predictions. The model predicts prostate and tumour volumes growth, along with serum PSA levels. This work represents a preliminary step towards the creation of digital twins for PCa patients, providing personalized insights into tumour growth.