A quinazoline-based HDAC inhibitor affects gene expression pathways involved in cholesterol biosynthesis and mevalonate in prostate cancer cells †

Chronic inflammation can lead to the development of cancers and resolution of inflammation is an ongoing challenge. Inflammation can result from dysregulation of the epigenome and a number of compounds that modify the epigenome are in clinical use. In this study the anti-inflammatory and anti-cancer effects of a quinazoline epigenetic-modulator compound were determined in prostate cancer cell lines using a non-hypothesis driven transcriptomics strategy utilising the Affymetrix PrimeView Human Gene Expression microarray. GATHER and IPA software were used to analyse the data and to provide information on significantly modified biological processes, pathways and networks. A number of genes were differentially expressed in both PC3 and DU145 prostate cancer cell lines. The top canonical pathways that frequently arose across both cell lines at a number of time points included cholesterol biosynthesis and metabolism, and the mevalonate pathway. Targeting of sterol and mevalonate pathways may be a powerful anticancer approach.


Introduction
Although inflammation is a necessary immune response, chronic inflammation can lead to the development of cancers as well as other chronic diseases, and there are ongoing efforts to either block pro-inflammatory mediators or to stimulate resolution of inflammation. 1 Deregulation of the epigenome can result in a range of chronic diseases associated with inflammation, and a number of compounds that modify the epigenome are in clinical use today. 2 Chromatin remodelling plays a key role in gene expression, and epigenetic modifications may be as important as mutations, insertions and deletions in tumour development and progression. 3 Unlike genetic mutations, epigenetic changes that have resulted in gene activation or silencing can sometimes be reversed by small molecules that modify the epigenome. One such group of small molecules are the histone deacetylase (HDAC) inhibitors. HDACs influence the expression of a number of key enzymes involved in pathways associated with apoptosis, cell cycle, tumour cell proliferation and inflammation, amongst others 4 and tumour progression is associated with an increase in HDAC activity. 5 However, although HDAC inhibitors have an impact on tumour and T cell lymphomas rather than non-malignant cells, their mechanism of action remains unclear. 4 HDACs and histone acetyltransferases (HATs) act in opposition to modify chromatin and thus control gene expression. 6 HDACs can repress transcription by bringing about chromatin condensation in response to the removal of acetyl groups from histone tails. 6 Not only have HDACs been found to be aberrantly recruited to ''inappropriate'' loci, but abnormal expression of HDACs 1, 2, 3 and 6 have been reported in numerous types of cancer e.g. gastric, breast, prostate, colorectal and cervical. 6 SN30028 is an HDAC inhibitor that was identified from an in-house compound library. 7 SN30028 (Fig. 1) is regarded as a quinazoline drug and was selected from the aforementioned compound library based on its anti-inflammatory activity and the strength of HDAC inhibition. 7,8  SN30028 decreased the activity of HDAC 1, 3 and 6 by 23%, 76% and 48% respectively. 9 HDACs 1 and 3, HDAC class I compounds, are restricted to the nucleus and are believed to play a key role in cell survival and proliferation. 10 Loss of HDAC 1 activity results in an overall reduction of deacetylase activity, reduced proliferation rates and increased levels of the cyclindependent kinase inhibitors p21 and p27. 11 HDAC 3 is important as it mediates gene expression of tumour necrosis factor (TNF) as well as the expression of other genes. 12 HDAC 6 belongs to the HDAC Class 2b group of compounds and is unique as it has two catalytic domains and a zinc finger. 6 HDAC 6 is of interest as it helps to protect against cellular stress by the regulation of heat shock protein 90 and alpha tubulin and down-regulation can bring about apoptosis and inhibition of metastasis. 13,14 Determining the effect of a particular compound on cancer cell lines can be challenging as effects on any one gene can be small, and these effects can also be broad. For this reason the anti-inflammatory and anti-cancer effects of a quinazoline epigenetic-modulator compound were determined, in prostate cancer cell lines, using a transcriptomics approach.
One advantage of transcriptomics is that experimental design is non-hypothesis driven, and provides sufficient sensitivity and breadth to examine the expression of thousands of genes simultaneously. 15 The Affymetrix PrimeView GeneChip Human Microarray was used as it is a ''perfect-match-only'' (probe to transcript) array 16 and therefore the false signal changes referred to by Li et al. 17 are less likely to arise.
In addition biological processes are likely to be represented by complex networks consisting of multiple signalling modules rather than a series of linear pathways 18 and thus a transcriptomics approach, followed by network analysis, was deemed preferable. The aim of this study was to determine the effect of SN30028 on differential gene expression in prostate cancer cell lines with a particular focus on inflammation and epigenetic modulation. Compounds, that restore histone acetylation, have potential as anti-inflammatory and anti-cancer drugs, 19 and we show evidence that the known 7 quinazoline-based HDAC inhibitor SN30028 ( Fig. 1) influenced cholesterol biosynthesis and mevalonate pathways.

Cell culture
The prostate cancer cell lines PC3, DU145 and LNCaP were obtained from the American Type Culture Collection (VA, USA) and maintained in Alpha 5 minimum essential media (Gibco, Invitrogen Corporation, New York, USA) with 10% fetal calf serum (Morgate Biotech, Hamilton, New Zealand) and 100 IU penicillin/streptomycin (Sigma Chemical Company, St Louis, MO, USA). Cells were cultured in a water jacketed incubator (Thermo Fisher Scientific, MA, USA) at 37 1C with 5% CO 2 , and passaged twice weekly. The 50% inhibitory concentrations (IC 50 s) for SN30028 and a number of other compounds were established and the anti-proliferative activity was determined relative to the reference drugs suberoylanilide hydroxamic acid (SAHA) (Sigma Chemical Company, St Louis, MO, USA) and 5 0 -aza-2-deoxycytidine (Sigma-Aldrich, St Louis, MO, USA). Dimethyl sulphoxide (DMSO) (Merck, Darmstadt, Germany) was used as a solvent for the compounds tested and hence was also used as a control. The IC 50 s were performed on PC3, DU145 and LNCaP using the sulforhodamine B colorimetric assay. 20 Three independent experiments were performed in quadruplicate.

HDAC activity
HDAC activity was assessed from total protein extracted from DU145, PC3 and LNCaP cells following treatment for 24 hours with a number of compounds (TSA, SAHA, SN30028, SN30029, SN30140, SN29887, SN29984, SN30711, 5-aza-2-deoxycitidine (5AZA)) (Sigma Aldrich, St Louis, USA) at their respective IC 50 concentrations. Total protein was extracted from each treated cell line grown in P100 plates. The cells were washed with PBS, lysed (lysis buffer containing 1% Triton X-100 (Sigma Aldrich, St Louis, USA)), frozen, then thawed and removed to a microcentrifuge tube, whereupon the cells were centrifuged and the supernatant collected. Protein concentration was assessed using the bicinchoninic acid (Sigma Aldrich, St Louis, USA) protein assay (as outlined in Lin et al.). 7 Using the manufacturer's protocol the HDAC Fluorometric Activity Assay (BIOMOL International -Cayman Chemical, Ann Arbor, USA) was used to measure the effect of compounds on HDAC activity from the extracted protein. Three independent experiments were performed in duplicate.

Transcriptomics
Two prostate cancer cell lines, namely DU145 and PC3 cells, were treated with SN30028 and harvested at 4, 24 and 96 hours after dosing (Fig. 2). In addition, cell lines were incubated with an equivalent amount of water for each of the three time points. These served as controls. DU145 and PC3 were chosen ahead of LNCaP as the aforementioned were shown to have greater sensitivity to HDAC inhibition than LNCaP 9 (see Section 3.1) A QIAshredder (Qiagen, MD, USA) was used to disrupt the cells and an RNeasy Mini kit (Qiagen, MD, USA) was used to isolate the RNA in triplicate, following the manufacturer's recommended protocol. The quantity and quality of RNA was determined using the NanoDrop ND 1000 spectrophotometer (NanoDrop Technologies, Wilmington, USA) and the Experion RNA StdSens Analysis Kit (Biorad, Hercules, USA) respectively. Approximately 500 ng of RNA with a RIN value of 48 were required for the PrimeView s Human Gene Expression arrays (Affymetrix, CA, USA). The RNA was reverse transcribed, biotin labelled, fragmented, hybridised to the gene chip and scanned 21 by the Centre for Genomics and Proteomics (New Zealand Genomics Ltd, University of Auckland, New Zealand).
The PrimeView s Human Gene Expression array uses 530 000 probes covering 36 000 transcripts and variants located in more than 20 000 genes. 22 Transcripts were measured independently by using multiple probes. The level of gene expression was associated with the probe/s targeting that specific gene and following adjustment for the solvent/media controls, differential gene expression was calculated. The workflow for the analysis of the gene expression array data is outlined in Fig. 2.
Confirmation of the microarray results was carried out using quantitative reverse transcription polymerase chain reaction (qRT-PCR) and 27 statistically significant differentially expressed genes were selected for validation. The 27 TaqMan probe sets were obtained from ThermoFisher (Pleasanton, USA). The aforementioned RNA (Section 2.2) was isolated from DU145 and PC3 following treatment for 4 and 24 h, converted to cDNA using a Quantitect Reverse Transcription Kit (Qiagen, Victoria, Australia) and PCR was performed (on the three biological repeats as well as non-template controls) on an Applied Biosystems 7900 thermocycler (Waltham, USA). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), hypoxanthine phosphor-ribosyltransferase 1 (HPRT1), and b-actin (ACTB) were tested for assessment as normalisation genes for RNA expression. A standard curve was calculated from technical triplicates using SDS2.3 and RQ Manager 2.2 software (Applied Biosystems, USA). The relative expression of each of the genes was calculated as fold change using a delta delta cycle threshold (Ct) method. 24 Thereafter fold changes from the Affymetrix and qRT-PCR experiments were compared.

Data analysis
To identify the differentially expressed genes, the data from the microarray experiments were subjected to Robust Multiarray Analysis 25 to assess for quantile normalisation and then filtered based on the following parameters: (adjusted p value o0.05; fold change 41.5). The ''Limma'' package within the statistical language R 26 was used to analyse the CEL files generated by the Affymetrix GeneChip Instrument System (Affymetrix, MD, USA).
Gene expression values were based on the absorbance values of the probes targeting each specific gene. Changes in average gene expression in response to treatment were calculated as fold changes. p-Values were calculated so as to determine whether the changes in gene expression were significant (Fig. 2). A multiple testing correction was applied using Benjamini-Hochberg's false-discovery rate 27 set at 5%. Functional, pathway and network analyses were performed by using GATHER (Gene Annotation Tool to Help Explain Relationships 28 ) and IPA (Ingenuity Pathway Analysis 29 ) software. Networks were ranked by a score that was assigned based on the probability that a collection of nodes within that network can occur by chance. A number of networks scored highly, but only the highest scoring network per cell line -treatment time combination were presented here so as to avoid biasing results through ''fishing'' (Fig. S1-S6, ESI †). The network score indicates the likelihood of Focus Genes in a network being found together, and it is clear from the high scores that one should be confident that these associations were not generated by chance. 30 IPA constructs networks using a stepwise process. Networks are constructed in such a way that both interconnectivity between Focus Genes and the number of Focus Genes in the network are optimised (using a network size of 35 or 70 genes/ proteins, as specified in the supplementary figure legends).

Results
The four day IC 50 concentrations were established for PC3, DU145 and LNCaP cell lines and the results are listed in Table 1.
With the exception of SN26855, DU145 cells showed the greatest tolerance to the novel compounds tested. In general, LNCaP cells appeared to be the most sensitive to the novel compounds.

HDAC inhibition
Following treatment of cells for 24 h, total protein was extracted. The effects of compounds of interest on HDAC activity are shown in Fig. 3-5. The HDAC activity of total protein extracted from the cell lines, varied from one cell line to another and varied amongst the novel compounds. HDAC activity of DU145 cells was markedly reduced following treatment with all the novel compounds, particularly SN30028 ( Fig. 4), whilst SN30140 had the largest reduction in HDAC activity in PC3 cells (Fig. 3). The novel compounds did not appear to reduce HDAC activity in LNCaP cells (Fig. 5).

Differential gene expression (Affymetrix PrimeViewt Microarrays)
Following the analysis of the Affymetrix PrimeViewt microarray data, it was found that a number of genes in PC3 and DU145 were differentially expressed ( p o 0.05; fold change 41.5) following treatment with SN30028. In total 59, 7 and 64 genes were up-regulated and 38, 6 and 34 genes were down-regulated in PC3 at 4, 24 and 96 hours respectively. In DU145 49, 108 and 796 genes were up-regulated and 3, 4 and 399 genes were down-regulated at 4, 24 and 96 hours respectively. A list of ten of the most differentially expressed genes for each cell linetime point combination, with statistically significant p-values and a fold change of 41.5 are shown in Table 2 (PC3) and Table 3 (DU145). The genes listed in the table were selected based on the absolute fold change and presented from largest to smallest, regardless of the direction of change.

Validation of the differential gene expression data generated by Affymetrix PrimeViewt microarray
To corroborate the differential gene expression data generated by the Affymetrix assay a validation was undertaken by comparing    the aforementioned data with that obtained from a select number of genes tested using qRT-PCR. qRT-PCR is widely used for accurate expression profiling of key genes 31 and is a suitable technique to validate differential gene expression data generated from microarray experiments. Genes were selected based on a high level of differential gene expression in one or more cell line -time point combinations, an important role in a biological pathway/ network relevant to cancer or inflammation, or effects on epigenetics.
Gene expression was measured using qRT-PCR in the following genes: AREG, ARID5B, CDKN2B, CYP1B1, CYP51A1, DHCR7, DUSP10, EGFR, EPGN, IFIT2, HMGCR, HMGCS1, GULP1, IDI1, INSIG1, KRT17, LDLR, MMP1, MMP3, NR4A3, PTGS2, RSG4, SCD, SGK2, SQLE, TM7SF2 and TP53INP1 (Table S1, ESI †). The aforementioned genes were selected based on the level of differential gene expression and relevance to cancer/epigenetic mechanisms. Although the magnitude of change varied between the two methods with qRT-PCR generating the higher value in general, the direction of change remained consistent, with the exception of the gene ARID5B. The gene expression level of ARID5B was down-regulated according to the results obtained from the Affymetrix array (fold change of À1.82), and up-regulated according to the results generated by qRT-PCR (fold change of 1.43).
GAPDH and HPRT1 were used as normalisation genes as they showed consistent results in both PC3 and DU145 cell lines across the Ct range of the 27 genes tested. Following normalisation, the fold changes in expression of the selected genes were compared between those generated from the Affymetrix microarray and those generated from qRT-PCR experiments (Table S1, ESI †).

GATHER
GATHER software was used to analyse the relationships between individual or groups of genes. Only those genes that were differentially expressed with a p-value r0.05 and a minimum fold change of 1.5 were uploaded to GATHER. SN30028 treatment modulated biological processes and the involved genes are summarised in Table S2 (ESI †). Those biological processes with Bayes factors Z3 and p-values r 0.05 were nominated. 32 The Gene Ontology (GO) database was used to obtain definitions of biological processes. 33 The biological processes with the highest Bayes factors following treatment of PC3 cells with SN30028 at 4, 24 and 96 hours were: the transforming growth factor beta receptor signalling process (3.53), cyclooxygenase process (4.72), and response to nutrients process (5.69), respectively. In DU145 cells the sterol biosynthesis and sterol metabolism processes had the highest Bayes factors at 4 (29.66 and 26.51 respectively) and 24 hours (51.26 and 49.75 respectively), and cell proliferation and cell cycle were the most affected biological processes at 96 hour (25.21 and 21.49 respectively) following SN30028 treatment.

Ingenuity pathways analysis
IPA software (IPA s , QIAGEN Redwood City, USA) was used to transform a list of genes into a set of relevant networks 29 and this software is widely used for biological pathway, function and disease, and molecular network analysis. 34,35 All results from each time-point-cell line combination were submitted to IPA with the filters and settings as follows: adjusted p-value of 0.05 and fold-change of 1.5. The most significantly affected canonical pathways and associated network functions in response to treatment by SN30028 at three time points in PC3 and DU145 cell lines generated through IPA are shown in Table 4. Using Fischer's exact test, the results were ranked. The ratios listed in the table represent the number of genes in the dataset that were in the canonical pathway stated, divided by the total number of genes in that particular pathway. Obtaining a list of genes and related pathways is informative, but it is the identification of the connections between the pathways that is important. The IPA networks are assembled based on connectivity between genes. Several networks were generated by IPA software from each microarray experiment using the IPA knowledge base, but only the strongest network for each cell line -time point combination are shown here (Fig. S1-S6, ESI †). The scores of the IPA network indicate how relevant the network is to the genes in the uploaded dataset (Score = Àlog 10 ( p-value)). It is evident from Table 4 that cancer-related pathways are highly affected by SN30028 treatment in both cell lines. The parameters were set at either 70 or 35 molecules per network and a direct interaction between the molecules. Only relevant genes with an adjusted p-value o0.05 and fold change 41.5 or oÀ1.5 were selected for network analysis.
IPA generated results regarding the most prominent diseases and functions associated with the network depicted in Fig. S1-S6 (ESI †). These diseases and functions included developmental disorders, cancer, inflammatory disorders, cell death and survival, cellular growth and proliferation, cellular function and maintenance and DNA replication amongst others. Using the Fisher's exact test IPA calculates a network score which is the log 10 of the p-value of the network. 29 The network scores ranged from 38 (Fig. S2, ESI †) to 84 (Fig. S6, ESI †) and between 13% (Fig. S2, ESI †) and 100% (Fig. S6, ESI †) of the genes in each network were differentially expressed.

Discussion
Regulation of gene expression is encoded by the genome and by the epigenome. The focus of this discussion is on the modulation of gene expression and the molecular interaction of differentially  expressed genes, represented in the form of networks, generated using transcriptomics and analysed using IPA. The most common pathways, as shown in Table 4, will be examined. Similarly, central nodes in networks will be discussed if they are evident in more than one network.
Novel compounds were used to treat prostate cancer cell lines, inhibition of cell proliferation was observed and HDAC activity was assessed to determine which compound-cell line combination to use in the transcriptomics experiments. DU145 was found to be more tolerant to the novel compounds than PC3 and LNCaP, and LNCaP was found to be the most sensitive to the novel compounds with respect to cell proliferation. In contrast, DU145 had the greatest response to the novel compounds with respect to HDAC inhibition, whereas LNCaP had the least or no inhibitory response to these compounds. In addition, there was a far greater response observed in DU145 than in PC3 when considering the number of genes differentially expressed, as well as the size of the response.
There are a number of phenotypic and genotypic differences amongst PC3, DU145 and LNCaP. LNCaP is androgen sensitive, whilst PC3 and DU145 are androgen independent. HDAC inhibitors interfere with androgen receptor activity 36 and therefore it is likely that the cell lines would respond differently. Seeing that PC3 and DU145 are androgen independent, it was not surprising that a change in the AR gene was not observed.
In addition to androgen sensitivity, the three cell lines of interest also differ with regards to TNFb. When treated with TNFb, cell proliferation was initially inhibited, whilst no effect was observed in LNCaP cells. 37 TGFb induces epigenetic changes to modulate cell proliferation, differentiation and migration, and TGFb may initiate cellular changes that facilitate its role as both a tumour suppressor during the early stages of tumour development, and as a tumour promoter in metastatic or later stage disease. 5 HDAC inhibitors may inhibit the activation of TGFb in epithelial cells 5 by blocking TGFb mediated epithelial-mesenchymal transition (EMT), which is essential for cell growth and invasion. 38 In our study we found that TGFb was down regulated in response to SN30028, with fold changes between 1.074 and 1.707 (adjusted p values were significant for DU145 96 h treatment only).

Biological processes and pathways
Based on the transcriptomics results and analysis using GATHER and IPA, a number of cancer related ''biological processes,'' and ''network associated functions'' were affected by SN300028 treatment of the prostate cancer cell lines, DU145 and PC3 at all three time points tested (Table 4). Using GATHER, the most commonly modulated biological processes in DU145 (with high Bayes factor values) were in relation to sterol and cholesterol metabolism and biosynthesis (Table S2, ESI †) at the 4 and 24 h time points, whilst at 96 h the more directly cancer related biological processes of cell proliferation, cell cycle and DNA replication were evident. This is consistent with data generated using IPA. The top canonical pathways generated using IPA included those related to cholesterol biosynthesis (PC3 -24 and 96 h, and DU145 -4 and 24 h) and the mevalonate pathway (PC3 -96 h and DU145 -4 h).
Other studies have been carried out using gene expression arrays, or targeted gene expression to assess levels of gene expression in normal versus adenocarcinoma or precursor adenocarcinoma tissues as well as in cell lines in response to HDAC inhibitors. [39][40][41][42][43] Variations were observed in different cancer cell lines in response to HDAC inhibitors, with, in some cases, non-overlapping cellular targets. 43 Using different cell lines and a different HDAC inhibitor, it is not surprising that  View Article Online different pathways were modulated, although some overlap was observed with respect to modification of gene expression. Cholesterol metabolism plays an important role in providing cells with compounds for growth and sterol biosynthesis is an essential metabolic component of cancers. 44 In addition, overexpression of cholesterol biosynthesis pathways has been previously detected in refractory breast cancers 45 and this is consistent with the data observed from the metastatic prostate cancer cell lines we tested. More specifically, cytochrome P450 1B1 (CYP1B1) is important for the synthesis of cholesterol steroids and lipids, and is well known for its role in drug metabolism. 46 CYP1B1 activity is inhibited by a number of anti-cancer agents and is commonly over-expressed in a variety of tumours. 47 In our study the expression of CYP1B1 was down regulated in PC3 cells treated for 4 and 24 h, and it is suggested that the inhibition of CYP1B1 is brought about through the inhibition of HDAC6 activity. 48 In contrast, inconsistencies have arisen, for example the HDAC inhibitors SAHA and TSA induced CYP1B1 expression in the human breast cancer MCF-7. 49 The mevalonate pathway has a broad influence and is associated with the cholesterol related biosynthesis pathways; is important in cellular metabolism; plays a role in the maintenance of cell membranes; is involved in steroid biosynthesis and can be disrupted by medication prescribed for bone-density disorders and high cholesterol levels. 50,51 In addition, the mevalonate pathway is an important target for anti-cancer therapy and inhibitors of this pathway target malignant cell growth 50,51 and are believed to act through the modification of methylation status of CpG sites in gene promoter regions involved in apoptosis and/or cell proliferation. 52 In addition, DU145 cells, after 96 hours of treatment, showed the top canonical pathway affected was ''molecular mechanisms of cancer'' with 46 of the 359 genes involved were differentially expressed.
All cell line -treatment time combinations showed cancer to be one of the top five diseases and disorders, with the exception of DU145 at 24 h with a closely related disease/disorder, namely inflammatory response listed in the top five. The number of cancer related genes that were differentially expressed for each cell line -treatment time combination ranged from 11 to 445 genes with p values ranging from 8.00 Â 10 À2 to 2.15 Â 10 À17 . It is clear that numerous cancer related genes were differentially expressed in prostate cancer cell lines in response to treatment and therefore the compound, SN30028 is of interest. Similarly, Chang et al. found that a large number of genes were differentially expressed in response to the HDAC inhibitor, TSA, in non-small cell lung cancer. 43 However, in the study by Chang et al., 43 the foldchanges observed were much higher than those reported here.
''Associated network functions'' were also listed as an output from IPA for each cell line -treatment time combination. Cell cycle, cellular growth and proliferation, DNA replication and repair, inflammatory response, and cancer all ranked highly in one or more of the cell line -treatment time combinations and thus it is evident that SN30028 has an impact on cancer related mechanisms.

Networks
The networks shown in Fig. S1 and S2 (ESI †) demonstrate that prostaglandin H synthase (PTGS2), otherwise known as COX2, connects with many genes in networks generated from PC3 treated for 4 and 24 hours and was one of the top ten differentially expressed genes following PC3 treatment for 4 h ( Table 2). The down regulation of COX2 was corroborated in qRT-PCR experiments (Table S1, ESI †). COX2 is down-regulated by treatment with SN30028 for 4 and 24 hours (PC3) by a fold change of À1.71 and À2.47 respectively. COX2 plays a key role in prostanoid production by catalysing the conversion of arachidonic acid to prostaglandin G2 and H2 resulting in an inflammatory response. 53,54 Numerous studies have been conducted to investigate COX2 inhibition and the reduction in inflammatory response 55 and its role in oncogenesis. 56 Non-steroidal anti-inflammatory drugs are a class of drugs that inhibit COX2 enzyme activity and thereby reduce the inflammatory response. 57,58 This supports the notion that SN30028 has antiinflammatory effects.
Matrix metalloproteinase 3 (MMP3) was initially up-regulated in PC3 cells treated with SN30028 (Table 2 and Fig. S1, ESI †), but this response was not maintained. In addition, although MMP3 acts on a number of different genes, the expression of these genes was not modified when MMP3 was up-regulated ( Fig. S1, ESI †). Although MMP3 is over-expressed in most human cancers, and is known to induce initial cancer cell-growth and differentiation, rather than act at a later stage in cancer progression, 59,60 it is also known to have many opposing functions 59 and to date MMP inhibitors have not been successful in the clinic. 60 Differential expression of numerous other central node or core genes is evident, but none of the central node genes arise in more than two top ranked networks representing each cell line -treatment time combination. Some of the differentially expressed central node genes include FOXO1 and 3, thought to be involved in triggering apoptosis or cell survival and are induced by oxidative stress; 61 amphiregulin (AREG) was downregulated in PC3 ( Fig. S2 and S3, ESI †) supporting the idea that SN30028 reduces inflammation and inhibits tumour development. 62 AREG is a ligand of epidermal growth factor receptor (EGFR) 62 and the fact that EGFR is up-regulated in DU145 (Fig. S4, ESI †) could either be a chance occurrence as this finding was only attained in DU145 treated for 4 h and therefore is an initial response that is not sustained, or it could be that SN30028 works through different mechanisms in the two cell lines. Up-regulation of EGFR is associated with prostate cancer progression and EGFR dysfunction induces cell survival, proliferation, invasion and metastasis and therefore was not an anticipated response. The AR gene is known to interact with EGFR and we wouldn't expect to see a change in expression as DU145 and PC3 are androgen insensitive. However, although these cell lines are regarded as androgen non-responsive, some authors have reported low level expression of AR mRNA, and treatment with interferon (IFN) resulted in up-regulation of AR protein levels. 63 This may be due to AR phosphorylation. In the experiments reported herein a change in AR expression in response to treatment with SN30028 was not noted.
Similar to the unexpected up-regulation of EGFR, the ETSrelated gene (ERG) was up-regulated (at low intensity) after 96 h (Fig. S6, ESI †) and is an unexpected response. ERG is a protooncogene, regulates cell proliferation, differentiation, angiogenesis, inflammation, apoptosis and can result in gene fusion products associated with prostate and other cancers. 64 ERG, when overexpressed, is able to regulate oncogenic pathways involving cMyc, AR and EZH2, 64 none of which were up-regulated in this study. ERG acts on neuronally expressed developmentally downregulated 4 (NEDD4), 64 which is a central node gene that was down-regulated in DU145 cells (Fig. S6, ESI †). NEDD4 is an oncoprotein that promotes degradation through the ubiquitination of its substrates and it is also thought to promote colon and lung carcinogenesis, be over-expressed in prostate, breast and bladder cancers, and promotes growth of colon cancer cells independently of PTEN and PI3K/AKT signalling. 65 Despite the fact that ERG is over-expressed in DU145 after 96 h of treatment and acts on NEDD4 (amongst other genes), SN30028 appears to inhibit NEDD4 which is desirable, although further work is required for the elucidation of the mechanisms involved in treating cancers in this way.

Conclusions
In conclusion SN30028, an HDAC inhibitor was synthesised and used in the treatment of prostate cancer cell lines for 4, 24 and 96 h, and assessed for differential gene expression. Using IPA software the data were analysed and pathways involved in cholesterol biosynthesis and mevalonate were affected most commonly in the cell line DU145. In addition to pathway and network analysis the most frequently affected diseases and disorders were assessed. Those related to inflammatory response and/or cancers were one of the top five diseases and disorders listed for each cell line-time point tested. SN30028 showed it could be a potential therapeutic agent in treating prostate cancer by targeting the sterol and mevalonate pathways, however further work is required to confirm these results.

Funding
This work was funded by the Auckland Cancer Society Research Centre.