Towards the contribution of the p38MAPK pathway to the dual role of TGFβ in cancer: A boolean model approach

Network dynamics investigation of omics-data-driven circadian-hypoxia crosstalk logical model in gallbladder cancer reveals key therapeutic target combinations

Recent findings in cancer research have pointed towards the bidirectional interaction between circadian and hypoxia pathways. However, little is known about their crosstalk mechanism. In this work, we aimed to investigate this crosstalk at a network level utilizing the omics information of gallbladder cancer. Differential gene expression and pathway enrichment analysis were used for selecting the crucial genes from both the pathways, followed by the construction of a logical crosstalk model using GINsim. Functional circuit identification and node perturbations were then performed. Significant node combinations were used to investigate the temporal behavior of the network through MaBoSS. Lastly, the model was validated using published in vitro experimentations. Four new positive circuits and a new axis viz. BMAL1/ HIF1αβ/ NANOG, responsible for stemness were identified. Through triple node perturbations viz.a. BMAL:CLOCK (KO or E1) + P53 (E1) + HIF1α (KO); b. P53 (E1) + HIF1α (KO) + MYC (E1); and c. HIF1α (KO) + MYC (E1) + EGFR (KO), the model was able to inhibit cancer growth and maintain a homeostatic condition. This work provides an architecture for drug simulation analysis to entrainment circadian rhythm and in vitro experiments for chronotherapy-related studies. Insight Box. Circadian rhythm and hypoxia are the key dysregulated processes which fuels-up the cancer growth. In the present work we have developed a gallbladder cancer (GBC) specific Boolean model, utilizing the RNASeq data from GBC dataset and tissue specific interactions. This work adequately models the bidirectional nature of interactions previously illustrated in experimental papers showing the effect of hypoxia on dysregulation of circadian rhythm and the influence of this disruption on progression towards metastasis. Through the dynamical study of the model and its response to different perturbations, we report novel triple node combinations that can be targeted to efficiently reduce GBC growth. This network can be used as a generalized framework to investigate different crosstalk pathways linked with cancer progression.

Common Attractors in Multiple Boolean Networks

Analyzing multiple networks is important to understand relevant features among different networks. Although many studies have been conducted for that purpose, not much attention has been paid to the analysis of attractors (i.e., steady states) in multiple networks. Therefore, we study common attractors and similar attractors in multiple networks to uncover hidden similarities and differences among networks using Boolean networks (BNs), where BNs have been used as a mathematical model of genetic networks and neural networks. We define three problems on detecting common attractors and similar attractors, and theoretically analyze the expected number of such objects for random BNs, where we assume that given networks have the same set of nodes (i.e., genes). We also present four methods for solving these problems. Computational experiments on randomly generated BNs are performed to demonstrate the efficiency of our proposed methods. In addition, experiments on a practical biological system, a BN model of the TGF- β signaling pathway, are performed. The result suggests that common attractors and similar attractors are useful for exploring tumor heterogeneity and homogeneity in eight cancers.

SNAIL driven by a feed forward loop motif promotes TGFβinduced epithelial to mesenchymal transition

Epithelial to Mesenchymal Transition (EMT) plays an important role in tissue regeneration, embryonic development, and cancer metastasis. Several signaling pathways are known to regulate EMT, among which the modulation of TGFβ(Transforming Growth Factor-β) induced EMT is crucial in several cancer types. Several mathematical models were built to explore the role of core regulatory circuit of ZEB/miR-200, SNAIL/miR-34 double negative feedback loops in modulating TGFβinduced EMT. Different emergent behavior including tristability, irreversible switching, existence of hybrid EMT states were inferred though these models. Some studies have explored the role of TGFβreceptor activation, SMADs nucleocytoplasmic shuttling and complex formation. Recent experiments have revealed that MDM2 along with SMAD complex regulates SNAIL expression driven EMT. Encouraged by this, in the present study we developed a mathematical model for p53/MDM2 dependent TGFβinduced EMT regulation. Inclusion of p53 brings in an additional mechanistic perspective in exploring the EM transition. The network formulated comprises a C1FFL moderating SNAIL expression involving MDM2 and SMAD complex, which functions as a noise filter and persistent detector. The C1FFL was also observed to operate as a coincidence detector driving the SNAIL dependent downstream signaling into phenotypic switching decision. Systems modelling and analysis of the devised network, displayed interesting dynamic behavior, systems response to various inputs stimulus, providing a better understanding of p53/MDM2 dependent TGF-βinduced Epithelial to Mesenchymal Transition.

In Silico Logical Modelling to Uncover Cooperative Interactions in Cancer

The multistep development of cancer involves the cooperation between multiple molecular lesions, as well as complex interactions between cancer cells and the surrounding tumour microenvironment. The search for these synergistic interactions using experimental models made tremendous contributions to our understanding of oncogenesis. Yet, these approaches remain labour-intensive and challenging. To tackle such a hurdle, an integrative, multidisciplinary effort is required. In this article, we highlight the use of logical computational models, combined with experimental validations, as an effective approach to identify cooperative mechanisms and therapeutic strategies in the context of cancer biology. In silico models overcome limitations of reductionist approaches by capturing tumour complexity and by generating powerful testable hypotheses. We review representative examples of logical models reported in the literature and their validation. We then provide further analyses of our logical model of Epithelium to Mesenchymal Transition (EMT), searching for additional cooperative interactions involving inputs from the tumour microenvironment and gain of function mutations in NOTCH.

Systems biology approach suggests new miRNAs as phenotypic stability factors in the epithelial-mesenchymal transition

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.