
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Cancer Therapy: Clinical |
Authors' Affiliations: 1 Department of Pharmacy and Pharmacology, The Netherlands Cancer Institute/Slotervaart Hospital; 2 Department of Medical Oncology, Antoni van Leeuwenhoek Hospital/The Netherlands Cancer Institute; 3 Department of Internal Medicine, Slotervaart Hospital, Amsterdam, the Netherlands; 4 Department of Internal Medicine, Section of Oncology, University Medical Center Utrecht; 5 Department of Biomedical Analysis, Faculty of Pharmaceutical Sciences, University of Utrecht, Utrecht, the Netherlands; 6 Department of Internal Medicine, Section of Oncology and Haematology, Medical Spectrum Twente, Enschede, the Netherlands; 7 Department of Internal Medicine, Section of Oncology, University Hospital, Maastricht, the Netherlands; and 8 Department of Medical Oncology, Radboud University Nijmegen Medical Centre, Nijmegen, the Netherlands
Requests for reprints: Tessa M. Bosch, Department of Pharmacy and Pharmacology, The Netherlands Cancer Institute/Slotervaart Hospital, Louwesweg 6, 1066 EC Amsterdam, the Netherlands. Phone: 31-20-512-4657; Fax: 31-20-512-4753; E-mail: aptbo{at}slz.nl.
| Abstract |
|---|
|
|
|---|
Experimental Design: Whole blood samples were obtained from patients with solid tumors and treated with docetaxel to quantify the exposure to docetaxel. DNA was collected to determine polymorphisms in the CYP3A and ABCB1 genes with DNA sequencing. A population pharmacokinetic analysis of docetaxel was done using nonlinear mixed-effect modeling.
Results: In total, 92 patients were assessable for pharmacokinetic analysis of docetaxel. A three-compartmental model adequately described the pharmacokinetics of docetaxel. Several polymorphisms in the CYP3A and ABCB1 genes were found, with allele frequencies of 0.54% to 48.4%. The homozygous C1236T polymorphism in the ABCB1 gene (ABCB1*8) was significantly correlated with a decreased docetaxel clearance (25%; P = 0.0039). No other relationships between polymorphisms and pharmacokinetic variables reached statistical significance. Furthermore, no relationship between haplotypes of CYP3A and ABCB1 and the pharmacokinetics could be identified.
Conclusions: The polymorphism C1236T in the ABCB1 gene was significantly related to docetaxel clearance. Our current finding may provide a meaningful tool to explain interindividual differences in docetaxel treatment in daily practice.
The metabolism of docetaxel consists of a CYP3A-mediated oxidation of the tert-butylpropionate side chain, which results in the formation of four metabolites with reduced cytotoxic activity (2). The elimination pathway is mediated by the membrane-localized, energy-dependent drug efflux ABC transporter, P-glycoprotein (ABCB1; MDR1; ref. 3). Several polymorphisms have been described in the CYP3A and ABCB1 genes. The exact functional significance of polymorphisms in the CYP3A4 gene is not yet known. CYP3A4*1B has been associated with an increased transcription in vitro (4). However, in vivo no effect of this polymorphism has been observed (5). CYP3A4*2 has been associated with a decreased nifedipine clearance (6); for CYP3A4*3, no effect of the genotype on metabolism of several substrates has been observed (7, 8); and CYP3A4*12 showed an altered enzyme activity for midazolam and testosterone (8). For CYP3A5, variant alleles have been described, such as CYP3A5*3, which causes alternative splicing and protein truncation, resulting in the absence of functional CYP3A5 from liver tissue (9). The polymorphism C3435T in the ABCB1 gene has been associated with decreased protein expression in homozygous mutant individuals (10). The C3435T polymorphism is often simultaneously found with C1236T and G2677T/A. These polymorphisms may also have an effect on the pharmacokinetics of substrates of ABCB1, but results of the in vivo relevance of ABCB1 polymorphisms have been contradictory (11).
In this study, we investigated the relationship between docetaxel pharmacokinetics and ABCB1 and CYP3A genotypes in more detail. For this, the presence of C1236T, G2677T/A, and C3435T in ABCB1; CYP3A4*1B, CYP3A4*2, CYP3A4*3, and CYP3A4*12; and CYP3A5*2 and CYP3A5*3 variant alleles was determined in cancer patients treated with docetaxel.
| Patients and Methods |
|---|
|
|
|---|
18 years; and have adequate hematopoietic, hepatic, and renal function (which was left to the discretion of the responsible oncologist). Treatment of the patient was done according to the local protocols for standard of care. Docetaxel was administered in doses of 50 to 100 mg/m2, infused over 1 hour. When at start liver function abnormalities were present, or when significant hematologic toxicity had developed during a preceding course, the dose of docetaxel was reduced according to standard clinical practice. Patients received routine supportive treatment (5'-hydroxytryptamine antagonists and dexamethasone) prophylactically. The study was approved by the Medical Ethical Committee of the Netherlands Cancer Institute and all patients provided written informed consent in accordance with institutional and governmental guidelines. All patients were treated between March 2002 and October 2005 at The Antoni van Leeuwenhoek Hospital/The Netherlands Cancer Institute, Amsterdam; UMC Utrecht, Utrecht; MST, Enschede; Academic Medical Hospital, Maastricht; UMC Nijmegen, Nijmegen; or Slotervaart Hospital, Amsterdam, all in the Netherlands.
Pharmacokinetics of docetaxel. For pharmacokinetic analysis, 5 mL of heparin blood samples were obtained before and at the end of infusion, 6 and 24 hours after start of infusion, according to the limited sampling strategy of Baille et al. (12). Blood samples were centrifuged, plasma was separated, and samples were immediately stored at 20°C until analysis. Docetaxel concentrations were determined with a validated liquid chromatography-tandem mass spectrometry assay earlier described (13). The validated range of docetaxel concentrations was 0.25 to 1,000 ng/mL and interassay accuracy and precision were between 10.2% and 1.02% and <12.8%, respectively.
Pharmacogenetic analysis. For pharmacogenetic screening, 5 mL of EDTA blood were taken before start of the docetaxel infusion. Genomic DNA was isolated according to the method of Boom et al. (14). PCR amplifications were done in 50-µL reactions with
100 ng of genomic DNA, 200 µmol/L deoxynucleotide triphosphates (Epicentre Technologies, Madison, WI), 10x PCR Buffer II (Applied Biosystems, Foster City, CA), MgCl2, 0.5 to 1 unit of AmpliTaq Gold (Applied Biosystems), and forward and reverse primers (Metabion, Planegg-Martinsried, Germany). The methods used for the amplification of the CYP genes have been described by Sata et al. (6) for CYP3A4 and by van Schaik et al. (15) for CYP3A5. Genetic polymorphisms in ABCB1 were all analyzed according to slightly modified methods previously described by Hoffmeyer et al. (10) and Kim et al. (16), respectively. PCR amplifications were done with a PTC-200 thermocycler system (MJ Research, Inc., Waltham, MA). Results of the PCR reaction were analyzed on a 2% agarose gel. DNA cycle sequencing was carried out essentially as described by the manufacturer (Applied Biosystems) in 20-µL reactions on a PTC-200 thermocycler (MJ Research). Residual dideoxy terminators were removed by ethanol/sodium acetate precipitation according to the protocol of the manufacturer (Applied Biosystems) and sequences analyzed on an Applied Biosystems 3100-Avant DNA sequencer. For sequence alignment, Seqscape v2.1 (Applied Biosystems) was used. Linkage disequilibrium between single-nucleotide polymorphisms on the same chromosome was done with the Graphical Overview of Linkage Disequilibrium (GOLD) software V1.1.0.0.9 Haplotype analysis was done with PHASE version 2.1 (17).
Population pharmacokinetic and pharmacogenetic data analysis. Population pharmacokinetic analysis was done using nonlinear mixed-effect modeling (NONMEM; double precision, version V, level 1.1; GloboMax LLC, Hanover, MD; ref. 18). PDx-Pop (version 1.1j Release 4; GloboMax LLC) was used as an interface for data and output processing and for modeling management. NONMEM uses a maximum likelihood criterion to simultaneously estimate population values of fixed effect (e.g., drug clearance or the influence of a certain polymorphism on drug clearance) and values of random effects (e.g., interindividual and residual variability).
The first-order conditional estimation method with interaction between interindividual, intraindividual, and residual variability was applied. Precision of variable estimates, as calculated by the Covariance Step of NONMEM, was evaluated.
A three-compartmental structural kinetic model previously described by Bruno et al. (19) was used as structural model. The relationship between polymorphisms in the CYP3A and ABCB1 genes and pharmacokinetic variables was tested by use of the log-likelihood test. For polymorphisms in the ABCB1 gene, the influence on docetaxel clearance and V1, V2 and V3 was determined; for CYP3A polymorphisms, only the influence on docetaxel clearance was determined. If the allele frequency of the polymorphisms was low (<15%), the assumption was made that individuals homozygous for this polymorphism had a two times larger effect than heterozygous individuals. For instance, a change in clearance (Cl) in CYP3A4*1B carriers was evaluated by use of the following equation:
![]() | (A) |
1 is the fractional change in clearance in CYP3A4*1B carriers.
If the allele frequencies were high for a polymorphism (>15%), a separate fixed effect was estimated for the different genotypes (wild-type, heterozygous, and homozygous mutant). For example, a change in docetaxel clearance in C3435T ABCB1 carriers was described as follows:
![]() | (B) |
2 and
3 are the fractional changes in clearance for heterozygous and homozygous carriers of the C3435T ABCB1 polymorphism, respectively. Furthermore, the influence of a haplotype on pharmacokinetic variables of docetaxel was determined. Statistical discrimination between hierarchical models was based on the log-likelihood ratio test using the minimal function of the objective value. P = 0.005 was considered statistically significant, which corresponds to a decrease in the minimal function of the objective value of 7.9 (degrees of freedom = 1) or 10.6 (degrees of freedom = 2).
To investigate whether the relationship between genotype and the pharmacokinetics remained significant after inclusion of known determinants of docetaxel clearance, the model previously described by Bruno et al. (19) was applied to our data set. In this model, the following relationship between clearance and covariates was found:
![]() | (C) |
1-acid glycoprotein (AAG), age (AGE), albumin (ALB), and hepatic function (HEP; i.e., aspartate aminotransferase and alkaline phosphatase, >60 IU and >300 IU, respectively) were the main predictors of docetaxel clearance. In our data set, not all covariates were available for all patients. Multiple imputation from the distribution of the covariates in the population was used for the missing data (20). Therefore, six data sets were created using this multiple imputation method and the models including the genotype relationship and without this relationship were applied.
Model validation. The bootstrap resampling technique was applied as an internal validation of the developed model. Basically, bootstrap replicates were generated by sampling randomly
65% from the original data set with replacement. The final model was fitted in the replicate data sets using the bootstrap option in the software package Wings for NONMEM (by N. Holford, version 408B, November 2004, Auckland, New Zealand; ref. 21) and variable estimates for each of the replicate data sets were obtained. The stability of the model was evaluated by visual inspection of the distribution of the model variables. Furthermore, the median variable values and 2.5 and 97.5 percentiles of the bootstrap replicates were compared with the estimated values of the original data set.
| Results |
|---|
|
|
|---|
|
|
|
|
|
|
The C1236T polymorphism in the ABCB1 gene was tested in the covariate model developed by Bruno et al. (19). Multiple imputation from the distribution of the covariates in our population was used to create values for the missing covariates, thereby developing six replicate data sets. The missing covariates ranged from 5% for BSA to 52% for AAG. The C1236T polymorphism remained statistically significant in all six replicate data sets obtained after multiple imputation (data not shown).
In this study, the following relationship between clearance and covariates was found:
![]() | (D) |
Model validation. From the original data set, 1,300 replicate data sets were generated and used for the evaluation of the stability of the final model. The mean values of the original data set were in the 2.5 to 97.5 percentiles of the bootstrap procedure (Table 2), indicating that all pharmacokinetic variables could be estimated with acceptable precision.
| Discussion |
|---|
|
|
|---|
The pharmacokinetic variables of docetaxel were well described by the previously published model of Bruno et al. (19). Overall, in the presented limited sampling strategy model, the pharmacokinetic variables of docetaxel could be estimated with acceptable precision. The mean clearance of docetaxel in this study was 54.2 L/h, which is in the range of docetaxel clearance values reported earlier (19, 23). The C1236T polymorphism is an additional determinant of docetaxel clearance, besides the known determinants of docetaxel clearance described by Bruno et al. (19). The AUC of docetaxel increased by 25% in homozygous mutant patients. Consequently, in these patients the docetaxel dose should be reduced by 25%.
The pharmacogenetic effect of C1236T on docetaxel clearance found in this study was explicit for homozygous mutant patients; heterozygous patients showed no difference in docetaxel clearance in comparison with wild-type patients. This might be explained by up-regulation of the wild-type allele in heterozygous patients. Mathijssen et al. (24) also showed that the polymorphism C1236T in the ABCB1 gene in homozygous mutant patients, and not in heterozygous mutant patients, caused a decreased clearance of irinotecan, an anticancer agent used in colorectal therapy.
Furthermore, Baker et al. (25) showed that the exposure to docetaxel was related to hematologic toxicity. Prospective identification of patients with a decreased docetaxel clearance, and thereby an increased risk of developing severe toxicity, may be important to decrease toxicity.
The homozygous synonymous C1236T polymorphism might probably have an indirect effect on the stability of mRNA (26). Shen et al. (27) suggested that allele-specific differences in RNA folding may influence downstream mRNA splicing, processing, or translational control and regulation. It is also possible that the glycine encoded by the synonymous single-nucleotide polymorphism has reduced translational activity. More complex mechanisms, such as gene-gene interactions, may also play a role. Recently, Wang et al. (28) identified the synonymous C3435T polymorphism as a main factor in allelic variation of mRNA expression of ABCB1. The wild-type allele of C3435T resulted in significantly higher mRNA expression than the mutant allele. The mutant alleles of C1236T and C3435T have different mRNA foldings, possibly causing a less efficient translation of the RNA. It was also shown that the C1236T polymorphism may not be the functional polymorphism causing allelic expression imbalance because this polymorphism did not affect mRNA expression (28). However, these results need to be confirmed in studies with larger number of samples. Evidently, more research is warranted to identify the molecular biological characteristics of C1236T that lead to altered P-glycoprotein function.
In our study, neither the C3435T nor the G2677T/A polymorphism was associated with altered pharmacokinetics of docetaxel. These polymorphisms were linked to the C1236T polymorphism, although the linkage was not very strong for the C1236T and C3435T polymorphisms. The linkage between C1236T and G2677T was strong, but the effect of the G2677A polymorphism could obscure the influence of the G2677T polymorphism on docetaxel pharmacokinetic variables. Consequently, haplotype analysis did not influence docetaxel clearance.
Other attempts to unravel the effect of polymorphisms in metabolizing enzymes and drug transporters on the pharmacokinetics of docetaxel have been undertaken. Goh et al. (29) could not detect an influence of the polymorphisms C3435T in ABCB1, CYP3A4*1B and CYP3A5*3, in an Asian population of 31 primarily nonsmall-cell lung cancer patients. It should, however, be noticed that in the study of Goh et al. (29), no patients were included who used substrates, inhibitors, or inducers for CYP3A. In our study, all patients received dexamethasone, which is a weak inducer of CYP3A.
No effect of polymorphisms in the CYP3A gene on the pharmacokinetics of docetaxel could be observed in this study. The allelic frequencies of mutations in the CYP3A gene in this population are too low to point out differences in docetaxel clearance. The clearance of docetaxel was related to the CYP3A activity measured by midazolam pharmacokinetics (29). Hirth et al. (30) and Yamamoto et al. (31) also showed a correlation between the CYP3A4 activity and docetaxel clearance using the erythromycin breath test and urinary metabolite of exogenous cortisol, respectively. Unfortunately, not all interindividual variation could be explained by measuring the CYP3A activity in patients.
It is important to screen whole genes, and not only coding regions, to identify possible relationships between polymorphisms or genes. Besides, to identify such effects, it is necessary to screen for polymorphisms in larger populations with defined ethnicities. Furthermore, details about the extent of the genomic region and allele frequencies of variants should be provided. A quantitative estimate of the prior probability of genes considered relevant should also be provided and a correction for multiple comparisons should be made due to multiple gene studies (32). First, the screening of polymorphisms and their effect on pharmacokinetics and pharmacodynamics of anticancer drugs have to be elucidated. Thereafter, high-throughput technology will find application in the simultaneous identification of polymorphisms in drug-metabolizing enzymes and drug transporters in large clinical settings.
In conclusion, the polymorphism C1236T in the ABCB1 gene was significantly related to docetaxel clearance and adds up to earlier defined determinants of docetaxel clearance. It is of interest to test in prospective trials whether dose-adaptation based on characterization of the C1236T status of ABCB1 will result in reduced interindividual variation of pharmacokinetics of docetaxel.
| 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.
9 http://www.sph.umich.edu/csg/abecasis/GOLD. ![]()
Received 12/ 2/05; revised 3/10/06; accepted 5/24/06.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
T. M. Sissung, C. E. Baum, J. Deeken, D. K. Price, J. Aragon-Ching, S. M. Steinberg, W. Dahut, A. Sparreboom, and W. D. Figg ABCB1 Genetic Variation Influences the Toxicity and Clinical Outcome of Patients with Androgen-Independent Prostate Cancer Treated with Docetaxel Clin. Cancer Res., July 15, 2008; 14(14): 4543 - 4549. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Puhalla, E. Mrozek, D. Young, S. Ottman, A. McVey, K. Kendra, N. J. Merriman, M. Knapp, T. Patel, M. E. Thompson, et al. Randomized Phase II Adjuvant Trial of Dose-Dense Docetaxel Before or After Doxorubicin Plus Cyclophosphamide in Axillary Node-Positive Breast Cancer J. Clin. Oncol., April 1, 2008; 26(10): 1691 - 1697. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Marsh, J. Paul, C. R. King, G. Gifford, H. L. McLeod, and R. Brown Pharmacogenetic Assessment of Toxicity and Outcome After Platinum Plus Taxane Chemotherapy in Ovarian Cancer: The Scottish Randomised Trial in Ovarian Cancer J. Clin. Oncol., October 10, 2007; 25(29): 4528 - 4535. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. D. Lewis, A. A. Miller, G. L. Rosner, J. E. Dowell, M. Valdivieso, M. V. Relling, M. J. Egorin, R. R. Bies, D. R. Hollis, E. G. Levine, et al. A Comparison of the Pharmacokinetics and Pharmacodynamics of Docetaxel between African-American and Caucasian Cancer Patients: CALGB 9871 Clin. Cancer Res., June 1, 2007; 13(11): 3302 - 3311. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |