Abstract
A genomics-based approach to identify pharmacodynamic biomarkers was used for a cyclin-dependent kinase inhibitory drug. R547 is a potent cyclin-dependent kinase inhibitor with a potent antiproliferative effect at pharmacologically relevant doses and is currently in phase I clinical trials. Using preclinical data derived from microarray experiments, we identified pharmacodynamic biomarkers to test in blood samples from patients in clinical trials. These candidate biomarkers were chosen based on several criteria: relevance to the mechanism of action of R547, dose responsiveness in preclinical models, and measurable expression in blood samples. We identified 26 potential biomarkers of R547 action and tested their clinical validity in patient blood samples by quantitative real-time PCR analysis. Based on the results, eight genes (FLJ44342, CD86, EGR1, MKI67, CCNB1, JUN, HEXIM1, and PFAAP5) were selected as dose-responsive pharmacodynamic biomarkers for phase II clinical trials. [Mol Cancer Ther 2009;8(9):2517–25]
- Pharmacodynamic biomarkers
- Preclinical research
- CDK inhibitors
- Translational medicine
- Biomarkers
Introduction
Dysregulation of the cell cycle is a hallmark of cancer (1). Cell cycle progression is regulated through the action of cyclin-dependent kinases (CDK) in conjunction with cyclins. CDKs typically remain at a constant concentration, whereas cyclins are present for only a short time in the cell cycle and then are degraded (2, 3). There have been efforts to develop small-molecule pharmacologic inhibitors of CDKs to block tumor cell growth through cell cycle arrest (4–11).
R547 is a potent CDK selective inhibitor (CDK1, CDK2, and CDK4; refs. 12, 13). It has a strong antiproliferative effect at pharmacologically relevant doses in 19 tested cell lines regardless of tissue type, p53, multidrug resistance, or retinoblastoma status. R547-treated cells undergo G1 and G2 arrest followed by apoptosis. Preclinical animal studies have shown acceptable toxicity and efficacy in tumor growth inhibition. R547 is currently in clinical phase I studies.
Using preclinical data, we endeavored to identify pharmacodynamic biomarkers with clinical utility. Blood cells are easily accessible and are used in clinical studies as a surrogate tissue to analyze drug activity. Therefore, a focus of this study was the identification of pharmacodynamic markers that would be particularly applicable to circulating WBC in a clinical setting. Here, we describe the identification and assessment of biomarker candidates from preclinical genomics-based studies and their testing in clinical samples. We report the identification of eight genes with potential utility as dose-responsive pharmacodynamic biomarkers for phase II clinical trials.
Materials and Methods
Reagents
R547 [4-amino-2-(1-methanesulfonylpiperidin-4-ylamino)pyrimidin-5-yl]-(2,3-difluoro-6-methoxyphenyl)methanone was synthesized by Discovery Chemistry, Hoffmann-La Roche. Its chemistry, synthesis, and crystallographic structures have been described (12).
Cell Culture
HCT116 (Rb+/p53+) and DU145 (Rb−/p53−) cells were maintained as described (13). MTT assays were used to determine dosing: for HCT116 cells, IC50 = 0.1 μmol/L R547, IC90 = 0.2 μmol/L, and 3× IC90 = 0.6 μmol/L; for DU145 cells, IC50 = 0.1 μmol/L R547, IC90 = 1.7 μmol/L, and 3× IC90 = 5.1 μmol/L. HCT116 cells were grown in vehicle (DMSO), IC50, IC90, or 3× IC90 dose of R547 for 0, 1, 2, 4, 6, and 24 h. DU145 cells were grown in vehicle (DMSO), IC50, IC90, or 3× IC90 dose of R547 for 0, 2, 6, and 24 h. Cells were harvested and lysed in Buffer RLT (Qiagen).
Human peripheral blood mononuclear cells (PBMC) were separated from whole blood of seven healthy donors and cultured in six T75 flasks per donor at 107 cells per flask in 10% heat-inactivated fetal bovine serum in RPMI. Cells from donors 1 to 6 were grown in vehicle, IC90 (0.2 μmol/L R547), and 3× IC90 (0.6 μmol/L R547) doses of R547 for 2 or 24 h. Donor 7 was used as the untreated 0 h time point. R547 concentrations were based on the HCT116 IC90 values. PBMC were scraped and centrifuged for 15 min at 1,500 rpm and the cell pellet was lysed by the addition of 1 mL Buffer RLT.
Microarray Processing
Total RNA was extracted from cells using the Qiagen RNeasy Mini Kit. Total RNA quality was assessed on the Agilent Bioanalyzer 2100.
RNA was converted into cDNA and cRNA according to the manufacturer's recommendation and using the manufacturer's kits (Affymetrix). Total RNA (15 μg) was converted to cDNA with the One-Cycle cDNA Synthesis Kit with the addition of reagents from the Poly-A RNA Control Kit. cRNA was produced with the IVT Labeling Kit. Both cDNA and cRNA were purified with the Sample Cleanup Module. Hybridization Mix contained, in addition to cRNA, reagents from the Hybridization Control Kit and Control Oligo B2.
Hybridization Mix was hybridized to Affymetrix Human Genome U133 Plus 2.0 Arrays. Staining and washing steps were done as suggested by the manufacturer (Affymetrix). Each hybridized Affymetrix GeneChip array was scanned with a GeneChip Scanner 3000 7G (Agilent/Affymetrix). Image analysis was done with the Affymetrix GCOS software.
Microarray Data Analysis
Affymetrix gene expression microarray data (GSE15396 SuperSeries, Gene Expression Omnibus, NIH) was first background-corrected and normalized via the RMA procedure followed by QQ normalization. Quality was assessed through a combination of methods including glyceraldehyde-3-phosphate dehydrogenase ratios, inspection of M/A plots, signal intensity distributions, and RMA residual plots. Quality was very high; only 3 of the 77 HCT116 samples and 7 of the 52 DU145 samples were omitted for quality.
Data for both cell line experiments (HCT116 and DU145) were analyzed in SAS (SAS Institute) via an area under the curve (AUC) method. The cell line experiments were modeled using a linear model including terms for dose, time, and their interaction. AUC was calculated from the linear model as a linear combination of terms using SAS estimate statements. Differences between the AUC values were determined via contrast. From the linear modeling results, "hits" were determined by filtering the P values and magnitudes of change for various forms of dose response. In general, the model P value was held at a threshold of <0.0001 to control for multiple comparisons; all following contrast P values were assessed at a 0.05 threshold.
For the PBMC study, a simpler linear model was used to compare the two doses to the control at the two time points. Hits were selected as probe sets with at least a 2-fold change and P < 0.05. The more relaxed P-value threshold on this data set reflects the greater variability and lower sample size available for the ex vivo samples compared with the cell line data.
Lists were then created of significant hits with specific patterns of dose response in the cell lines or significant changes reflective of dose response in the PBMC (e.g., “responsive at IC90 and above” or “responding with dose to a ceiling of the IC90”). These lists formed the basis for further analysis and candidate biomarker selection.
Process of Gene Ranking
Genes were selected for evaluation in clinical samples after consideration of their cellular activity, gene regulation, and biological pathway classification using Ingenuity Pathways Analysis. In addition, the Roche Expression Database of Affymetrix profiling experiments was examined to determine expected expression levels in the blood of cancer patients (Asterand; Proteogenix).
Clinical Trial
Patients were required to be ages ≥18 years with histologically or cytologically confirmed locally advanced solid tumors or chronic lymphocytic leukemia or refractory lymphomas. Patients with solid tumors must have had measurable or evaluable disease by Response Evaluation Criteria in Solid Tumors (14). Other key eligibility criteria include Eastern Cooperative Oncology Group performance status 0 to 2 and adequate bone marrow, hepatic, and renal function, defined as absolute neutrophil count 1,500/μL, platelet count 75,000/μL, hemoglobin 8.0 g/dL, total bilirubin level ≤1.5 mg/dL, aspartate aminotransferase or alanine aminotransferase levels ≤2 times the upper limit of normal or ≤5 times the upper limit of normal if due to liver metastases, and serum creatinine ≤1.5 mg/dL or creatinine clearance ≤60 mL/min. Patients with unstable central nervous system metastases (symptomatic disease requiring steroids or progressive disease by computed tomography/magnetic resonance imaging) were ineligible as were patients with neuropathy grade ≥2, recent myocardial infarction or congestive heart failure, history of cerebrovascular accident, or treatment with chemotherapy, radiotherapy, or immunotherapy within 3 weeks of start of study medication. Clinical characteristics are shown in Table 1.
Clinical characteristics
The clinical study was conducted at four institutions within the United States: University of Colorado Health Sciences Center, University of Texas M.D. Anderson Cancer Center, Cancer Institute of New Jersey, and Carolinas Hematology/Oncology. Patients were required to provide written informed consent before enrollment into the study. Institutional review board approvals were obtained from the participating centers before patient enrollment.
Drug was dosed on days 1 and 8 by 90-min injections (cohorts 1-7) or 180-min injections (cohorts 6A and 7A) every 21 days (Table 2). Two PAXgene tubes (PreAnalytiX, Hombrechtikon) of blood were drawn from each patient on day 1 just before injection of drug (time point 301), 3 h after the completion of injection (time point 302), and on day 2, 24 h after the start of the injection (time point 303). Samples were drawn for the second dose at the same time points 0, 3, and 24 h (304, 305, and 306) on days 8 and 9. There were some missing data points due to patients missing a blood draw or poor sample quality (degraded RNA or low yield).
Patient cohort assignment
Quantitative Real-time PCR
RNA extraction was done from the PAXgene blood on a BioRobot MDx following the manufacturer's protocol (Qiagen), with omission of the 80°C heating step. RNA was quantified with the Quant-IT RiboGreen kit (Molecular Probes, Invitrogen) at 50-fold dilution against standard curves of Escherichia coli ribosomal RNA (Roche Diagnostics).
cDNA synthesis was done with the SuperScript II First-Strand Synthesis SuperMix for quantitative real-time PCR (qRT-PCR; Invitrogen) on 400 ng total RNA following the manufacturer's protocol but with omission of the RNase H digest. For each reverse transcription reaction, Universal Human Reference total RNA (Stratagene) was run as a positive and negative control (nonenzyme control) on the same plate. Controls were assayed by qRT-PCR.
Predesigned gene expression assays were obtained from Applied Biosystems. Where possible, exon-spanning assays were selected to ensure cDNA specificity. In addition, assays were selected to lie either within the Affymetrix probe sequence of interest or as 3′ as possible within the coding sequence of the reference mRNAs of interest. One assay was custom designed. All gene expression assays were done on an ABI PRISM 7900HT Sequence Detection System (Applied Biosystems) with the recommended standard settings. All assays were prepared with the TaqMan Universal PCR Master Mix (Applied Biosystems) following the manufacturer's recommendations. All runs included a standard curve dilution of cDNA, a nonenzyme control, a nontemplate control, and a calibrator sample. Standard curve cDNA was synthesized from Universal Human Reference total RNA (Stratagene) and the calibrator sample from a pool of total blood RNA from healthy donors. cDNA samples were diluted 10-fold in molecular-grade water, and 2 μL were added to 18 μL predistributed assay Master Mix. This corresponds to cDNA from 4 ng total RNA. All samples on a plate were assayed with one assay for a gene of interest and the endogenous control gene assay, GUSB (β-glucuronidase, TaqMan Gene Expression Assays; Applied Biosystems). Each measurement was done in triplicate.
Statistical Analysis of qRT-PCR Data
ΔCt values were used for statistical analysis. The nonparametric Jonckheere-Tepstra back (J-T test) trend test was used to test for a consistently increasing or decreasing relationship between dose and gene expression level (monotonic dose-response trend with expression). Additionally, AUC- and Cmax-based tests were done, testing for differences across the doses (data not shown).
For graphical representation for data analysis, both −ΔΔCt values and fold changes were used. Fold changes are calculated from ΔΔCt with the following formula: 2−ΔΔCt. ΔΔCt is defined as the normalized value of ΔCt using time point 301 or 304 as a reference.
Results
To identify candidate biomarkers, we chose two proliferating cancer cell lines (the HCT116 colon cancer cell line and the DU145 prostate cancer cell line) and nonproliferating PBMC as a surrogate for normal tissue. We conducted a series of expression profiling experiments using Affymetrix GeneChip arrays (for detailed diagram, see ref. 15), designed to determine the dose responsiveness (vehicle, IC50, IC90, and 3× IC90) of these cells to the CDK inhibitor R547 over varying time points (0, 1, 2, 4, 6, and 24 h). Numbers of genes found to be differentially expressed in a dose-responsive manner in the different comparisons are shown in Fig. 1. These genes were narrowed to 88 by several criteria. First, all genes dysregulated on treatment with R547 common to the PBMC and either or both cell lines were included based on the hypothesis that PBMC, being nondividing primary human cells, are most similar to patient blood samples. Next, we selected genes with the highest fold changes from the genes significantly (P ≤ 0.05) dysregulated on treatment in the PBMC alone. Top differentially expressed genes (both up-regulated and down-regulated) in the overlap of the two cell lines were also included because these genes may better represent cancer cells.
Diagram showing method of microarray analysis to isolate possible biomarkers. Each experimental model is represented by a different colored circle: blue circle, HCT116 cell line; green circle, DU145 cell line; red circle, PBMC. All genes that were differentially expressed in either HCT116 or DU145 cell line and PBMC plus selected genes differentially expressed in both HCT116 and DU145 or in the PBMC comprised our candidate biomarkers (yellow).
We ranked the 88 genes based on personal expertise (15) and the following four criteria: (a) The gene was mechanistically related to the action of the compound based on literature and pathway analysis. (b) The gene was expected to be detectable in blood based on its expression level in the PBMC data set and in the Roche Expression Database containing donor whole-blood samples collected in RNA PAXgene tubes. (c) The expression levels of the gene showed a dose response. (d) Avoid genes with expression that changes in response to many oncology compounds.
The top 26 genes from this set were chosen for analysis in the human phase I patient samples (Table 3). The assays for all genes, except CCL7, showed measurable amplification in preliminary blood samples and were subsequently run on the clinical blood samples. A variety of measures were used to assess the robustness and quality of the measurements and technologies, including median Ct level, the range of technical triplicates, the directionality of the expression change as measured by gene expression microarray and qRT-PCR, and the level of significance of any changes in expression of the endogenous control gene compared with that of the target gene.
Median fold changes for genes selected preclinically for RT-PCR follow-up in patient samples from the clinical trial
After reviewing the qRT-PCR fold changes (Tables 3 and 4) and robustness of PCR assays, the genes were ranked as candidate biomarkers. The top four genes were FLJ44342, CD86, EGR1, and MK167 (Fig. 2, top). The second tier genes CCNB1, JUN, HEXIM1, and PFAAP5 are shown in Fig. 2 (bottom). These genes all showed a statistically significant dose response to R547. Nine genes showed a borderline or incomplete pattern of significance in the dose response (CALM2, CES1, CHN2, KITLG, KRAS, NME1, PRDM2, SAP30, and STAT3). Eight genes showed no correlation between treatment and expression level (BCCIP, BUB3, CDC14B, MAF, MAP3K5, PAIP1, PHDLA1, and WEE1) despite the response in the preclinical studies.
P values resulting from statistical testing of qRT-PCR results
Mean fold change within each cohort by time point for the top four biomarkers (top) and the second tier (bottom). X axis, different time points (see Materials and Methods for further discussion) at which clinical samples were taken; Y axis, fold change (calculated by −ΔΔCt) of each sample compared with its baseline time point.
Table 3 details the fold changes for each gene in the three preclinical microarray studies as well as the qRT-PCR results from the clinical blood samples. Table 4 shows the results of statistical testing for genes showing a significant dose response (boldfaced genes were significant).
The qRT-PCR results (Table 3) for cDNA FLJ44342 validated the microarray results for PBMC. CD86 (B-lymphocyte activation antigen b7-2) had a highly significant down-regulation in the higher dose groups at the 3 h time point, which returned to baseline after 24 h; this pattern was repeated in the second dose. EGR1 (early growth response 1) was up at the 3 h time point and returned to baseline by the 24 h time point as measured by qRT-PCR, although there were mixed results with microarray. MKI67 (proliferation-related Ki-67 antigen) first increased and then decreased below baseline by qRT-PCR and showed a decrease in all microarray experiments. CCNB1 (cyclin B1) was down-regulated in a dose-dependent fashion at both the 3 h time point and the 24 h time point by qRT-PCR and by microarray (data not shown). A decrease in JUN (V-JUN avian sarcoma virus 17 oncogene homologue) was seen in most conditions by both qRT-PCR and Affymetrix. HEXIM1 (hexamethylene bisacetamide-inducible protein 1) was up by qRT-PCR in these studies. PFAAP5 (N4BP2L2, NEDD4 binding protein 2-like 2) was down at the 3 h time point by qRT-PCR, although the results were mixed by microarray.
Discussion
Preclinical genomics models of the CDK inhibitor R547 have proven to be predictive of drug-responsive genes in oncology patients. A multidisciplinary process of gene selection from several preclinical Affymetrix GeneChip experiments identified a panel of potential pharmacodynamic biomarkers. Eight of these genes changed clinically in a dose-responsive fashion, indicating that the drug is reaching its target. More detailed exploration of these top eight genes revealed seven with known links to cancer and one that was unknown (FLJ44342).
cDNA FLJ44342 is interrogated by the Affymetrix probe set 226419_s_at. At the time of gene selection for this study, this probe set was annotated as belonging to a 5′ exon of SFRS1. When selecting qRT-PCR assays for each gene, we discovered that it was no longer annotated as such. We therefore designed primers and probes to the same region as the Affymetrix probe set. The transcript for this gene, Q6ZTQ9_Human, is an Ensembl Known Protein Coding region, located between the genes for VEZF1 and SFRS1 on chromosome 17q23.2 (16). The protein is not a member of any known family but does contain a low complexity region and a signal peptide cleavage site.
CD86 is expressed on the surface of antigen-presenting cells (17). It was found to be expressed at higher levels in thyroid carcinoma than in normal thyroid and there was a trend for higher tumor recurrence in CD86-expressing tissues (18). Given these findings, a cancer therapeutic agent should act to lower the levels of this gene in tissues as seen in this study both preclinically and in phase I patients.
The literature on EGR1 is mixed. Microvascular endothelial cell growth, neovascularization, tumor angiogenesis, and tumor growth are all critically dependent on EGR1 (19). Apoptosis requires intact EGR1 as well as p53 (20), but EGR1 confers resistance to apoptosis signals through inhibition of Fas expression leading to insensitivity to FasL (21). EGR1−/− mice have decreased progression of prostate cancer (22). However, studies suggest that it is a tumor suppressor gene (23). EGR1-dependent expression of TGFB1 inhibited cancer growth in cell lines (24) and some human tumor lines have little or no EGR1 compared with normals (25). The qRT-PCR and Affymetrix results from our studies fit with the expectation that an anticancer agent would increase tumor suppressor genes.
MK167 is a marker of proliferating cells (26). Its expression might be necessary for proliferation. Frameshift mutations have been found in gastric carcinoma (27) and its up-regulation has been associated with renal cell carcinoma (28). Our results showing a decrease in the expression of this gene is suggestive of slowing cell proliferation.
Up-regulation of CCNB1 is associated with numerous forms of cancer including lung cancer (29) and breast cancer (30). This gene product forms a complex with p34 (cdc2) creating maturation promoting factor driving the G2-M transition (31). CDC2, also known as CDK1, is a direct target of R547 (13). Interestingly, although there was no change in cdc2 at the mRNA level, there was a decrease in transcription of CCNB1 both in the preclinical Affymetrix studies and in the phase I qRT-PCR analyses. It has been shown that decreased CCNB1 levels via small interfering RNA led to delay in formation of the metaphase plate and of mitosis progression but not cell cycle arrest due to the loss of the spindle checkpoint (through the loss of CCNB1; refs. 32, 33).
Jun is a known oncogene. Mutant Jun is associated with breast cancer in humans (34). Chemically induced hepatocellular carcinoma produces fewer tumors of smaller size on inactivation of Jun (35). Our studies showed a decrease in this transcript in most conditions.
HEXIM1 competes with C2TA for binding to CCNT1. When C2TA binds CCNT1, it forms positive transcription elongation factor b (P-TEFb) complex that turns on transcription of MHCII genes. If HEXIM1 were up, it would block MHCII gene transcription by sequestration of CCNT1 (36). CDK9 and CCNT1 are binding partners with approximately 80% of CDK9 complexed with CCNT1 (37). Each CDK9/CCNT1 complex is an active P-TEFb molecule that can phosphorylate the COOH-terminal domain of the largest subunit of RNA polymerase II, facilitating the transition from abortive elongation to productive elongation (38). About half of P-TEFb exists in an inactive complex with 7SK (which inhibits CDK9 kinase activity) and HEXIM1 protein. The other half is competitively associated with BRD4 and is active (39). Given that HEXIM1 was up by qRT-PCR, we can hypothesize that more P-TEFb was inactive, thus decreasing the levels of overall transcription.
The function of the PFAAP5 protein is not understood, although it contains a polynucleotide kinase domain [ref. 40; CNPase domain (41)] and has been included in a 19-gene pancreatic cancer–specific expression profile (42). It is up-regulated by TP73 (43). TP73 is downstream of ATM/cABL, blocks G1-S transition, and induces apoptosis (44). Also downstream of ATM/cABL is Jun (45), which has the alternate effect on the G1-S transition.
The genes selected for qRT-PCR follow-up covered several different interesting patterns and were not always the most obvious choices from either a single statistical or biological perspective. For example, CD86 was chosen based on high blood expression level and adequate fold change in PBMC with limited biological evidence for its selection. However, in patients, this gene had a near ideal result, that is, the highest fold change and most significant response to R547 in patients. This gene was the only one of the top eight to come exclusively from the PBMC data, whereas another gene with more biological rationale for its selection, MAF, had a similar preclinical pattern but did not show a response in patients. Three of the most responsive genes in the clinical samples (FLJ44342, HEXIM1, and PFAAP5) were top hits in all three of the preclinical experiments, whereas the four remaining top hits were only derived from the two cell line experiments. This suggests that use of a wide variety of preclinical models is most effective for the identification of predictors of response in patient blood samples. Additionally, as shown for FLJ44342, even if a probe set interrogates an unnamed gene, it can still be biologically responsive and therefore, of interest.
In conclusion, our rate of success was ∼8 genes confirmed of 26 tested suggesting that, in the early phases of development, inclusion of ≥20 genes (within hospital guidelines) may be necessary to detect a pharmacodynamic signal and be most effective for narrowing a list down to a manageable number of pharmacodynamic biomarkers. In addition, these data also suggest that the CDK inhibitor R547 affects physiologic pathways in patients well enough to predict these expression changes in whole-blood samples. Finally, the data suggest that these eight genes could be used to further understand a potentially variable response to our drug among patients in phase II (efficacy) clinical trails.
Disclosure of Potential Conflicts of Interest
All authors are current or former employees of Hoffmann-La Roche.
Acknowledgments
We thank Hongjin Bian and Lyubomir Vassilev for critical discussions and the Cooperative Human Tissue Network for the tissues queried in the Roche Expression Database.
Footnotes
-
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.
- Received February 5, 2009.
- Revision received June 17, 2009.
- Accepted July 10, 2009.
- © 2009 American Association for Cancer Research.