Bayesian neural networks, decision boundary, explainable artificial intelligence, machine learning, model uncertainty, posterior predictive check, prior probability, probabilistic programming, scatter plot, variational posterior distribution

The emergence and growth of drug-resistant cancer cell subpopulations during anti-cancer treatment is a major challenge for cancer therapies. Combination therapies are usually applied for overcoming drug resistance. In the present paper, we explored the evolution outcome of tumor cell populations under different combination schedules of chemotherapy and p53 vaccine, by construction of replicator dynamical model for sensitive cells, chemotherapy-resistant cells and p53 vaccine-resistant cells. The local asymptotic stability analysis of the evolutionary stable points revealed that cancer population could evolve to the population with single subpopulation, or coexistence of sensitive cells and p53 vaccine-resistant cells, or coexistence of chemotherapy-resistant cells and p53 vaccine-resistant cells under different monotherapy or combination schedules. The design of adaptive therapy schedules that maintain the subpopulations under control is also demonstrated by sequential and periodic application of combination treatment strategies based on the evolutionary velocity and evolutionary absorbing regions. Applying a new replicator dynamical model, we further explored the supportive effects of sensitive cancer cells on targeted therapy-resistant cells revealed in mice experiments. It was shown that the supportive effects of sensitive cells could drive the evolution of cell population from sensitive cells to coexistence of sensitive cells and one type of targeted therapy-resistant cells.