The epithelial–mesenchymal transition (EMT) is a cellular programme on which epithelial cells undergo a phenotypic transition to mesenchymal ones acquiring metastatic properties such as mobility and invasion. TGF-β signalling can promote the EMT process. However, the dynamics of the concentration response of TGF-β-induced EMT is not well explained. In this work, we propose a logical model of TGF-β dose dependence of EMT in MCF10A breast cells. The model outcomes agree with experimentally observed phenotypes for the wild-type and perturbed/mutated cases. As important findings of the model, it predicts the coexistence of more than one hybrid state and that the circuit between TWIST1 and miR-129 is involved in their stabilization. Thus, miR-129 should be considered as a phenotypic stability factor. The circuit TWIST1/miR-129 associates with ZEB1-mediated circuits involving miRNAs 200, 1199, 340, and the protein GRHL2 to stabilize the hybrid state. Additionally, we found a possible new autocrine mechanism composed of the circuit involving TGF-β, miR-200, and SNAIL1 that contributes to the stabilization of the mesenchymal state. Therefore, our work can extend our comprehension of TGF-β-induced EMT in MCF10A cells to potentially improve the strategies for breast cancer treatment.

It is well established that microRNA-21 (miR-21) targets phosphatase and tensin homolog (PTEN), facilitating epithelial-to-mesenchymal transition (EMT) and drug resistance in cancer. Recent evidence indicates that PTEN activates its pseudogene-derived long non-coding RNA, PTENP1, which in turn inhibits miR-21. However, the dynamics of PTEN, miR-21, and PTENP1 in the DNA damage response (DDR) remain unclear. Thus, we propose a dynamic Boolean network model by integrating the published literature from various cancers. Our model shows good agreement with the experimental findings from breast cancer, hepatocellular carcinoma (HCC), and oral squamous cell carcinoma (OSCC), elucidating how DDR activation transitions from the intra-S phase to the G2 checkpoint, leading to a cascade of cellular responses such as cell cycle arrest, senescence, autophagy, apoptosis, drug resistance, and EMT. Model validation underscores the roles of PTENP1, miR-21, and PTEN in modulating EMT and drug resistance. Furthermore, our analysis reveals nine novel feedback loops, eight positive and one negative, mediated by PTEN and implicated in DDR cell fate determination, including pathways related to drug resistance and EMT. Our work presents a comprehensive framework for investigating cellular responses following DDR, underscoring the therapeutic potential of targeting PTEN, miR-21, and PTENP1 in cancer treatment.

Recent studies showed that induced microRNA-449a (miR-449a) enhances a G2/M cell cycle checkpoint arrest in prostate cancer (LNCaP) and lung adenocarcinoma cell lines. In the case of LNCaP cells, upregulated miR-449a directly downregulates c-Myc that is required to induce the cell cycle regulators Cdc25A and Cdc2/CyclinB whose inactivation blocks G2 to M phase transition. However, the molecular mechanisms involved are yet unclear, although in other prostate cancer cells the interactions among p53, miR-449a and Sirt-1 can affect the induction of the G2/M arrest. In order to clarify these molecular mechanisms, in this work we propose a boolean model of the G2/M checkpoint arrest regulation contemplating the influence of miR-449a. The model shows that the cell fate determination between two cellular phenotypes: G2/M-Arrest for DNA repair and G2/M-induced apoptosis is stochastic and influenced by miR-449a state of activation. The results were compared with experimental data available presenting agreement. We also found that several feedback loops are involved in this cell fate regulation and we indicate, through in silico gain or loss of function perturbations of genes, which of these feedback loops are more efficient to favor a specific phenotype.

The tumor suppressor p53, in its wild-type form, plays a central role in cellular homeostasis by regulating senescence, apoptosis, and autophagy within the DNA damage response (DDR). Recent findings suggest that wild-type p53 also governs ferroptosis, an iron-dependent cell death process driven by lipid peroxidation. Post-translational modifications of p53 generate proteoforms that significantly enhance its functional diversity in regulating these mechanisms. A key target in this process is the cystine/glutamate transporter (xCT), which is essential for redox balance and ferroptosis resistance. Additionally, p53-induced miR-34c-5p suppresses cancer cell proliferation and drug resistance by modulating Myc, an oncogene further influenced by non-coding RNAs like circular RNA NOTCH1 (CricNOTCH1) and long non-coding RNA MALAT1. However, the exact role of these molecules in ferroptosis remains unclear. To address this, we introduce the first dynamic Boolean model that delineates the influence of these ncRNAs and p53 on ferroptosis, apoptosis, and senescence within the DDR context. Validated through gain- and loss-of-function perturbations, our model closely aligns with experimental observations in cancers such as oral squamous cell carcinoma, nasopharyngeal carcinoma, and osteosarcoma. The model identifies crucial positive feedback loops (CricNOTCH1/miR-34c/Myc, MALAT1/miR-34c/Myc, and Myc/xCT) and highlights the therapeutic potential of using p53 proteoforms and ncRNAs to combat drug resistance and induce cancer cell death.

Cell fate determination is a complex process that is frequently described as cells traveling on rugged pathways, beginning with DNA damage response (DDR). Tumor protein p53 (p53) and phosphatase and tensin homolog (PTEN) are two critical players in this process. Although both of these proteins are known to be key cell fate regulators, the exact mechanism by which they collaborate in the DDR remains unknown. Thus, we propose a dynamic Boolean network. Our model incorporates experimental data obtained from NSCLC cells and is the first of its kind. Our network’s wild-type system shows that DDR activates the G2/M checkpoint, and this triggers a cascade of events, involving p53 and PTEN, that ultimately lead to the four potential phenotypes: cell cycle arrest, senescence, autophagy, and apoptosis (quadra-stable dynamics). The network predictions correspond with the gain-and-loss of function investigations in the additional two cell lines (HeLa and MCF-7). Our findings imply that p53 and PTEN act as molecular switches that activate or deactivate specific pathways to govern cell fate decisions. Thus, our network facilitates the direct investigation of quadruplicate cell fate decisions in DDR. Therefore, we concluded that concurrently controlling PTEN and p53 dynamics may be a viable strategy for enhancing clinical outcomes.

The long non-coding RNA metastasis-associated lung adenocarcinoma transcript 1 (lncRNA MALAT1) and microRNA-145-5p (miR-145) axis play a pivotal role in regulating drug resistance, apoptosis and senescence in non-small cell lung cancer (NSCLC). MALAT1 drives drug resistance by suppressing miR-145 and activating MUC1, thereby inhibiting ferroptosis; however, its precise role in regulating ferroptosis in NSCLC remains unclear. Therefore, we propose a computational modelling approach to unravel the impact of the MALAT1/miR-145 axis on ferroptosis and drug resistance, to identify potential therapeutic strategies that promote ferroptosis. Using Boolean logic and a stochastic updating scheme, we developed and validated a robust regulatory model that encompasses ferroptosis, apoptosis, senescence and drug resistance pathways. The model, based on extensive literature and validated through gain- and loss-of-function perturbations, demonstrated strong alignment with observed clinical data that were not included in its construction. Our analysis identified three previously unreported feedback loops, miR-145/Wip1/p53, miR-145/Myc/MALAT1 and miR-145/MUC1/BMI1, establishing miR-145 as a central regulator in NSCLC. Perturbations targeting MALAT1 and wild-type p53-induced phosphatase 1 (Wip1) revealed potential therapeutic opportunities, with miR-145 activation emerging as a promising strategy to induce ferroptosis and overcome drug resistance. These findings highlight the MALAT1/miR-145 axis as a transformative therapeutic target, presenting a computational foundation to advance NSCLC treatment strategies.