Epidermal growth factor receptor (EGFR) and COX-2 inhibitors synergistically inhibit head and neck squamous cell carcinoma tumorigenesis in preclinical studies. We conducted a phase I and pharmacokinetic study with the erlotinib and celecoxib combination in patients with advanced premalignant lesions. Thirty-six subjects with oral leukoplakia, mild, moderate, or severe dysplasia, or carcinoma in situ were screened for study participation; 12 consented and received therapy for a median of 5.38 months. Erlotinib was escalated following a standard 3+3 design at 50, 75, and 100 mg orally daily and celecoxib was fixed at 400 mg twice daily for 6 months. Biopsy of lesions and cytobrush of normal mucosa were performed at baseline, 3, 6, and 12 months. Erlotinib pharmacokinetics were analyzed in 10 subjects. The maximum tolerated dose of erlotinib with celecoxib 400 mg BID was 50 mg per day with skin rash being the main observed toxicity. Overall histologic response rate was 63% (complete response, 43%; partial response, 14%; stable disease, 29%; and disease progression, 14%). With median follow-up of 36 months, mean time to progression to higher-grade dysplasia or carcinoma was 25.4 months. Downregulation of EGFR and p-ERK in follow-up biopsies correlated with response to treatment. Larger average erlotinib V/F (approximately 308 L) and CL/F (8.3 L/h) compared with previous studies may be related to relatively large average bodyweights. Average erlotinib t1/2 was 25.6 hours. Encouraging responses to the celecoxib and erlotinib combination correlated with EGFR pathway inhibition. Although erlotinib-related rash was the main limitation to dose escalation, the intervention was well tolerated. Cancer Prev Res; 7(3); 283–91. Β©2013 AACR.

See related commentary by Wirth, p. 279

The epithelial growth factor receptor (EGFR) is expressed in a wide variety of malignant tumors, including head and neck, colon, pancreatic, non–small cell lung, breast, kidney, ovarian, bladder carcinomas and gliomas (1–3). The incidence of EGFR expression in head and neck squamous cell carcinoma (HNSCC) is over 90%, suggesting that EGFR inhibition may be effective in HNSCC (4–6). One approach to block EGFR activity involves the use of small molecule tyrosine kinase inhibitors (TKI) that target the intracellular domain of EGFR. This ability to specifically inhibit intracellular tyrosine kinase activity has been observed over precise dose ranges (7, 8).

A broad range of laboratory investigations, animal models, and epidemiologic studies provide evidence that inhibition of COX-2 pathways may contribute to cancer treatment in general (9–11) and HNSCC in particular (12, 13). In HNSCC, COX-2 is expressed in both tumor tissue and adjacent epithelium, with increased expression in invasive carcinoma compared with normal epithelium. COX-2 inhibition has been shown to result in cell growth inhibition in HNSCC cell lines (14). Furthermore, we have previously elucidated the differential expression pattern of COX-2 in stages of head and neck premalignant lesions and invasive carcinoma (15). Our findings indicated that COX-2 is involved in early and intermediate stages of carcinogenesis in HNSCC. COX-2 levels increased progressively throughout all stages of carcinogenesis. This may reflect a role for COX-2 in this process, further supporting the rationale for COX-2 inhibition as a valid strategy for cancer chemoprevention. Based on this evidence, a number of chemoprevention and therapeutic trials in HNSCC using COX-2 inhibitors (COX-2I) are under way (16). Evidence suggests that COX inhibitors, including nonsteroidal anti-inflammatory drugs (NSAID), protect against a variety of tumors (9, 17, 18). In patients with familial adenomatous polyposis, celecoxib, a selective COX-2I, caused a 28% reduction in the number of colorectal polyps compared with 4.5% reduction for placebo (19). In HNSCC, Western blot analysis showed expression of COX-2 in mucosa of subjects at different stages of carcinogenesis and not in normal mucosa, suggesting a possible role for COX-2 inhibition in HNSCC chemoprevention (20). Previous studies suggested that NSAIDs might have a similar effect in delaying the growth of head and neck tumor cell lines (21–23). Furthermore, prior studies using the COX-2I as a single agent in oral premalignant lesions revealed evidence of improvement in the degree of dysplasia (24).

Preclinical and clinical studies have described the interaction between EGFR and COX-2 and shown that targeting these two pathways can synergistically or additively block progression of HNSCC growth in vitro and in vivo (14, 25, 26), providing a scientific basis for using the combination of EGFR-TKI and COX-2I as a chemopreventive approach in HNSCC. We report here the results of a phase I clinical trial and pharmacokinetic studies of this combination in subjects with premalignant lesions (i.e., leukoplakia, erythroplakia, and/or erythroleukoplakia) of the oral cavity oropharynx or larynx. The correlation of response with biomarker modulation was reported in a separate publication (27).

Patient accrual

Between October 24, 2006, and June 28, 2012, 36 subjects with documented premalignant lesions, including mild (mild-D), moderate, or severe oral leukoplakia, and carcinoma in situ (CIS) were screened. Lesion sites included oral cavity oropharynx, and the larynx. Participants were eligible regardless of their smoking status. Eligibility requirements included Eastern Cooperative Oncology Group (ECOG) performance status of 0 to 1, age β‰₯18 years, adequate bone marrow, liver and renal function, signed written informed consent, negative serum pregnancy test (Ξ²-HCG) within 72 hours of receiving treatment, willingness to use appropriate contraception during study participation for women of child-bearing potential, adequate pulmonary function (FEV1 and Forced vital capacity β‰₯60% by spirometry), adequate cardiac function (echocardiogram with normal Left Ventricular ejection fraction).

Subjects were excluded from participation if they had documented hyperplasia only, acute intercurrent illness or recent major surgery, history of previous malignancies excluding stage I or II cancers rendered disease-free more than 1 year from time of consent, documented pregnancy, were breast feeding, had active cardiovascular events, including angina, unstable angina, palpitation, tachycardia, arrhythmia, a recent cerebrovascular accident or myocardial infarction within 6 months from enrollment, documented history of coagulopathy and/or were taking warfarin or warfarin-derivative anticoagulants, hypertension not adequately controlled by medication, history of congestive heart failure greater than New York Heart Association (NYHA) grade 2, confusion, disorientation, or a history of major psychiatric illness impairing their understanding of the informed consent, history of intake of COX-2I or EGFR-TKI within 3 months of study entry, documented history of interstitial lung disease or known connective tissue disease, history of NSAID-induced ulcers or participation in a clinical trial using an investigational drug within 12 months of enrollment.

Study participants were required to have a complete inspection of the oral cavity oropharynx and larynx. A baseline biopsy for initial diagnosis and grading was mandatory. All evaluable patients were required to have repeated biopsies at 3, 6, and 12 months from initiation of therapy. Biopsies of suspicious lesions in the oral cavity or laryngeal lesions were performed using standard β€œpunch biopsy” procedure with approximately a 3-mm punch. If a follow-up biopsy was deemed difficult to obtain at the suspicious site, a brushing or wash was performed instead. Histologic assessments were performed by a head and neck pathologist. Buccal scrapings (cytobrush) for target lesions and normal buccal mucosa were performed at baseline, 3, 6, and 12 months as surrogates for biomarker modulation.

Drug treatment

Participants received a fixed dose of celecoxib 400 mg orally BID continuously for 6 months. Erlotinib was dose escalated at 3 dose levels of 50, 75, and 100 mg orally every day for 6 months. Dose escalation followed a standard 3+3 escalation design.

Definition of response

Response evaluation was based on pathologic examination of the degree of dysplasia observed and recorded by an expert head and neck pathologist. Pathologic complete response was defined as complete disappearance of dysplasia from the epithelium. Pathologic partial response was defined as improvement of dysplasia by at least one degree (i.e., severe dysplasia becomes moderate dysplasia). Pathologic minor response or stable disease was defined as minor focal improvement without change of degree of dysplasia (i.e., focal improvement from moderate to mild dysplasia with still moderate dysplasia overall) or no pathologic changes after treatment. Pathologic progressive disease was defined as worsening by at least one degree of dysplasia (i.e., mild to moderate dysplasia) or development of invasive cancer on or following treatment.

Pharmacokinetic studies

Blood sampling for measurement of erlotinib pharmacokinetics was performed before drug administration (t = 0), and 0.5, 1, 3, 6, 8, 24, and 48 hours after administration. Blood was centrifuged using a refrigerated centrifuge and plasma collected and frozen before being assayed for erlotinib.

Assays of erlotinib and active metabolites

Erlotinib was extracted from plasma samples and plasma concentrations were measured by liquid chromatography/mass spectrometry (LC/MS), as detailed in Supplementary Materials. A standard curve was constructed using weighted (1/X2) linear regressions of the peak area (Y) of erlotinib against the corresponding nominal concentrations of erlotinib (X, ng/mL) in blank plasma. The liquid chromatography/tandem mass spectrometry (LC/MS-MS) assay was validated with specificity, precision (coefficients of variation <15%), accuracy (>85%), matrix effects, and linearity (0.1 to 500 ng/mL; r > 0.99). The limits of detection and quantitation for erlotinib were 0.5 and 2 pg on the column. Recoveries of erlotinib from human plasma were 80% to 82%.

Pharmacokinetic data analysis and modeling

Plasma concentrations of erlotinib were used for pharmacokinetic modeling. Assay results (confirmed by patient drug diaries) indicated variations with respect to timings of erlotinib administration, and certain individuals were missing blood samples at critical times, making the data unsuitable for compartmental modeling on a per individual basis. However, modeling was feasible using population pharmacokinetic methods with the NONMEM modeling program (28). The population pharmacokinetic of erlotinib was previously modeled in patients with cancer enrolled in phase II studies, using an open one-compartmental model with first-order absorption into plasma (29–31). Here, a similar modeling approach was used to study the pharmacokinetics of erlotinib in comparatively healthy subjects with premalignant lesions, using the NONMEM software (ver 7.1; ICON Development Solutions; ref. 32) run with PLT tools (version 2.6; PLTsoft). Goodness of fit was analyzed using commonly used graphical methods for population pharmacokinetic, and the influence of body weight on V/F and CL/F was explored using previously reported covariate model structures (29–31). Methods are detailed in Supplementary Materials.

Patient characteristics

A total of 36 subjects with pathologic dysplasia (mild-D, moderate, and severe) or CIS were screened. Nineteen screened subjects were not enrolled due to social or personal reasons mostly related to reasons of convenience such as transportation or other personal commitments, in addition to concerns about toxicity, comorbidities, or ineligibility. One patient with mild-D was enrolled after the protocol was amended to include mild-D. Seventeen subjects were enrolled on the study, 3 of whom withdrew consent (one male of 69 years of age, and two females of 43 and 42 years of age). Two patients who signed informed consent were deemed to be screen failures, one secondary to prior history of oral squamous cell carcinoma and the other secondary to a history of cardiac arrhythmias. A total of 12 subjects received therapy. None of the 12 patients who received therapy had a history of prior treated malignancy. Their characteristics with baseline pathology, response, and duration of response are described in Table 1. Two subjects were taken off treatment before response evaluation due to the following reasons: grade 3 rash, elevated serum creatinine, or urosepsis. Two subjects chose to withdraw from therapy before first response evaluation. One patient was excluded from the efficacy analysis as she was found to have microinvasive carcinoma on her resected pretherapy tissue biopsy. A total of 7 subjects were evaluable for response.

Table 1.

Evaluable patient characteristics, baseline pathology, pathologic response, and time to progression

PatientAge (y)GenderSiteLesionPathologic responseProgression/discontinuationTime to progression or last evaluation (months)/cause of withdrawal
1Β 50Β FΒ TongueΒ MDΒ CRΒ YΒ 4Β 
2Β 42Β MΒ TongueΒ MDΒ PRΒ YΒ 26Β 
3Β 41Β FΒ TongueΒ SDΒ SDiΒ YΒ 26Β 
4Β 65Β FΒ TongueΒ SDΒ PDΒ YΒ 6Β 
5Β 62Β MΒ Buccal mucosaΒ SDΒ CRΒ YΒ 55Β 
6Β 64Β MΒ LarynxΒ Mild-DΒ SDiΒ YΒ 6Β 
7Β 40Β MΒ LarynxΒ MDΒ CRΒ NΒ 36Β 
8Β 69Β FΒ TongueΒ HGΒ NEΒ D (elevated creatinine)Β Renal insufficiencyΒ 
9Β 74Β FΒ TongueΒ MDΒ NEΒ D (grade 3 rash)Β RashΒ 
10Β 87Β FΒ BOTΒ SDΒ NEΒ D (withdrew)Β Patient chose to withdrawΒ 
11Β 46Β FΒ TongueΒ HGΒ NEΒ D (withdrew)Β Withdrew after signing consentΒ 
12Β 57Β FΒ TongueΒ HGΒ N/AΒ N/AΒ Excluded from efficacy analysisΒ 
PatientAge (y)GenderSiteLesionPathologic responseProgression/discontinuationTime to progression or last evaluation (months)/cause of withdrawal
1Β 50Β FΒ TongueΒ MDΒ CRΒ YΒ 4Β 
2Β 42Β MΒ TongueΒ MDΒ PRΒ YΒ 26Β 
3Β 41Β FΒ TongueΒ SDΒ SDiΒ YΒ 26Β 
4Β 65Β FΒ TongueΒ SDΒ PDΒ YΒ 6Β 
5Β 62Β MΒ Buccal mucosaΒ SDΒ CRΒ YΒ 55Β 
6Β 64Β MΒ LarynxΒ Mild-DΒ SDiΒ YΒ 6Β 
7Β 40Β MΒ LarynxΒ MDΒ CRΒ NΒ 36Β 
8Β 69Β FΒ TongueΒ HGΒ NEΒ D (elevated creatinine)Β Renal insufficiencyΒ 
9Β 74Β FΒ TongueΒ MDΒ NEΒ D (grade 3 rash)Β RashΒ 
10Β 87Β FΒ BOTΒ SDΒ NEΒ D (withdrew)Β Patient chose to withdrawΒ 
11Β 46Β FΒ TongueΒ HGΒ NEΒ D (withdrew)Β Withdrew after signing consentΒ 
12Β 57Β FΒ TongueΒ HGΒ N/AΒ N/AΒ Excluded from efficacy analysisΒ 

Abbreviations: BOT, base of tongue; D, discontinued therapy before the first evaluation of response; F, female; HG, high-grade dysplasia; M, male; MD, moderate dysplasia; N, no progression as of the last follow-up; NE, nonevaluable; SD, severe dysplasia; SDi, stable disease; Y, documented progression.

Dose escalation and toxicity

Three subjects were enrolled on cohort 1 (erlotinib 50 mg daily) with no dose limiting toxicities (DLT) observed, allowing escalation to cohort 2 (erlotinib 75 mg daily). One patient had a grade 3 rash at dose level 2 and cohort 2 was expanded to 6 subjects. As a second subject had grade 3 rash, the maximum allowable dose of erlotinib was deemed to be 75 mg. One patient was enrolled at a dose level 3 (erlotinib 100 mg daily) before toxicity analysis of the patients who received 75 mg/day and received continued therapy for the entire duration of the study with no documented grade 3 toxicities. Following dose deescalation, 2 additional patients received erlotinib 50 mg/day. None of the 5 subjects treated at the 50 mg/day dose experienced a DLT. The maximum tolerated dose (MTD) of erlotinib in combination with celecoxib 400 mg BID was therefore determined to be 50 mg/day. The observed grade 2–3 toxicities were as follows for the 12 treated patients: rash 2 of 12 (17%), mucositis 1 of 12 (17%), mouth sores 3 of 12 (25%), leukopenia 1 of 12 (8%), anemia 1 of 12 (8%), hyperglycemia 2 of 12 (17%), hypoalbuminemia 1 of 12 (8%), hypoglycemia 1 of 12 (8%), throat infection 1 of 12 (8%), urosepsis 1 of 12 (8%), and elevated creatinine 1 of 12 (8%; Table 2).

Table 2.

Toxicities including grades 1 to 4

ToxicityGrade 1Grade 2Grade 3Grade 4
Abdominal crampingΒ 2Β Β Β Β 
AlopeciaΒ 2Β Β Β Β 
AnemiaΒ 2Β 1Β Β Β 
AnxietyΒ 2Β Β Β Β 
Decreased proteinΒ 2Β Β Β Β 
LeukopeniaΒ 1Β 1Β Β Β 
DepressionΒ 3Β Β Β Β 
DiarrheaΒ 5Β Β Β Β 
Dry EyesΒ 4Β Β Β Β 
Dry SkinΒ 6Β Β Β Β 
Elevated LDHΒ 3Β Β Β Β 
Elevated serum creatinineΒ 4Β 1Β Β Β 
Elevated alkaline phosphataseΒ 3Β Β Β Β 
Elevated ALTΒ 5Β Β Β Β 
Elevated ASTΒ 4Β Β Β Β 
FatigueΒ 6Β Β Β Β 
HyperbilirubinemiaΒ 2Β Β Β Β 
HypercholesterolemiaΒ 2Β Β Β Β 
HyperglycemiaΒ 7Β 2Β Β Β 
HypoalbuminemiaΒ 3Β 1Β Β Β 
HypocalcemiaΒ 4Β Β Β Β 
HypoglycemiaΒ 1Β 1Β Β Β 
HypokalemiaΒ 2Β Β Β Β 
HyponatremiaΒ 3Β Β Β Β 
Mouth soresΒ 9Β 3Β Β Β 
MucositisΒ 3Β Β 1Β Β 
NauseaΒ 4Β Β Β Β 
NeuropathyΒ 3Β Β Β Β 
PruritisΒ 2Β Β Β Β 
RashΒ 8Β Β 2Β Β 
Shortness of breathΒ 3Β Β Β Β 
Strep throatΒ Β 1Β Β Β 
UrosepsisΒ 0Β 0Β 1Β Β 
VomitingΒ 2Β Β Β Β 
ToxicityGrade 1Grade 2Grade 3Grade 4
Abdominal crampingΒ 2Β Β Β Β 
AlopeciaΒ 2Β Β Β Β 
AnemiaΒ 2Β 1Β Β Β 
AnxietyΒ 2Β Β Β Β 
Decreased proteinΒ 2Β Β Β Β 
LeukopeniaΒ 1Β 1Β Β Β 
DepressionΒ 3Β Β Β Β 
DiarrheaΒ 5Β Β Β Β 
Dry EyesΒ 4Β Β Β Β 
Dry SkinΒ 6Β Β Β Β 
Elevated LDHΒ 3Β Β Β Β 
Elevated serum creatinineΒ 4Β 1Β Β Β 
Elevated alkaline phosphataseΒ 3Β Β Β Β 
Elevated ALTΒ 5Β Β Β Β 
Elevated ASTΒ 4Β Β Β Β 
FatigueΒ 6Β Β Β Β 
HyperbilirubinemiaΒ 2Β Β Β Β 
HypercholesterolemiaΒ 2Β Β Β Β 
HyperglycemiaΒ 7Β 2Β Β Β 
HypoalbuminemiaΒ 3Β 1Β Β Β 
HypocalcemiaΒ 4Β Β Β Β 
HypoglycemiaΒ 1Β 1Β Β Β 
HypokalemiaΒ 2Β Β Β Β 
HyponatremiaΒ 3Β Β Β Β 
Mouth soresΒ 9Β 3Β Β Β 
MucositisΒ 3Β Β 1Β Β 
NauseaΒ 4Β Β Β Β 
NeuropathyΒ 3Β Β Β Β 
PruritisΒ 2Β Β Β Β 
RashΒ 8Β Β 2Β Β 
Shortness of breathΒ 3Β Β Β Β 
Strep throatΒ Β 1Β Β Β 
UrosepsisΒ 0Β 0Β 1Β Β 
VomitingΒ 2Β Β Β Β 

Clinical outcome

Responses were evaluated for 7 subjects using the last documented histologic response at 3, 6, or 12 months. At baseline, 3 of 7 (43%) had moderate dysplasia, 3 of 7 had severe dysplasia (43%), and 1 of 7 had mild-D (14%). Three of seven achieved a complete remission (CR, 43%), 1 of 7 partial remission (PR, 14%), 1 of 7 progressive disease (PD, 14%), and 2 of 7 had stable disease(29%) with an overall Response Rate of 57% and overall clinical benefit of 86% (histologic responses to therapy are shown in Fig. 1). The median time to achieving a documented response from the time of enrollment was 5.6 months. All 3 patients who achieved a CR also had complete disappearance of their lesions by visual inspection. At the time of the last analysis, 6 of 7 patients (85%) had documented progression: 1 to stage I invasive carcinoma 6 months after completion of therapy for moderate dysplasia of the lateral tongue, 1 to stage II oral cavity carcinoma 4 months after completion of therapy, 1 to invasive squamous cell carcinoma after having stable severe dysplasia for 26 months, 2 with recurrent moderate dysplasia and severe dysplasia at 26 and 6 months, respectively, and 1 patient with recurrent high-grade dysplasia 55 months after achieving a CR for severe dysplasia of the buccal mucosa. One patient continues to be in CR 36 months after treatment of a high-grade dysplasia of the vocal cord. The median duration of follow-up for the 7 evaluable subjects was 36 months. The mean time to progression to a higher-grade dysplasia or carcinoma was 25.4 months (Fig. 2). No patients had died by the time of data analysis.

Figure 1.

Histologic responses from biopsies of 2 patients treated with celecoxib and erlotinib. A, patient 1, baseline severe dysplasia; B, patient 1, hyperplasia and no dysplasia at 3 months; C, patient 2, baseline moderate to severe dysplasia; D, patient 2, mild dysplasia at 3 months; E, patient 2, hyperplasia and no dysplasia at 6 months; F, patient 2, hyperkeratosis at 12 months.

Figure 1.

Histologic responses from biopsies of 2 patients treated with celecoxib and erlotinib. A, patient 1, baseline severe dysplasia; B, patient 1, hyperplasia and no dysplasia at 3 months; C, patient 2, baseline moderate to severe dysplasia; D, patient 2, mild dysplasia at 3 months; E, patient 2, hyperplasia and no dysplasia at 6 months; F, patient 2, hyperkeratosis at 12 months.

Close modal
Figure 2.

Time to progression to higher-grade dysplasia or carcinoma. Mean time to progression was 25.4 months. Median time for follow-up was 36 months.

Figure 2.

Time to progression to higher-grade dysplasia or carcinoma. Mean time to progression was 25.4 months. Median time for follow-up was 36 months.

Close modal

Correlative biomarker studies

Tissue analysis comparing baseline biopsies with follow-up samples after treatment revealed decreases in EGFR and p-ERK levels in the patients' last biopsies. These changes were correlated significantly with clinical responses to therapy (P = 0.019 and P = 0.006 for EGFR and p-ERK, respectively). The correlative biomarker studies were described in detail elsewhere (27).

Pharmacokinetic studies

Plasma concentrations of erlotinib and the demethyl metabolite (OSI-420) from 10 subjects administered were plotted versus time (Fig. 3A). Erlotinib concentrations only were modeled, as OSI-420 concentrations were <10% of erlotinib and curves were parallel. Four subjects administered erlotinib after the 24-hour blood draw, whereas the rest did so before their 24-hour dose. Also, 6 subjects had plasma concentrations similar to Cmax before the 48-hour sampling time. Based on the previously reported erlotinib Tmax (2–3 hours) and long elimination half-life (t1/2 > 24 hour), these undocumented timings were inputted as 3 hours before the 48-hour blood draw (31).

Figure 3.

A, plasma concentrations of erlotinib and OSI-420 (demethyl erlotinib) versus time. Subjects self-administered 50 mg (⁠⁠, n = 3), 75 mg (O, n = 6), or 100 mg erlotinib (Ξ”, n = 1) every day. Closed symbols and solid lines represent erlotinib concentrations, open symbols and dashed lines represent the OSI-420 metabolite. Plasma concentrations of erlotinib were used for pharmacokinetic modeling. B, observed versus final model predicted plasma concentrations of erlotinib. The left plot is based on the β€œtypical” prediction without including model variability (IIV), whereas the right plot is based on the complete model including IIV (post hoc fit). The linear regression coefficients (r2, observed vs. predicted) for the respective plots were 0.508 and 0.783, respectively. Inclusion of the error model (post hoc model) decreased the standard errors of prediction of the model from 22% to 15%. C, CWRES (39) were plotted versus time and post hoc predicted concentrations. The curvilinear plots depict local regression trend smoother fits to the data. The curvilinear plots in C and D depict local regression trend smoother (R Supersmoother) fits to the data. D, model validation by β€œprecision corrected visual predictive check” (40), in which observations are normalized to the median (75 mg every day) dosing scheme. One thousand simulations were performed per subject and used to compute percentile ranges (P5, P25, P50, P75, and P95) normalized to median dose and body weight. Actual observations were also normalized, and superimposed over these ranges for comparison. The shaded area indicates 90% confidence interval; solid lines indicate percentiles: 2.5, 97.5 (red); 5, 95 (blue); 25, 75 (green); 50 (black). Dashed lines indicate percentiles 5, 50, and 95 of observations.

Figure 3.

A, plasma concentrations of erlotinib and OSI-420 (demethyl erlotinib) versus time. Subjects self-administered 50 mg (⁠⁠, n = 3), 75 mg (O, n = 6), or 100 mg erlotinib (Ξ”, n = 1) every day. Closed symbols and solid lines represent erlotinib concentrations, open symbols and dashed lines represent the OSI-420 metabolite. Plasma concentrations of erlotinib were used for pharmacokinetic modeling. B, observed versus final model predicted plasma concentrations of erlotinib. The left plot is based on the β€œtypical” prediction without including model variability (IIV), whereas the right plot is based on the complete model including IIV (post hoc fit). The linear regression coefficients (r2, observed vs. predicted) for the respective plots were 0.508 and 0.783, respectively. Inclusion of the error model (post hoc model) decreased the standard errors of prediction of the model from 22% to 15%. C, CWRES (39) were plotted versus time and post hoc predicted concentrations. The curvilinear plots depict local regression trend smoother fits to the data. The curvilinear plots in C and D depict local regression trend smoother (R Supersmoother) fits to the data. D, model validation by β€œprecision corrected visual predictive check” (40), in which observations are normalized to the median (75 mg every day) dosing scheme. One thousand simulations were performed per subject and used to compute percentile ranges (P5, P25, P50, P75, and P95) normalized to median dose and body weight. Actual observations were also normalized, and superimposed over these ranges for comparison. The shaded area indicates 90% confidence interval; solid lines indicate percentiles: 2.5, 97.5 (red); 5, 95 (blue); 25, 75 (green); 50 (black). Dashed lines indicate percentiles 5, 50, and 95 of observations.

Close modal

The pharmacokinetic parameters of a basic model with no covariates (model 1), and a model relating V/F as proportional with subject body weight (model 2), are summarized in Table 3. Both estimated V/F, CL/F, and absorption rate constant (Ka) with reasonable precision. Including body weight as a covariate of V/F improved the model fit (difference in objective function = 5.09), without the need for additional fitted parameters. Unlike previous reports, inclusion of body weight as a covariate of CL/F did not improve model fit further, possibly due to the limited range of body weights in this study (29). There was >20% shrinkage for Ka in model 1 and for Ka and V/F in model 2, indicative of underestimation of interindividual parameter variability (IIV) for these parameters. Neither model allowed estimation of the intersubject variability (IIV) of CL/F. Model 2 was used for all further analyses.

Table 3.

Population pharmacokinetic parameters of erlotinib

OF (βˆ’2LL)βˆ’45.72βˆ’50.81
Parameters (RSE, %)Β Basic modelΒ kg on VΒ 
 CL/F (L/h)Β 8.34 (4.99)Β 8.29 (5.10)Β 
 V1/F (L)Β 295.2 (12.69)Β 308.35 (10.03)Β 
 Ka (h βˆ’ 1)Β 1.00 (29.55)Β 1.02 (31.09)Β 
IIVa, Οƒ2 (% shrinkage)Β 
 IIV of V/FΒ 33.02 (14.92)Β 21.19 (28.01)Β 
 IIV of KaΒ 83.87 (20.18)Β 66.09 (22.36)Β 
 σ2 (%)Β 31.97 (11.63)Β 32.7 (9.48)Β 
Number significant digitsΒ 3.8Β 4.9Β 
OF (βˆ’2LL)βˆ’45.72βˆ’50.81
Parameters (RSE, %)Β Basic modelΒ kg on VΒ 
 CL/F (L/h)Β 8.34 (4.99)Β 8.29 (5.10)Β 
 V1/F (L)Β 295.2 (12.69)Β 308.35 (10.03)Β 
 Ka (h βˆ’ 1)Β 1.00 (29.55)Β 1.02 (31.09)Β 
IIVa, Οƒ2 (% shrinkage)Β 
 IIV of V/FΒ 33.02 (14.92)Β 21.19 (28.01)Β 
 IIV of KaΒ 83.87 (20.18)Β 66.09 (22.36)Β 
 σ2 (%)Β 31.97 (11.63)Β 32.7 (9.48)Β 
Number significant digitsΒ 3.8Β 4.9Β 

NOTE: No reference intravenous dose was administered so that the fraction (F) of orally absorbed drug is unknown, as no reference intravenous dose was administered (assumed 1 for calculations). CL/F, oral clearance, V is the volume of distribution (normalized for an 81.7-kg individual in the case of the covariate model); Ka, first-order oral absorption rate constant. IIV are interindividual variance estimates, Οƒ2 = residual (intraindividual) variance and OF is the value of NONMEM objective function. Relative standard errors (RSE) of parameters were calculated as % CV.

aIIV, Οƒ2, and RSE were reported as percentage coefficient of variations (% CV), noting that the formula for % CV for a log-normal distributed parameter = √(eβˆ§Ο‰2 βˆ’ 1) Γ— 100, in which Ο‰2 is the variance of a log-normally distributed parameter. The percentage of shrinkage estimates was obtained from the NONMEM output.

There was reasonable symmetry of predicted population and post hoc concentrations versus observed data (Fig. 3B), with a slight tendency to overestimate certain lower concentrations. As expected, the individual post hoc predictions (with IIV) had greater correlations with observations (r2 = 0.508 and 0.783, respectively). Conditioned weighted residuals (CWRES) versus time and post hoc predictions (Fig. 3C) were confirmatory of an adequately fitted model. Internal model validation using a precision corrected visual predictive check plot (Fig. 3D) for model 2 resulted in 12% and 14% of observations greater than or less than the predicted 95th and 5th percentiles, respectively, which could indicate some variance shrinkage. However, the performance may be adequate given only 51 plasma concentrations. Due to data limitations, a study of the influence of celecoxib on the pharmacokinetics of erlotinib was not feasible in this study.

The combination of erlotinib and celecoxib has been shown to synergistically inhibit head and neck cancer cell growth in preclinical studies performed by our group (14, 25). Prior phase I trials have examined the clinical activity of this combination with radiation therapy in patients with recurrent HNSCC and advanced lung cancer (33). Unlike our study, the dose of erlotinib was fixed in these studies and the celecoxib dose was escalated. Grade 3 toxicities, including dermatitis and mucositis, were also reported. Our study focused on primary chemoprevention in healthy subjects with premalignancies; therefore, our threshold of defining DLTs was lower, given the chemopreventive nature of our investigation. Although both celecoxib and erlotinib have been used as single agents in chemoprevention studies for different tumors, including colorectal cancer and HNSCC (19, 34, 35), the combination has not been extensively explored in the chemoprevention setting.

The results of this phase I clinical trial support our hypothesis that EGFR-TKI and COX-2I can inhibit tumorigenesis in premalignant lesions of the head and neck. Based on our conservative approach in assessing toxicities, we determined the MTD of erlotinib in this combination to be 50 mg/day. This regimen was fairly well tolerated with erlotinib-induced rash being the most commonly observed grade 3 toxicity preventing dose escalation. We believe erlotinib at 50 mg per day would be a well-tolerated dose using this combination.

A major pitfall of our study is the small number of patients who were evaluable for response. Despite this, our correlative biomarker studies showing a downregulation in EGFR and p-ERK levels, which correlated with response to therapy (27), as well as an observed pathologic response rate of 63% with a mean duration of 20.9 months support the clinical benefit of this intervention. This will, however, require further confirmation in a larger phase II study.

There are multiple challenges and obstacles in implementing any chemopreventive intervention in a large population of healthy subjects. These include difficulty in accrual, toxicity, and quality-of-life implications in an otherwise healthy patient population, among others. As chemopreventive agents need to be delivered continuously over a rather extended period of time, it is important to select agents with a very favorable toxicity profile. Erlotinib-induced skin rash has become an accepted consequence in patients with advanced metastatic non–small cell lung cancer when compared with the dangers and side effects of cytotoxic chemotherapy. This is clearly not the case in a chemoprevention patient population such as ours. It could be argued that this is a rather high price to pay in the chemopreventive setting. In our experience, the toxicity of our intervention has been one major limitation to accrual and to dose escalation on the study. Other chemopreventive trials in head and neck cancer are incorporating erlotinib or celecoxib as a single agent or in combination with other agents. In a study using single-agent celecoxib in oral premalignant lesions, 12 of 18 biopsies showed improvement in the degree of dysplasia after 12 weeks (24). Whereas the majority of our treated subjects had severe or high-grade dysplasia, 50% of subjects in the Wirth study had mild dysplasia. Still these results raise the question of whether there is an advantage in combining the two agents versus using single agents. The results of trials using single-agent erlotinib are eagerly awaited. Our rationale for combining celecoxib and erlotinib stems from our preclinical data showing a clear synergism between EGFR and COX-2 blockade as far as inhibition of HNSCC progression and growth in vitro and in vivo (14, 25, 26).

The toxicity concern clearly opens the door to exploring better tolerated agents such as natural compounds as future chemopreventive agents, which is currently the focus of several groups including ours (36–38).

Our pharmacokinetic results indicated a larger average V/F (308.35 L) in our study compared with previous studies (approximately 220 L), which is roughly proportional to the average body weights in the two studies (approximately 69 and 82 kg, respectively). Furthermore, the average CL/F in the present study (8.29 L/h) was larger than that reported (approximately 3.95 L/h). Because CL/F of erlotinib increases less than proportionally with body weight (29), it is likely that factors in addition to body weight were responsible for the higher CL/F, e.g., the lower dose used in this study (50 to 100 mg every day) compared with 150 mg every day in the previous phase II trials. Furthermore, Thomas and colleagues reported that erlotinib CL/F was influenced by drug transporter polymorphisms (e.g., ABCB1, CYP3A5, and ABCG2; refs. 29, 31), which were not characterized in this study. The average t1/2 of erlotinib observed in this study (ln2 Γ— V/CL = 25.6 h) was shorter than in the study with patients with cancer (approximately 37.5 hours).

In summary, despite concerns about toxicity, this phase I study demonstrated clinical responses that correlated with downregulation of activated protein levels of the EGFR pathway. This is further supported by our preclinical findings (27). By the time of the last evaluation most patients had progressed to higher-grade dysplasia or invasive carcinoma; however, in some instances progression to dysplasia took more than 4 years from the time of documented CR and 1 patient continues to have a durable CR after 3 years. The preliminary efficacy signal is therefore encouraging and supports future evaluation of this regimen in a phase II efficacy trial. Our study remains the first to combine an EGFR-TKI and COX-2I with pharmacokinetic analysis for patients with premalignancies of the head and neck. The half-life of erlotinib in our study was shorter than that in other studies (approximately 37.5 hours). Larger studies are necessary to determine whether coadministration of celecoxib influences the pharmacokinetics of erlotinib.

No potential conflicts of interest were disclosed.

Conception and design: N.F. Saba, F.R. Khuri, Z.G. Chen, D.M. Shin

Development of methodology: C.S. Yang, Y. Zhao, Z.G. Chen

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): N.F. Saba, C.S. Yang, Y. Zhao, G. Sica, S. MΓΌller, R. Moreno-Williams, M. Lewis, W. Grist, A.Y. Chen, C.E. Moore, T.K. Owonikoko, S. Ramalingam, J.J. Beitler, S. Nannapaneni, H.J.C. Shin, J.R. Grandis, F.R. Khuri

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): N.F. Saba, S.J. Hurwitz, C.S. Yang, Y. Zhao, A.Y. Chen, T.K. Owonikoko, S. Ramalingam, J.J. Beitler, S. Nannapaneni, H.J.C. Shin, F.R. Khuri, Z.G. Chen, D.M. Shin

Writing, review, and/or revision of the manuscript: N.F. Saba, S.A. Kono, Y. Zhao, G. Sica, S. MΓΌller, A.Y. Chen, T.K. Owonikoko, J.J. Beitler, S. Nannapaneni, J.R. Grandis, F.R. Khuri, Z.G. Chen, D.M. Shin

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): N.F. Saba, R. Moreno-Williams

Study supervision: N.F. Saba, S. Ramalingam, F.R. Khuri, Z.G. Chen, D.M. Shin

The authors thank Genentech for providing erlotinib for this study, Dr. Anthea Hammond for her critical reading of the article, and Dana Ray for her administrative support of the study.

This study was supported by grants from the NIH (U01 CA101244 and Specialized Programs of Research Excellence in Head and Neck Cancer, P50 CA128613) to D.M. Shin and R01CA0983722, P50CA097190, and the American Cancer Society (to J.R. Grandis).

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.

1.
Rusch
β€ˆ
V
,
Klimstra
β€ˆ
D
,
Venkatraman
β€ˆ
E
,
Pisters
β€ˆ
PW
,
Langenfeld
β€ˆ
J
,
Dmitrovsky
β€ˆ
E
.β€ˆ
Overexpression of the epidermal growth factor receptor and its ligand transforming growth factor alpha is frequent in resectable non-small cell lung cancer but does not predict tumor progression
.
Clin Cancer Res
β€ˆ
1997
;
3
:
515
–
22
.
2.
Klijn
β€ˆ
JG
,
Berns
β€ˆ
PM
,
Schmitz
β€ˆ
PI
,
Foekens
β€ˆ
JA
.β€ˆ
The clinical significance of epidermal growth factor receptor (EGF-R) in human breast cancer: a review on 5232 patients
.
Endocr Rev
β€ˆ
1992
;
13
:
3
–
17
.
3.
Yao
β€ˆ
M
,
Shuin
β€ˆ
T
,
Misaki
β€ˆ
H
,
Kubota
β€ˆ
Y
.β€ˆ
Enhanced expression of c-myc and epidermal growth factor receptor (C-erbB-1) genes in primary human renal cancer
.
Cancer Res
β€ˆ
1988
;
48
:
6753
–
7
.
4.
Shin
β€ˆ
DM
,
Ro
β€ˆ
JY
,
Hong
β€ˆ
WK
,
Hittelman
β€ˆ
WN
.β€ˆ
Dysregulation of epidermal growth factor receptor expression in premalignant lesions during head and neck tumorigenesis
.
Cancer Res
β€ˆ
1994
;
54
:
3153
–
9
.
5.
Rubin Grandis
β€ˆ
J
,
Melhem
β€ˆ
MF
,
Barnes
β€ˆ
EL
,
Tweardy
β€ˆ
DJ
.β€ˆ
Quantitative immunohistochemical analysis of transforming growth factor-alpha and epidermal growth factor receptor in patients with squamous cell carcinoma of the head and neck
.
Cancer
β€ˆ
1996
;
78
:
1284
–
92
.
6.
Haddad
β€ˆ
RI
,
Shin
β€ˆ
DM
.β€ˆ
Recent advances in head and neck cancer
.
N Engl J Med
β€ˆ
2008
;
359
:
1143
–
54
.
7.
Noonberg
β€ˆ
SB
,
Benz
β€ˆ
CC
.β€ˆ
Tyrosine kinase inhibitors targeted to the epidermal growth factor receptor subfamily: role as anticancer agents
.
Drugs
β€ˆ
2000
;
59
:
753
–
67
.
8.
Saba
β€ˆ
NF
,
Khuri
β€ˆ
FR
,
Shin
β€ˆ
DM
.β€ˆ
Targeting the epidermal growth factor receptor. Trials in head and neck and lung cancer
.
Oncology (Williston Park)
β€ˆ
2006
;
20
:
153
–
61
.
9.
Wang
β€ˆ
D
,
Dubois
β€ˆ
RN
.β€ˆ
Prostaglandins and cancer
.
Gut
β€ˆ
2006
;
55
:
115
–
22
.
10.
Dannenberg
β€ˆ
AJ
,
Subbaramaiah
β€ˆ
K
.β€ˆ
Targeting cyclooxygenase-2 in human neoplasia: rationale and promise
.
Cancer Cell
β€ˆ
2003
;
4
:
431
–
6
.
11.
Dannenberg
β€ˆ
AJ
,
Lippman
β€ˆ
SM
,
Mann
β€ˆ
JR
,
Subbaramaiah
β€ˆ
K
,
DuBois
β€ˆ
RN
.β€ˆ
Cyclooxygenase-2 and epidermal growth factor receptor: pharmacologic targets for chemoprevention
.
J Clin Oncol
β€ˆ
2005
;
23
:
254
–
66
.
12.
Cohen
β€ˆ
EG
,
Almahmeed
β€ˆ
T
,
Du
β€ˆ
B
,
Golijanin
β€ˆ
D
,
Boyle
β€ˆ
JO
,
Soslow
β€ˆ
RA
, et al
β€ˆ
Microsomal prostaglandin E synthase-1 is overexpressed in head and neck squamous cell carcinoma
.
Clin Cancer Res
β€ˆ
2003
;
9
:
3425
–
30
.
13.
Chan
β€ˆ
G
,
Boyle
β€ˆ
JO
,
Yang
β€ˆ
EK
,
Zhang
β€ˆ
F
,
Sacks
β€ˆ
PG
,
Shah
β€ˆ
JP
, et al
β€ˆ
Cyclooxygenase-2 expression is up-regulated in squamous cell carcinoma of the head and neck
.
Cancer Res
β€ˆ
1999
;
59
:
991
–
4
.
14.
Chen
β€ˆ
Z
,
Zhang
β€ˆ
X
,
Li
β€ˆ
M
,
Wang
β€ˆ
Z
,
Wieand
β€ˆ
HS
,
Grandis
β€ˆ
JR
, et al
β€ˆ
Simultaneously targeting epidermal growth factor receptor tyrosine kinase and cyclooxygenase-2, an efficient approach to inhibition of squamous cell carcinoma of the head and neck
.
Clin Cancer Res
β€ˆ
2004
;
10
:
5930
–
9
.
15.
Saba
β€ˆ
NF
,
Choi
β€ˆ
M
,
Muller
β€ˆ
S
,
Shin
β€ˆ
HJ
,
Tighiouart
β€ˆ
M
,
Papadimitrakopoulou
β€ˆ
VA
, et al
β€ˆ
Role of cyclooxygenase-2 in tumor progression and survival of head and neck squamous cell carcinoma
.
Cancer Prev Res
β€ˆ
2009
;
2
:
823
–
9
.
16.
Papadimitrakopoulou
β€ˆ
VA
,
William
β€ˆ
WN
β€ˆJr.
,
Dannenberg
β€ˆ
AJ
,
Lippman
β€ˆ
SM
,
Lee
β€ˆ
JJ
,
Ondrey
β€ˆ
FG
, et al
β€ˆ
Pilot randomized phase II study of celecoxib in oral premalignant lesions
.
Clin Cancer Res
β€ˆ
2008
;
14
:
2095
–
101
.
17.
Backlund
β€ˆ
MG
,
Mann
β€ˆ
JR
,
Dubois
β€ˆ
RN
.β€ˆ
Mechanisms for the prevention of gastrointestinal cancer: the role of prostaglandin E2
.
Oncology
β€ˆ
2005
;
69
Suppl 1
:
28
–
32
.
18.
Williams
β€ˆ
CS
,
Watson
β€ˆ
AJ
,
Sheng
β€ˆ
H
,
Helou
β€ˆ
R
,
Shao
β€ˆ
J
,
DuBois
β€ˆ
RN
.β€ˆ
Celecoxib prevents tumor growth in vivo without toxicity to normal gut: lack of correlation between in vitro and in vivo models
.
Cancer Res
β€ˆ
2000
;
60
:
6045
–
51
.
19.
Steinbach
β€ˆ
G
,
Lynch
β€ˆ
PM
,
Phillips
β€ˆ
RK
,
Wallace
β€ˆ
MH
,
Hawk
β€ˆ
E
,
Gordon
β€ˆ
GB
, et al
β€ˆ
The effect of celecoxib, a cyclooxygenase-2 inhibitor, in familial adenomatous polyposis
.
N Engl J Med
β€ˆ
2000
;
342
:
1946
–
52
.
20.
Nathan
β€ˆ
CO
,
Leskov
β€ˆ
IL
,
Lin
β€ˆ
M
,
Abreo
β€ˆ
FW
,
Shi
β€ˆ
R
,
Hartman
β€ˆ
GH
, et al
β€ˆ
COX-2 expression in dysplasia of the head and neck: correlation with elF4E
.
Cancer
β€ˆ
2001
;
92
:
1888
–
95
.
21.
Ondrey
β€ˆ
FG
,
Juhn
β€ˆ
SK
,
Adams
β€ˆ
GL
.β€ˆ
Inhibition of head and neck tumor cell growth with arachidonic acid metabolism inhibition
.
Laryngoscope
β€ˆ
1996
;
106
:
129
–
34
.
22.
Scioscia
β€ˆ
KA
,
Snyderman
β€ˆ
CH
,
Rueger
β€ˆ
R
,
Reddy
β€ˆ
J
,
D'Amico
β€ˆ
F
,
Comsa
β€ˆ
S
, et al
β€ˆ
Role of arachidonic acid metabolites in tumor growth inhibition by nonsteroidal antiinflammatory drugs
.
Am J Otolaryngol
β€ˆ
1997
;
18
:
1
–
8
.
23.
Panje
β€ˆ
WR
.β€ˆ
Regression of head and neck carcinoma with a prostaglandin-synthesis inhibitor
.
Arch Otolaryngol
β€ˆ
1981
;
107
:
658
–
63
.
24.
Wirth
β€ˆ
LJ
,
Krane
β€ˆ
JF
,
Li
β€ˆ
Y
,
Othus
β€ˆ
M
,
Moran
β€ˆ
AE
,
Dorfman
β€ˆ
DM
,
et al
.β€ˆ
A pilot surrogate endpoint biomarker study of celecoxib in oral premalignant lesions
.
Cancer Prev Res
β€ˆ
2008
;
1
:
339
–
48
.
25.
Zhang
β€ˆ
X
,
Chen
β€ˆ
ZG
,
Choe
β€ˆ
MS
,
Lin
β€ˆ
Y
,
Sun
β€ˆ
SY
,
Wieand
β€ˆ
HS
, et al
β€ˆ
Tumor growth inhibition by simultaneously blocking epidermal growth factor receptor and cyclooxygenase-2 in a xenograft model
.
Clin Cancer Res
β€ˆ
2005
;
11
:
6261
–
9
.
26.
Choe
β€ˆ
MS
,
Zhang
β€ˆ
X
,
Shin
β€ˆ
HJ
,
Shin
β€ˆ
DM
,
Chen
β€ˆ
ZG
.β€ˆ
Interaction between epidermal growth factor receptor- and cyclooxygenase 2-mediated pathways and its implications for the chemoprevention of head and neck cancer
.
Mol Cancer Ther
β€ˆ
2005
;
4
:
1448
–
55
.
27.
Shin
β€ˆ
DM
,
Zhang
β€ˆ
H
,
Saba
β€ˆ
NF
,
Chen
β€ˆ
AY
,
Nannapaneni
β€ˆ
S
,
Amin
β€ˆ
AR
, et al
β€ˆ
Chemoprevention of head and neck cancer by simultaneous blocking of epidermal growth factor receptor and cyclooxygenase-2 signaling pathways: preclinical and clinical studies
.
Clin Cancer Res
β€ˆ
2013
;
19
:
1244
–
56
.
28.
Sheiner
β€ˆ
LB
,
Beal
β€ˆ
SL
.β€ˆ
Evaluation of methods for estimating population pharmacokinetic parameters. III. Monoexponential model: routine clinical pharmacokinetic data
.
J Pharmacokinet Biopharm
β€ˆ
1983
;
11
:
303
–
19
.
29.
White-Koning
β€ˆ
M
,
Civade
β€ˆ
E
,
Geoerger
β€ˆ
B
,
Thomas
β€ˆ
F
,
Le Deley
β€ˆ
MC
,
Hennebelle
β€ˆ
I
, et al
β€ˆ
Population analysis of erlotinib in adults and children reveals pharmacokinetic characteristics as the main factor explaining tolerance particularities in children
.
Clin Cancer Res
β€ˆ
2011
;
17
:
4862
–
71
.
30.
Lu
β€ˆ
JF
,
Eppler
β€ˆ
SM
,
Wolf
β€ˆ
J
,
Hamilton
β€ˆ
M
,
Rakhit
β€ˆ
A
,
Bruno
β€ˆ
R
, et al
β€ˆ
Clinical pharmacokinetics of erlotinib in patients with solid tumors and exposure-safety relationship in patients with non-small cell lung cancer
.
Clin Pharmacol Ther
β€ˆ
2006
;
80
:
136
–
45
.
31.
Thomas
β€ˆ
F
,
Rochaix
β€ˆ
P
,
White-Koning
β€ˆ
M
,
Hennebelle
β€ˆ
I
,
Sarini
β€ˆ
J
,
Benlyazid
β€ˆ
A
, et al
β€ˆ
Population pharmacokinetics of erlotinib and its pharmacokinetic/pharmacodynamic relationships in head and neck squamous cell carcinoma
.
Eur J Cancer
β€ˆ
2009
;
45
:
2316
–
23
.
32.
Beal
β€ˆ
S SL
,
Bookman
β€ˆ
A
,
Bauer
β€ˆ
RJ
.β€ˆ
NONMEM user's guides. 1989–2009
.
Ellicott City, MD
:
ICON Development Solutions
;β€ˆ
2009
.
33.
Kao
β€ˆ
J
,
Genden
β€ˆ
EM
,
Chen
β€ˆ
CT
,
Rivera
β€ˆ
M
,
Tong
β€ˆ
CC
,
Misiukiewicz
β€ˆ
K
, et al
β€ˆ
Phase 1 trial of concurrent erlotinib, celecoxib, and reirradiation for recurrent head and neck cancer
.
Cancer
β€ˆ
2011
;
117
:
3173
–
81
.
34.
Dogne
β€ˆ
JM
,
Hanson
β€ˆ
J
,
Supuran
β€ˆ
C
,
Pratico
β€ˆ
D
.β€ˆ
Coxibs and cardiovascular side-effects: from light to shadow
.
Curr Pharm Des
β€ˆ
2006
;
12
:
971
–
5
.
35.
Lin
β€ˆ
DT
,
Subbaramaiah
β€ˆ
K
,
Shah
β€ˆ
JP
,
Dannenberg
β€ˆ
AJ
,
Boyle
β€ˆ
JO
.β€ˆ
Cyclooxygenase-2: a novel molecular target for the prevention and treatment of head and neck cancer
.
Head Neck
β€ˆ
2002
;
24
:
792
–
9
.
36.
Amin
β€ˆ
AR
,
Kucuk
β€ˆ
O
,
Khuri
β€ˆ
FR
,
Shin
β€ˆ
DM
.β€ˆ
Perspectives for cancer prevention with natural compounds
.
J Clin Oncol
β€ˆ
2009
;
27
:
2712
–
25
.
37.
Kim
β€ˆ
JW
,
Amin
β€ˆ
AR
,
Shin
β€ˆ
DM
.β€ˆ
Chemoprevention of head and neck cancer with green tea polyphenols
.
Cancer Prev Res
β€ˆ
2010
;
3
:
900
–
9
.
38.
Amin
β€ˆ
AR
,
Wang
β€ˆ
D
,
Zhang
β€ˆ
H
,
Peng
β€ˆ
S
,
Shin
β€ˆ
HJ
,
Brandes
β€ˆ
JC
, et al
β€ˆ
Enhanced anti-tumor activity by the combination of the natural compounds (-)-epigallocatechin-3-gallate and luteolin: potential role of p53
.
J Biol Chem
β€ˆ
2010
;
285
:
34557
–
65
.
39.
Hooker
β€ˆ
AC
,
Staatz
β€ˆ
CE
,
Karlsson
β€ˆ
MO
.β€ˆ
Conditional weighted residuals (CWRES): a model diagnostic for the FOCE method
.
Pharm Res
β€ˆ
2007
;
24
:
2187
–
97
.
40.
Bergstrand
β€ˆ
M
,
Hooker
β€ˆ
AC
,
Wallin
β€ˆ
JE
,
Karlsson
β€ˆ
MO
.β€ˆ
Prediction-corrected visual predictive checks for diagnosing nonlinear mixed-effects models
.
AAPS J
β€ˆ
2011
;
13
:
143
–
51
.

Supplementary data