Insight into glucocorticoid receptor signalling through interactome model analysis

Glucocorticoid hormones (GCs) are used to treat a variety of diseases because of their potent anti-inflammatory effect and their ability to induce apoptosis in lymphoid malignancies through the glucocorticoid receptor (GR). Despite ongoing research, high glucocorticoid efficacy and widespread usage in medicine, resistance, disease relapse and toxicity remain factors that need addressing. Understanding the mechanisms of glucocorticoid signalling and how resistance may arise is highly important towards improving therapy. To gain insight into this we undertook a systems biology approach, aiming to generate a Boolean model of the glucocorticoid receptor protein interaction network that encapsulates functional relationships between the GR, its target genes or genes that target GR, and the interactions between the genes that interact with the GR. This model named GEB052 consists of 52 nodes representing genes or proteins, the model input (GC) and model outputs (cell death and inflammation), connected by 241 logical interactions of activation or inhibition. 323 changes in the relationships between model constituents following in silico knockouts were uncovered, and steady-state analysis followed by cell-based microarray genome-wide model validation led to an average of 57% correct predictions, which was taken further by assessment of model predictions against patient microarray data. Lastly, semi-quantitative model analysis via microarray data superimposed onto the model with a score flow algorithm has also been performed, which demonstrated significantly higher correct prediction ratios (average of 80%), and the model has been assessed as a predictive clinical tool using published patient microarray data. In summary we present an in silico simulation of the glucocorticoid receptor interaction network, linked to downstream biological processes that can be analysed to uncover relationships between GR and its interactants. Ultimately the model provides a platform for future development both by directing laboratory research and allowing for incorporation of further components, encapsulating more interactions/genes involved in glucocorticoid receptor signalling

p53 modeling as a route to mesothelioma patients stratification and novel therapeutic identification

Background Malignant pleural mesothelioma (MPM) is an orphan disease that is difficult to treat using traditional chemotherapy, an approach which has been effective in other types of cancer. Most chemotherapeutics cause DNA damage leading to cell death. Recent discoveries have highlighted a potential role for the p53 tumor suppressor in this disease. Given the pivotal role of p53 in the DNA damage response, here we investigated the predictive power of the p53 interactome model for MPM patients’ stratification. Methods We used bioinformatics approaches including omics type analysis of data from MPM cells and from MPM patients in order to predict which pathways are crucial for patients’ survival. Analysis of the PKT206 model of the p53 network was validated by microarrays from the Mero-14 MPM cell line and RNA-seq data from 71 MPM patients, whilst statistical analysis was used to identify the deregulated pathways and predict therapeutic schemes by linking the affected pathway with the patients’ clinical state. Results In silico simulations demonstrated successful predictions ranging from 52 to 85% depending on the drug, algorithm or sample used for validation. Clinical outcomes of individual patients stratified in three groups and simulation comparisons identified 30 genes that correlated with survival. In patients carrying wild-type p53 either treated or not treated with chemotherapy, FEN1 and MMP2 exhibited the highest inverse correlation, whereas in untreated patients bearing mutated p53, SIAH1 negatively correlated with survival. Numerous repositioned and experimental drugs targeting FEN1 and MMP2 were identified and selected drugs tested. Epinephrine and myricetin, which target FEN1, have shown cytotoxic effect on Mero-14 cells whereas marimastat and batimastat, which target MMP2 demonstrated a modest but significant inhibitory effect on MPM cell migration. Finally, 8 genes displayed correlation with disease stage, which may have diagnostic implications. Conclusions Clinical decisions related to MPM personalized therapy based on individual patients’ genetic profile and previous chemotherapeutic treatment could be reached using computational tools and the predictions reported in this study upon further testing in animal models

Differential regulation of cell death pathways by the microenvironment correlates with chemoresistance and survival in leukaemia

Glucocorticoids (GCs) and topoisomerase II inhibitors are used to treat acute lymphoblastic leukaemia (ALL) as they induce death in lymphoid cells through the glucocorticoid receptor (GR) and p53 respectively. Mechanisms underlying ALL cell death and the contribution of the bone marrow microenvironment to drug response/resistance remain unclear. The role of the microenvironment and the identification of chemoresistance determinants were studied by transcriptomic analysis in ALL cells treated with Dexamethasone (Dex), and Etoposide (Etop) grown in the presence or absence of bone marrow conditioned media (CM). The necroptotic (RIPK1) and the apoptotic (caspase-8/3) markers were downregulated by CM, whereas the inhibitory effects of chemotherapy on the autophagy marker Beclin-1 (BECN1) were reduced suggesting CM exerts cytoprotective effects. GCs upregulated the RIPK1 ubiquitinating factor BIRC3 (cIAP2), in GC-sensitive (CEM-C7-14) but not in resistant (CEM-C1-15) cells. In addition, CM selectively affected GR phosphorylation in a site and cell-specific manner. GR is recruited to RIPK1, BECN1 and BIRC3 promoters in the sensitive but not in the resistant cells with phosphorylated GR forms being generally less recruited in the presence of hormone. FACS analysis and caspase-8 assays demonstrated that CM promoted a pro-survival trend. High molecular weight proteins reacting with the RIPK1 antibody were modified upon incubation with the BIRC3 inhibitor AT406 in CEM-C7-14 cells suggesting that they represent ubiquitinated forms of RIPK1. Our data suggest that there is a correlation between microenvironment-induced ALL proliferation and altered response to chemotherapy