Construction and Experimental Validation of a Petri Net Model of Wnt/β-Catenin Signaling

The Axin scaffold protects the kinase GSK3β from cross-pathway inhibition

Multiple signaling pathways regulate the kinase GSK3β by inhibitory phosphorylation at Ser9, which then occupies the GSK3β priming pocket and blocks substrate binding. Since this mechanism should affect GSK3β activity toward all primed substrates, it is unclear why Ser9 phosphorylation does not affect other GSK3β-dependent pathways, such as Wnt signaling. We used biochemical reconstitution and cell culture assays to evaluate how Wnt-associated GSK3β is insulated from cross-activation by other signals. We found that the Wnt-specific scaffold protein Axin allosterically protects GSK3β from phosphorylation at Ser9 by upstream kinases, which prevents accumulation of pS9-GSK3β in the Axin•GSK3β complex. Scaffold proteins that protect bound proteins from alternative pathway reactions could provide a general mechanism to insulate signaling pathways from improper crosstalk.

Increased Expression of LEF1 and β-Catenin in Invasive Micropapillary Carcinoma of the Breast is Associated With Lymphovascular Invasion and Lymph Node Metastasis

Invasive micropapillary breast carcinoma (IMPC) is a rare breast cancer subtype characterized by small tumor cell clusters with loss of stromal attachment, an inside-out growth appearance, and lymphotropism. IMPC is associated with high incidence of lymphovascular invasion (LVI) and lymph node metastasis. Activated Wnt signaling has been implicated in the metastasis of other aggressive breast tumors, including triple-negative and basal-like carcinomas. In this study, we tested whether activated Wnt signaling could be detected in IMPC. Upon ligand binding, the central mediator of the Wnt pathway, β-catenin, accumulates in the cytosol and translocates to the nucleus where it forms a complex with lymphoid enhancer-binding factor 1 (LEF1) to regulate transcription. We performed immunostaining for β-catenin and LEF1 on a well-annotated cohort of 40 breast tumors and nodal metastases displaying micropapillary histopathology. Strong nuclear accumulation of β-catenin was not observed, however a dim cytosolic and/or nuclear accumulation of β-catenin was sometimes seen in IMPC and this expression pattern was significantly associated with nodal metastasis. β-catenin expression correlated with the upregulation of LEF1 in IMPC. LEF1 expression was detected in 26 of 40 (65%) cases and was specifically enriched at the invasive front of the tumor and in tumor clusters undergoing LVI. Detection of LEF1 expression in the primary tumor was associated with an increased rate of LVI, lymph node metastasis, and disease relapse. LEF1 and β-catenin expression levels were significantly higher in metastases compared with primary tumors. In summary, this study demonstrates an association between the upregulation of β-catenin/LEF1 and the metastasis of IMPC.

Quantitative live-cell imaging and computational modeling shed new light on endogenous WNT/CTNNB1 signaling dynamics

WNT/CTNNB1 signaling regulates tissue development and homeostasis in all multicellular animals, but the underlying molecular mechanism remains incompletely understood. Specifically, quantitative insight into endogenous protein behavior is missing. Here, we combine CRISPR/Cas9-mediated genome editing and quantitative live-cell microscopy to measure the dynamics, diffusion characteristics and absolute concentrations of fluorescently tagged, endogenous CTNNB1 in human cells under both physiological and oncogenic conditions. State-of-the-art imaging reveals that a substantial fraction of CTNNB1 resides in slow-diffusing cytoplasmic complexes, irrespective of the activation status of the pathway. This cytoplasmic CTNNB1 complex undergoes a major reduction in size when WNT/CTNNB1 is (hyper)activated. Based on our biophysical measurements, we build a computational model of WNT/CTNNB1 signaling. Our integrated experimental and computational approach reveals that WNT pathway activation regulates the dynamic distribution of free and complexed CTNNB1 across different subcellular compartments through three regulatory nodes: the destruction complex, nucleocytoplasmic shuttling, and nuclear retention.

The Scaffold Protein Axin Promotes Signaling Specificity within the Wnt Pathway by Suppressing Competing Kinase Reactions

Scaffold proteins are thought to promote signaling specificity by accelerating reactions between bound kinase and substrate proteins. To test the long-standing hypothesis that the scaffold protein Axin accelerates glycogen synthase kinase 3β (GSK3β)-mediated phosphorylation of β-catenin in the Wnt signaling network, we measured GSK3β reaction rates with multiple substrates in a minimal, biochemically reconstituted system. We observed an unexpectedly small, ∼2-fold Axin-mediated rate increase for the β-catenin reaction when measured in isolation. In contrast, when both β-catenin and non-Wnt pathway substrates are present, Axin accelerates the β-catenin reaction by preventing competition with alternative substrates. At high competitor concentrations, Axin produces >10-fold rate effects. Thus, while Axin alone does not markedly accelerate the β-catenin reaction, in physiological settings where multiple GSK3β substrates are present, Axin may promote signaling specificity by suppressing interactions with competing, non-Wnt pathway targets. This mechanism for scaffold-mediated control of competition enables a shared kinase to perform distinct functions in multiple signaling networks.

A framework for exhaustive modelling of genetic interaction patterns using Petri nets

MOTIVATION

Genetic interaction (GI) patterns are characterized by the phenotypes of interacting single and double mutated gene pairs. Uncovering the regulatory mechanisms of GIs would provide a better understanding of their role in biological processes, diseases and drug response. Computational analyses can provide insights into the underpinning mechanisms of GIs.

RESULTS

In this study, we present a framework for exhaustive modelling of GI patterns using Petri nets (PN). Four-node models were defined and generated on three levels with restrictions, to enable an exhaustive approach. Simulations suggest ∼5 million models of GIs. Generalizing these we propose putative mechanisms for the GI patterns, inversion and suppression. We demonstrate that exhaustive PN modelling enables reasoning about mechanisms of GIs when only the phenotypes of gene pairs are known. The framework can be applied to other GI or genetic regulatory datasets.

AVAILABILITY AND IMPLEMENTATION

The framework is available at http://www.ibi.vu.nl/programs/ExhMod.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Aberrant WNT/CTNNB1 Signaling as a Therapeutic Target in Human Breast Cancer: Weighing the Evidence

WNT signaling is crucial for tissue morphogenesis during development in all multicellular animals. After birth, WNT/CTNNB1 responsive stem cells are responsible for tissue homeostasis in various organs and hyperactive WNT/CTNNB1 signaling is observed in many different human cancers. The first link between WNT signaling and breast cancer was established almost 40 years ago, when Wnt1 was identified as a proto-oncogene capable of driving mammary tumor formation in mice. Since that discovery, there has been a dedicated search for aberrant WNT signaling in human breast cancer. However, much debate and controversy persist regarding the importance of WNT signaling for the initiation, progression or maintenance of different breast cancer subtypes. As the first drugs designed to block functional WNT signaling have entered clinical trials, many questions about the role of aberrant WNT signaling in human breast cancer remain. Here, we discuss three major research gaps in this area. First, we still lack a basic understanding of the function of WNT signaling in normal human breast development and physiology. Second, the overall extent and precise effect of (epi)genetic changes affecting the WNT pathway in different breast cancer subtypes are still unknown. Which underlying molecular and cell biological mechanisms are disrupted as a result also awaits further scrutiny. Third, we survey the current status of targeted therapeutics that are aimed at interfering with the WNT pathway in breast cancer patients and highlight the importance and complexity of selecting the subset of patients that may benefit from treatment.

The ability of transcription factors to differentially regulate gene expression is a crucial component of the mechanism underlying inversion, a frequently observed genetic interaction pattern

Genetic interactions, a phenomenon whereby combinations of mutations lead to unexpected effects, reflect how cellular processes are wired and play an important role in complex genetic diseases. Understanding the molecular basis of genetic interactions is crucial for deciphering pathway organization as well as understanding the relationship between genetic variation and disease. Several hypothetical molecular mechanisms have been linked to different genetic interaction types. However, differences in genetic interaction patterns and their underlying mechanisms have not yet been compared systematically between different functional gene classes. Here, differences in the occurrence and types of genetic interactions are compared for two classes, gene-specific transcription factors (GSTFs) and signaling genes (kinases and phosphatases). Genome-wide gene expression data for 63 single and double deletion mutants in baker's yeast reveals that the two most common genetic interaction patterns are buffering and inversion. Buffering is typically associated with redundancy and is well understood. In inversion, genes show opposite behavior in the double mutant compared to the corresponding single mutants. The underlying mechanism is poorly understood. Although both classes show buffering and inversion patterns, the prevalence of inversion is much stronger in GSTFs. To decipher potential mechanisms, a Petri Net modeling approach was employed, where genes are represented as nodes and relationships between genes as edges. This allowed over 9 million possible three and four node models to be exhaustively enumerated. The models show that a quantitative difference in interaction strength is a strict requirement for obtaining inversion. In addition, this difference is frequently accompanied with a second gene that shows buffering. Taken together, these results provide a mechanistic explanation for inversion. Furthermore, the ability of transcription factors to differentially regulate expression of their targets provides a likely explanation why inversion is more prevalent for GSTFs compared to kinases and phosphatases.

Aurora kinase A (AURKA) interaction with Wnt and Ras-MAPK signalling pathways in colorectal cancer

Hyperactivation of Wnt and Ras-MAPK signalling are common events in development of colorectal adenomas. Further progression from adenoma-to-carcinoma is frequently associated with 20q gain and overexpression of Aurora kinase A (AURKA). Interestingly, AURKA has been shown to further enhance Wnt and Ras-MAPK signalling. However, the molecular details of these interactions in driving colorectal carcinogenesis remain poorly understood. Here we first performed differential expression analysis (DEA) of AURKA knockdown in two colorectal cancer (CRC) cell lines with 20q gain and AURKA overexpression. Next, using an exact algorithm, Heinz, we computed the largest connected protein-protein interaction (PPI) network module of significantly deregulated genes in the two CRC cell lines. The DEA and the Heinz analyses suggest 20 Wnt and Ras-MAPK signalling genes being deregulated by AURKA, whereof β-catenin and KRAS occurred in both cell lines. Finally, shortest path analysis over the PPI network revealed eight 'connecting genes' between AURKA and these Wnt and Ras-MAPK signalling genes, of which UBE2D1, DICER1, CDK6 and RACGAP1 occurred in both cell lines. This study, first, confirms that AURKA influences deregulation of Wnt and Ras-MAPK signalling genes, and second, suggests mechanisms in CRC cell lines describing these interactions.

MicroRNA pharmacogenomics based integrated model of miR-17-92 cluster in sorafenib resistant HCC cells reveals a strategy to forestall drug resistance

Among solid tumors, hepatocellular carcinoma (HCC) emerges as a prototypical therapy-resistant tumor. Considering the emerging sorafenib resistance crisis in HCC, future studies are urgently required to overcome resistance. Recently noncoding RNAs (ncRNAs) have emerged as significant regulators in signalling pathways involved in cancer drug resistance and pharmacologically targeting these ncRNAs might be a novel stratagem to reverse drug resistance. In the current study, using a hybrid Petri net based computational model, we have investigated the harmonious effect of miR-17-92 cluster inhibitors/mimics and circular RNAs on sorafenib resistant HCC cells in order to explore potential resistance mechanisms and to identify putative targets for sorafenib-resistant HCC cells. An integrated model was developed that incorporates seven miRNAs belonging to miR-17-92 cluster (hsa-miR-17-5p, hsa-miR-17-3p, hsa-miR-19a, hsa-miR-19b, hsa-miR-18a, hsa-miR-20a and hsa-miR-92) and crosstalk of two signaling pathways (EGFR and IL-6) that are differentially regulated by these miRNAs. The mechanistic connection was proposed by the correlation between members belonging to miR-17-92 cluster and corresponding changes in the protein levels of their targets in HCC, specifically those targets that have verified importance in sorafenib resistance. Current findings uncovered potential pathway features, underlining the significance of developing modulators of this cluster to combat drug resistance in HCC.

SATB2/β-catenin/TCF-LEF pathway induces cellular transformation by generating cancer stem cells in colorectal cancer

Recent studies have demonstrated the involvement of colorectal cancer (CRC) stem cells (CSC) in transformation, cancer progression and metastasis. The main goal of this paper was to examine the molecular mechanisms by which SATB2 induced malignant transformation of colorectal epithelial cells. SATB2 induced malignant transformation and these transformed cells gained the characteristics of CSCs by expressing stem cell markers (CD44, CD133, LGR5 and DCLK1) and transcription factors (c-Myc, Nanog and Sox2). Overexpression of SATB2 in normal colorectal epithelial cells increased cell motility, migration and invasion, which were associated with an increase in N-cadherin and Zeb1, and decrease in E-cadherin expression. SATB2 overexpression also upregulated XIAP and cyclin D1, suggesting its role in cell survival and cell cycle. Furthermore, the expression of SATB2 was positively correlated with β-catenin expression in CRC. In contrary, depletion of SATB2 inhibited cell proliferation, colony formation, cell motility and expression of β-catenin, Snail, Slug, Zeb1 and N-cadherin, and upregulated E-cadherin. Furthermore, SATB2 silencing inhibited the expression of stem cell markers, pluripotency maintaining transcription factors, cell cycle and cell proliferation/survival genes and TCF/LEF targets. Finally, β-catenin/TCF-LEF pathway mediated the biological effects of SATB2 in CSCs. These studies support the role of SATB2/β-catenin/TCF-LEF pathway in transformation and carcinogenesis.

SATB2/β-catenin/TCF-LEF pathway induces cellular transformation by generating cancer stem cells in colorectal cancer

Recent studies have demonstrated the involvement of colorectal cancer (CRC) stem cells (CSC) in transformation, cancer progression and metastasis. The main goal of this paper was to examine the molecular mechanisms by which SATB2 induced malignant transformation of colorectal epithelial cells. SATB2 induced malignant transformation and these transformed cells gained the characteristics of CSCs by expressing stem cell markers (CD44, CD133, LGR5 and DCLK1) and transcription factors (c-Myc, Nanog and Sox2). Overexpression of SATB2 in normal colorectal epithelial cells increased cell motility, migration and invasion, which were associated with an increase in N-cadherin and Zeb1, and decrease in E-cadherin expression. SATB2 overexpression also upregulated XIAP and cyclin D1, suggesting its role in cell survival and cell cycle. Furthermore, the expression of SATB2 was positively correlated with β-catenin expression in CRC. In contrary, depletion of SATB2 inhibited cell proliferation, colony formation, cell motility and expression of β-catenin, Snail, Slug, Zeb1 and N-cadherin, and upregulated E-cadherin. Furthermore, SATB2 silencing inhibited the expression of stem cell markers, pluripotency maintaining transcription factors, cell cycle and cell proliferation/survival genes and TCF/LEF targets. Finally, β-catenin/TCF-LEF pathway mediated the biological effects of SATB2 in CSCs. These studies support the role of SATB2/β-catenin/TCF-LEF pathway in transformation and carcinogenesis.

BioASF: a framework for automatically generating executable pathway models specified in BioPAX

Motivation: Biological pathways play a key role in most cellular functions. To better understand these functions, diverse computational and cell biology researchers use biological pathway data for various analysis and modeling purposes. For specifying these biological pathways, a community of researchers has defined BioPAX and provided various tools for creating, validating and visualizing BioPAX models. However, a generic software framework for simulating BioPAX models is missing. Here, we attempt to fill this gap by introducing a generic simulation framework for BioPAX. The framework explicitly separates the execution model from the model structure as provided by BioPAX, with the advantage that the modelling process becomes more reproducible and intrinsically more modular; this ensures natural biological constraints are satisfied upon execution. The framework is based on the principles of discrete event systems and multi-agent systems, and is capable of automatically generating a hierarchical multi-agent system for a given BioPAX model.

Results: To demonstrate the applicability of the framework, we simulated two types of biological network models: a gene regulatory network modeling the haematopoietic stem cell regulators and a signal transduction network modeling the Wnt/β-catenin signaling pathway. We observed that the results of the simulations performed using our framework were entirely consistent with the simulation results reported by the researchers who developed the original models in a proprietary language.

Availability and Implementation: The framework, implemented in Java, is open source and its source code, documentation and tutorial are available at http://www.ibi.vu.nl/programs/BioASF.

Contact:j.heringa@vu.nl

Value of the level of methylation of RASSF1A and WIF-1 in tissue and serum in neoadjuvant chemotherapeutic assessment for advanced breast cancer

This study assessed the clinical efficacy of the neoadjuvant chemotherapy TAC scheme in treatment of patients with locally advanced breast cancer, and the value of the level of Ras association domain family 1A (RASSF1A) gene methylation and the Wnt inhibitory factor (WIF)-1 gene in tissue and serum of patients in clinical outcome prediction. In total, 126 patients were consecutively selected to receive TAC scheme (docetaxel, pirarubicin/epirubicin and cyclophosphamide) for at least four cycles with the total effective rate. The incidence of complications, progression-free survival and survival rate were recorded. Tumor tissues and peripheral blood samples collected in this study was used to detect methylation positive rate of RASSF1A and WIF-1 by methylation-specific PCR method and the relative level of expression of RASSF1A and WIF-1 mRNA by reverse transcription PCR method. Of the 126 patients, there were 18 cases with complete response (CR), 32 cases with partial response (PR), 50 cases with stable disease (SD), and 26 cases with disease progression (PD) with a total effective rate of 79.37%. Comparison on baseline data of effective group and ineffective group showed no difference (P>0.05), and comparison on adverse reactions occurrence showed no difference (P>0.05). Progression-free survival of the effective group was prolonged with a significant increase in survival rate (P<0.05). Positive rates of RASSF1A methylation and WIF-1 in tissue and serum of the patients in the effective group were significantly lower than those in the ineffective group, but the mRNA of RASSF1A and WIF-mRNA was significantly higher than the ineffective group (P<0.05). The sensitivity of clinical outcome prediction using tissue RASSF1A methylation was 67.0%, the specificity 15.4%, positive predictive value 69.0% and negative predictive value 31.0%. The above-mentioned indexes of tissue WIF-1 were 76.0, 31.4, 72.2 and 27.8, respectively. The indexes of serum RASSF1A were 85.0, 50.0, 76.2 and 23.8%, respectively, and the indexes of serum WIF-1 were 94.0, 75.0, 81.0 and 19.0%, respectively. The receiver operating characteristic curve analysis suggested that the accuracy of clinical outcome prediction using tissue RASSF1A mRNA level was 0.812. The sensitivity 85.2%, the specificity 76.3% and the critical value 0.4256. These indexes of tissue WIF-1 were 0.833, 86.7%, 75.4% and 0.3562 for CR, PR, SD and PD, respectively. These indexes of serum RASSF1A were 0.864, 88.3%, 77.4% and 0.2564, respectively, and for serum WIF-1 were 0.882, 89.4%, 73.5% and 0.1562, respectively. In conclusion, the detection of RASSF1A and WIF-1 gene methylation and level of mRNA expression in tissue and serum of patients with locally advanced breast cancer has an important application value in predicting clinical efficacy of neoadjuvant chemotherapy of the TAC scheme.