Mathematical Modeling and Dynamical Analysis for Tumor Cells and Tumor Propagating Cells Controlled by G9a Inhibitors

It is of critical importance to comprehend the biological environment and core tumor populations when trying to design successful therapeutic solutions for fighting cancers. In several diseases, G9a has been recognized as a novel epigenetic therapeutic target, and its blockage can shift tumor cells (TCs) toward tumor propagating cells (TPCs). This study combines mathematical modeling based on ordinary differential equations and dynamical analysis to quantitatively and qualitatively understand the interactions among G9a, TCs, and TPCs, denoted as G9a-TC-TPC. We propose four different dynamical systems with the impact of the strong Allee effect, named the Hill-Hill system, Logistic-Logistic system, Hill-Logistic system and Logistic-Hill system, to simulate different biological processes through the Hill functions and the Logistic functions that are often used in the models of biological systems. Based on theoretical analysis of these models, including the positivity, boundedness and stability of equilibria, we find that the Hill-Logistic system can display bistable states that correspond to the wild-type tumors and the aggressive tumors. Consequently, we use bifurcation analysis and numerical simulations to illustrate the complicated dynamical behavior of this system. It has been shown that under a specific therapy that changes the relative apoptotic rate of TCs (G9a suppresses the apoptosis of TCs), which can affect the bistability and instability of the system, the wide-type state can be obtained. We also discover that the relative handling time of TCs and TPCs can cooperatively enhance bistability, whereas the cooperative coefficient of feedback can contribute to all tumor cells moving from high-level monostability to bistable states in a restricted region, then to low-level stable states. These results offer new insights for more precisely understanding epigenetic therapy treatments with G9a.