Translational Predictive Biomarker Analysis of the Phase 1b Sorafenib and Bevacizumab Study Expansion Cohort*

Predictive biomarkers are needed to triage patients to the best therapy. We prospectively planned examination of sequential blood, biopsy, and functional imaging with which to confirm the mechanism and to identify potential predictive biomarkers in a phase Ib clinical trial expansion of patients with solid tumors receiving sorafenib/bevacizumab. The maximally tolerated doses of sorafenib at 200 mg twice daily with bevacizumab at 5 mg/kg every other week were given to biopsiable patients. Patients were randomized to receive either sorafenib or bevacizumab monotherapy for the first 28-day cycle with the second drug added with cycle 2. Biopsies, dynamic contrast-enhanced MRI, and fluorodeoxyglucose-proton emission tomography were done pre-therapy and at 2 and 6 weeks (2 weeks into combination therapy). Tumor and serum proteomics, Ras/Raf mutational analysis, and functional imaging results were examined individually and across the dataset to identify potential changes predictive of response to therapy and those that confirm the biochemical drug mechanism(s). Therapy with sorafenib/bevacizumab resulted in clinical benefit in 45% of this mixed solid tumor group. ERK activation and microvessel density were decreased with monotherapy treatment with sorafenib or bevacizumab, respectively; whereas a decreased signal over the group of total AKT, phospho(p)-VEGF receptor2, p-endothelial nitric-oxide synthase, b-RAF, and cleaved poly(ADP-ribose) polymerase was associated with earlier progression of disease. Tumor metabolic activity decreased in those patients with clinical benefits lasting longer than 4 months, and activity increased with progression of disease. Cleavage of caspase 3 and poly(ADP-ribose) polymerase was increased, and Ki67 expression decreased in patients with prolonged clinical benefits, consistent with decreased proliferation and increased apoptosis. The conglomerate analysis, incorporating pharmacodynamic and tumor biochemistry, demonstrated sorafenib/bevacizumab-targeted vascular activity in the tumor. Results suggest potential biomarkers for which changes, as a group, during early therapeutic exposure may predict clinical benefit.

Sorafenib and bevacizumab have demonstrated clinical utility as single agents or in combination with chemotherapy for solid tumors. Sorafenib, initially developed as a c-Raf kinase inhibitor, also has potent inhibitory activity against the vascular endothelial growth factor receptor-2 (VEGFR2). 1 Clinical activity has been shown for bevacizumab, the humanized neutralizing monoclonal antibody against VEGF, also alone and in chemotherapy combinations (1)(2)(3)(4)(5). The role of combining two agents with overlapping target biology had not yet been studied.
We tested the clinical hypothesis that signal interruption at collaborative pathway points, both vertical and horizontal interactions, may yield equal or greater effect than the agents in isolation in a phase I trial combining bevacizumab and sorafenib (NCT00095459), and we now report the translational analyses (6). Sorafenib was selected for its ability to target both receptor and cytosolic kinases important in a variety of activated cells in the tumor microenvironment, including stromal, endothelial, and malignant cells. Because such kinase inhibitor treatment has been shown to up-regulate production of proangiogenic cytokines, we added bevacizumab to reduce VEGF ligand availability and augment inhibition of endothelial cells. We observed the clinical benefit, including partial response and prolonged disease stabilization, using attenuated doses of the individual agents as determined by safety assessments during the trial; partial response or disease stabilization of at least 4 months occurred in 59% of the daily sorafenib cohort and in 55% of those on the intermittent, 5 of 7 days, sorafenib schedule (6,7). These benefits lasted up to 37ϩ months with over 25% of patients receiving 12 or more months of therapy. The trial prospectively planned comprehensive translational assessment using a randomized drug addition design (Fig. 1A) to evaluate individual drug target specificity and combination drug effects to iden-tify potential predictive biomarkers to examine in the ongoing phase II study of sorafenib/bevacizumab in ovarian cancer.
Predictive biomarkers are increasingly important for the advancement of targeted therapies. Such knowledge should allow more effective triage of patients to interventions more likely to provide clinical benefit. Biomarkers that predict drug response may consist of direct measures of activity, such as modulation of biochemical signals in the tumor (8) or those that yield pharmacodynamic measures, such as functional imaging (9,10). Changes in metabolic activity and/or blood flow using dynamic imaging may fall into both categories with decreased glucose uptake due to reduced glucose delivery and/or reduced glucose metabolism and altered vascular permeability in response to attenuation of the VEGF drive. Aggregate analysis of these varied translational measures may yield a more detailed view of the cancer and the drug combination, allowing broader dissection into potential predictive biomarkers. Linking modulation of activity with clinical benefit is a first step in validating prospective biomarkers.
We designed a novel drug administration schema from which to examine the contribution of both sorafenib and bevacizumab on the modulation of tumor and the tumor microenvironment behavior. The biochemical and imaging data demonstrate changes consistent with alteration of tumor vascularity, demonstrate direct association of target effect with clinical outcome across solid tumor types, and confirm the benefit of complementary pathway targeting. Reduction in blood flow, up-regulation of cytokine production, and inhibition of a set of anti-apoptotic anti-proliferative signaling events together may define potentially predictive changes to examine in early drug administration in subsequent trials.

Patients and Samples
-Details of the NCI IRB-approved phase I study and expansion cohort have been published previously (6). The maximally tolerated dose, 200 mg of sorafenib twice daily and 5 mg/kg bevacizumab every 2 weeks in 28-day cycles, was administered in this expansion cohort in a novel randomization schema ( Fig. 1A; one additional patient was enrolled after submission of the phase I report). Blood samples were obtained, processed, and stored within 4 h of sampling (6). Peripheral blood mononuclear cells were collected once for mutational analyses, and serum and plasma were collected monthly. Patients underwent elective image-guided percutaneous 18-gauge core needle biopsy, under separate informed consent and local anesthesia, by an interventional radiologist. Tissue samples were frozen immediately in OCT in the interventional radiologist suite and stored at Ϫ80°C until analysis. If the quality of the initial biopsy was poor or tumor was lacking, subsequent biopsies were aborted per protocol.

Functional Imaging
Dynamic Contrast Enhanced Magnetic Resonance (MR)-All MR images were obtained at base line and at protocol-defined points. Region-of-interest MR measurements were obtained from one selected target lesion, independent of the biopsy site, from scans within 48 h prior to each of the three biopsies. Imaging was performed on a 1.5-T MR system (GE Healthcare or Philips Achieva, Best, The Netherlands) using dedicated receive-only phased array coils. Continuous 30-s imaging data sets were obtained before, during, and after administration of the contrast medium for a total of 8 min, resulting in 23 repeated datasets. T2-weighted images (time of repetition (TR)/time of echo (TE) 4600/100 ms, a section thickness of 6 mm, 400-mm field of view, and a matrix of 320 ϫ 320) were used to locate the target tumor. Next, unenhanced T1-weighted images (TR/TE 9/3.6 ms, a 5°f lip angle, 5-mm-thick sections, 400-mm field of view, and a matrix of 256 ϫ 256) were obtained with a three-dimensional spoiled gradientecho sequence to determine the tissue T1 map. Finally, DCE-MR images (TR/TE 9/3.6 ms, a 15°flip angle, 5-mm-thick sections through the entire target lesion, 400-mm field of view, an acquisition time of 30 s per data set, and a matrix of 256 ϫ 256) were obtained with a three-dimensional spoiled gradient-echo sequence. After three base-line unenhanced image acquisitions, an automatic injector (Medrad Spectris, Indianola, PA) was used to intravenously infuse gadopentetate dimeglumine (Magnevist; Bayer Healthcare Pharmaceuticals, Wayne, NJ) at 0.3 ml/s, for a total of 0.1 mmol/kg of body weight (typically 15-20 ml), followed by a 50-ml normal saline flush.
DCE-MR Image Analysis-MR data were analyzed using a twocompartment model based on the general kinetic (GKM) Kety model (11,12). Three parameters derived from the curve-fitting GKM algorithm were used to generate quantitative parameters as follows: K trans , the forward contrast transfer rate; k ep , the reverse contrast transfer rate; V e , extravascular, extracellular volume fraction of the tumor. The DCE-MRI model incorporates an arterial input function derived from large arteries (e.g. aorta) and a T1 map for converting signal intensity to gadolinium concentration. It is based on a twocompartment model that assumes that the vascular space is in rapid equilibrium with the extravascular, extracellular space, and it further assumes a rapid water exchange between intra-and extracellular water. The GKM model was programmed in an IDL-based (Interactive Data Language; Research Systems Inc., Boulder, CO) research tool (Cine Tool, GE Healthcare).
FDG-PET-Patients fasted at least 6 h prior to the intravenous injection of [ 18 F]fluorodeoxyglucose (15 mCi). Emission images (8 min) were obtained in two-dimensional mode from the upper thigh to the base of skull starting ϳ60 min after injection. Transmission scans (3 min) were obtained for attenuation correction. Scans were performed using a GE Advance scanner (General Electric Medical Systems) with a 15-cm axial field of view. PET images were reconstructed on a 256 ϫ 256 matrix using an iterative algorithm provided by the manufacturer. SUV values corrected for lean body mass were obtained using the maximum pixel activity value within a region of interest drawn over index lesions. Index lesions assessed on MR were identified and SUV measurements performed.

Protein Analysis and Quantitation
Immunohistochemistry-Serial frozen sections for the set of three time points were fixed in acetone and incubated with primary antibody (Table I). CD31-stained slides were pre-heated in a microwave oven for 5 min; no other antigen retrieval was done. Visualization was achieved using the Dako EnVision TM ϩ peroxidase system. Appropriate positive and negative controls were included in each staining experiment, and all specimens containing Ͻ100 tumor cells were excluded. Staining was considered positive when localized to the nucleus (Ki-67 and p85-PARP) or cytoplasm (VEGF, CD31, and p85-PARP). Staining extent was scored blinded by a gynecopathologist (BD Biosciences) for the number of positive cells on a scale of 0 -4 for Ki-67, p85-PARP, and VEGF as follows: 0 ϭ no staining, 1 ϭ 1-5%, 2 ϭ 6 -25%, 3 ϭ 26 -75%, and 4 ϭ 76 -100% of tumor cells (13). CD31-positive vessels are presented as the average number of vessels per 10 high power fields (ϫ400 magnification).
RPPA-Tissue quantity assessment and quantitation, lysis, serial dilution, and antibody analysis was done according to our optimized procedure (14). The core from each biopsy set with the greatest percentage of tumor, and the least necrosis was selected for use after biopsy pathology was reviewed (Fig. 2), and tissue area was measured using the Veritas laser capture microdissection system (Arcturus, Sunnyvale, CA). Final tissue quantity was 8 m ϫ 30 mm 2 , which approximated 30 g of total tissue. Two independent replicate tissue sample sets from each set of biopsies were used. Lysates were randomly distributed onto the RPPA slides in triplicate, undiluted and serially diluted 1:2, on 40 replicate nitrocellulose membrane-coated slides (Whatman) by an Aushon 2470 arrayer (Aushon, Billerica, MA).
Internal controls included lysates from human microvascular endothelial cells, HeLa cells, and A431 cells treated with growth factors or apoptosis-inducing agents (15) as follows: human microvascular endothelial cells Ϯ VEGF (50 ng/ml for 2 min); HeLa cells Ϯ etoposide (25 M for 5 h), or EGF (100 ng/ml for 2 h), and untreated A431 cells for vascular cell activation, apoptosis, and receptor tyrosine kinase pathway activation, respectively. A sample buffer-only negative control was printed with each dilution replicate for ambient background.
Target proteins were detected by specific validated antibodies from the same lot tested for validation and optimization (Table I), using an avidin-biotin amplification system and stained with diaminobenzidine (DAKO, Carpinteria, CA). Colloidal gold stain (Bio-Rad) was used to quantitate total protein load per spot. Stained slides were scanned; spot intensity was measured, and intensity values were normalized to total protein as described (14).

Mutational Analysis
For DNA extraction, 8-m tissue recuts, providing a total area of 10 -20 mm 2 , were provided to the Clinical Molecular Profiling Core, Genetics Branch, Center for Cancer Research, NCI, National Institutes of Health (DA). OCT was removed by washing samples in 500 l of TE and centrifuging the tissue for 10 min (500 ϫ g). DNA was then extracted from the tissues with the QIAamp DNA Micro Kit (Qiagen, Valencia, CA). Mutational analysis was done using standard Sanger DNA sequencing protocols after amplification and applied M13tagged primers against all exons in BRAF, KRAS, NRAS, and HRAS.
Amplicons were evaluated by gel electrophoresis, cleaned (Exo-Sap It; United States Biochemical Corp., Cleveland, OH), and then sequenced using BigDyeா terminator version 3.1 cycle sequencing kit (Applied Biosystems). GenBank TM reference sequences NM_ 004333.4, NM_033360.2, NM_002524.3, and NM_005343.2, respectively, were used in the DNA sequence analysis (Mutation Surveyor software, Softgenetics, State College, PA) with default settings. Single nucleotide polymorphisms without commensurate amino acid change did not require repeat sequencing, although base pair changes resulting in a mutant protein underwent repeat sequencing for confirmation.

Statistical Considerations
Ten patients per arm, each with three sequential usable biopsies, were needed to provide 80% power to detect a difference equal to 1.5 S.D. of the change from base line for eight primary parameters with an  2. Example of core needle biopsy specimens. Representative 18-gauge core needle biopsy specimens pre-treatment, after 2 weeks of monotherapy, and at 6 weeks is shown. Patient E (Pt E) received sorafenib for cycle 1 monotherapy and patient R (Pt R) received bevacizumab monotherapy on cycle 1. Slides are stained with H&E, and this is a Veritas laser capture microdissection digital picture taken during tissue assessment and measurement. Arrows point to areas of viable tumor.
overall ␣ ϭ 0.05 per arm. All results are considered exploratory and all tests used nonparametric statistics, with a Hochberg adjustment for multiple comparisons. Relative changes from base line or from week 2 were used in the analyses, as these were less dependent on the earlier time point values than were absolute changes. Continuous data between the two groups were examined with an exact Wilcoxon rank sum test. Spearman correlation analysis determined the correlation between continuous parameters. The strength of the correlation was considered to be more important in interpretation than were the p values, which test if r ϭ 0. Strong correlation is indicated by ͉r͉ Ͼ 0.70, moderate if 0.5 Ͻ ͉r͉ Ͻ0.7, and of decreasing strength if below 0.5.
Protein intensity fold-change ratios were calculated from the RPPA intensity values. Ratios were calculated between tumors sampled prior to initiating single agent administration, after 2 weeks of single agent (2 weeks/0 weeks ratio), and after 2 weeks of dual agent administration (6 weeks/2 weeks for second drug addition and 6 weeks/0 weeks for total change over time). These ratios were mediancentered and clustered using Cluster 2.0 software. Student's t test for differences between average protein intensity ratios at these intervals between patients in cluster 1 versus cluster 2 was calculated using Microsoft Excel by grouping all the ratios for individual clusters. A second analysis using Student's t test evaluated the ratio of signal at 6 weeks versus on-study to identify a core group of putative biomarkers of treatment interval. Time-to-event end points, such as progression-free survival, were computed using Kaplan-Meier statistics, with comparisons made using a log-rank test.

RESULTS
Patients and Specimens-Twenty eight patients were enrolled in the expansion translational cohort to get 10 patients with triplet biopsies (Consort diagram, Fig. 1B). Patient demographics, tumor characteristics, and clinical activity are shown in Table II. After final analysis, only 19 patients had three serial biopsies with at least one core per time point of adequate tissue quality and quantity for use. Reasons for incomplete biopsy sets included poor quality, fluid only, or Յ50% tumor cells (1 patient each), refusal (1), safety (2), or study removal due to toxicity or disease progression (3). Twenty six patients had adequate peripheral blood mononu-clear cell ascertainment and serial monthly blood sampling (aggregate VEGF concentration results were reported (6)). Twenty three patients had all three planned [ 18 F]FDG-PET scans, and 16 patients had all planned MRs; most common causes for incomplete imaging series were obesity and claustrophobia for MR, and scheduling for out-of-town patients.
Biochemical Illustration of Mechanism-RPPA was used to confirm the presence, activation, and inhibition of select putative biochemical/predictive targets. ERK activation, the Raf kinase downstream event, was significantly reduced with 2 weeks of sorafenib but not in patients receiving bevacizumab (p ϭ 0.02; Fig. 3A). Patients randomized to single agent bevacizumab had a lower CD31 vessel count at 2 weeks (p ϭ 0.05; Fig. 3B), consistent with the greater activity of ligand neutralization. Increased circulating VEGF concentrations have been observed with multiple inhibitors of angiogenic signaling pathways (5,16), and were confirmed in this patient cohort (6). Consistent with this, we found decreased tissue pERK at 2 weeks correlated with higher circulating VEGF concentration at both 2 and 6 weeks (r ϭ Ϫ0.60, p ϭ 0.024; r ϭ Ϫ0.63, p ϭ 0.0501, respectively; Table III).
Combination Therapy Results in Reduction in Vascular Flow and Metabolism-Functional imaging was used to examine metabolism and vascular flow. A trend for reduced index lesion FDG-PET SUVs at 6 weeks was seen in patients who attained a PR or had SD for Ն4 cycles (p ϭ 0.047; Table III). Vascular perfusion by MR-K trans increased early (2 weeks; r ϭ 0.72, p ϭ 0.042) in patients in whom increased VEGF staining was observed in the 6-week biopsy, a pharmacodynamic validation of the biochemical observation. Vascular permeability, k ep , increased at 6 weeks in patients who had increased CD31-positive microvessels at 2 weeks (r ϭ 0.75, p ϭ 0.02). The increase in k ep in these tumor areas could be due to increased microvessels, increased permeability of those microvessels, or both. k ep was lower at the 2-week mark in patients randomized to receive sorafenib first (p ϭ 0.01), suggesting that it had a detectable anti-vascular effect. Similar effects were not seen with bevacizumab; this may be in part due to imaging being done after a single dose of 5 mg/kg bevacizumab, below the doses for which functional imaging changes have been reported. These individual comparisons show the mechanism of the combination and suggest end points for further analysis.
Association with Disease Behavior-Biomarker value is based upon linkage to a pertinent clinical event. We next examined the relationship between the biochemical and pharmacodynamic findings and clinical benefit (Table III). A strong relationship between decreased p-ERK (p ϭ 0.011) and p-AKT (p ϭ 0.015) was seen with combination treatment in cases where significant tumor regression or necrosis was present in post-treatment biopsies. Tumors of patients who received more cycles of therapy expressed less cleaved PARP at base line (r ϭ Ϫ0.45, p ϭ 0.053); this might be interpreted to mean there was more viable tumor for re- sponse. Changes in cleaved PARP correlated with patient receipt of a greater number of cycles (r ϭ 0.59; p ϭ 0.0072), suggesting increased apoptosis was associated with longer progression-free intervals. Similarly, a moderate correlation was found between the early apoptosis event of change in caspase 3 cleavage (r ϭ 0.50, p ϭ 0.03) and decreased Ki67 staining with longer time on treatment (r ϭ Ϫ0.70, p ϭ 0.016), both of which would infer net tumor loss. In addition, patients receiving more cycles of therapy also had decreased levels of pVEGFR2 at baseline (r ϭ 0.45, p ϭ 0.053); this has also been observed in other studies suggesting that anti-angiogenic activity may be more active with moderate angiogenesis rather than robust activity. A trend was observed between decreased PET SUV values on monotherapy in patients who received more cycles (r ϭ Ϫ0.41, p ϭ 0.044). This suggests that PET changes at 2 weeks might be a predictive biomarker. Continued reduction in SUV values on the 6-week PET scan with combination therapy occurred in patients who had a partial response or prolonged S.D. (p ϭ 0.047). Together, these findings suggest an overall inhibition of angiogenesis as the driving mechanism.
The value of a comprehensive translational data collection is the ability to examine the contribution of molecular and functional results in aggregate. We next explored the target protein changes in the patients as a grouped analysis. Unsupervised hierarchical clustering of both proteins and patients segregated into two major clusters (Fig. 4A). Next, the relative change in protein intensity with a single agent was compared with dual agent administration in patients from the two different clusters (Fig. 4B). The overall protein changes between the two groups of patients was statistically significantly different (p ϭ 0.00006), with tumors from cluster 1 patients having an average decrease in protein intensity values with addition of the second agent, and an overall increase in tumors from cluster 2 patients. This difference correlated with clinical benefit, as indicated by the number of cycles each patient received (Fig. 4C). Patients in cluster 1, whose tumors responded to the addition of the second agent with a decrease in protein intensity, had a statistically significantly longer time on the study (p ϭ 0.01) and thus progression-free survival. Within this analysis, relative to cluster 1, reduced apoptosis (decreased cleaved PARP) and reduced total BRAF, pVEGFR2, AKT, and p-endothelial nitric-oxide synthase grouped together as a potential poor outcome biomarker (p ϭ 0.011; Fig. 4D).
Other Planned Analyses-Angiogenic potential of patient sera in a rat aortic ring assay was unrevealing (17). Tumor KRAS and BRAF mutational analysis revealed one patient with a codon 12 KRAS mutation; this patient could not be assessed for clinical benefit due to removal from the study for uncontrolled hypertension. CYP3A4 SNP analyses yielded no significant differences between patients (18 -20). DISCUSSION Increasingly, the therapeutic potential for targeted agent combinations is being recognized. How best to combine agents is still a quandary, as is how to select the best targeted therapies for patients with recurrent solid tumors. Illustration of target affected by molecular therapeutics, linkage of target modulation with clinical benefit, and determination of predictive utility are needed. Following our ongoing interest in targeting both the tumor and its microenvironment, we elected to combine the VEGF-neutralizing monoclonal antibody bevacizumab with the c-RAF/VEGFR2-targeted sorafenib, hypothesizing that these agents together could affect tumor cells, stromal cells, and vascular elements. To address those objectives, we incorporated serial tumor biopsies, functional imaging, and blood collections to provide biochemical and clinical measures of target effect and anti-tumor activity. The trial design, novel at the time, randomized an initial round of monotherapy for target validation, followed by continued combination therapy with which to examine potentially predictive changes. We show here that the expected targets were hit and correlated with clinical benefit, as defined by continued study treatment (treatment interval). Aggregate proteomic analysis described a panel of biochemical changes in the tumor, occurring between pretreatment and 6 weeks that statistically and significantly correlated with treatment interval. The set of end points identified in Fig. 4D will be examined for potential validation in prospectively collected samples in the phase II trial of sorafenib and bevacizumab in ovarian cancer patients that have recently completed accrual.
An ideal predictive biomarker is one that can be applied prior to treatment and will accurately forecast response to therapy (8,21). They can be used for trial enrichment, patient triage, and to describe biosimilars. Incorporation in our trial of a monotherapy element and multiple translational end points, coupled with the clinical benefits observed, provided the template for predictive biomarker discovery. We validated both FDG-PET and DCE-MRI as measures of tumor vascularity in the context of sorafenib and bevacizumab treatment. Change in vascular flow as estimated by functional imaging has been shown to be predictive of the benefit of imatinib when used in gastrointestinal tumors (22), and there are reports where FDG-PET changes can be predictive of response to bevacizumabbased therapy for glioma, colorectal, and head and neck cancers (23)(24)(25). In our study, both functional imaging modalities contributed to the demonstration that the combination of sorafenib and bevacizumab affected their benefit through an anti-angiogenic drive.
Incorporation of a discovery of the proteomics component with assessment of multiple pharmacodynamic end points provided a dataset from which we could look for changes that would predict prolonged treatment intervals. Aggregate anal- ysis of the proteomic changes against treatment duration segregated two groups in an unsupervised clustering. Biochemical changes between pretreatment and 2-week monotherapy that clustered in group 2, the short treatment interval group, include reduction in B-RAF, pVEGFR2, total EGF receptor, and total AKT. The most striking changes, shown in  (1) and addition of the second agent (3) identifies statistically significant (p ϭ 0.00006) differences between clusters 1 and 2 from (A). C, Kaplan-Meier statistics demonstrate different clinical outcome for clusters 1 and 2 (p ϭ 0.01) indicating that an increase in signaling events with addition of the second agent predicts poor performance to therapy. D, exploratory predictive biomarker. Further analysis of changes occurring at 6 weeks relative to on-study yielded a subset of proteomic changes that together and individually were attenuated by patients in cluster 1 who had clinical benefit and were increased in cluster 2 patients who did not respond to therapy. This identifies a potential biomarker subset for further validation. this potential negative biomarker group could be equally important to moving agents such as these forward in combinations.
It is unclear whether we can improve clinical benefit by using targeted agents in series (vertically) or in parallel (horizontally). This study incorporated putative interactions in both directions, VEGF/VEGFR2 and VEGF/c-RAF, that were then queried as to target inhibition and for predictive value in the aggregate analysis (Fig. 5). Biochemical parameters differed with single agent therapy; those differences were lost and broader overlapping signal inhibition was observed upon the combination of the two agents. This breadth of altered regulation by the combination could, in part, explain the increase in toxicity and also the clinical benefit seen with the addition of bevacizumab to sorafenib. Hand-foot syndrome, a recognized sorafenib toxicity, was increased in frequency with addition of bevacizumab and then further increased when bevacizumab dose increased. Furthermore, there is a possibility of increasing value of an otherwise minimally active drug when used in a rational biochemical combination. The 19 ovarian cancer patients treated on this phase 1 trial had a 47% response rate with a progression free survival median of 6 months (mean, 11 months; range 4 -37) (7). Other studies incorporating single agent bevacizumab in recurrent ovarian cancer patients had response rates ranging from 16 to 21% with median duration of 4.4 -4.7 months (26,27), and single agent sorafenib has limited activity in ovarian cancer (28,29).
Identification of predictive biomarkers for anti-angiogenic therapy is a critical need. Approaches taken to date have included genomic and genetic studies, and proteomic studies at the serologic and tissue levels, as in our study. Several investigators have suggested that SNPs in the VEGF pathway genes may have predictive value. Recently, Lambrechts et al. (30) evaluated SNPs in the VEGF pathway in a subset of samples from the randomized placebo-controlled trials, AViTA (BO17706) and AVOREN (BO17705), randomizing Ϯ bevacizumab for gemcitabine/erlotinib in pancreatic cancer and interferon ␣2a for renal cell carcinoma, respectively. Their discovery phase identified a synonymous SNP affecting Tyr 1213 -VEGFR1 in its tyrosine kinase domain as associated with progression-free survival in the bevacizumab group (HR 2.1; p ϭ 0.00014); this was validated independently as predictive in the progression-free survival of the bevacizumab group in AVOREN (HR 1.81, p ϭ 0.033). This SNP increased VEGFR1 expression and downstream VEGFR1 signaling. Serologic proteomics has been applied by several groups. Prespecified serologic end points were analyzed in the AVAGAST randomized placebo-controlled trial of the role of bevacizumab with cisplatin/capecitabine therapy for advanced gastric cancer. They show that high pretreatment circulating concentrations of VEGF-A and low pretreatment circulating neuropilin-1 were prognostic for improved survival (HR 0.7) (31). A similar survival prognostic, but not predictive, value was observed for VEGF-A concentrations in a series of other studies of bevacizumab Ϯ chemotherapy in solid tumors (32,33). Angiome profiling was also done in a gemcitabine Ϯ bevacizumab-randomized trial (CALGB 80303) identifying a different survival prognostic signature in the bevacizumab versus placebo arms, both of which included IGFBP1 and PDGF-AA (34). Our study is the first to look at a broad array of activated protein end points over time in tumor tissue from which to develop hypothesis-generating predictive angio-biomarkers.
Illustration of biochemical mechanism and identification of predictive biomarkers remain a challenge. Biomarkers can be genetic, epigenetic, and/or proteomic, and investigators must balance advantages and disadvantages of the multiple analysis platforms for each. We focused on the RPPA format to assess biomarkers at the protein level, where the effects of genetic, epigenetic, and post-translational modifications coalesce to truly affect cellular function. We applied our previously reported RPPA optimization methods and controls (15). Using an approach containing monotherapy and serially sampled translational end points, we confirmed on-target activity of activated ERK for sorafenib, and microvessel density as a marker of tumor vascular response for bevacizumab. Although neither result was unexpected, the analysis indicated that the biological effect of the combination was not broadly against tumor, stroma, and vascular elements as hypothesized, but prominently against the vasculature defining this as a true anti-angiogenic combination. End point analysis identified a panel of biochemical changes measured in the tumor tissue that correlated significantly with clinical benefit to the combination. This biomarker panel will be examined prospectively as a potential predictive biomarker change in the ongoing phase II study of sorafenib and bevacizumab in ovarian cancer NCT00436215. The ever-increasing classes of targeted agents available for single and combination studies preclude examination of every possible combination. Currently, rational combina- tions must be designed based upon cellular target, interactive biochemical and/or biological targets, and toxicity profiles. This is an example where target, inhibition, and outcome harmonized in a tolerable and active regimen.