This paper describes a simulation study which aims at optimizing the therapy for the control of Chronic Myeloid Leukemia according to the following objectives: the reduction of the administered drug and vaccine amounts, the establishment of a auto-immune response and the long-term control of disease without reducing the effective of therapy with respect to the full treatment. A therapy optimization method is developed defining and solving a Model Predictive Control algorithm, preceded by an accurate Initial Guess search based on Monte-Carlo like approach. Simulation results show that the suggested procedure achieves the proposed goals.