A cancer sub-population competition model reveals optimal levels of immune response that minimize tumor size

Breast cancer is a complex disease with significant phenotypic heterogeneity of cells, even within a single breast tumor. Emerging evidence underscores the significance of intratumoral competition, which can serve as a key contributor to cancer drug resistance, imparting substantial clinical implications. Understanding the competitive dynamics is paramount as it can significantly influence disease progression and treatment outcomes. In the present work, a mathematical model was developed using a system of differential equations to describe the dynamic interactions between two cancer subtypes (each further classified into cancer stem cells and tumor cells) and innate immune cells. The purpose of the model is to comprehensively understand the competitive interactions between the heterogeneous subpopulations. The equilibrium points and stability analysis for each equilibrium point were established. Model simulations showed that the competition between two cancer subtypes directly affects the number of both species. When competition between two cancer subtypes is strong, increasing the immune response rate specific to themore competitive species effectively reduces the tumor size. However, if the competition is relatively weak, an optimal immune response rate is required to minimize the total number of tumor cells. Rates below the optimal level fail to reduce the population of the stronger species, whereas rates above the optimal level can lead to the recurrence of the weaker species. Overall, this model provides insights into breast cancer dynamics and guides the development of effective treatment strategies.