Abstract
The serine/threonine kinase Akt is a promising target in cancer. We previously identified five phosphatidylinositol ether lipid analogues (PIA) that inhibited Akt activation and selectively killed lung and breast cancer cells with high levels of Akt activity. To assess the spectrum of activity in other cell types and to compare PIAs with other inhibitors of the phosphatidylinositol 3-kinase (PI3K)/Akt/mammalian target of rapamycin (mTOR) pathway, we compared growth inhibition by PIAs against the PI3K inhibitors LY294002 and wortmannin and the mTOR inhibitor rapamycin in the NCI60 cell line panel. Although each of these compounds inhibited the growth of all the cell lines, distinct patterns were observed. The PIAs were the least potent but the most cytotoxic. The broad spectrum of activity of PIAs was confirmed in vivo in hollow fiber assays. The response to PIAs was significantly correlated with levels of active but not total Akt in the NCI60, as assessed using COMPARE analysis. However, a number of molecular targets were identified whose expression was more highly correlated with sensitivity to PIAs than active Akt. Expression of these molecular targets did not overlap with those that correlated with sensitivity to LY294002, wortmannin, or rapamycin. A COMPARE analysis of the National Cancer Institute chemical screening database revealed that the patterns of activity of PIAs correlated best with patterns of activity of other lipid-based compounds. These studies show that although PIAs are widely active in cancer cells, which correlates with the presence of its intended target, active Akt, PIAs are biologically distinct from other known inhibitors of the PI3K/Akt/mTOR pathway. [Mol Cancer Ther 2006;5(3):713–22]
Introduction
Akt is a serine/threonine kinase that controls cellular growth, migration, metabolism,5
J.J. Gills, S.S. Castillo, C. Zhang, P.A. Petukhov, R. Memmott, N. Warfel, J. Han, A.P. Kozikowski, P.A. Dennis, unpublished data.
The development of PIAs is one example of the many efforts in academia and industry to develop inhibitors of the PI3K/Akt/mammalian target of rapamycin (mTOR) pathway. Because many components in the pathway have been targeted, including PI3K, PDK-1, Akt, and mTOR, comparison studies will be required to reveal which inhibitors are truly targeted and which are most effective. In that vein, we assessed the spectrum of activity of PIAs in vitro and in vivo and compared responses to PIAs in the NCI60 cell line panel with that of two PI3K inhibitors (LY294002 and wortmannin) and an mTOR inhibitor (rapamycin). We show that cellular responses to PIAs indeed correlate with levels of active Akt. Using COMPARE analysis, we identify novel molecular targets that are highly correlated with response to PIAs and identify structurally related compounds that have similar patterns of activity.
Materials and Methods
NCI60 Cell Line Screen
Methods for evaluation of cell growth inhibition in the NCI60 cell line panel were published previously (3, 4). Briefly, PIAs were solubilized in DMSO, diluted into RPMI 1640 + 5% fetal bovine serum, and added to 96-well plates containing cell lines that were previously cultured for 24 hours. After a 48-hour incubation with the PIAs, the media were removed, and the cells were fixed and stained with sulforhodamine B. Unbound dye was removed with five washes of 1% acetic acid, and the plates were allowed to air dry. The dye was resolubilized in Tris buffer, and the absorbance at 515 nm was measured. The concentration that produced 50% growth inhibition compared with a DMSO control (GI50), total growth inhibition, or 0% growth, compared with a DMSO control (TGI), and the concentration that produced death of 50% of the cells present at the start of the experiment (LC50) were determined.
Hollow Fiber Assays
Cells were grown inside polyvinylidene difluoride fibers and inserted into the i.p. and s.c. compartments of nude mice as described previously (5). PIA23 was dissolved in 10% DMSO in saline and given i.p. every day for 4 days, with dosing beginning on the third day after fiber implantation at doses of 25 and 37.5 mg/kg.
Correlation of Sensitivity of NCI60 Cell Line Panel to PIAs with Levels of Phosphorylated and Total Akt
Levels of phosphorylated Akt (phospho-Akt) at S473 and T308, as well as total Akt protein were measured previously in the NCI60 cell line panel by immunoblotting, and these data are available through the Developmental Therapeutics Program web site.6
The COMPARE algorithm (6) was used to determine Pearson correlation coefficients between the −log(GI50) of the PIAs and the measured levels of total Akt and phospho-Akt in the cell lines.Correlation of Sensitivity of the PIAs to Other Molecular Targets in the 60 Cell Line Panel
To date the Developmental Therapeutics Program has assembled over 220,000 molecular target measurements in the NCI60 cell line panel. The COMPARE algorithm and the JMP statistical software package (SAS Institute, Cary, NC) was used to determine Pearson correlation coefficients between the −log(GI50) of the PIAs and the levels of molecular targets that have been measured in the NCI60 cell line panel. Given the uncertainty in an individual measurement from a single microarray, we first identified a subset of genes whose measurement on multiple microarray platforms were well correlated (Pearson correlation coefficient ≥0.5). This was done using microarray data sets available through the Developmental Therapeutics Program web site and included four arrays: a cDNA array (7, 8), an Affymetrix HUM6000 array (from Millenium Pharmaceuticals, Cambridge, MA), an Affymetrix U95A array done in triplicate (Novartis, Cambridge, MA), and an Affymetrix U95A-E array (Gene Logic, Gaithersburg, MD). Using this subset of 5,796 genes (∼20,000 total measurements), we did correlation analyses with the growth inhibition patterns of these compounds. The results were further filtered to include only those genes where measurements on multiple microarrays gave significant correlations with growth inhibition by the PIAs.
Correlation of Sensitivity of the PIAs to Other Compounds That Have Been Tested in the 60 Cell Line Panel
As of 2005, the National Cancer Institute has evaluated compounds for anticancer activity, with publicly available data for ∼43,000 compounds. The COMPARE algorithm was used to search for compounds that have been tested in the NCI60 cell line panel that have a similar sensitivity profile to the PIAs.
Results
Spectrum of Activity
To assess the spectrum of cancer cell types that are sensitive to PIAs, we screened PIA5, PIA6, PIA23, PIA24, and PIA25 (NSC726850, NSC726851, NSC726852, NSC726853, and NSC726854) in the NCI60 cell line panel. All PIAs were effective in inhibiting the growth of all of the cell lines, with the GI50s for most falling between 1 and 20 μmol/L. The composite dose-response curves were similar for all five PIAs (Fig. 1, left and middle). For the majority of the cell lines, PIAs begin to inhibit growth at 1 μmol/L. Above 10 μmol/L, the growth of all cell lines is sharply inhibited. PIAs are uniformly cytotoxic at 100 μmol/L. To determine if this was characteristic of other inhibitors of the pathway, we compared the sensitivity of the PIAs that inhibit Akt with three other inhibitors of the pathway: LY294002 (NSC697286) and wortmannin (NSC221019), both of which are PI3K inhibitors, and rapamycin (NSC226080), which is an mTOR inhibitor (Fig. 1, right). LY294002 and wortmannin have more gradually shaped dose-response growth curves than the PIAs. Rapamycin, unlike PIAs, LY294002, or wortmannin, is a very potent cytostatic agent (10 nmol/L to 10 μmol/L). At doses between 10 and 100 μmol/L, however, rapamycin is cytotoxic. These studies show that PIAs are less potent than the PI3K inhibitors or rapamycin. At the highest doses tested, PIAs are widely cytotoxic like rapamycin and more cytotoxic than LY294002 or wortmannin.
Growth curves of PIAs, LY294004, wortmannin, and rapamycin in the NCI60 cell line panel. Overlays of the percent growth of all 60 cell lines following a 48-h incubation with varying doses of PIAs, LY294004, wortmannin, and rapamycin are shown. Dose-response curves were generated as described in Materials and Methods.
Growth curves of PIAs, LY294004, wortmannin, and rapamycin in the NCI60 cell line panel. Overlays of the percent growth of all 60 cell lines following a 48-h incubation with varying doses of PIAs, LY294004, wortmannin, and rapamycin are shown. Dose-response curves were generated as described in Materials and Methods.
To determine if differences could be discerned between individual PIAs, we compared the efficacy of each PIA in each cell line (Fig. 2). Overall, the mean graphs of the GI50 concentrations of the five PIAs are similar, except that PIA24 seemed less efficacious. (The TGI and LC50 data for PIAs are shown in Supplementary Figs. S1 and S2, respectively.)7
Supplementary material for this article are available at Molecular Cancer Therapeutics Online (http://mct.aacrjournals.org/).
Mean graphs displaying the GI50 of PIA5, PIA6, PIA23, PIA24, and PIA25 in the NCI60 cell line panel. The central line is the arithmetic mean of the GI50 of all the cells lines. Bars to the right or left indicate greater or lesser sensitivity to PIAs than the average of all the cell lines.
Mean graphs displaying the GI50 of PIA5, PIA6, PIA23, PIA24, and PIA25 in the NCI60 cell line panel. The central line is the arithmetic mean of the GI50 of all the cells lines. Bars to the right or left indicate greater or lesser sensitivity to PIAs than the average of all the cell lines.
To assess the spectrum of activity of PIAs in vivo, we assayed an active PIA (PIA23) in hollow fiber assays that use 12 of the 60 cell lines in the NCI60. PIA23 was effective in inhibiting the growth of 8 of 12 of the cell lines grown inside hollow fibers placed in either the i.p. or the s.c. compartments, and cytotoxicity was noted (Table 1). PIA23 inhibited growth in an equal number of fibers located in the s.c. and i.p. compartments. This is notable because PIA23 was given by i.p. injection, which shows that PIA23 was able to reach effective circulating concentrations.
In vivo efficacy of PIA23 in hollow fiber assays
. | i.p. . | s.c. . |
---|---|---|
H522 | — | 2 |
UACC-62 | — | — |
U251 | — | — |
H23 | 2 | 4 |
MDA-MB-231 | — | — |
SW620 | 2 | 2 |
MDA-MB-435 | 4 | 4 |
OVCAR-5 | — | 2 |
SF-295 | 4 | 2 |
LOX | 2 | — |
COLO205 | — | — |
OVCAR3 | 2 | — |
Total | 16 | 16 |
. | i.p. . | s.c. . |
---|---|---|
H522 | — | 2 |
UACC-62 | — | — |
U251 | — | — |
H23 | 2 | 4 |
MDA-MB-231 | — | — |
SW620 | 2 | 2 |
MDA-MB-435 | 4 | 4 |
OVCAR-5 | — | 2 |
SF-295 | 4 | 2 |
LOX | 2 | — |
COLO205 | — | — |
OVCAR3 | 2 | — |
Total | 16 | 16 |
NOTE: A score of 2 indicates that growth was inhibited ≥50% within the fiber for a particular dose. The number 4 indicates it decreased growth ≥50% in hollow fibers in mice that received both the 25 and 37.5 mg/kg dose levels.
We then compared responsiveness in the hollow fiber assays with levels of Akt activation. The most sensitive cell line in the assay was the cell line MDA-MB-435, formerly believed to be a breast cancer cell line but recently was shown to be virtually identical to the melanoma line M14 (9). MDA-MB-435 cells have been previously shown to contain a high level of endogenous Akt activity (10). Conversely, MDA-MB-231 cells, which did not respond to PIA23, have low levels of active Akt and do not respond to LY294002 (10). The correlation with active Akt was not absolute, because the cell line with the most phospho-Akt on the panel, U251, did not respond to PIA23, whereas a cell line with a relatively low amount of phospho-Akt, NCI-H522, did respond. This shows that response to PIA treatment is not solely dependent on Akt and may rely on expression of other targets.
Correlation of PIA Activity with Akt Status
To further evaluate whether the cytotoxicity of the PIAs is related to modulation of endogenous Akt activity, which was suggested by our earlier studies, the sensitivity to PIAs (GI50) was compared with the levels of active Akt (as defined by phosphorylation of S473 or T308), and total Akt in the NCI60 cell line panel (Table 2). The Pearson correlation coefficients (PCC) that compare response to PIAs with S473 phosphorylation were marginally statistically significant for all PIAs (PCC ≥0.28). The PCC for T308 was marginally statistically significant for all but PIA24. In contrast, the PCCs for total Akt were below zero. This confirms that sensitivity of the PIAs correlates positively with the presence of active Akt but not total Akt.
Comparison of sensitivity of PIAs to levels of phosphorylated and total Akt in the NCI60 cell line panel
. | PS473 . | PT308 . | Akt . |
---|---|---|---|
PIA5 | 0.31 | 0.34 | −0.26 |
PIA6 | 0.28 | 0.31 | −0.31 |
PIA23 | 0.31 | 0.29 | −0.09 |
PIA24 | 0.3 | 0.24 | −0.1 |
PIA25 | 0.32 | 0.29 | −0.29 |
. | PS473 . | PT308 . | Akt . |
---|---|---|---|
PIA5 | 0.31 | 0.34 | −0.26 |
PIA6 | 0.28 | 0.31 | −0.31 |
PIA23 | 0.31 | 0.29 | −0.09 |
PIA24 | 0.3 | 0.24 | −0.1 |
PIA25 | 0.32 | 0.29 | −0.29 |
NOTE: The COMPARE program was used to generate Pearson correlation coefficients between the −log(GI50) of the PIAs and the levels of phospho-Akt (S473 and T308) and total Akt present in the cell lines as measured by immunoblot analyses. PCCs range from a value of 1.0 for a perfect match, to −1.0 for a perfect inverse match, with a value of 0, indicating no correlation between measurements.
Correlation of PIA Activity with Expression of Other Molecular Targets
To determine whether expression of other molecular targets might correlate with sensitivity to PIAs, we did analyses using a subset of the public microarray data. Multiple molecular targets were identified whose expression correlated positively with response to PIAs, indicating that cell lines with higher levels of expression of these genes tended to be more sensitive to the compounds. The 10 highest scoring genes with PCCs ≥ 0.4 are rank listed in Table 3. The highest correlations for each PIA were higher than that observed with active Akt. The highest PCC for any target was 0.67 for the association between the −log(GI50) of PIA24 and syntaxin 1A. Similar analysis was done for molecular targets whose expression negatively correlated with sensitivity to PIAs, and multiple molecular targets were identified whose PCCs were less than or equal to −0.4 (Table 4). Overall, there were more positively correlated targets that were identified than negatively correlated ones. For example, the only negative correlating gene for PIA24 was STAT3. Although no overarching categorization was obvious, several genes that regulate chromatin (HDAC1-positively correlated, SMARCA3-negatively correlated) and transcription (KLF4, NSEP1, c-myc, SOX9, and ETS2-positively correlated, AP-2α-negatively correlated) were identified.
Molecular targets whose levels in the NCI60 cell line panel correlate most positively with PIA sensitivity
Rank | PCC | Gene | Description | |
PIA5 | 1 | 0.60 | KLF4 | Kruppel-like factor 4 (gut) |
2 | 0.57 | DPP4 | Dipeptidylpeptidase 4 (CD26, adenosine deaminase complexing protein 2) | |
3 | 0.54 | ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 (avian) | |
4 | 0.53 | SLCO4A1 | Organic anion transporter SLC21A12 | |
5 | 0.53 | MAN2A1 | Mannosidase, alpha, class 2A, member 1 | |
6 | 0.49 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
7 | 0.49 | HDAC1 | Histone deacetylase 1 | |
8 | 0.48 | CLDN4 | Claudin 4 | |
9 | 0.47 | RGS2 | Regulator of G-protein signaling 2, 24 kDa | |
10 | 0.47 | AIM1 | Absent in melanoma 1 | |
PIA6 | 1 | 0.62 | KLF4 | Kruppel-like factor 4 (gut) |
2 | 0.55 | SLCO4A1 | Organic anion transporter SLC21A12 | |
3 | 0.54 | MAN2A1 | Mannosidase, alpha, class 2A, member 1 | |
4 | 0.52 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
5 | 0.50 | EBP | Emopamil binding protein (sterol isomerase) | |
6 | 0.50 | HDAC1 | Histone deacetylase 1 | |
7 | 0.49 | DPP4 | Dipeptidylpeptidase 4 (CD26, adenosine deaminase complexing protein 2) | |
8 | 0.49 | PPARG | Peroxisome proliferative activated receptor, γ | |
9 | 0.49 | HSPE1 | Heat shock 10-kDa protein 1 (chaperonin 10) | |
10 | 0.48 | SOX9 | SRY (sex-determining region Y)-box 9 (campomelic dysplasia, autosomal sex-reversal) | |
PIA23 | 1 | 0.49 | LRP8 | Low-density lipoprotein receptor-related protein 8, apolipoprotein e receptor |
2 | 0.48 | ATRX | Alpha thalassemia/mental retardation syndrome X-linked (RAD54 homologue) | |
3 | 0.47 | GUK1 | Guanylate kinase 1 | |
4 | 0.47 | NUPL2 | Nucleoporin like 2 | |
5 | 0.46 | PMAIP1 | Phorbol-12-myristate-13-acetate-induced protein 1 | |
6 | 0.46 | PTK2B | PYK2 (protein tyrosine kinase) | |
7 | 0.45 | EPHA3 | EphA3 (ephrin receptor, tyrosine kinase) | |
8 | 0.45 | TXK | TXK (protein tyrosine kinase) | |
9 | 0.45 | HAPLN1 | Hyaluronan and proteoglycan link protein 1 | |
10 | 0.45 | AIM1 | Absent in melanoma 1 | |
PIA24 | 1 | 0.67 | STX1A | Syntaxin 1A (brain) |
2 | 0.59 | HAPLN1 | Hyaluronan and proteoglycan link protein 1 | |
3 | 0.58 | PPIF | Peptidylprolyl isomerase F (cyclophilin F) | |
4 | 0.58 | ANXA7 | Annexin A7 | |
5 | 0.57 | GUK1 | Guanylate kinase 1 | |
6 | 0.56 | CDKN2B | Cyclin-dependent kinase inhibitor 2B (p15, inhibits CDK4) | |
7 | 0.55 | HOXB13 | Homeobox B13 | |
8 | 0.55 | RGS2 | Regulator of G-protein signaling 2, 24 kDa | |
9 | 0.55 | NEFL | Neurofilament, light polypeptide 68 kDa | |
10 | 0.55 | KLF4 | Kruppel-like factor 4 (gut) | |
PIA25 | 1 | 0.51 | NUPL2 | Nucleoporin like 2 |
2 | 0.51 | HDAC1 | Histone deacetylase 1 | |
3 | 0.50 | AIM1 | Absent in melanoma 1 | |
4 | 0.50 | ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 (avian) | |
5 | 0.49 | SLCO4A1 | Organic anion transporter SLC21A12 | |
6 | 0.48 | RPS24 | Ribosomal protein S24 | |
7 | 0.48 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
8 | 0.48 | KLF4 | Kruppel-like factor 4 (gut) | |
9 | 0.47 | MYC | c-Myc (basic helix-loop-helix transcription factor family) | |
10 | 0.47 | AKAP1 | A kinase (PRKA) anchor protein 1 |
Rank | PCC | Gene | Description | |
PIA5 | 1 | 0.60 | KLF4 | Kruppel-like factor 4 (gut) |
2 | 0.57 | DPP4 | Dipeptidylpeptidase 4 (CD26, adenosine deaminase complexing protein 2) | |
3 | 0.54 | ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 (avian) | |
4 | 0.53 | SLCO4A1 | Organic anion transporter SLC21A12 | |
5 | 0.53 | MAN2A1 | Mannosidase, alpha, class 2A, member 1 | |
6 | 0.49 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
7 | 0.49 | HDAC1 | Histone deacetylase 1 | |
8 | 0.48 | CLDN4 | Claudin 4 | |
9 | 0.47 | RGS2 | Regulator of G-protein signaling 2, 24 kDa | |
10 | 0.47 | AIM1 | Absent in melanoma 1 | |
PIA6 | 1 | 0.62 | KLF4 | Kruppel-like factor 4 (gut) |
2 | 0.55 | SLCO4A1 | Organic anion transporter SLC21A12 | |
3 | 0.54 | MAN2A1 | Mannosidase, alpha, class 2A, member 1 | |
4 | 0.52 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
5 | 0.50 | EBP | Emopamil binding protein (sterol isomerase) | |
6 | 0.50 | HDAC1 | Histone deacetylase 1 | |
7 | 0.49 | DPP4 | Dipeptidylpeptidase 4 (CD26, adenosine deaminase complexing protein 2) | |
8 | 0.49 | PPARG | Peroxisome proliferative activated receptor, γ | |
9 | 0.49 | HSPE1 | Heat shock 10-kDa protein 1 (chaperonin 10) | |
10 | 0.48 | SOX9 | SRY (sex-determining region Y)-box 9 (campomelic dysplasia, autosomal sex-reversal) | |
PIA23 | 1 | 0.49 | LRP8 | Low-density lipoprotein receptor-related protein 8, apolipoprotein e receptor |
2 | 0.48 | ATRX | Alpha thalassemia/mental retardation syndrome X-linked (RAD54 homologue) | |
3 | 0.47 | GUK1 | Guanylate kinase 1 | |
4 | 0.47 | NUPL2 | Nucleoporin like 2 | |
5 | 0.46 | PMAIP1 | Phorbol-12-myristate-13-acetate-induced protein 1 | |
6 | 0.46 | PTK2B | PYK2 (protein tyrosine kinase) | |
7 | 0.45 | EPHA3 | EphA3 (ephrin receptor, tyrosine kinase) | |
8 | 0.45 | TXK | TXK (protein tyrosine kinase) | |
9 | 0.45 | HAPLN1 | Hyaluronan and proteoglycan link protein 1 | |
10 | 0.45 | AIM1 | Absent in melanoma 1 | |
PIA24 | 1 | 0.67 | STX1A | Syntaxin 1A (brain) |
2 | 0.59 | HAPLN1 | Hyaluronan and proteoglycan link protein 1 | |
3 | 0.58 | PPIF | Peptidylprolyl isomerase F (cyclophilin F) | |
4 | 0.58 | ANXA7 | Annexin A7 | |
5 | 0.57 | GUK1 | Guanylate kinase 1 | |
6 | 0.56 | CDKN2B | Cyclin-dependent kinase inhibitor 2B (p15, inhibits CDK4) | |
7 | 0.55 | HOXB13 | Homeobox B13 | |
8 | 0.55 | RGS2 | Regulator of G-protein signaling 2, 24 kDa | |
9 | 0.55 | NEFL | Neurofilament, light polypeptide 68 kDa | |
10 | 0.55 | KLF4 | Kruppel-like factor 4 (gut) | |
PIA25 | 1 | 0.51 | NUPL2 | Nucleoporin like 2 |
2 | 0.51 | HDAC1 | Histone deacetylase 1 | |
3 | 0.50 | AIM1 | Absent in melanoma 1 | |
4 | 0.50 | ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 (avian) | |
5 | 0.49 | SLCO4A1 | Organic anion transporter SLC21A12 | |
6 | 0.48 | RPS24 | Ribosomal protein S24 | |
7 | 0.48 | PRSS3 | Protease, serine, 3 (mesotrypsin) | |
8 | 0.48 | KLF4 | Kruppel-like factor 4 (gut) | |
9 | 0.47 | MYC | c-Myc (basic helix-loop-helix transcription factor family) | |
10 | 0.47 | AKAP1 | A kinase (PRKA) anchor protein 1 |
NOTE: The COMPARE program was used to generate a list of molecular targets whose levels correlate with the −log(GI50) of the PIAs. Molecular targets with the top 10 highest positive PCCs are shown for each PIA.
Molecular targets whose levels in the NCI60 cell line panel correlate most negatively with PIA sensitivity
. | Rank . | PCC . | Gene . | Description . |
---|---|---|---|---|
PIA5 | 1 | −0.50 | SMARCA3 | SWI/SNF–related, actin-dependent regulator of chromatin, subfamily a, member 3 |
2 | −0.50 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
3 | −0.49 | HRMT1L1 | HMT1 hnRNP methyltransferase-like 1 (Saccharomyces cerevisiae) | |
4 | −0.48 | MPDZ | Multiple PDZ domain protein | |
5 | −0.46 | CCNE | Cyclin E | |
6 | −0.46 | PRAME | Preferentially expressed antigen in melanoma | |
7 | −0.46 | PHF10 | PHD finger protein 10 | |
8 | −0.46 | MYH10 | Myosin, heavy polypeptide 10, nonmuscle | |
9 | −0.46 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
10 | −0.45 | CHST10 | Carbohydrate sulfotransferase 10 | |
PIA6 | 1 | −0.50 | CCNE | Cyclin E |
2 | −0.49 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) | |
3 | −0.47 | HRMT1L1 | HMT1 hnRNP methyltransferase-like 1 (S. cerevisiae) | |
4 | −0.47 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
5 | −0.47 | SMARCA3 | SWI/SNF–related, actin-dependent regulator of chromatin, subfamily a, member 3 | |
6 | −0.44 | APLP1 | Amyloid beta (A4) precursor-like protein 1 | |
7 | −0.43 | USP6 | Ubiquitin-specific protease 6 (Tre-2 oncogene) | |
8 | −0.43 | VAT1 | Vesicle amine transport protein 1 homologue (Torreya californica) | |
9 | −0.43 | USP32 | Ubiquitin-specific protease 32 | |
10 | −0.43 | PRAME | Preferentially expressed antigen in melanoma | |
PIA23 | 1 | −0.56 | UBL3 | Ubiquitin-like 3 |
2 | −0.54 | CCNE | Cyclin E | |
3 | −0.53 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
4 | −0.50 | RIPK1 | Receptor (TNFRSF)-interacting serine-threonine kinase 1 | |
5 | −0.49 | PTTG1IP | Pituitary tumor-transforming 1 interacting protein | |
6 | −0.47 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) | |
7 | −0.44 | LAPTM4B | Lysosomal-associated protein transmembrane 4 beta | |
8 | −0.43 | ERBB3 | V-erb-b2 erythroblastic leukemia viral oncogene homologue 3 (avian) | |
9 | −0.43 | PRAME | Preferentially expressed antigen in melanoma | |
10 | −0.43 | ELOVL2 | Elongation of very long chain fatty acids (FEN1/Elo2, SUR4/Elo3, yeast)-like 2 | |
PIA24 | 1 | −0.44 | STAT3 | Signal transducer and activator of transcription 3 |
PIA25 | 1 | −0.61 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) |
2 | −0.57 | CCNE | Cyclin E | |
3 | −0.54 | MPDZ | Multiple PDZ domain protein | |
4 | −0.53 | PRAME | Preferentially expressed antigen in melanoma | |
5 | −0.51 | MYH10 | Myosin, heavy polypeptide 10, nonmuscle | |
6 | −0.51 | DRPLA | Dentatorubral-pallidoluysian atrophy (atrophin-1) | |
7 | −0.50 | PPT2 | Palmitoyl-protein thioesterase 2 | |
8 | −0.50 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
9 | −0.50 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
10 | −0.48 | UBL3 | Ubiquitin-like 3 |
. | Rank . | PCC . | Gene . | Description . |
---|---|---|---|---|
PIA5 | 1 | −0.50 | SMARCA3 | SWI/SNF–related, actin-dependent regulator of chromatin, subfamily a, member 3 |
2 | −0.50 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
3 | −0.49 | HRMT1L1 | HMT1 hnRNP methyltransferase-like 1 (Saccharomyces cerevisiae) | |
4 | −0.48 | MPDZ | Multiple PDZ domain protein | |
5 | −0.46 | CCNE | Cyclin E | |
6 | −0.46 | PRAME | Preferentially expressed antigen in melanoma | |
7 | −0.46 | PHF10 | PHD finger protein 10 | |
8 | −0.46 | MYH10 | Myosin, heavy polypeptide 10, nonmuscle | |
9 | −0.46 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
10 | −0.45 | CHST10 | Carbohydrate sulfotransferase 10 | |
PIA6 | 1 | −0.50 | CCNE | Cyclin E |
2 | −0.49 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) | |
3 | −0.47 | HRMT1L1 | HMT1 hnRNP methyltransferase-like 1 (S. cerevisiae) | |
4 | −0.47 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
5 | −0.47 | SMARCA3 | SWI/SNF–related, actin-dependent regulator of chromatin, subfamily a, member 3 | |
6 | −0.44 | APLP1 | Amyloid beta (A4) precursor-like protein 1 | |
7 | −0.43 | USP6 | Ubiquitin-specific protease 6 (Tre-2 oncogene) | |
8 | −0.43 | VAT1 | Vesicle amine transport protein 1 homologue (Torreya californica) | |
9 | −0.43 | USP32 | Ubiquitin-specific protease 32 | |
10 | −0.43 | PRAME | Preferentially expressed antigen in melanoma | |
PIA23 | 1 | −0.56 | UBL3 | Ubiquitin-like 3 |
2 | −0.54 | CCNE | Cyclin E | |
3 | −0.53 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
4 | −0.50 | RIPK1 | Receptor (TNFRSF)-interacting serine-threonine kinase 1 | |
5 | −0.49 | PTTG1IP | Pituitary tumor-transforming 1 interacting protein | |
6 | −0.47 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) | |
7 | −0.44 | LAPTM4B | Lysosomal-associated protein transmembrane 4 beta | |
8 | −0.43 | ERBB3 | V-erb-b2 erythroblastic leukemia viral oncogene homologue 3 (avian) | |
9 | −0.43 | PRAME | Preferentially expressed antigen in melanoma | |
10 | −0.43 | ELOVL2 | Elongation of very long chain fatty acids (FEN1/Elo2, SUR4/Elo3, yeast)-like 2 | |
PIA24 | 1 | −0.44 | STAT3 | Signal transducer and activator of transcription 3 |
PIA25 | 1 | −0.61 | TFAP2A | Transcription factor AP-2 alpha (activating enhancer binding protein 2 alpha) |
2 | −0.57 | CCNE | Cyclin E | |
3 | −0.54 | MPDZ | Multiple PDZ domain protein | |
4 | −0.53 | PRAME | Preferentially expressed antigen in melanoma | |
5 | −0.51 | MYH10 | Myosin, heavy polypeptide 10, nonmuscle | |
6 | −0.51 | DRPLA | Dentatorubral-pallidoluysian atrophy (atrophin-1) | |
7 | −0.50 | PPT2 | Palmitoyl-protein thioesterase 2 | |
8 | −0.50 | DZIP3 | Zinc finger DAZ interacting protein 3 | |
9 | −0.50 | ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | |
10 | −0.48 | UBL3 | Ubiquitin-like 3 |
NOTE: The COMPARE program was used to generate a list of molecular targets whose levels correlate with the −log(GI50) of the PIAs. Shown are the molecular targets with 10 most negative PCCs for each PIA.
Because the active PIAs are structurally related and kill cancer cells similarly, we sought to identify common targets that correlated with response to PIAs. We compared molecular targets that had PCC ≥ 0.4 and identified those that were common to three or more PIAs (Table 5). Fourteen targets positively correlated with sensitivity to multiple PIAs. Three of 14 positively correlated targets (mesotrypsin, c-myc, and guanylate kinase 1) were common across all five PIAs. Although there were no negatively correlated targets that were in common with PIA24, five targets were common across the other four PIAs (cyclin E, ErbB3, palmitoyl-protein thioesterase 2, preferentially expressed antigen in melanoma, and myosin heavy polypeptide 10; Table 6). Of the group of five PIAs, PIA5, PIA6, and PIA25 shared the most targets in common, whereas PIA24 had the least number of targets in common with any of the other PIAs. These data show that other targets correlate better with responsiveness than active Akt and suggest that PIA24 has different markers for responsiveness than other PIAs.
Common molecular targets whose levels positively correlate with sensitivity to PIAs
Gene . | . | Pearson correlation coefficient . | . | . | . | . | ||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | PIA5 . | PIA6 . | PIA23 . | PIA24 . | PIA25 . | ||||
KLF4 | Kruppel-like factor 4 (gut) | 0.60 | 0.62 | 0.55 | 0.48 | |||||
DPP4 | Dipeptidylpeptidase 4 | 0.57 | 0.49 | 0.46 | ||||||
SLCO4A1 | Organic anion transporter SLC21A12 | 0.53 | 0.55 | 0.51 | ||||||
PRSS3 | Protease, serine, 3 (mesotrypsin) | 0.49 | 0.52 | 0.41 | 0.52 | 0.48 | ||||
RGS2 | Regulator of G-protein signaling 2, 24 kDa | 0.47 | 0.47 | 0.55 | 0.41 | |||||
AIM1 | Absent in melanoma 1 | 0.47 | 0.44 | 0.45 | 0.50 | |||||
AKAP1 | A kinase (PRKA) anchor protein 1 | 0.44 | 0.48 | 0.47 | ||||||
MYC | c-Myc | 0.42 | 0.44 | 0.45 | 0.47 | 0.47 | ||||
HDAC1 | Histone deacetylase 1 | 0.50 | 0.40 | 0.51 | ||||||
ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 | 0.45 | 0.48 | 0.50 | ||||||
MAN2A1 | Mannosidase, alpha, class 2A, member 1 | 0.53 | 0.54 | 0.47 | ||||||
SOX9 | SRY (sex determining region Y)-box 9 | 0.45 | 0.48 | 0.45 | ||||||
GUK1 | Guanylate kinase 1 | 0.44 | 0.41 | 0.47 | 0.57 | 0.41 | ||||
NSEP1 | Nuclease sensitive element binding protein 1 | 0.42 | 0.41 | 0.46 |
Gene . | . | Pearson correlation coefficient . | . | . | . | . | ||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | PIA5 . | PIA6 . | PIA23 . | PIA24 . | PIA25 . | ||||
KLF4 | Kruppel-like factor 4 (gut) | 0.60 | 0.62 | 0.55 | 0.48 | |||||
DPP4 | Dipeptidylpeptidase 4 | 0.57 | 0.49 | 0.46 | ||||||
SLCO4A1 | Organic anion transporter SLC21A12 | 0.53 | 0.55 | 0.51 | ||||||
PRSS3 | Protease, serine, 3 (mesotrypsin) | 0.49 | 0.52 | 0.41 | 0.52 | 0.48 | ||||
RGS2 | Regulator of G-protein signaling 2, 24 kDa | 0.47 | 0.47 | 0.55 | 0.41 | |||||
AIM1 | Absent in melanoma 1 | 0.47 | 0.44 | 0.45 | 0.50 | |||||
AKAP1 | A kinase (PRKA) anchor protein 1 | 0.44 | 0.48 | 0.47 | ||||||
MYC | c-Myc | 0.42 | 0.44 | 0.45 | 0.47 | 0.47 | ||||
HDAC1 | Histone deacetylase 1 | 0.50 | 0.40 | 0.51 | ||||||
ETS2 | V-ets erythroblastosis virus E26 oncogene homologue 2 | 0.45 | 0.48 | 0.50 | ||||||
MAN2A1 | Mannosidase, alpha, class 2A, member 1 | 0.53 | 0.54 | 0.47 | ||||||
SOX9 | SRY (sex determining region Y)-box 9 | 0.45 | 0.48 | 0.45 | ||||||
GUK1 | Guanylate kinase 1 | 0.44 | 0.41 | 0.47 | 0.57 | 0.41 | ||||
NSEP1 | Nuclease sensitive element binding protein 1 | 0.42 | 0.41 | 0.46 |
NOTE: Included are positively correlated targets that were common across three or more PIAs (PCC more than or equal to 0.4).
Common molecular targets whose levels negatively correlate with sensitivity to PIAs
Gene . | . | Pearson correlation coefficient . | . | . | . | . | ||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | PIA5 . | PIA6 . | PIA23 . | PIA24 . | PIA25 . | ||||
CCNE | Cyclin E | −0.46 | −0.50 | −0.54 | −0.57 | |||||
ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | −0.46 | −0.42 | −0.53 | −0.50 | |||||
TFAP2A | Transcription factor AP-2 alpha | −0.49 | −0.47 | −0.61 | ||||||
DZIP3 | Zinc finger DAZ interacting protein 3 | −0.50 | −0.47 | −0.50 | ||||||
MPDZ | Multiple PDZ domain protein | −0.48 | −0.42 | −0.54 | ||||||
PHF10 | PHD finger protein 10 | −0.46 | −0.40 | −0.40 | ||||||
UBL3 | Ubiquitin-like 3 | −0.40 | −0.56 | −0.48 | ||||||
PPT2 | Palmitoyl-protein thioesterase 2 | −0.41 | −0.42 | −0.41 | −0.50 | |||||
RHOQ | Ras homologue gene family, member Q | −0.41 | −0.40 | −0.43 | ||||||
PRAME | Preferentially expressed antigen in melanoma | −0.46 | −0.43 | −0.43 | −0.53 | |||||
CHST10 | Carbohydrate sulfotransferase 10 | −0.45 | −0.41 | −0.44 | ||||||
MYH10 | Myosin, heavy polypeptide 10, nonmuscle | −0.46 | −0.43 | −0.41 | −0.51 | |||||
VAT1 | Vesicle amine transport protein 1 homologue (T. californica) | −0.44 | −0.43 | −0.47 | ||||||
APLP1 | Amyloid beta (A4) precursor-like protein 1 | −0.43 | −0.44 | −0.41 | ||||||
USP6 | Ubiquitin-specific protease 6 (Tre-2 oncogene) | −0.42 | −0.43 | −0.47 | ||||||
LAPTM4B | Lysosomal-associated protein transmembrane 4 beta | −0.40 | −0.44 | −0.47 |
Gene . | . | Pearson correlation coefficient . | . | . | . | . | ||||
---|---|---|---|---|---|---|---|---|---|---|
. | . | PIA5 . | PIA6 . | PIA23 . | PIA24 . | PIA25 . | ||||
CCNE | Cyclin E | −0.46 | −0.50 | −0.54 | −0.57 | |||||
ERBB3 | ErbB3 (receptor protein-tyrosine kinase) | −0.46 | −0.42 | −0.53 | −0.50 | |||||
TFAP2A | Transcription factor AP-2 alpha | −0.49 | −0.47 | −0.61 | ||||||
DZIP3 | Zinc finger DAZ interacting protein 3 | −0.50 | −0.47 | −0.50 | ||||||
MPDZ | Multiple PDZ domain protein | −0.48 | −0.42 | −0.54 | ||||||
PHF10 | PHD finger protein 10 | −0.46 | −0.40 | −0.40 | ||||||
UBL3 | Ubiquitin-like 3 | −0.40 | −0.56 | −0.48 | ||||||
PPT2 | Palmitoyl-protein thioesterase 2 | −0.41 | −0.42 | −0.41 | −0.50 | |||||
RHOQ | Ras homologue gene family, member Q | −0.41 | −0.40 | −0.43 | ||||||
PRAME | Preferentially expressed antigen in melanoma | −0.46 | −0.43 | −0.43 | −0.53 | |||||
CHST10 | Carbohydrate sulfotransferase 10 | −0.45 | −0.41 | −0.44 | ||||||
MYH10 | Myosin, heavy polypeptide 10, nonmuscle | −0.46 | −0.43 | −0.41 | −0.51 | |||||
VAT1 | Vesicle amine transport protein 1 homologue (T. californica) | −0.44 | −0.43 | −0.47 | ||||||
APLP1 | Amyloid beta (A4) precursor-like protein 1 | −0.43 | −0.44 | −0.41 | ||||||
USP6 | Ubiquitin-specific protease 6 (Tre-2 oncogene) | −0.42 | −0.43 | −0.47 | ||||||
LAPTM4B | Lysosomal-associated protein transmembrane 4 beta | −0.40 | −0.44 | −0.47 |
NOTE: Included are negatively correlated targets that were common to at least three PIAs (PCC less than or equal to −0.4).
Correlation of Activity of PIAs with Other Compounds
Because of the apparent distinction between PIAs and other inhibitors of the PI3K/Akt/mTOR pathway, we used COMPARE to search the NCI60 database to determine if other compounds have shown similar sensitivity profiles in the NCI60 cell line panel. When PIAs were compared against each other, PIA5, PIA6, PIA23, and PIA25 exhibited high correlation coefficients (0.74 < PCC < 0.91), indicating their similarity (data not shown). When PIA24 was compared against the other PIAs, the highest PCC was 0.67 for the correlation with PIA6 (data not shown). When PIA5, PIA6, PIA23, or PIA25 was used as a seed against the entire database, structurally similar compounds were identified (Fig. 3A-E). The compound whose GI50 pattern showed the most similarity to the PIAs was NSC643826. This had a coefficient of 0.78 compared with PIA5, 0.70 with PIA6, 0.69 with PIA23, and 0.76 with PIA25. This and many of the other top matching compounds are phospholipids, like the PIAs. Some of these compounds have an amino group attaching the phosphate group with the lipid side chain, whereas the PIAs have an ether linkage attaching the side chain. Unlike the PIAs, which have an inositol head group, all of the phospholipids that were identified contained nitrogen in the head group (many within a choline moiety). Two of the PIAs (PIA5 and PIA6) had similarity to the alkylphospholipids (miltefosine or perifosine) in their top matches. PIA24 did not show similarity to any phospholipids in its top five matches (Fig. 3D). Although PIA24 correlated with two structures that had hydrocarbon tails, these did not contain a phosphate group. Overall, probing the database with PIA24 yielded the most diverse structures. To confirm the validity of this approach, we compared the activity profiles of LY294002, wortmannin, and rapamycin and obtained values of (0.31 < PCC < 0.49) for all combinations, indicating that this type of analysis is capable of associating inhibitors that target the same pathway. These studies show that PIAs have the greatest similarity to other lipid-based molecules and validate the biological distinction between PIAs and other inhibitors of the pathway, such as LY294002, wortmannin, or rapamycin.
Compounds with similar sensitivity profiles to the PIAs in the NCI60 cell line panel. The COMPARE program was used to compare the sensitivity of PIAs to other compounds that have been screened (and for which data is publicly available) in the National Cancer Institute database. Correlations between PIAs are excluded from these tables (see text). Structures of the top five compounds that had the highest Pearson correlation coefficient scores with (A) PIA5, (B) PIA6, (C) PIA23, (D) PIA24, and (E) PIA25 used as a seed.
Compounds with similar sensitivity profiles to the PIAs in the NCI60 cell line panel. The COMPARE program was used to compare the sensitivity of PIAs to other compounds that have been screened (and for which data is publicly available) in the National Cancer Institute database. Correlations between PIAs are excluded from these tables (see text). Structures of the top five compounds that had the highest Pearson correlation coefficient scores with (A) PIA5, (B) PIA6, (C) PIA23, (D) PIA24, and (E) PIA25 used as a seed.
Discussion
Our studies indicate that PIAs can kill a variety of cancer cells in a manner that is dependent upon Akt activity yet is distinct from other inhibitors of the pathway. In the NCI60 cell line panel, PIAs were able to inhibit the growth of all the cell lines within the dose range tested (up to 100 μmol/L), with the GI50 for the majority of the cell lines falling between 1 and 20 μmol/L. Although broadly effective, these compounds showed only moderate potency. A possible reason for this might be the affinity of the lipid side chain of PIAs for serum, a phenomenon that has been observed previously (2). Despite this moderate potency, activity in 8 of 12 cell lines in hollow fiber assays was observed, indicating that effective circulating levels of PIAs can be achieved.
How do PIAs compare to other inhibitors of the pathway, such as LY294002 and rapamycin? The PI3K inhibitors LY294002 and wortmannin caused growth inhibition in a similar dose range as the PIAs; however, PIAs caused more cytotoxicity than LY294002 or wortmannin at the higher doses. Although the sensitivity patterns of rapamycin, LY294002, and wortmannin correlated with each other, a COMPARE analysis of the PIAs with these other inhibitors did not yield any significant correlation. This indicates that the pattern of growth inhibition in the different cell types to PIAs versus these other agents is dissimilar. Therefore, one might predict that clinical activity and/or toxicities of these compounds might also be dissimilar. Clinical trials that assess tolerability and efficacy will determine which inhibitors are the most promising drugs.
Establishing the relationship between target modulation by an agent and anticancer effects is an important step in targeted drug development. Sensitivity to PIAs weakly correlated with levels of phospho-Akt but not total Akt in the NCI60 cell line panel. However, we also found a number of other targets whose expression was more highly correlated to PIA sensitivity than active Akt. Several of these were common across three or more PIAs. It is unknown whether PIAs modulate any of these targets directly, but these targets that positively or negatively correlate with sensitivity to PIAs could potentially serve as predictive biomarkers for responsiveness to PIAs.
The COMPARE algorithm was also used to identify compounds that have been screened in the NCI60 that have similar sensitivity profiles to the PIAs. This yielded phospholipids that differed in the head group and linkage to the lipid side chain, as well as length of the side chain. One group of lipid-based agents that was identified (alkylphospholipids) has diverse functions and is currently under development as anticancer therapies. Interestingly, the alkylphospholipids that correlated with PIA5 and PIA6 (miltefosine and perifosine) also inhibit the translocation and activation of Akt (11, 12).
These studies show a wide spectrum of activity for PIAs, provide evidence of in vivo activity, validate active Akt as a predictive factor for PIAs, and identify new molecular targets that are associated with response to PIAs. Whether levels of expression of these novel targets contribute to the cytotoxic effects of PIAs, or whether they may serve as useful predictors of PIA sensitivity, will be the subjects of future research. Regardless, the wide spectrum of activity, the in vivo efficacy, and the novel combination of biological activities of PIAs that includes inhibition of Akt and activation of p38α highlight the potential of these compounds as cancer therapeutics.
Grant support: Intramural Research Program of the NIH/National Cancer Institute, Center for Cancer Research and federal funds from the National Cancer Institute/NIH under NCI contract NO1-CO-12400.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
Note: The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does mention of trade names, commercial products, or organizations imply endorsement by the U.S. Government.
Acknowledgments
We thank the members of the Dennis lab for helpful discussions.