- Open Access
Gene expression-based screening for inhibitors of PDGFR signaling
Genome Biology volume 9, Article number: R47 (2008)
Here we describe a proof-of-concept experiment designed to explore the possibility of using gene expression-based high-throughput screening (GE-HTS) to find inhibitors of a signaling cascade, using platelet derived growth factor receptor (PDGFR) signaling as the example. The previously unrecognized ability of aurintricarboxylic acid to inhibit PDGFR signaling, discovered through a screen of 1,739 compounds, demonstrates the feasibility and generalizability of GE-HTS for the discovery of small molecule modulators of any signaling pathway of interest.
High throughput screening of small-molecule libraries is a well-established and highly productive tool for the identification of chemical compounds targeting a specific protein function of interest. Traditionally, the high-throughput screening for modulators of molecular pathways involves cell-free biochemical assays, or in some cases, highly specialized cell-based phenotypic assays . However, in many cases the optimal target for therapeutic intervention is not known, or the development of a suitable phenotypic read-out is not technically feasible. For example, it is becoming increasingly of interest to modulate the activity of particular signal transduction pathways, but the components of such pathways are in many cases only partially known. It would therefore be of interest to develop a screening approach that could identify inhibitors of such pathways without first defining the biochemical target of candidate small molecules. Here we demonstrate that it is possible to use mRNA expression levels as a read-out to infer activity of a signal transduction pathway, thus establishing a general approach to screening for modulators of signal transduction pathways.
Endogenous mRNA expression has been previously successfully used as a surrogate of cellular states in high-throughput screening for compounds inducing differentiation of acute myeloid leukemia cells, and for identifying inhibitors of androgen receptor-mediated transcriptional activation in prostate cancer [2–5]. It is not obvious, however, that gene expression signatures could be used to identify inhibitors of signal transduction pathways that are regulated at the level of post-translational modification (phosphorylation), as opposed to transcriptional regulation.
To test the feasibility of using gene expression-based high-throughput screening (GE-HTS) to identify inhibitors of a signaling pathway, we chose platelet derived growth factor receptor (PDGFR) signaling for a proof-of-concept study, with particular emphasis on downstream activation of the extracellular regulated kinase (ERK) pathway (also known as the p42/p44 mitogen activated protein (MAP) kinase pathway) as a target pathway for the screen. The ERK pathway plays a major role in the control of cell growth, cell differentiation and cell survival . The protein kinase cascade Raf>mitogen/extracellular signal-regulated kinase (MEK)>ERK, also referred to as the MAP kinase module, is activated in mammalian cells through receptor tyrosine kinases, G-protein coupled receptors and integrins . Activated ERKs translocate to the nucleus where they phosphorylate transcription factors. The ERK pathway is often upregulated in human tumors , and as such is an attractive target for anticancer therapy. Furthermore, because the pathway has been extensively studied, many experimental tools are available with which to interrogate the pathway. We demonstrate here that indeed small molecule inhibitors of the PDGFR/ERK pathway can be discovered using the GE-HTS approach.
Identification of a signature of PDGFR/ERK activity
In GE-HTS, a gene expression signature is used as a surrogate of a biological state. In the present context, we sought to define a signature of ERK activation mediated by PDGFR stimulation. Specifically, we treated SH-SY5Y neuroblastoma cells with the BB homodimer of PDGF (PDGF-BB), which resulted in PDGFRβ phosphorylation and subsequent ERK activation. We selected PDGFRβ over PDGFRα for our studies because of previous observations that PDGFRα might mediate functions of other PDGF isoforms in addition to PDGF-A [7, 8]. The activation state of the members of the PDGFβ pathway can be traced by increase in their phosphorylation levels shortly after introduction of the growth factor . In particular, ERK phosphorylation peaks at about 15-20 minutes after induction, and then decreases to background levels some 20-30 minutes later . Accordingly, we performed gene expression profiling using Affymetrix U133A arrays 30 minutes following PDGF stimulation, thereby identifying those genes whose expression is correlated with PDGFR activity. In order to identify the component of the gene expression signature that was attributable to ERK activation by PDGFR (as opposed to other pathways downstream of PDGFR), we also pretreated the cells with the MEK inhibitor U0126 and the ERK inhibitor apigenin, and repeated the gene expression profiling studies (Figure 1a).
To define the signature of ERK activation, we developed and applied a rank-pairwise comparison algorithm as described in Materials and methods. We note that the genes identified in this manner are chosen because of their ability to reflect the PDGF-stimulated state - not because of their necessarily being critical effectors of PDGFR signaling. The top three genes identified in this fashion were those for c-fos, early growth response 1 (EGR1), and activity-regulated cytoskeleton-associated protein (ARC). All three genes were previously shown to be upregulated by activation of ERK, and we further confirmed their regulation by reverse transcriptase (RT)-PCR (Figure 1b) [11–13]. Two additional genes, ribosomal protein RPL23A and ATP5B, were selected as internal controls, because their expression was not significantly altered by PDGFR activation.
High-throughput screening to find inhibitors of the PDGFR/ERK pathway
Having defined a gene expression signature of PDGFR/ERK activation, we next sought to screen a library of small molecules to find compounds that would reverse the signature (for primary screen data, see Additional data file 1). We chose TIP5 fibroblast cells for the high-throughput screen instead of SH-SY5Y neuroblastoma cells used to define the gene expression signature. Both TIP5 and SH-SY5Y cells have wild-type PDGFR/ERK signaling, which makes it unnecessary to employ mutant and/or constitutively activated PDGFR cascades. TIP5 cells, however, were more adherent to 384-well plates, making them more amenable to the screening setting.
The screen was performed as follows. TIP5 cells were plated in 384-well plates, serum-starved overnight and compounds then added by pin transfer. The compound library, previously described in , consisted of 1,739 chemicals with previously established biological functions. Some of the compounds have been approved for use in humans by the Food and Drug Administration. After a 30 minute compound-incubation period, PDGF-BB was added. 45 minutes later, the growth medium was discarded, and cells were lysed. RNA was then extracted, the signature genes amplified by RT-PCR, and the PCR amplicons quantified by single-base extension mass spectrometry, as we previously described  (Figure 1c). Cells were treated in triplicate at two concentrations (approximately 10 μM and 50 μM). Compounds were defined as hits if the expression of two marker genes, c-fos and EGR1, normalized by expression of control genes was significantly (more than one standard deviation) lower than average expression in all positive control wells. Compounds that inhibited the signature of the activated PDGFR/ERK pathway in four out of six replicas were selected as hits for further characterization.
Validation of hit compounds
Three wells met the hit selection criteria: aurintricarboxylic acid (ATA; free acid), aurintricarboxylic acid triammonium salt (aluminon), and quinacrine dihydrochloride (mepacrine) (Figure 2a,b); all three were therefore selected for further studies. Western analysis of total lysates from cells treated with these compounds demonstrated that both ATA and its salt (which in solution is identical to ATA), but not quinacrine dihydrochloride, abrogated PDGF-mediated phosphorylation of ERK (Figure 3a), thereby identifying ATA as an inhibitor of the ERK pathway. Quinacrine dihydrochloride did not inhibit ERK phosphorylation, but it has been previously shown to be a non-specific inhibitor of phospholipase A2 . Activated ERK phosphorylates phospholipase A2 , and as a result stimulates transcription of the c-fos and EGR1 genes, two components of our ERK signature .
We then relaxed hit selection criteria, and identified nine more potential candidates. However, further study indicated that none of these nine additional compounds affected activation of the PDGFR/ERK pathway.
Disruption of phosphorylation of ERK by ATA was an indication that ATA inhibited the PDGFR/ERK pathway upstream of ERK. Subsequent analysis indicated that phosphorylation of both MEK (Figure 3b) and PDGFR (Figure 3c) was abrogated by ATA, thus pointing to PDGFR as a possible ATA target.
To address the possibility that ATA might in some fashion deplete PDGF ligand from the growth medium, TIP5 cells were first incubated with ATA for 30 minutes. Next, the cells were washed thrice with serum-free medium and then stimulated with PDGF. As shown in Figure 3d, PDGFR phosphorylation remained inhibited, suggesting that PDGF ligand was unlikely to be the target of ATA.
The experiments described so far indicated that ATA inhibits PDGF-mediated ERK phosphorylation by inhibiting PDGFR phosphorylation. To localize the portion of PDGFR targeted by ATA, we utilized a series of chimeric receptor constructs (Figure 4a). The first chimera, TEL/PDGFR, is a naturally occurring, leukemia-associated fusion of the oligomerization domain of the transcription factor TEL (ETV6) to the transmembrane and cytoplasmic domains of PDGFR, resulting in constitutive activation of PDGFR . As shown in Figure 4b, ATA was unable to inhibit TEL/PDGFR phosphorylation at concentrations as high as 100 μM, indicating that ATA does not target the transmembrane or cytoplasmic portions of PDGFR present in the TEL/PDGFR chimera.
The next chimera, termed PER, is composed of the extracellular domain of PDGFR and the transmembrane and cytoplasmic domains of epidermal growth factor receptor (EGFR) . ATA inhibited PER phosphorylation in PER-PC12 cells (Figure 4c), thus mapping the site of ATA action to the extracellular domain of PDGFR. To exclude the possibility of ATA inhibiting any receptor tyrosine kinase extracellular domain, we tested ATA against a third chimera, EKR, consisting of the extracellular domain of EGFR and the transmembrane and cytoplasmic domains of c-KIT . ATA failed to inhibit EKR (Figure 5a), indicating that ATA exhibits some specificity for the PDGFR extracellular domain. Similarly, ATA failed to inhibit insulin-like growth factor (IGF)-induced phosphorylation of IGF1 receptor (IGF1R; Figure 5b), or EGF-induced phosphorylation of EGFR (Figure 5c) . Interestingly, ATA did inhibit stem cell factor (SCF)-mediated activation of cKIT (Figure 5d). The cKIT and PDGFR extracellular domains have 41% sequence similarity (26% identity), whereas no significant homology is seen between the extracellular domains of PDGFR and EGFR or IGF1R.
We note that whereas phosphorylation of the PER chimera is PDGF-dependent (and ATA inhibitable) in PER-PC12 cells, PER is constitutively active in 501 MEL and MCF7 cells, and in those contexts PER phosphorylation is not fully abrogated by ATA (Figure 6a,b). These experiments further point to the possibility of ATA inhibiting PDGF binding to the extracellular domain of PDGFR and disrupting ligand-mediated activation of the receptor.
Structure-activity relationships in the series of ATA analogues
In order to characterize the features of the ATA molecule required for biological activity, we analyzed a diverse set of ATA structural analogs (Figure S1 in Additional data file 2) available from the Available Chemicals Directory . We split compounds into three groups to test three different hypotheses on the structure-activity relationship in the series. The activities of methylenedisalicylic acid, salicylic acid and 3-methylsalicylic acid (Figure S1a in Additional data file 2) were analyzed to examine if the skeletal-triphenylmethane structure of ATA was essential to its activity. Aurin, uranine and phenolphthalein sodium salt (Figure S1b in Additional data file 2) were tested to evaluate the roles the carboxyl and hydroxyl groups on the triphenylmethane scaffold play in the inhibitory potency of ATA. Compounds in the third group (Figure S1c in Additional data file 2) were evaluated to test the effect of various modifications of the phenyl rings on the inhibitory properties of ATA. No compounds in the series inhibited PDGFR at concentrations sufficient for ATA inhibition (less than 5 μM). In the first group, methylenedisalicylic acid (Figure 7a), but not methylsalicylic or salicylic acids inhibited PDGFR phosphorylation at 50 μM, suggesting that increasing the number of substituted salicylic acid moieties from one to three boosts the inhibitory potency of ATA. The positions and number of carboxyl and hydroxyl groups were essential for PDGFR inhibition, as indicated by the fact that no compounds in the second group inhibited PDGFR at 100 μM concentration. These results corroborate earlier reports that both the aurin triphenyl methane ring system and the carboxylic acid groups are necessary for ATA inhibitory properties .
In the third group, Basic Violet 3, Ethyl Violet and Victoria Pure Blue BO inhibited PDGFR in the 5-10 μM range (Figure 7b-d). Interestingly, these three compounds exhibited less specific patterns of receptor inhibition than ATA, inhibiting not only cKIT, but also EGFR and IGF1R at 10-100 μM (Figure 8). Moreover, different from ATA, Ethyl Violet and Victoria Pure Blue BO readily translocated across the cell membrane, as indicated by their inhibition of cytoplasmic TEL/PDGFR in Ba/F3 cells at 10 μM (Figure 9). Taken together, these results suggest that the inhibitory mechanism of Basic Violet 3, Ethyl Violet and Victoria Pure Blue BO is different from the extracellular receptor inhibition mechanism of ATA.
In this report, we describe the proof-of-concept efforts to approach the discovery of inhibitors of signal transduction using a novel chemical genomic approach. We discovered a previously unknown property of the triphenylmethane derivative ATA, using GE-HTS. Having defined a signature of PDGFR activation, we screened a library of bioactive small molecules for compounds capable of turning off the signature. Importantly, the screen required neither a highly specialized signal transduction assay, nor prior knowledge of the protein to be targeted. In principle, small molecules acting upstream, downstream or at the level of PDGFR itself would be captured by the screen.
Two compounds in the library met pre-established criteria for hits abrogating the PDGFR/ERK activation signature. The hit compounds reproducibly inhibited the signature in follow-up studies, indicating that the false positive rate of the screen was quite low. One of the hits, quinacrine dihydrochloride, is a known inhibitor of phospholipase A2, a known regulator of ERK signaling [14–16]. The other compound, ATA, was a novel discovery, and was therefore followed up in greater detail.
ATA is a polymeric carboxylated triphenylmethane derivate with a molecular weight range of 422-6,500 , that has displayed a wide range of biological activity in in vitro biochemical assays. For example, ATA has been reported to inhibit enzymes involved in protein-nucleic acid interactions, including DNA and RNA polymerases, reverse transcriptase, nucleases, primases, topoisomerases, ribonucleotide reductases, aminoacyl-tRNA synthetase, nuclear factor-kappaB and HIV-1 integration protein . In addition, ATA has also been shown to inhibit other classes of proteins in vitro, including phosphatases , NAD(H)/NADP(H)-requiring enzymes , aminopropyltransferases , mu- and m-calpain , delta aminolevulinic acid dehydratase , glucose-6-phosphate dehydrogenase , phenylalanine:tRNA ligase  and kinases, such as phosphofructokinase , ERK, p38 MAP kinase, IkappaB kinase , inositol-1,4,5-trisphosphate 3-kinase and inositol polypohosphate multikinase . In vitro inhibition of protein synthesis has also been described .
Biological activity of ATA has also been observed in vivo, although in most cases only at rather high concentrations. For example, ATA is reported to obviate binding of interferon-alpha to its receptor in the 12-50 μM range , to prevent von Willebrand factor binding to platelet receptor glycoprotein Ib , and to block binding of gp120, the HIV coat protein, to its receptor, CD4 . Similarly, ATA has been shown to be a N-methyl-D-aspartate (NMDA) receptor antagonist with an IC50 of 26.9 μM and was reported to antagonize excitotoxicity at both NMDA and non-NMDA glutamate receptors in the 50-100 μM range . ATA inhibited progesterone receptor at 100-500 μM , estradiol receptor at 100-200 μM , and glucocorticoid receptor complex at 50-200 μM . ATA also was reported to activate IGF1R (25-100 μM)  and erbB4 (10 μM) . These studies suggest that ATA has a range of biological activities, most of which, however, are observable only at quite high concentrations, in many cases as high as 100 μM.
More limited activity has been reported at lower concentrations of ATA. For example, at 1-5 μM, ATA was reported to reverse the transformed phenotype of cells transfected with basic fibroblast growth factor fused to a signal peptide sequence (spbFGF cells) . It was suggested, on the basis of ATA fluorescence studies, that ATA binds to acidic fibroblast growth factor, altering its physicochemical properties and decreasing its mitogenic activity , although these results were not confirmed by more direct biochemical methods. The observed ATA interactions in this setting take place at the cell surface, consistent with the finding that ATA does not readily penetrate cellular membranes. ATA is not taken up by HeLa cells, VERO cells, rabbit reticulocytes, or a variety of bacterial cells . Accordingly, ATA did not inhibit intracellular proteins, even at concentrations hundreds of times higher than those required for inhibition in vitro . Only at high concentrations (500 μM) was intracellular ATA fluorescence detectable . It seems most likely, therefore, that our observed effects of ATA on PDGFR activity occur at the cell surface.
Consistent with this notion, our analysis indicated that all signaling downstream of PDGFR was inhibited by ATA, and ATA wash-out experiments suggested that ATA did not abrogate the signaling by binding and inactivating PDGF. Furthermore, analysis of chimeric PDGFR constructs localized the ATA effect to the PDGFR extracellular domain. Interestingly, modest concentrations of ATA (2-5 μM) also inhibited activity of the related receptor tyrosine kinase cKIT, which shares sequence homology with PDGFR in the extracellular domain, whereas kinases lacking such homology (for example, IGFR and EGFR) were inhibited only at concentrations of 100 μM. It is possible that the previously described inhibition of JAK/STAT signaling by ATA [32, 44] is attributable to its inhibition of PDGFR family receptor tyrosine kinases, known to be upstream activators of the JAK/STAT pathway [45, 46].
The polymeric nature of ATA may make it unattractive as a therapeutic agent and, moreover, multiple highly potent PDGFR kinase inhibitors have been previously reported . Our work establishes proof of concept, however, for the notion that mRNA expression signatures can be effectively used as a read-out for the identification of inhibitors of signal transduction, often thought approachable only through the direct examination of protein phosphorylation states. We note that indeed antibody-based high-throughput screens have been reported , but such assays obviously require the availability of a sufficiently sensitive and specific antibody for this purpose. For many, if not most, proteins of interest, such high quality antibodies are not available. The ability to convert any biological process or cell state into a completely generic gene expression signature that can be monitored in high throughput and at low cost is therefore attractive.
The implementation of the GE-HTS concept described here involves the detection of multiplexed RT-PCR signature genes by a single-base-extension reaction followed by MALDI-TOF (matrix assisted laser desorption ionization-time of flight) mass spectrometry . While this method was effective in the study described here, it has several limitations. For example, conventional RT-PCR amplification is not easily multiplexed, and the ability to simultaneously detect multiple amplicons by the mass spectrometric method is limited. Lastly, the approach can become expensive if extended to the ultra-high-throughput setting. We have therefore modified the approach to allow for the efficient amplification of up to 100 transcripts using a ligation-mediated amplification method, followed by detection on polystyrene beads via flow cytometry, as we recently described [3, 5, 49]. The present study, however, establishes that the GE-HTS concept can be applied to screening for modulators of signal transduction, representing a general approach to the discovery of compounds that affect any signaling pathway of interest.
Materials and methods
The EKR construct  was kindly provided by Dr. Ullrich, Department of Molecular Biology, Max-Planck-Institut fur Biochemie. PER chimera  was a gift of Dr. Tyson and Dr. Bradshaw, Department of Physiology and Biophysics, University of California, Irvine.
Chemical compounds apigenin, U0126, quinacrine dihydrochloride and ATA were obtained from Calbiochem ; Methyl Violet B base, Rhodamine 6 G tetrafluoroborate, sulforhodamine, Ethyl Violet, Victoria Pure Blue BO, Rhodamine B, Lissamine Green B, Methyl Violet 2B, Rhodamine 6G, (L-Asp)2Rhodamine 110 TFA, Rhodamine 110 chloride, Eosin B, Rhodamine 123 hydrate, Rhodamine 19 perchlorate, Acid Fuchsin calcium salt, p-Rosolic acid, Basic Violet2, Gentian Violet, pararosaniline hydrochloride, and salicylic acid were purchased from Sigma ; and 3-methylsalicylic acid, 5,5'-methylenedisalicylic acid, phenolphthalein sodium salt, and Uranine K were obtained from ABCR .
Cell culture reagents RPMI 1640, Dulbecco's modified Eagle's medium (DMEM), and HAM's F-10 were purchased from Mediatech , penicillin and streptomycin from Invitrogen , and fetal bovine serum from Sigma. p44/42 MAP kinase, phospho-p44/42 MAP kinase (Thr202/Tyr204), MEK1/2, phospho-MEK1/2 (Ser217/221), PDGF-BB, phospho-PDGFRβ (Tyr751), phospho-EGFR (Tyr1068), cKIT, Phospho-cKIT (Tyr719), IGF-IαR, and Phospho-IGF-IR (Tyr1131)/insulin receptor (Tyr1146) antibodies were obtained from Cell Signaling. EGFR and mouse cKIT antibodies were purchased from Santa Cruz Biotechnology . Alfa-tubulin antibody was obtained from Sigma.
SH-SY5Y neuroblastoma cells were purchased from American Type Culture Collection . The IL3-dependent pro-B lymphoid cell line Ba/F3 and Ba/F3 cells expressing TEL/PDGFRβ [17, 59] were obtained from Dr. Gary Gilliland. TIP5 primary fibroblasts  were a gift from Dr. Stephen Lessnick. We thank Dr. Ruth Halaban for 501 MEL human melanoma cells. PER-expressing PC12 cells were generously provided by Dr. Darren Tyson. SH-SY5Y, PC12, TIP5 and MCF7 cells were cultured in DMEM, BaF3 cells and BaF3 cells expressing TEL/PDGFRβ were maintained in RPMI 1640 medium, and 501 MEL cells were grown in Ham's 10 medium. Medium for IL3-dependent Baf3 cells was supplemented with 0.05 ng/ml IL3. Media for all cell lines except PC12 contained 10% fetal bovine serum, 10 U/ml penicillin, and 10 μg/ml streptomycin. PC12 cells were grown in DMEM with 15% horse serum, 5% fetal bovine serum, 10 U/ml penicillin, and 10 μg/ml streptomycin. All cells were grown at 37°C in 5% CO2.
Characterization of the activation signature for ERK/PDGFR pathway
SH-SY5Y cells were grown to confluence and starved overnight in serum-free medium in order to silence any sustained effects from growth factor signaling. Prior to induction with 50 ng/ml PDGF, cells were treated with pathway inhibitors 74 μM apigenin or 50 μM U0126, or with dimethyl sulfoxide (DMSO) as carrier for 60 minutes. Total RNA was isolated 30 minutes after PDGF addition. Experiments were performed in duplicate. The RNA was processed and hybridized with Affymetrix U133A GeneChips as described in .
To define the ERK/PDGFR activation signature, a pair-ranking algorithm was used. Genes were ranked according to their raw expression values on each chip. Ten genes with maximum change in ranking were selected for each one of three pairs of conditions: cells with PDGF versus cells without PDGF, cells with PDGF versus cells with PDGF and apigenin, and cells with PDGF versus cells with PDGF and U0126. Three genes common to all three conditions were selected as a signature of the activated ERK/PDGFR pathway. The signature was then trimmed from three to two genes based on their relative strength of expression in TIP5 cells.
TIP5 cells were grown to confluence and starved overnight with 20 μl serum-free medium per well of 384-well plates. We added 20 μl of compounds diluted in media so that the final concentration of compounds would be approximately 10 μM in three out of six replicas, and 50 μM in three remaining replicas. Media containing carrier (DMSO) was added to control wells instead of compounds. The compound library was composed of 1,739 chemicals either approved for use in humans by the Food and Drug Administration or extensively biologically characterized [2, 62] (the full list of tested compounds is available in these publications). After 30 minutes of compound treatment, cells were induced with 40 μl of PDGF diluted in media (final PDGF concentration 40 ng/ml). PDGF was added to half of control wells to measure PDGF response; only media was added to the remaining control wells. After 40 minutes of PDGF induction, media was discarded, cells were lysed and RNA was extracted and quantified as described in .
Briefly, 15 μl of lysis solution containing a hypotonic detergent, dithiothreitol and RNAse inhibitor were added to medium-free cells for 15 minutes. The lysates were transferred to a 384-well oligo-dT-coated plate and incubated with 6 μl of 2.5× binding buffer. After 30 minutes of incubation lysates were discarded and reverse transcription was carried out in a 5 μl Moloney murine leukemia virus (M-MuLV) RT reaction at 37°C for 2 h.
After incubation, the RT mixture was discarded and multiplex PCR was carried in a 5 μl volume. The resulting mixture was treated with shrimp alkaline phosphatase and the single-base extension reactions were carried out in 9 μl reaction volumes with 1× Thermosequenase buffer, 2.7 μM of each primer, 0.2 mM of each ddNTP and 0.58 units per reaction of Thermosequenase as described in .
The lysis buffers, 384-well custom-coated oligo-dT plates, and M-MuLV were purchased from Pierce  and were used according to a modified version of the Express Direct mRNA Capture and RT-PCR system. The shrimp alkaline phosphatase, Thermosequenase buffer, ddNTP and Thermosequenase were obtained from Sequenom .
The primers used for multiplex PCR reactions were: EGR1, 5'-AGC GGA TAA CAC CTC ATA CCC ATC CCC TGT-3' and 5'-AGC GGA TAA CTG TCC TGG GAG AAA AGG TTG-3'; c-fos, 5'-AGC GGA TAA CGC TTC CCT TGA TCT GAC TGG-3' and 5'-AGC GGA TAA CAT GAT GCT GGG AAC AGG AAG-3'; ATP5B, 5'-AGC GGA TAA CCA AAG CCC ATG GTG GTT ACT-3' and 5'-AGC GGA TAA CGC CCA ATA ATG CAG ACA CCT-3'; RPL23A, 5'-AGC GGA TAA CAA GAA GAA GAT CCG CAC GTC-3' and 5'-AGC GGA TAA CCG AAT CAG GGT GTT GAC CTT-3'.
The following primers were used for single-base extension reactions: EGR1, 5'-TTC CCC CTG CTT TCC CG-3'; c-fos, 5'-TGC CTC TCC TCA ATG ACC CT-3'; ATP5B, 5'-GAC TGT GGC TGA ATA CTT CA-3'; RPL23A, 5'-GTC TGC CAT GAA GAA GAT AGA A-3'.
To select compounds that inhibited expression of the pathway signature, the following procedure was performed. For each compound on each plate four ratios were determined: expression of EGR1 versus expression of ATP5B(VEGR1/ATP5B), EGR1 versus RPL23A (VEGR1/RPL23A), c-fos versus ATP5B (Vc-fos/ATP5B), and c-fos versus RPL23A (Vc-fos/RPL23A). For each plate a median (μ) and standard deviation (σ) were determined for each of four ratios. A compound was considered a plate hit if each of the four ratios for this compound were at least one standard deviation smaller than the median (V < (μ - σ): VEGR1/ATP5B < μEGR1/ATP5B - σEGR1/ATP5B; VEGR1/RPL23A < μEGR>1/RPL23A - σEGR1/RPL23A; Vc-fos/ATP5B < μc-fos/ATP5B - σc-fos/ATP5B; and Vc-fos/RPL23A < μc-fos/RPL23A - σc-fos/RPL23A). Compounds that were plate hits in four out of six replicas were selected for further consideration.
Western blotting and transfection
For transfection experiments, 501 MEL cells or TIP5 cells were grown overnight to 50% confluence and transfected using Fugene 6 transfection reagent (Roche ) as recommended by the manufacturer. Then 24 h after transfection, medium was exchanged for a serum free one, and cells were serum starved overnight.
Otherwise, adherent cells (TIP5, MEL 501) were grown to confluence, serum starved overnight, and treated with compounds and growth factors as indicated. Cells growing in suspension (BaF3 cells and BaF3 cells expressing TEL/PDGFR protein) were grown to 106 cells/ml and treated with compounds as indicated. After treatment media was removed, adherent cells were scraped with Sample Buffer from Cell Signaling, and suspension cells were pelleted and resuspended in Sample Buffer. The resulting lysates of approximately 1 × 105 cells were boiled, chilled, run on 4-15% gradient gels from BioRad , transferred to a polyvinylidene difluoride membrane from Millipore , blocked, probed and visualized as recommended by the antibody manufacturers.
Comparative sequence analysis between PDGF (UniProtKB/Swiss-Prot entry P09619), cKIT (UniProtKB/Swiss-Prot entry P10721), EGFR (UniProtKB/Swiss-Prot entry P00533), and IGFR (UniProtKB/Swiss-Prot entry P08069) was performed with BLAST 2 SEQUENCES .
Averaging and normalization of high-throughput screen expression levels of marker genes c-fos and EGR1 for Figure 2b
Each primary screen replica plate contained 16 wells with PDGFβ and carrier DMSO as a positive control for PDGF activation (called PDGF in Figure 2b), and 16 wells with carrier only as a negative control for PDGF activation (called No PDGF in Figure 2b). The expression levels of marker genes normalized by expression of control gene (ratios c-fos/ATP5B and EGR1/ATP5B) were averaged for 16 PDGF wells to have a single value for the positive PDGF control per plate, and for 16 No PDGF wells to have a single value for the negative No PDGF control per replica plate. Only one well was allocated for each hit compound on a single replica plate.
To compare data between replica plates in Figure 2b, the ratios c-fos/ATP5B and EGR1/ATP5B were adjusted to be equal to 1 for positive PDGF control. This means that on each replica plate the marker/control ratios in all wells were divided by the corresponding value for the positive PDGF control for this plate. The procedure was performed independently for both c-fos/ATP5B and EGR1/ATP5B ratios. As a result of this procedure, the c-fos/ATP5B and EGR1/ATP5B ratios for the hit compounds and for the No PDGF control on each plate were divided by PDGF control c-fos/ATP5B and EGR1/ATP5B ratios for this plate. The resulting adjusted values were then averaged between three replica plates.
The data have been deposited in the Gene Expression Omnibus  with accession number GSE7403.
Dulbecco's modified Eagle's medium
epidermal growth factor
extracellular regulated kinase
gene expression-based high-throughput screening
insulin-like growth factor
mitogen activated protein
mitogen/extracellular signal-regulated kinase
Moloney murine leukemia virus
platelet derived growth factor
receptor tyrosine kinase
stem cell factor.
Johnston P: Cellular assays in HTS. Methods Mol Biol. 2002, 190: 107-116.
Stegmaier K, Ross KN, Colavito SA, O'Malley S, Stockwell BR, Golub TR: Gene expression-based high-throughput screening(GE-HTS) and application to leukemia differentiation. Nat Genet. 2004, 36: 257-263. 10.1038/ng1305.
Hieronymus H, Lamb J, Ross KN, Peng XP, Clement C, Rodina A, Nieto M, Du J, Stegmaier K, Raj SM, Maloney KN, Clardy J, Hahn WC, Chiosis G, Golub TR: Gene expression signature-based chemical genomic prediction identifies a novel class of HSP90 pathway modulators. Cancer Cell. 2006, 10: 321-330. 10.1016/j.ccr.2006.09.005.
Wei G, Twomey D, Lamb J, Schlis K, Agarwal J, Stam RW, Opferman JT, Sallan SE, den Boer ML, Pieters R, Golub TR, Armstrong SA: Gene expression-based chemical genomics identifies rapamycin as a modulator of MCL1 and glucocorticoid resistance. Cancer Cell. 2006, 10: 331-342. 10.1016/j.ccr.2006.09.006.
Stegmaier K, Wong JS, Ross KN, Chow KT, Peck D, Wright RD, Lessnick SL, Kung AL, Golub TR: Signature-based small molecule screening identifies cytosine arabinoside as an EWS/FLI modulator in Ewing sarcoma. PLoS Med. 2007, 4: e122-10.1371/journal.pmed.0040122.
Steelman LS, Pohnert SC, Shelton JG, Franklin RA, Bertrand FE, McCubrey JA: JAK/STAT, Raf/MEK/ERK, PI3K/Akt and BCR-ABL in cell cycle progression and leukemogenesis. Leukemia. 2004, 18: 189-218. 10.1038/sj.leu.2403241.
Soriano P: The PDGF alpha receptor is required for neural crest cell development and for normal patterning of the somites. Development. 1997, 124: 2691-2700.
Boström H, Willetts K, Pekny M, Levéen P, Lindahl P, Hedstrand H, Pekna M, Hellström M, Gebre-Medhin S, Schalling M, Nilsson M, Kurland S, Törnell J, Heath JK, Betsholtz C: PDGF-A signaling is a critical event in lung alveolar myofibroblast development and alveogenesis. Cell. 1996, 85: 863-873. 10.1016/S0092-8674(00)81270-2.
Heldin CH, Westermark B: Mechanism of action and in vivo role of platelet-derived growth factor. Physiol Rev. 1999, 79: 1283-1316.
Kolch W: Meaningful relationships: the regulation of the Ras/Raf/MEK/ERK pathway by protein interactions. Biochem J. 2000, 351: 289-305. 10.1042/0264-6021:3510289.
Westwick JK, Cox AD, Der CJ, Cobb MH, Hibi M, Karin M, Brenner DA: Oncogenic Ras activates c-Jun via a separate pathway from the activation of extracellular signal-regulated kinases. Proc Natl Acad Sci USA. 1994, 91: 6030-6034. 10.1073/pnas.91.13.6030.
Harada T, Morooka T, Ogawa S, Nishida E: ERK induces p35, a neuron-specific activator of Cdk5, through induction of Egr1. Nat Cell Biol. 2001, 3: 453-459. 10.1038/35074516.
Ying SW, Futter M, Rosenblum K, Webber MJ, Hunt SP, Bliss TV, Bramham CR: Brain-derived neurotrophic factor induces long-term potentiation in intact adult hippocampus: requirement for ERK activation coupled to CREB and upregulation of Arc synthesis. J Neurosci. 2002, 22: 1532-1540.
Martin TW, Wysolmerski RB, Lagunoff D: Phosphatidylcholine metabolism in endothelial cells: evidence for phospholipase A and a novel Ca2+-independent phospholipase C. Biochim Biophys Acta. 1987, 917: 296-307.
Croxtall JD, Choudhury Q, Flower RJ: Glucocorticoids act within minutes to inhibit recruitment of signalling factors to activated EGF receptors through a receptor-dependent, transcription-independent mechanism. Br J Pharmacol. 2000, 130: 289-298. 10.1038/sj.bjp.0703272.
Lerea LS, Carlson NG, McNamara JO: N-methyl-D-aspartate receptors activate transcription of c-fos and NGFI-A by distinct phospholipase A2-requiring intracellular signaling pathways. Mol Pharmacol. 1995, 47: 1119-1125.
Golub TR, Barker GF, Lovett M, Gilliland DG: Fusion of PDGF receptor beta to a novel ets-like gene, tel, in chronic myelomonocytic leukemia with t(5;12) chromosomal translocation. Cell. 1994, 77: 307-316. 10.1016/0092-8674(94)90322-0.
Tyson DR, Larkin S, Hamai Y, Bradshaw RA: PC12 cell activation by epidermal growth factor receptor: role of autophosphorylation sites. Int J Dev Neurosci. 2003, 21: 63-74. 10.1016/S0736-5748(02)00139-9.
Kubota Y, Angelotti T, Niederfellner G, Herbst R, Ullrich A: Activation of phosphatidylinositol 3-kinase is necessary for differentiation of FDC-P1 cells following stimulation of type III receptor tyrosine kinases. Cell Growth Differ. 1998, 9: 247-256.
Blume-Jensen P, Hunter T: Oncogenic kinase signalling. Nature. 2001, 411: 355-365. 10.1038/35077225.
Available Chemicals Directory. [http://www.mdli.com/products/experiment/available_chem_dir/index.jsp]
Beery R, Haimsohn M, Wertheim N, Hemi R, Nir U, Karasik A, Kanety H, Geiger A: Activation of the insulin-like growth factor 1 signaling pathway by the antiapoptotic agents aurintricarboxylic acid and Evans blue. Endocrinology. 2001, 142: 3098-3107. 10.1210/en.142.7.3098.
Wang P, Kozlowski J, Cushman M: Isolation and structure elucidation of low molecular weight components of aurintricarboxylic acid (ATA). J Org Chem. 1992, 57: 3861-3866. 10.1021/jo00040a027.
Cho H, Lee DY, Shrestha S, Shim YS, Kim KC, Kim MK, Lee KH, Won J, Kang JS: Aurintricarboxylic acid translocates across the plasma membrane, inhibits protein tyrosine phosphatase and prevents apoptosis in PC12 cells. Mol Cells. 2004, 18: 46-52.
Thompson DC, Reed M: Inhibition of NAD(H)/NADP(H)-requiring enzymes by aurintricarboxylic acid. Toxicol Lett. 1995, 81: 141-149. 10.1016/0378-4274(95)03418-8.
Hibasami H, Tanaka M, Nagai J: Kinetic studies on inhibition of aminopropyltransferases by aurintricarboxylic acid in vitro. Chem Biol Interact. 1982, 40: 319-323. 10.1016/0009-2797(82)90154-5.
Posner A, Raser KJ, Hajimohammadreza I, Yuen PW, Wang KK: Aurintricarboxylic acid is an inhibitor of mu- and m-calpain. Biochem Mol Biol Int. 1995, 36: 291-299.
Allain P, Leblondel G, Streichenberger G: Anti-inflammatory effect of aluminum or aurin tricarboxylic acid. C R Seances Soc Biol Fil. 1980, 174: 68-73.
Hsu CL, Joshi JG: Aurinetricarboxylic acid: a potent inhibitor of glucose-6-phosphate dehydrogenase. Chem Biol Interact. 1977, 17: 61-67. 10.1016/0009-2797(77)90072-2.
Carias JR, Julien R: Competitive inhibition of Phe: tRNA ligase of sheep embryos by aurintricarboxylic acid. FEBS Lett. 1975, 56: 303-306. 10.1016/0014-5793(75)81114-8.
McCune SA, Foe LG, Kemp RG, Jurin RR: Aurintricarboxylic acid is a potent inhibitor of phosphofructokinase. Biochem J. 1989, 259: 925-927.
Tsi CJ, Chao Y, Chen CW, Lin WW: Aurintricarboxylic acid protects against cell death caused by lipopolysaccharide in macrophages by decreasing inducible nitric-oxide synthase induction via IkappaB kinase, extracellular signal-regulated kinase, and p38 mitogen-activated protein kinase inhibition. Mol Pharmacol. 2002, 62: 90-101. 10.1124/mol.62.1.90.
Mayr GW, Windhorst S, Hillemeier K: Antiproliferative plant and synthetic polyphenolics are specific inhibitors of vertebrate inositol-1,4,5-trisphosphate 3-kinases and inositol polyphosphate multikinase. J Biol Chem. 2005, 280: 13229-13240. 10.1074/jbc.M500545200.
Stewart ML, Grollman AP, Huang MT: Aurintricarboxylic acid: inhibitor of initiation of protein synthesis. Proc Natl Acad Sci USA. 1971, 68: 97-101. 10.1073/pnas.68.1.97.
Gan YX, Weaver JL, Pine PS, Zoon KC, Aszalos A: Aurin tricarboxylic acid, the anti-AIDS compound, prevents the binding of interferon-alpha to its receptor. Biochem Biophys Res Commun. 1990, 172: 1298-1303. 10.1016/0006-291X(90)91590-O.
Weinstein M, Vosburgh E, Philips M, Turner N, Chute-Rose L, Moake J: Isolation from commercial aurintricarboxylic acid of the most effective polymeric inhibitors of von Willebrand factor interaction with platelet glycoprotein Ib. Comparison with other polyanionic and polyaromatic polymers. Blood. 1991, 78: 2291-2298.
Zeevalk GD, Schoepp D, Nicklas WJ: Excitotoxicity at both NMDA and non-NMDA glutamate receptors is antagonized by aurintricarboxylic acid: evidence for differing mechanisms of action. J Neurochem. 1995, 64: 1749-1758.
Moudgil VK, Eessalu TE: Sensitivity of progesterone receptors to aurintricarboxylic acid. Inhibition of nuclear uptake, ATP and DNA binding. Biochim Biophys Acta. 1980, 627: 301-312.
Moudgil VK, Weekes GA: Inhibition of hen oviduct estradiol receptor by aurintricarboxylic acid. FEBS Lett. 1978, 94: 324-326. 10.1016/0014-5793(78)80967-3.
Okada N, Koizumi S: Tyrosine phosphorylation of ErbB4 is stimulated by aurintricarboxylic acid in human neuroblastoma SH-SY5Y cells. Biochem Biophys Res Commun. 1997, 230: 266-269. 10.1006/bbrc.1996.5934.
Benezra M, Vlodavsky I, Yayon A, Bar-Shavit R, Regan J, Chang M, Ben-Sasson S: Reversal of basic fibroblast growth factor-mediated autocrine cell transformation by aromatic anionic compounds. Cancer Res. 1992, 52: 5656-5662.
Lozano RM, Rivas G, Giménez-Gallego G: Destabilization, oligomerization and inhibition of the mitogenic activity of acidic fibroblast-growth factor by aurintricarboxylic acid. Eur J Biochem. 1997, 248: 30-36. 10.1111/j.1432-1033.1997.t01-2-00030.x.
Grollman AP: Inhibition of messenger ribonucleic acid attachment to ribosomes. II. Proposed mechanism for the design of novel antiviral agents. Antimicrobial Agents and Chemotherapy. Edited by: Hobby GL. 1968, New York: American Society of Microbiology, 36-40.
Chen CW, Chao Y, Chang YH, Hsu MJ, Lin WW: Inhibition of cytokine-induced JAK-STAT signalling pathways by an endonuclease inhibitor aurintricarboxylic acid. Br J Pharmacol. 2002, 137: 1011-1020. 10.1038/sj.bjp.0704955.
Novak U, Harpur AG, Paradiso L, Kanagasundaram V, Jaworowski A, Wilks AF, Hamilton JA: Colony-stimulating factor 1-induced STAT1 and STAT3 activation is accompanied by phosphorylation of Tyk2 in macrophages and Tyk2 and JAK1 in fibroblasts. Blood. 1995, 86: 2948-2956.
Mui AL, Wakao H, O'Farrell AM, Harada N, Miyajima A: Interleukin-3, granulocyte-macrophage colony stimulating factor and interleukin-5 transduce signals through two STAT5 homologs. EMBO J. 1995, 14: 1166-1175.
Arslan MA, Kutuk O, Basaga H: Protein kinases as drug targets in cancer. Curr Cancer Drug Targets. 2006, 6: 623-634. 10.2174/156800906778742479.
Bailey SN, Sabatini DM, Stockwell BR: Microarrays of small molecules embedded in biodegradable polymers for use in mammalian cell-based screens. Proc Natl Acad Sci USA. 2004, 46: 16144-16149. 10.1073/pnas.0404425101.
Peck D, Crawford ED, Ross KN, Stegmaier K, Golub TR, Lamb J: A method for high-throughput gene expression signature analysis. Genome Biol. 2006, 7: R61-10.1186/gb-2006-7-7-r61.
Cell Signaling. [http://www.cellsignal.com]
R&D Systems. [http://www.rndsystems.com]
Santa Cruz Biotechnology. [http://www.scbt.com]
American Type Culture Collection. [http://www.atcc.org]
Carroll M, Tomasson MH, Barker GF, Golub TR, Gilliland DG: The TEL/platelet-derived growth factor beta receptor (PDGF beta R) fusion in chronic myelomonocytic leukemia is a transforming protein that self-associates and activates PDGF beta R kinase-dependent signaling pathways. Proc Natl Acad Sci USA. 1996, 93: 14845-14850. 10.1073/pnas.93.25.14845.
Lessnick SL, Dacwag CS, Golub TR: The Ewing's sarcoma oncoprotein EWS/FLI induces a p53-dependent growth arrest in primary human fibroblasts. Cancer Cell. 2002, 1: 393-401. 10.1016/S1535-6108(02)00056-9.
Shipp MA, Ross KN, Tamayo P, Weng AP, Kutok JL, Aguiar RC, Gaasenbeek M, Angelo M, Reich M, Pinkus GS, Ray TS, Koval MA, Last KW, Norton A, Lister TA, Mesirov J, Neuberg DS, Lander ES, Aster JC, Golub TR: Diffuse large B-cell lymphoma outcome prediction by gene-expression profiling and supervised machine learning. Nat Med. 2002, 8: 68-74. 10.1038/nm0102-68.
Root DE, Flaherty SP, Kelley BP, Stockwell BR: Biological mechanism profiling using an annotated compound library. Chem Biol. 2003, 10: 881-892. 10.1016/j.chembiol.2003.08.009.
Tatusova TA, Madden TL: BLAST 2 Sequences, a new tool for comparing protein and nucleotide sequences. FEMS Microbiol Lett. 1999, 174: 247-250. 10.1111/j.1574-6968.1999.tb13575.x.
Gene Expression Omnibus. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/geo/]
We thank members of the Golub laboratory, the Broad Institute Cancer Program, and, in particular, K. Stegmaier, D. Peck, C. Ladd-Acosta, S. Gabriel, S. Lessnick, J. Grabarek and M. Reich for their contributions, H. Widlund and V. Igras for cell culture advice and Eric Lander for helpful discussions. AA and TG were supported by the Howard Hughes Medical Institute. BS was supported in part by a Career Award at the Scientific Interface from the Burroughs Wellcome Fund and by the National Cancer Institute (CA097061). These funding agencies played no role in the design of experiments, interpretation of data, writing of the manuscript, or planning of its submission for publication.
AA designed, performed and analyzed experiments, and drafted the manuscript. BS guided experimental design, provided chemical libraries, and assisted in writing of the manuscript. TG supervised the design and execution of experiments, and participated in writing of the manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Antipova, A.A., Stockwell, B.R. & Golub, T.R. Gene expression-based screening for inhibitors of PDGFR signaling. Genome Biol 9, R47 (2008) doi:10.1186/gb-2008-9-3-r47
- Epidermal Growth Factor Receptor
- Additional Data File
- Platelet Derive Growth Factor Receptor
- Extracellular Regulate Kinase
- Methyl Violet