Modeling the Mechanism of GR/c-Jun/Erg Crosstalk in Apoptosis of Acute Lymphoblastic Leukemia

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.

Here we present modelling of the glucocorticoid receptor (GR) signalling network. The GR is the effector for a class of drugs known as corticosteroids, which are widely used in medicine for their anti-inflammatory effects and ability to induce apoptosis in leukaemic cells. However, side effects, treatment-related toxicity and glucocorticoid resistance remain and therefore increased understanding of the glucocorticoid receptor mechanism of action may improve therapeutic outcomes. The GEB052 model presented herein has been used to generate predictions for how the network is altered between glucocorticoid-sensitive and glucocorticoid-resistant scenarios, and these predictions have been verified using published gene expression data from established cell lines (for both qualitative and semi-quantitative analysis). The model has also been preliminarily assessed as a predictive clinical tool by correlating model predictions with clinical outcomes of thirteen leukaemia patients. Thus, the GEB052 model demonstrates successful modelling to understand GR function. GEB052 provides accurate predictions and has indicated potential routes through which glucocorticoid resistance may arise. The work presented herein thus demonstrates a proof-of-principle of this modelling approach to furthering GR research, and provides insight into potential mechanisms of corticosteroids resistance.