Analysis of IL6 Signal Transduction Model using Reduced Rank Regression

This work presents the application of reduced rank regression to the field of systems biology. A computational approach is used to investigate the mechanisms of the janus-associated kinases/signal transducers and transcription factors (JAK/STAT) and mitogen activated protein kinases (MAPK) signal transduction pathways in hepatic cells stimulated by interleukin-6. The results obtained identify the contribution of individual reactions to the dynamics of the model. These findings are compared to previously available results from sensitivity analysis of the model which focused on the parameters involved and their effect. This application of reduced rank regression allows for an understanding of the individual reaction terms involved in the modelled signal transduction pathways and has the benefit of being computationally inexpensive. The obtained results complement existing findings and also confirm the importance of several protein complexes in the MAPK pathway which hints at benefits that can be achieved by further refining the model.

S100A2 is a BRCA1/p63 coregulated tumour suppressor gene with roles in the regulation of mutant p53 stability

Here, we show for the first time that the familial breast/ovarian cancer susceptibility gene, BRCA1, along with interacting ΔNp63 proteins, transcriptionally upregulate the putative tumour suppressor protein, S100A2. Both BRCA1 and ΔNp63 proteins are required for S100A2 expression. BRCA1 requires ΔNp63 proteins for recruitment to the S100A2 proximal promoter region, while exogenous expression of individual ΔNp63 proteins cannot activate S100A2 transcription in the absence of a functional BRCA1. Consequently, mutation of the ΔNp63/p53 response element within the S100A2 promoter completely abrogates the ability of BRCA1 to upregulate S100A2. S100A2 shows growth control features in a range of cell models. Transient or stable exogenous S100A2 expression inhibits the growth of BRCA1 mutant and basal-like breast cancer cell lines, while short interfering RNA (siRNA) knockdown of S100A2 in non-tumorigenic cells results in enhanced proliferation. S100A2 modulates binding of mutant p53 to HSP90, which is required for efficient folding of mutant p53 proteins, by competing for binding to HSP70/HSP90 organising protein (HOP). HOP is a cochaperone that is required for the efficient transfer of proteins from HSP70 to HSP90. Loss of S100A2 leads to an HSP90-dependent stabilisation of mutant p53 with a concomitant loss of p63. Accordingly, S100A2-deficient cells are more sensitive to the HSP-90 inhibitor, 17-N-allylamino-17-demethoxygeldanamycin, potentially representing a novel therapeutic strategy for S100A2- and BRCA1-deficient cancers. Taken together, these data demonstrate the importance of S100A2 downstream of the BRCA1/ΔNp63 signalling axis in modulating transcriptional responses and enforcing growth control mechanisms through destabilisation of mutant p53.

Genome-wide characterization reveals complex interplay between TP53 and TP63 in response to genotoxic stress

In response to genotoxic stress the TP53 tumour suppressor activates target gene expression to induce cell cycle arrest or apoptosis depending on the extent of DNA damage. These canonical activities can be repressed by TP63 in normal stratifying epithelia to maintain proliferative capacity or drive proliferation of squamous cell carcinomas, where TP63 is frequently overexpressed/amplified. Here we use ChIP-sequencing, integrated with microarray analysis, to define the genome-wide interplay between TP53 and TP63 in response to genotoxic stress in normal cells. We reveal that TP53 and TP63 bind to overlapping, but distinct cistromes of sites through utilization of distinctive consensus motifs and that TP53 is constitutively bound to a number of sites. We demonstrate that cisplatin and adriamycin elicit distinct effects on TP53 and TP63 binding events, through which TP53 can induce or repress transcription of an extensive network of genes by direct binding and/or modulation of TP63 activity. Collectively, this results in a global TP53-dependent repression of cell cycle progression, mitosis and DNA damage repair concomitant with activation of anti-proliferative and pro-apoptotic canonical target genes. Further analyses reveal that in the absence of genotoxic stress TP63 plays an important role in maintaining expression of DNA repair genes, loss of which results in defective repair.

Bioinformatic evaluation of transcriptional regulation of WNT pathway genes with reference to diabetic nephropathy

Objective: Diabetic nephropathy (DN) is a microvascular complication of diabetes. Members of the WNT/ β-catenin pathways have been implicated in interstitial fibrosis and glomerular sclerosis, characteristic hallmarks of DN. These processes are controlled, in part, by transcription factors (TFs), proteins which bind to gene promoter regions attenuating their regulation. We sought to identify predicted cis-acting transcription factor binding sites (TFBS) over-represented within the promoter regions of WNT pathway members compared to genes across the genome.Methods: We assessed the frequency of 62 TFBS motifs from the JASPAR databases on 65 WNT pathway genes. P-values were estimated on the hypergeometric distribution for each TF. Gene expression profiles of enriched motifs were examined from DN-related datasets to assess clinical significance.Results: TFBS motifs transcription factor AP-2 alpha (TFAP2A), myeloid zinc finger 1 (MZF1), and specificity protein 1 (SP1) were significantly enriched within WNT pathway genes (P-values<6.83x10-29, 1.34x10-11 and 3.01x10-6 respectively). MZF1 gene expression was significantly increased in DN in a whole kidney dataset (fold change = 1.16; 16% increase; P = 0.03). TFAP2A gene expression was decreased in an independent dataset (fold change = -1.02; P = 0.03). SP1 was not differentially expressed in any datasets examined.Conclusions: Three TFBS profiles are significantly enriched within the WNT pathway genes examined highlighting the use of in silico analyses for identifying key regulators of this pathway. Modification of TF binding to gene promoter regions involved in DN pathology may limit progression, making refinement of targeted therapeutic strategies possible through clearer delineation of their role.

Role played by BRCA1 in transcriptional regulation in response to therapy

BRCA1 (breast-cancer susceptibility gene 1) is a tumour suppressor, implicated in the hereditary predisposition to breast and ovarian cancer. BRCA1 has been implicated in a number of cellular processes including DNA repair and recombination, cell cycle checkpoint control, chromatin remodelling and ubiquitination. In addition, substantial data now exist to suggest a role for BRCA1 in transcriptional regulation; BRCA1 has been shown to interact with the Pol II holoenzyme complex and to interact with multiple transcription factors, such as p53 and c-Myc. We have previously identified a range of BRCA1 transcriptional targets and have linked these to specific cellular pathways, including cell cycle checkpoint activation and apoptosis. Current research is focused on the transcriptional mechanisms that underpin the association of BRCA1 deficiency with increased sensitivity to DNA damage-based chemotherapy and resistance to spindle poisons.

Inhibition of cellular proliferation by the Wilms tumor suppressor WT1 requires association with the inducible chaperone Hsp70

The Wilms tumor suppressor WT1 encodes a zinc finger transcription factor that is expressed in glomerular podocytes during a narrow window in kidney development. By immunoprecipitation and protein microsequencing analysis, we have identified a major cellular protein associated with endogenous WT1 to be the inducible chaperone Hsp70. WT1 and Hsp70 are physically associated in embryonic rat kidney cells, in primary Wilms tumor specimens and in cultured cells with inducible expression of WT1. Colocalization of WT1 and Hsp70 is evident within podocytes of the developing kidney, and Hsp70 is recruited to the characteristic subnuclear clusters that contain WT1. The amino-terminal transactivation domain of WT1 is required for binding to Hsp70, and expression of that domain itself is sufficient to induce expression of Hsp70 through the heat shock element (HSE). Substitution of a heterologous Hsp70-binding domain derived from human DNAJ is sufficient to restore the functional properties of a WT1 protein with an amino-terminal deletion, an effect that is abrogated by a point mutation in DNAJ that reduces binding to Hsp70. These observations indicate that Hsp70 is an important cofactor for the function of WT1, and suggest a potential role for this chaperone during kidney differentiation.

Activation of the c-Jun N-terminal kinase (JNK) signaling pathway is essential during PBOX-6-induced apoptosis in chronic myelogenous leukemia (CML) cells

The mitogen-activated protein (MAP) kinase family is activated in response to a wide variety of external stress signals such as UV irradiation, heat shock, and many chemotherapeutic drugs and leads to the induction of apoptosis. A novel series of pyrrolo-1,5-benzoxazepines have been shown to potently induce apoptosis in chronic myelogenous leukemia (CML) cells, which are resistant to many chemotherapeutic agents. In this study we have delineated part of the mechanism by which a representative compound known as PBOX-6 induces apoptosis. We have investigated whether PBOX-6 induces activation of MAP kinase signaling pathways in CML cells. Treatment of K562 cells with PBOX-6 resulted in the transient activation of two JNK isoforms, JNK1 and JNK2. In contrast, PBOX-6 did not activate the extracellular signal-regulated kinase (ERK) or p38. Apoptosis was found to occur independently of the small GTPases Ras, Rac, and Cdc42 but involved phosphorylation of the JNK substrates, c-Jun and ATF-2. Pretreatment of K562 cells with the JNK inhibitor, dicoumarol, abolished PBOX-6-induced phosphorylation of c-Jun and ATF-2 and inhibited the induced apoptosis, suggesting that JNK activation is an essential component of the apoptotic pathway induced by PBOX-6. Consistent with this finding, transfection of K562 cells with the JNK scaffold protein, JIP-1, inhibited JNK activity and apoptosis induced by PBOX-6. JIP-1 specifically scaffolds JNK, MKK7, and members of the mixed-lineage kinase (MLK) family, implicating these kinases upstream of JNK in the apoptotic pathway induced by PBOX-6 in K562 cells.

Radiation to control transgene expression in tumors

Significant evidence has accumulated indicating that certain genes are induced by ionising radiation. An implication of this observation is that their promoter regions include radiation-responsive sequences. These sequences have been isolated in the promoter of several genes including Erg-1, p21/WAF-1, GADD45alpha and t-PA. The mechanism by which radiation induces gene expression remains unclear but involves putative binding sites for selected transcription factors and/or p53. Consensus CC(A/T)6GG sequences have been localized in the Erg-1 promoter and are referred to as serum response elements or CArG elements. The tandem combination of CArG elements has been shown to improve gene expression levels, with a 9-copy motif conferring maximum inducibility. The response of these genes to ionising radiation appears to follow a sigmoid relationship with time and dose. Therapeutic induction of suicide genes and significant cytotoxicity can be achieved at clinically relevant x-rays doses both in vitro and in vivo but was found to be cell-type dependent. Radiation-inducible gene therapy can be potentially enhanced by exploiting hypoxia through the inclusion of hypoxia-response element motifs in the expression cassette, the use of the anaerobic bacteria or the use of neutron irradiation. These results are encouraging and provide significant evidence that gene therapy targeted to the radiation field is a reasonably attractive therapeutic option and could help overcome hypoxic radioresistant tumors.

Glucocorticoid receptor and Klf4 co-regulate anti-inflammatory genes in keratinocytes

The glucocorticoid (GC) receptor (GR) and Kruppel-like factor Klf4 are transcription factors that play major roles in skin homeostasis. However, whether these transcription factors cooperate in binding genomic regulatory regions in epidermal keratinocytes was not known. Here, we show that in dexamethasone-treated keratinocytes GR and Klf4 are recruited to genomic regions containing adjacent GR and KLF binding motifs to control transcription of the anti-inflammatory genes Tsc22d3 and Zfp36. GR- and Klf4 loss of function experiments showed total GR but partial Klf4 requirement for full gene induction in response to dexamethasone. In wild type keratinocytes induced to differentiate, GR and Klf4 protein expression increased concomitant with Tsc22d3 and Zfp36 up-regulation. In contrast, GR-deficient cells failed to differentiate or fully induce Klf4, Tsc22d3 and Zfp36 correlating with increased expression of the epithelium-specific Trp63, a known transcriptional repressor of Klf4. The identified transcriptional cooperation between GR and Klf4 may determine cell-type specific regulation and have implications for developing therapies for skin diseases.

Myc-dependent purine biosynthesis affects nucleolar stress and therapy response in prostate cancer

The androgen receptor is a key transcription factor contributing to the development of all stages of prostate cancer (PCa). In addition, other transcription factors have been associated with poor prognosis in PCa, amongst which c-Myc (MYC) is a well-established oncogene in many other cancers. We have previously reported that the AR promotes glycolysis and anabolic metabolism; many of these metabolic pathways are also MYC-regulated in other cancers. In this study, we report that in PCa cells de novo purine biosynthesis and the subsequent conversion to XMP is tightly regulated by MYC and independent of AR activity. We characterized two enzymes, PAICS and IMPDH2, within the pathway as PCa biomarkers in tissue samples and report increased efficacy of established anti-androgens in combination with a clinically approved IMPDH inhibitor, mycophenolic acid (MPA). Treatment with MPA led to a significant reduction in cellular guanosine triphosphate (GTP) levels accompanied by nucleolar stress and p53 stabilization. In conclusion, targeting purine biosynthesis provides an opportunity to perturb PCa metabolism and enhance tumour suppressive stress responses.

HES6 drives a critical AR transcriptional programme to induce castration-resistant prostate cancer through activation of an E2F1-mediated cell cycle network

Castrate-resistant prostate cancer (CRPC) is poorly characterized and heterogeneous and while the androgen receptor (AR) is of singular importance, other factors such as c-Myc and the E2F family also play a role in later stage disease. HES6 is a transcription co-factor associated with stem cell characteristics in neural tissue. Here we show that HES6 is up-regulated in aggressive human prostate cancer and drives castration-resistant tumour growth in the absence of ligand binding by enhancing the transcriptional activity of the AR, which is preferentially directed to a regulatory network enriched for transcription factors such as E2F1. In the clinical setting, we have uncovered a HES6-associated signature that predicts poor outcome in prostate cancer, which can be pharmacologically targeted by inhibition of PLK1 with restoration of sensitivity to castration. We have therefore shown for the first time the critical role of HES6 in the development of CRPC and identified its potential in patient-specific therapeutic strategies.

HOXB13, RFX6 and prostate cancer risk

A new study shows that HOXB13 is preferentially recruited to the risk allele of a prostate cancer-associated SNP, enhancing the expression of RFX6, a driver of prostate cancer cell migration and predictor of disease progression. The work illustrates how a single risk locus contributes both to prostate cancer incidence and, through functional follow-up, to disease progression.

Molecular subtyping of primary prostate cancer reveals specific and shared target genes of different ETS rearrangements

This work aimed to evaluate whether ETS transcription factors frequently involved in rearrangements in prostate carcinomas (PCa), namely ERG and ETV1, regulate specific or shared target genes. We performed differential expression analysis on nine normal prostate tissues and 50 PCa enriched for different ETS rearrangements using exon-level expression microarrays, followed by in vitro validation using cell line models. We found specific deregulation of 57 genes in ERG-positive PCa and 15 genes in ETV1-positive PCa, whereas deregulation of 27 genes was shared in both tumor subtypes. We further showed that the expression of seven tumor-associated ERG target genes (PLA1A, CACNA1D, ATP8A2, HLA-DMB, PDE3B, TDRD1, and TMBIM1) and two tumor-associated ETV1 target genes (FKBP10 and GLYATL2) was significantly affected by specific ETS silencing in VCaP and LNCaP cell line models, respectively, whereas the expression of three candidate ERG and ETV1 shared targets (GRPR, KCNH8, and TMEM45B) was significantly affected by silencing of either ETS. Interestingly, we demonstrate that the expression of TDRD1, the topmost overexpressed gene of our list of ERG-specific candidate targets, is inversely correlated with the methylation levels of a CpG island found at -66 bp of the transcription start site in PCa and that TDRD1 expression is regulated by direct binding of ERG to the CpG island in VCaP cells. We conclude that ETS transcription factors regulate specific and shared target genes and that TDRD1, FKBP10, and GRPR are promising therapeutic targets and can serve as diagnostic markers for molecular subtypes of PCa harboring specific fusion gene rearrangements.

Mapping protein-DNA interactions using ChIP-sequencing

Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.

Nuclear translocation and functions of growth factor receptors

Cellular signal transduction in response to environmental signals involves a relay of precisely regulated signal amplifying and damping events. A prototypical signaling relay involves ligands binding to cell surface receptors and triggering the activation of downstream enzymes to ultimately affect the subcellular distribution and activity of DNA-binding proteins that regulate gene expression. These so-called signal transduction cascades have dominated our view of signaling for decades. More recently evidence has accumulated that components of these cascades can be multifunctional, in effect playing a conventional role for example as a cell surface receptor for a ligand whilst also having alternative functions for example as transcriptional regulators in the nucleus. This raises new challenges for researchers. What are the cues/triggers that determine which role such proteins play? What are the trafficking pathways which regulate the spatial distribution of such proteins so that they can perform nuclear functions and under what circumstances are these alternative functions most relevant?