Abstract
Quantitative modeling of microscopic genes regulatory mechanisms in an individual cell is a crucial step towards understanding various macroscopic physiological phenomena of cell populations. Based on the regulatory mechanisms of genes zeb1 and cdh1 in the growth and development of breast cancer cells, we propose a kinetic model at the level of single cell. By constructing the effective landscape of underlying stationary probability for the genes expressions, it is found that (i) each breast cancer cell has three phenotypic states (i.e., the stem-like, basal, and luminal states) which correspond to three attractions of the probability landscape. (ii) The interconversions between phenotypic states can be induced by the noise intensity and the property of phenotypic switching is quantified by the mean first-passage time. (iii) Under certain conditions, the probabilities of each cancer cell appearing in the three states are consistent with the macroscopic phenotypic equilibrium proportions in the breast cancer SUM159 cell line. (iv) Our kinetic model involving the TGF-β signal can also qualitatively explain several macroscopic physiological phenomena of breast cancer cells, such as the "TGF-β paradox" in tumor therapy, the five clinical subtypes of breast cancer cells, and the effects of transient TGF-β on breast cancer metastasis.