A mathematical-computational model explains phenotypic plasticity changes by mutations underlying hepatic cancer stem cells emergence

In several carcinomas, including hepatocellular carcinoma, it has been demonstrated that cancer stem cells (CSCs) have enhanced invasiveness and therapy resistance compared to differentiated cancer cells. Mathematical-computational tools could be valuable for integrating experimental results and understanding the phenotypic plasticity mechanisms for CSCs emergence. Based on the literature review, we constructed a Boolean model that recovers eight stable states (attractors) corresponding to the gene expression profile of hepatocytes and mesenchymal cells in senescent, quiescent, proliferative, and stem-like states. The epigenetic landscape associated with the regulatory network was analyzed. We observed that the loss of p53, p16, RB, or the constitutive activation of β-catenin and YAP1 increases the robustness of the proliferative stem-like phenotypes. Additionally, we found that p53 inactivation facilitates the transition of proliferative hepatocytes into stem-like mesenchymal phenotype. Thus, phenotypic plasticity may be altered, and stem-like phenotypes related to CSCs may be easier to attain following the mutation acquisition.