Assessment of activity levels for CYP2D6*1, CYP2D6*2, and CYP2D6*41 genes by population pharmacokinetics of dextromethorphan

Semi-physiological Pharmacokinetic Model of Clozapine and Norclozapine in Healthy, Non-smoking Volunteers: The Impact of Race and Genetics

BACKGROUND AND OBJECTIVES: Clozapine is the medication of choice for treatment-resistant schizophrenia. However, it has a complex metabolism and unexplained interindividual variability. The current work aims to develop a pharmacokinetic model of clozapine and norclozapine in non-smokers and assess the impact of demographic and genetic predictors.

METHODS

Healthy volunteers were recruited in a population pharmacokinetic study. Blood samples were collected at 30 min and 1, 2, 3, 5 and 8 h following a single flat dose of clozapine (12.5 mg). The clozapine and norclozapine concentrations were measured via high-performance liquid chromatography-ultraviolet method. A semi-physiological pharmacokinetic model of clozapine and norclozapine was developed using nonlinear mixed-effects modeling. Clinical and genetic predictors were evaluated, including CYP1A2 (rs762551) and ABCB1 (rs2032582), using restriction fragment length polymorphism.

RESULTS

A total of 270 samples were collected from 33 participants. The data were best described using a two-compartment model for clozapine and a two-compartment model for norclozapine with first-order absorption and elimination and pre-systemic metabolism. The estimated (relative standard error) clearance of clozapine and norclozapine were 27 L h-1 (31.5 %) and 19.6 L h-1 (30%), respectively. Clozapine clearance was lower in sub-Saharan Africans (n = 4) and higher in Caucasians (n = 9) than Asians (n = 20). Participants with CYP1A2 (rs762551) (n = 18) and ABCB1 (rs2032582) (n = 12) mutant alleles had lower clozapine clearance in the univariate analysis.

CONCLUSIONS

This is the first study to develop a semi-physiological pharmacokinetic model of clozapine and norclozapine accounting for the pre-systemic metabolism. Asians required lower doses of clozapine as compared with Caucasians, while clozapine pharmacokinetics in sub-Saharan Africans should be further investigated in larger trials.

Circulating Biomarkers Instead of Genotyping to Establish Metabolizer Phenotypes

Pharmacogenomics (PGx) enables personalized treatment for the prediction of drug response and to avoid adverse drug reactions. Currently, PGx mainly relies on the genetic information of absorption, distribution, metabolism, and excretion (ADME) targets such as drug-metabolizing enzymes or transporters to predict differences in the patient's phenotype. However, there is evidence that the phenotype-genotype concordance is limited. Thus, we discuss different phenotyping strategies using exogenous xenobiotics (e.g., drug cocktails) or endogenous compounds for phenotype prediction. In particular, minimally invasive approaches focusing on liquid biopsies offer great potential to preemptively determine metabolic and transport capacities. Early studies indicate that ADME phenotyping using exosomes released from the liver is reliable. In addition, pharmacometric modeling and artificial intelligence improve phenotype prediction. However, further prospective studies are needed to demonstrate the clinical utility of individualized treatment based on phenotyping strategies, not only relying on genetics. The present review summarizes current knowledge and limitations.

Interethnic differences in drug response: projected impact of genetic variations in the Saudi population

Ethnicity is known to have an impact on drug responses. This is particularly important for drugs that have a narrow therapeutic window, nonlinearity in pharmacokinetics and are metabolized by enzymes that demonstrate genetic polymorphisms. However, most clinical trials are conducted among Caucasians, which might limit the usefulness of the findings of such studies for other ethnicities. The representation of participants from Saudi Arabia in global clinical trials is low. Therefore, there is a paucity of evidence to assess the impact of ethnic variability in the Saudi population on drug response. In this article, the authors assess the projected impact of genetic polymorphisms in drug-metabolizing enzymes and drug targets on drug response in the Saudi population.

A frequent CYP2D6 variant promotes skipping of exon 3 and reduces CYP2D6 protein expression in human liver samples

CYP2D6 is one of the most polymorphic drug-metabolizing enzymes in the liver. While genetic CYP2D6 variants serve as clinical biomarkers to predict CYP2D6 activity, large inter-person variability in CYP2D6 expression remains unaccounted for. Previous results suggest that there is variable expression of a CYP2D6 splice isoform with an in-frame deletion of exon 3 (CYP2D6ΔE3) encoding a protein lacking numerous active site residues. Here, using fragment analysis and RT-qPCR, we revealed that rs1058164 G (MAF = 27%-43%) is associated with increased formation of CYP2D6∆E3 in human liver samples (1.4-2.5-fold) and transfected cells. Furthermore, western blots showed that rs1058164 G was associated with a 50% decrease in full-length hepatic CYP2D6 protein expression. In addition, by studying a larger liver cohort, we confirmed our previous results that rs16947 (CYP2D6*2) reduces full-length CYP2D6 mRNA by increasing the production of an unstable splice isoform lacking exon 6 (CYP2D6ΔE6) and that the impact of CYP2D6ΔE6 is offset in carriers of the downstream enhancer variant rs5758550. The three frequent SNPs (rs1058164, rs16947, and rs5758550) form various 3-SNP-haplotypes, each with distinct CYP2D6 expression characteristics. Using an expression score (ES) system, we tested the impact of the 3-SNP-haplotype on improving the standard model to predict hepatic CYP2D6 protein expression based on genotype. A model that incorporates the 3-SNP-haplotype provided the best fit for CYP2D6 expression and also accounted for more variability in CYP2D6 protein levels (59%) than a model based on the accepted standard (36%) or one that only adds rs16947 and rs5758550 (42%). Clinical studies are needed to determine whether including the 3-SNP-haplotype alongside current standard CYP2D6 models improves the predictive value of CYP2D6 panels.

Determination of CYP450 activities in diabetes mellitus rats by a UHPLC-MS/MS method

In this study, we investigated the effect of type 1 diabetes mellitus on the modulation of the activities of CYP450s in dynamics by a UHPLC-MS/MS method. The diabetic rat model was constructed by an intraperitoneal single injection of streptozotocin. Fasting blood glucose levels > 16.7 mmol/L were considered as diabetic. The rats were given a cocktail of four probe drugs (10 mg/kg phenacetin, 1 mg/kg tolbutamide, 10 mg/kg metoprolol, and 10 mg/kg midazolam) by oral administration for the pharmacokinetic study. Thereafter, the metabolic ratio (MR) of the metabolites to probe substrates were determined. The results indicated that two weeks after diabetes was induced, diabetes increased the MRs of acetaminophen/phenacetin (CYP1A2) and 4-hydroxyl tolbutamide/tolbutamide (CYP2C9); however, it decreased the MRs of α-hydroxy metoprolol/metoprolol (CYP2D6) and 1-hydroxy midazolam/midazolam (CYP3A4). Two months after diabetes was induced, diabetes increased the MRs of acetaminophen/phenacetin and 4-hydroxyl tolbutamide/tolbutamide. The MR of α-hydroxy metoprolol/metoprolol was decreased and the MR of 1-hydroxy midazolam/midazolam was increased but the difference was not significant. According to the results, CYP1A2 and CYP2C9 activities were enhanced in the diabetic rats. and CYP2D6 activity was inhibited in a short period of diabetes; however, the decrease in CYP2D6 activity was not significant in the long period. CYP3A4 activity was decreased in a short period of diabetes and increased in a long period of diabetes but was not significant in the two periods. This study suggests the activity change rule of the CYP450 enzyme system in diabetes mellitus, which can provide a reference for precise clinical medication.

Estimating the In Vivo Function of CYP2D6 Alleles through Population Pharmacokinetic Modeling of Brexpiprazole

Accurate prediction of CYP2D6 phenotype from genotype information is important to support safe and efficacious pharmacotherapy with CYP2D6 substrates. To facilitate accurate CYP2D6 genotype-phenotype translation, there remains a need to investigate the enzyme activity associated with individual CYP2D6 alleles using large clinical data sets. This study aimed to quantify and compare the in vivo function of different CYP2D6 alleles through population pharmacokinetic (PopPK) modeling of brexpiprazole using data from 13 clinical studies. A PopPK model of brexpiprazole and its two metabolites, DM-3411 and DM-3412, was developed based on plasma concentration samples from 826 individuals. As the minor metabolite, DM-3412, is formed via CYP2D6, the metabolic ratio of DM-3412:brexpiprazole calculated from the PopPK parameter estimates was used as a surrogate measure of CYP2D6 activity. A CYP2D6 genotype-phenotype analysis based on 496 subjects showed that the CYP2D6*2 allele (n = 183) was associated with only 10% enzyme activity relative to the wild-type allele (CYP2D6*1) and a low enzyme activity was consistently observed across genotypes containing CYP2D6*2. Among the decreased function alleles, the following enzyme activities relative to CYP2D6*1 were estimated: 23% for CYP2D6*9 (n = 20), 32% for CYP2D6*10 (n = 62), 64% for CYP2D6*14 (n = 1), 4% for CYP2D6*17 (n = 37), 4% for CYP2D6*29 (n = 13), and 9% for CYP2D6*41 (n = 64). These findings imply that a lower functional value would more accurately reflect the in vivo function of many reduced function CYP2D6 alleles in the metabolism of brexpiprazole. The low enzyme activity observed for CYP2D6*2, which has also been reported by others, suggests that the allele exhibits substrate-specific enzyme activity.

Physiologically based pharmacokinetic (PBPK) modeling of the role of CYP2D6 polymorphism for metabolic phenotyping with dextromethorphan

The cytochrome P450 2D6 (CYP2D6) is a key xenobiotic-metabolizing enzyme involved in the clearance of many drugs. Genetic polymorphisms in CYP2D6 contribute to the large inter-individual variability in drug metabolism and could affect metabolic phenotyping of CYP2D6 probe substances such as dextromethorphan (DXM). To study this question, we (i) established an extensive pharmacokinetics dataset for DXM; and (ii) developed and validated a physiologically based pharmacokinetic (PBPK) model of DXM and its metabolites dextrorphan (DXO) and dextrorphan O-glucuronide (DXO-Glu) based on the data. Drug-gene interactions (DGI) were introduced by accounting for changes in CYP2D6 enzyme kinetics depending on activity score (AS), which in combination with AS for individual polymorphisms allowed us to model CYP2D6 gene variants. Variability in CYP3A4 and CYP2D6 activity was modeled based on in vitro data from human liver microsomes. Model predictions are in very good agreement with pharmacokinetics data for CYP2D6 polymorphisms, CYP2D6 activity as described by the AS system, and CYP2D6 metabolic phenotypes (UM, EM, IM, PM). The model was applied to investigate the genotype-phenotype association and the role of CYP2D6 polymorphisms for metabolic phenotyping using the urinary cumulative metabolic ratio (UCMR), DXM/(DXO + DXO-Glu). The effect of parameters on UCMR was studied via sensitivity analysis. Model predictions indicate very good robustness against the intervention protocol (i.e. application form, dosing amount, dissolution rate, and sampling time) and good robustness against physiological variation. The model is capable of estimating the UCMR dispersion within and across populations depending on activity scores. Moreover, the distribution of UCMR and the risk of genotype-phenotype mismatch could be estimated for populations with known CYP2D6 genotype frequencies. The model can be applied for individual prediction of UCMR and metabolic phenotype based on CYP2D6 genotype. Both, model and database are freely available for reuse.

The regioselectivity of the interaction between dextromethorphan and CYP2D6

CYP2D6 is an important enzyme of the cytochrome P450 superfamily, and catalyzes nearly 25% of the drugs sold in the market. For decades, the interactions and metabolism between CYP2D6 and substrates have been a hot topic. However, the key factors of the catalytic regioselectivity for CYP2D6 still remain controversial. Here, we construct four systems to explore the interaction between dextromethorphan (DM) and CYP2D6. A new binding mode of CYP2D6 is defined, and two key residues (residue Asp301 and residue Glu216) are discovered working simultaneously to stabilize the DM at the reactive site by forming water bridge hydrogen bonds when CYP2D6 binds DM. Our results also indicate that the substrate concentration could mediate the binding mode between the substrate and CYP2D6 by decreasing the volume of the catalytic pocket, which is not conducive to the O-demethylation of DM but benefits the N-demethylation of DM. These results could shed light on the process of CYP2D6 binding to the substrate, and help to better understand the regioselectivity of CYP2D6 catalyzing the substrates.

CYP2D6 Reduced Function Variants and Genotype/Phenotype Translations of CYP2D6 Intermediate Metabolizers: Implications for Personalized Drug Dosing in Psychiatry

Genetic differences in cytochrome P450 (CYP)-mediated metabolism have been known for several decades. The clinically most important polymorphic CYP enzyme is CYP2D6, which plays a key role in the metabolism of many antidepressants and antipsychotics, along with a range of non-psychiatric medications. Dose individualization based on CYP2D6 genotype to improve the effect and safety of drug treatment has been an ambition for a long time. Clinical use of CYP2D6 genotyping is steadily increasing; however, for pre-emptive genotyping to be successful in predicting individual dose requirements, high precision of genotype-to-phenotype translations are required. Recently, guidelines for assigning CYP2D6 enzyme activity scores of CYP2D6 variant alleles, and subsequent diplotype-to-phenotype translations, were published by the Clinical Pharmacogenetics Implementation Consortium (CPIC) and the Dutch Pharmacogenetics Working Group. Consensus on assigning activity scores of CYP2D6 variant alleles and translating diplotype scores into CYP2D6 poor, intermediate, normal, or ultrarapid metabolizer groups were obtained by consulting 37 international experts. While assigning enzyme activities of non-functional (score 0) and fully functional (score 1) alleles are straightforward, reduced function variant alleles are more complex. In this article, we present data showing that the assigned activity scores of reduced function variant alleles in current guidelines are not of sufficient precision; especially not for CYP2D6*41, where the guideline activity score is 0.5 compared to 0.05–0.15 in pharmacogenetic studies. Due to these discrepancies, CYP2D6 genotypes with similar guidelinediplotype scores exhibit substantial differences in CYP2D6 metabolizer phenotypes. Thus, it is important that the guidelines are updated to be valid in predicting individual dose requirements of psychiatric drugs and others metabolized by CYP2D6.

Quantification of In Vivo Metabolic Activity of CYP2D6 Genotypes and Alleles Through Population Pharmacokinetic Analysis of Vortioxetine

Assignment of CYP2D6 phenotype from genotype data can be challenging and despite efforts to standardize translation, there is currently no universally accepted method. To facilitate standardization, there remains a need to precisely quantify the in vivo function of different CYP2D6 genotypes. Vortioxetine is metabolized to its major metabolite, Lu AA34443, primarily via CYP2D6. The aim of this study was to quantify the in vivo CYP2D6 activity of different CYP2D6 alleles and genotypes through population pharmacokinetic (PopPK) modeling of vortioxetine and Lu AA34443. Plasma concentration data of vortioxetine and Lu AA34443 from 1,140 subjects originating from 29 clinical pharmacology studies were pooled for the analysis. A joint PopPK model described the pharmacokinetics of vortioxetine and Lu AA34443 simultaneously and provided estimates of the CYP2D6-mediated metabolism for each subject. Subjects normally classified as CYP2D6 intermediate metabolizers (IMs) showed different levels of CYP2D6 activity with carriers of one fully functional allele and one null function allele having 77% higher CYP2D6 activity compared with carriers of two decreased function alleles (P < 0.0001). The decreased function alleles were associated with different levels of reduction of CYP2D6 activity. Fixing the activity of fully functional alleles to 1.0, the relative activities of CYP2D6*9, CYP2D6*10, CYP2D6*17, and CYP2D6*41 were 0.22, 0.37, 0.17, and 0.21, respectively. The activity of CYP2D6*10 was shown to be significantly greater than that of CYP2D6*17 (P = 0.01) and CYP2D6*41 (P = 0.02). These results warrant further discussion of current CYP2D6 genotype-phenotype classification systems particularly regarding decreased function alleles and the IM phenotype.

Tri-Allelic Haplotypes Determine and Differentiate Functionally Normal Allele CYP2D6*2 and Impaired Allele CYP2D6*41

CYP2D6 metabolizes 20-25% of all clinically used drugs and its complex genetic polymorphism is a major determinant of drug safety and efficacy. We investigated the basis for the functional difference between the two common alleles *2 (g.2851C>T + g.4181G>C, normal function) and *41 (additional intronic g.2989G>A, reduced function). A recently reported far-distant enhancer polymorphism rs5758550A/G linked to *2 has been suggested to play a decisive role. Genotyping of two white cohorts confirmed strong linkage of rs5758550G to *2, whereas no influence was found on metabolic ratio of sparteine or hepatic expression. Genomic plasmid constructs carrying individual variants or combinations thereof were expressed in COS1 and Huh7 cells. Both g.2851C>T(R296C) and g.2989G>A reduced enzyme activity and protein levels similarly by ~ 50-65% compared to reference (*1), whereas the double variant had only ~ 20% activity. Although the unexpected loss of function caused by g.2851C>T was compensated by g.4181G>C (mimicking the EM-phenotype of *2), the additional loss of function due to intronic g.2989G>A in the triple variant was not compensated (mimicking the IM-phenotype of *41). We also confirmed increased erroneous splicing in carriers of g.2989G>A but not of g.2851C>T as a likely explanation for the impaired function of *41. In conclusion, our data demonstrate g.2989G>A as causal variant of impaired allele CYP2D6*41 whereas triple-haplotypes have to be considered to explain the functional difference between *2 and *41. These data are important for genotyping strategies and clinical implementation of CYP2D6 pharmacogenetics.

Long-Distance Phasing of a Tentative "Enhancer" Single-Nucleotide Polymorphism With CYP2D6 Star Allele Definitions

BACKGROUND

The CYP2D6 gene locus has been extensively studied over decades, yet a portion of variability in CYP2D6 activity cannot be explained by known sequence variations within the gene, copy number variation, or structural rearrangements. It was proposed that rs5758550, located 116 kb downstream of the CYP2D6 gene locus, increases gene expression and thus contributes to variability in CYP2D6 activity. This finding has, however, not been validated. The purpose of the study was to address a major technological barrier, i.e., experimentally linking rs5758550, also referred to as the "enhancer" single-nucleotide polymorphism (SNP), to CYP2D6 haplotypes >100 kb away. To overcome this challenge is essential to ultimately determine the contribution of the "enhancer" SNP to interindividual variability in CYP2D6 activity.

METHODS

A large ethnically mixed population sample (n=3,162) was computationally phased to determine linkage between the "enhancer" SNP and CYP2D6 haplotypes (or star alleles). To experimentally validate predicted linkages, DropPhase2D6, a digital droplet PCR (ddPCR)-based method was developed. 10X Genomics Linked-Reads were utilized as a proof of concept.

RESULTS

Phasing predicted that the "enhancer" SNP can occur on numerous CYP2D6 haplotypes including CYP2D6*1, *2, *5, and *41 and suggested that linkage is incomplete, i.e., a portion of these alleles do not have the "enhancer" SNP. Phasing also revealed differences among the European and African ancestry data sets regarding the proportion of alleles with and without the "enhancer" SNP. DropPhase2D6 was utilized to confirm or refute the predicted "enhancer" SNP location for individual samples, e.g., of n=3 samples genotyped as *1/*41, rs5758550 was on the *41 allele of two samples and on the *1 allele of one sample. Our findings highlight that the location of the "enhancer" SNP must not be assigned by "default." Furthermore, linkage between the "enhancer" SNP and CYP2D6 star allele haplotypes was confirmed with 10X Genomics technology.

CONCLUSIONS

Since the "enhancer" SNP can be present on a portion of normal, decreased, or no function alleles, the phase of the "enhancer" SNP must be considered when investigating the impact of the "enhancer" SNP on CYP2D6 activity.

Simultaneous Semi-Mechanistic Population Pharmacokinetic Modeling Analysis of Enalapril and Enalaprilat Serum and Urine Concentrations From Child Appropriate Orodispersible Minitablets

Enalapril is recommended as the first line of therapy and is proven to improve survival rates for treatment of Pediatric Heart Failure; however, an approved drug and child appropriate dosage formulation is still absent. The present analysis was conducted to perform a detailed model informed population pharmacokinetic analysis of prodrug enalapril and its active metabolite enalaprilat in serum and urine. Further, a model informed dosage form population-pharmacokinetic analysis was conducted to evaluate differences in pharmacokinetics of enalapril and its active metabolite enalaprilat when prodrug was administered to 24 healthy adults in a crossover, two periods, two treatments, phase I clinical trial using child-appropriate orodispersible mini-tablets (ODMT) and reference (Renitec®) dosage formulation. A simultaneous semi-mechanistic population-pharmacokinetic model was developed using NONMEM software, which predicted full profile serum and urine concentrations of enalapril and enalaprilat. First-order conditional estimation with interaction was used for parameter estimation. Transit compartments added using Erlang distribution method to predicted enalapril absorption and enalaprilat formation phases. Normalized body weight was identified as covariate related to enalapril volume of distribution. Visual predictive check (VPC) plots and conducted bootstrap analysis validated the model. The data from the two formulations were pooled for population-pharmacokinetic analysis and covariate effect of the formulation was found on mean transit time (MTT1) of enalapril absorption. In addition, data of each formulation were modeled separately and the estimated parameters of each individual administered both formulations were correlated using paired samples Wilcoxon rank test (p < 0.05 = significant) which also showed only a significant difference (p = 0.03) in MTT1 i.e., 5 min early appearance of enalapril from ODMT compared to reference tablets. No difference in the pharmacokinetics of active enalaprilat was found from the ODMT compared to the reference formulation. The population pharmacokinetic analysis provided detailed information about the pharmacokinetics of enalapril and enalaprilat, which showed that the ODMT formulation might have similar pharmacodynamic response compared to the reference formulation.

Significantly lower CYP2D6 metabolism measured as the O/N-desmethylvenlafaxine metabolic ratio in carriers of CYP2D6*41 versus CYP2D6*9 or CYP2D6*10: a study on therapeutic drug monitoring data from 1003 genotyped Scandinavian patients

AIMS

CYP2D6*9, CYP2D6*10 and CYP2D6*41 are the most frequent reduced-function CYP2D6 alleles in Caucasians. Despite lacking in vivo evidence, they are collectively classified with an enzyme activity score of 0.5. Thus, the aim of this study was to compare the functional impact of CYP2D6*9, CYP2D6*10 and CYP2D6*41 on CYP2D6 metabolism in a large patient population.

METHODS

A total of 1003 patients (mainly Caucasians) with data on CYP2D6 genotype and serum concentrations of venlafaxine and metabolites were included from a therapeutic drug monitoring service in Oslo, Norway. The O-desmethyl-to-N-desmethyl-venlafaxine metabolic ratio (MR) was applied as CYP2D6 biomarker and compared (Mann-Whitney) between carriers of CYP2D6*9-10 (merged) and CYP2D6*41, either combined with CYP2D6*1 or non-coding (null) alleles. MR subgroup estimates were obtained by multiple linear regression for calculations of CYP2D6*9-10 and CYP2D6*41 activity scores.

RESULTS

MR was significantly lower in carriers of CYP2D6*41 than CYP2D6*9-10 (P < 0.002). The majority of CYP2D6*41/null carriers (86.7%) had MR in the observed range of CYP2D6null/null carriers compared with the minority of CYP2D6*9-10/null carriers (17.4%). CYP2D6 genotype explained 60.7% of MR variability in the multivariate analysis providing subgroup estimates of 9.54 (95% CI; 7.45-12.20), 3.55 (2.06-6.10), 1.33 (0.87-2.05) and 0.47 (0.35-0.61) in carriers of CYP2D6*1/null (n = 269), CYP2D6*9-10/null (n = 17), CYP2D6*41/null (n = 30) and CYP2D6null/null (n = 95), respectively. Based on these estimates, the calculated activity score of CYP2D6*41 was 0.095 compared to 0.34 for CYP2D6*9-10.

CONCLUSIONS

CYP2D6 metabolism measured as the O/N-desmethylvenlafaxine ratio is significantly lower in Scandinavian carriers of CYP2D6*41 vs. CYP2D6*9-10. Thus, these alleles should be differentiated when classifying CYP2D6 phenotype from genotype.

Microsecond MD simulations of human CYP2D6 wild-type and five allelic variants reveal mechanistic insights on the function

Characterization of cytochrome P450 2D6 (CYP2D6) and the impact of the major identified allelic variants on the activity of one of the most dominating drug-metabolising enzymes is essential to increase drug safety and avoid adverse reactions. Microsecond molecular dynamics simulations have been performed to capture the dynamic signatures of this complex enzyme and five allelic variants with diverse enzymatic activity. In addition to the apo simulations, three substrates (bufuralol, veliparib and tamoxifen) and two inhibitors (prinomastat and quinidine) were included to explore their influence on the structure and dynamical features of the enzyme. Our results indicate that the altered enzyme activity can be attributed to changes in the hydrogen bonding network within the active site, and local structural differences in flexibility, position and shape of the binding pocket. In particular, the increased (CYP2D6*53) or the decreased (CYP2D6*17) activity seems to be related to a change in dynamics of mainly the BC loop due to a modified hydrogen bonding network around this region. In addition, the smallest active site volume was found for CYP2D6*4 (no activity). CYP2D6*2 (normal activity) showed no major differences in dynamic behaviour compared to the wild-type.

Association of Impaired Cytochrome P450 2D6 Activity Genotype and Phenotype With Therapeutic Efficacy of Primaquine Treatment for Latent Plasmodium vivax Malaria

IMPORTANCE

Latent hepatic Plasmodium vivax hypnozoites provoke repeated clinical attacks called relapses. Only primaquine phosphate kills hypnozoites, and its therapeutic activity may depend on naturally polymorphic cytochrome P450 2D6 isotype (CYP2D6) activity.

OBJECTIVE

To examine the association of impaired CYP2D6 genotypes and CYP2D6 metabolic phenotypes with therapeutic failure of directly observed high-dose primaquine treatment for P vivax malaria relapse.

DESIGN, SETTING, AND PARTICIPANTS

Nested case-control study of patients who, in July 2014, completed a randomized clinical trial of directly observed primaquine treatment for radical cure of acute P vivax malaria in an area of Indonesia where reinfection during 1 year of posttreatment follow-up was improbable. A total of 177 of 180 patients with P vivax malaria completed the clinical trial of primaquine treatment to prevent relapse; 151 were eligible for recruitment as controls. After screening, 59 potential control individuals (no relapse) and 26 potential case patients (relapse) were considered, and 36 controls and 21 cases were enrolled.

EXPOSURES

Cases and controls were exposed to P vivax malaria and primaquine therapy but had variable exposure to the enzymatic activity of CYP2D6, classified as impaired by a genotype-determined qualitative phenotype (poor or intermediate), genotype-determined activity score less than 1.5, or a log of the 24-hour pooled urine dextromethorphan-dextrorphan metabolic ratio greater than -1.0.

MAIN OUTCOMES AND MEASURES

Unadjusted odds ratios (ORs) of relapse with impaired CYP2D6 metabolism determined by genotype or measured by urinary dextromethorphan-dextrorphan metabolic ratio.

RESULTS

Among the 21 cases (mean [SD] age, 30.5 [6.3] years; all male) and 36 controls (mean [SD] age, 29.0 [3.6] years; all male), 6 CYP2D6 alleles (*1, *2, *4, *5, *10, and *41) occurred as 12 distinct genotypes, with model activity scores ranging from 0.0 to 2.0. Among 32 patients with genotypic activity scores of 1.0 or less, 18 had experienced relapse, whereas among the 25 with scores higher than 1.0, 3 had experienced relapse (OR, 9.4; 95% CI, 2.1-57.0; P = .001). When the log of the metabolic ratio of dextromethorphan-dextrorphan was -1.0 or less, only 1 of 18 patients experienced relapse, whereas above that threshold (consistent with low metabolic activity), 20 of 39 patients experienced relapse (OR, 18; 95% CI, 2.2-148.0; P = .007).

CONCLUSIONS AND RELEVANCE

Genotype-determined and directly measured impaired levels of CYP2D6 activity were associated with elevated risk of therapeutic failure. These findings suggest a natural variability in CYP2D6-dependent metabolism of primaquine as a key determinant of therapeutic efficacy against latent P vivax malaria.

Ten Years' Experience with the CYP2D6 Activity Score: A Perspective on Future Investigations to Improve Clinical Predictions for Precision Therapeutics

The seminal paper on the CYP2D6 Activity Score (AS) was first published ten years ago and, since its introduction in 2008, it has been widely accepted in the field of pharmacogenetics. This scoring system facilitates the translation of highly complex CYP2D6 diplotype data into a patient’s phenotype to guide drug therapy and is at the core of all CYP2D6 gene/drug pair guidelines issued by the Clinical Pharmacogenetics Implementation Consortium (CPIC). The AS, however, only explains a portion of the variability observed among individuals and ethnicities. In this review, we provide an overview of sources in addition to CYP2D6 genotype that contribute to the variability in CYP2D6-mediated drug metabolism and discuss other factors, genetic and non-genetic, that likely contribute to the observed variability in CYP2D6 enzymatic activity.

Pharmacokinetics and Pharmacodynamics of Meloxicam in East Asian Populations: The Role of Ethnicity on Drug Response

We aimed to reanalyze the differences in the pharmacokinetics (PKs) of meloxicam in East Asian populations based on a population approach using previously published data and to investigate the factors found in population PK analysis that affect the pharmacodynamics (PDs) of meloxicam. Population PK analysis was performed in 119 healthy male subjects (30 Japanese, 30 Chinese, 29 Korean, and 30 white) under strictly controlled trial conditions with regulated meals and a single lot of the drug. We found that CYP2C9 genotype and lean body mass were statistically significant predictors of clearance and volume of distribution, respectively. A statistical significant difference in the PK parameters between ethnic groups could not be identified. Simulations using PK/PD models showed that CYP2C9 genotype is the factor that affects the PDs of meloxicam. The genetic polymorphisms highlighted in this study would be beneficial for conducting clinical trials in East Asians with similar genetic backgrounds.

CYP2D6 phenotypes are associated with adverse outcomes related to opioid medications

BACKGROUND

Variation in the CYP2D6 gene may affect response to opioids in both poor and ultrarapid metabolizers, but data demonstrating such associations have been mixed, and the impact of variants on toxicity-related symptoms (e.g., nausea) is unclear. Therefore, we examined the association between CYP2D6 phenotype and poor pain control or other adverse symptoms related to the use of opioids in a sample of primary care patients.

MATERIALS AND METHODS

We identified all patients in the Mayo Clinic RIGHT Protocol who were prescribed an opioid medication between July 01, 2013 and June 30, 2015, and categorized patients into three phenotypes: poor, intermediate to extensive, or ultrarapid CYP2D6 metabolizers. We reviewed the electronic health record of these patients for indications of poor pain control or adverse symptoms related to medication use. Associations between phenotype and outcomes were assessed using Chi-square tests and logistic regression.

RESULTS

Overall, 257 (25% of RIGHT Protocol participants) patients received at least one opioid prescription; of these, 40 (15%) were poor metabolizers, 146 (57%) were intermediate to extensive metabolizers, and 71 (28%) were ultrarapid metabolizers. We removed patients that were prescribed a CYP2D6 inhibitor medication (n=38). After adjusting for age and sex, patients with a poor or ultrarapid phenotype were 2.7 times more likely to experience either poor pain control or an adverse symptom related to the prescription compared to patients with an intermediate to extensive phenotype (odds ratio: 2.68; 95% CI: 1.39, 5.17; p=0.003).

CONCLUSION

Our results suggest that >30% of patients with a poor or ultrarapid CYP2D6 phenotype may experience an adverse outcome after being prescribed codeine, tramadol, oxycodone, or hydrocodone. These medications are frequently prescribed for pain relief, and ~39% of the US population is expected to carry one of these phenotypes, suggesting that the population-level impact of these gene-drug interactions could be substantial.

CYP2D6 phenotypes are associated with adverse outcomes related to opioid medications

Background

Variation in the CYP2D6 gene may affect response to opioids in both poor and ultrarapid metabolizers, but data demonstrating such associations have been mixed, and the impact of variants on toxicity-related symptoms (e.g., nausea) is unclear. Therefore, we examined the association between CYP2D6 phenotype and poor pain control or other adverse symptoms related to the use of opioids in a sample of primary care patients.

Materials and methods

We identified all patients in the Mayo Clinic RIGHT Protocol who were prescribed an opioid medication between July 01, 2013 and June 30, 2015, and categorized patients into three phenotypes: poor, intermediate to extensive, or ultrarapid CYP2D6 metabolizers. We reviewed the electronic health record of these patients for indications of poor pain control or adverse symptoms related to medication use. Associations between phenotype and outcomes were assessed using Chi-square tests and logistic regression.

Results

Overall, 257 (25% of RIGHT Protocol participants) patients received at least one opioid prescription; of these, 40 (15%) were poor metabolizers, 146 (57%) were intermediate to extensive metabolizers, and 71 (28%) were ultrarapid metabolizers. We removed patients that were prescribed a CYP2D6 inhibitor medication (n=38). After adjusting for age and sex, patients with a poor or ultrarapid phenotype were 2.7 times more likely to experience either poor pain control or an adverse symptom related to the prescription compared to patients with an intermediate to extensive phenotype (odds ratio: 2.68; 95% CI: 1.39, 5.17; p=0.003).

Conclusion

Our results suggest that >30% of patients with a poor or ultrarapid CYP2D6 phenotype may experience an adverse outcome after being prescribed codeine, tramadol, oxycodone, or hydrocodone. These medications are frequently prescribed for pain relief, and ~39% of the US population is expected to carry one of these phenotypes, suggesting that the population-level impact of these gene–drug interactions could be substantial.

Application of a physiologically based pharmacokinetic model for the evaluation of single-point plasma phenotyping method of CYP2D6

PURPOSE

Determining metabolic ratio from single-point plasma is potentially a good phenotyping method of CYP2D6 to reduce the required time interval and increase the reliability of data. It is difficult to conduct large sample size clinical trials to evaluate this phenotyping method for multiple plasma points. A physiologically based pharmacokinetic (PBPK) model can be developed to do simulations based on the large virtual Chinese population and evaluate single-point plasma phenotyping method of CYP2D6.

METHODS

Pharmacokinetic data of dextromethorphan (DM) and its metabolite dextrorphan (DX) after oral administration were used for model development. The SimCYP® model incorporating Chinese demographic, physiological, and enzyme data was used to simulate DM and DX pharmacokinetics in different phenotype groups.

RESULTS

The ratios of the simulated to the observed mean AUC and Cmax of DM were 1.01 and 0.81 for extensive metabolizers (EMs), 0.90 and 0.81 for intermediate metabolizers (IMs), and 1.12 and 0.84 for poor metabolizers (PMs). The ratios of the simulated to the observed mean AUC and Cmax of DX were 1.12 and 0.89 for EMs, 0.66 and 0.62 for IMs. All ratios were within the predefined criterion of 0.5-2. The simulations of DM and DX pharmacokinetic profiles in 1000 virtual Chinese subjects with reported frequencies of different phenotypes indicated that statistically significant correlations were found between metabolic ratio of DM to DX (MRDM/DX) from AUC and from single-point plasma from 1 to 30h (all p-values <0.001).

CONCLUSION

MRDM/DX from single-point plasma from 1 to 30h after the administration of 30mg controlled-release DM could predict the MRDM/DX from AUC well and could be used as the phenotyping method of CYP2D6 for EMs, IMs, and PMs.

Influence of the cytochrome P450 2D6 *10/*10 genotype on the pharmacokinetics of paroxetine in Japanese patients with major depressive disorder: a population pharmacokinetic analysis

OBJECTIVE

Although the reduced function of the cytochrome P450 2D6*10 (CYP2D6*10) allele is common among Asian populations, existing evidence does not support paroxetine therapy adjustments for patients who have the CYP2D6*10 allele. In this study, we attempted to evaluate the degree of the impact of different CYP2D6 genotypes on the pharmacokinetic (PK) variability of paroxetine in a Japanese population using a population PK approach.

METHODS

This retrospective study included 179 Japanese patients with major depressive disorder who were being treated with paroxetine. CYP2D6*1, *2, *5, *10, and *41 polymorphisms were observed. A total of 306 steady-state concentrations for paroxetine were collected from the patients. A nonlinear mixed-effects model identified the apparent Michaelis-Menten constant (Km) and the maximum velocity (Vmax) of paroxetine; the covariates included CYP2D6 genotypes, patient age, body weight, sex, and daily paroxetine dose.

RESULTS

The allele frequencies of CYP2D6*1, *2, *5, *10, and *41 were 39.4, 14.5, 4.5, 41.1, and 0.6%, respectively. There was no poor metabolizer who had two nonfunctional CYP2D6*5 alleles. A one-compartment model showed that the apparent Km value was decreased by 20.6% in patients with the CYP2D6*10/*10 genotype in comparison with the other CYP2D6 genotypes. Female sex also influenced the apparent Km values. No PK parameters were affected by the presence of one CYP2D6*5 allele.

CONCLUSION

Unexpectedly, elimination was accelerated in individuals with the CYP2D6*10/*10 genotype. Our results show that the presence of one CYP2D6*5 allele or that of any CYP2D6*10 allele may have no major effect on paroxetine PKs in the steady state.

Integrated semi-physiological pharmacokinetic model for both sunitinib and its active metabolite SU12662

AIMS

Previously published pharmacokinetic (PK) models for sunitinib and its active metabolite SU12662 were based on a limited dataset or lacked important elements such as correlations between sunitinib and its metabolite. The current study aimed to develop an improved PK model that circumvented these limitations and to prove the utility of the PK model in treatment optimization in clinical practice.

METHODS

One thousand two hundred and five plasma samples from 70 cancer patients were collected from three PK studies with sunitinib and SU12662. A semi-physiological PK model for sunitinib and SU12662 was developed incorporating pre-systemic metabolism using non-linear mixed effects modelling (nonmem). Allometric scaling based on body weight was applied. The final model was used for simulation of the PK of different treatment regimens.

RESULTS

Sunitinib and SU12662 PK were best described by a one and two compartment model, respectively. Introduction of pre-systemic formation of SU12662 strongly improved model fit, compared with solely systemic metabolism. The clearance of sunitinib and SU12662 was estimated at 35.7 (relative standard error (RSE) 5.7%) l h(-1) and 17.1 (RSE 7.4%) l h(-1), respectively for 70 kg patients. Correlation coefficients were estimated between inter-individual variability of both clearances, both volumes of distribution and between clearance and volume of distribution of SU12662 as 0.53, 0.48 and 0.45, respectively. Simulation of the PK model predicted correctly the ratio of patients who did not reach proposed PK targets for efficacy.

CONCLUSIONS

A semi-physiological PK model for sunitinib and SU12662 in cancer patients was presented including pre-systemic metabolism. The model was superior to previous PK models in many aspects.

Higher than expected clozapine serum level and clozapine/norclozapine ratio due to CYP450 gene polymorphisms

Aims: We investigated a case of the unexpectedly high clozapine/norclozapine serum levels in a non-smoker, female patient. Method: The RFLP genotype method was used for the loss-of-function CYP450 gene polymorphisms detection. Results: The genotype revealed the presence of CYP2D6*41/*41, rs28371725, NM_000106.5:c.985+39G>A and CYP1A2*1/*1D rs35694136, NM_000761.4:c.-1635delT. Conclusion: An appropriate clozapine serum level is required for both optimization of clinical response and avoidance of troublesome, or even life-threatening side effects. Although not the only factor, genetic polymorphisms responsible for the metabolism of clozapine play an important role towards such goals. Pharmacogenetic tests can guide dose adjustment in order to obtain an optimum serum level, which can lead to a better therapeutic response and minimization of harmful side effects.

Detection of an endogenous urinary biomarker associated with CYP2D6 activity using global metabolomics

AIM

We sought to discover endogenous urinary biomarkers of human CYP2D6 activity.

PATIENTS & METHODS

Healthy pediatric subjects (n = 189) were phenotyped using dextromethorphan and randomized for candidate biomarker selection and validation. Global urinary metabolomics was performed using liquid chromatography quadrupole time-of-flight mass spectrometry. Candidate biomarkers were tested in adults receiving fluoxetine, a CYP2D6 inhibitor.

RESULTS

A biomarker, M1 (m/z 444.3102) was correlated with CYP2D6 activity in both the pediatric training and validation sets. Poor metabolizers had undetectable levels of M1, whereas it was present in subjects with other phenotypes. In adult subjects, a 9.56-fold decrease in M1 abundance was observed during CYP2D6 inhibition.

CONCLUSION

Identification and validation of M1 may provide a noninvasive means of CYP2D6 phenotyping.

In vivo assessment of the metabolic activity of CYP2D6 diplotypes and alleles

AIMS

A prospectively enrolled patient cohort was used to assess whether the prediction of CYP2D6 phenotype activity from genotype data could be improved by reclassification of diplotypes or alleles.

METHODS

Three hundred and fifty-five patients receiving tamoxifen 20 mg were genotyped for CYP2D6 and tamoxifen metabolite concentrations were measured. The endoxifen : N-desmethly-tamoxifen metabolic ratio, as a surrogate of CYP2D6 activity, was compared across four diplotypes (EM/IM, EM/PM, IM/IM, IM/PM) that are typically collapsed into an intermediate metabolizer (IM) phenotype. The relative metabolic activity of each allele type (UM, EM, IM, and PM) and each EM and IM allele was estimated for comparison with the activity scores typically assigned, 2, 1, 0.5 and 0, respectively.

RESULTS

Each of the four IM diplotypes have distinct CYP2D6 activity from each other and from the EM and PM phenotype groups (each P < 0.05). Setting the activity of an EM allele at 1.0, the relative activities of a UM, IM and PM allele were 0.85, 0.67 and 0.52, respectively. The activity of the EM alleles were statistically different (P < 0.0001), with the CYP2D6*2 allele (scaled activity = 0.63) closer in activity to an IM than an EM allele. The activity of the IM alleles were also statistically different (P = 0.014).

CONCLUSION

The current systems for translating CYP2D6 genotype into phenotype are not optimally calibrated, particularly in regards to IM diplotypes and the *2 allele. Additional research is needed to improve the prediction of CYP2D6 activity from genetic data for individualized dosing of CYP2D6 dependent drugs.

Population pharmacokinetic modelling to assess the impact of CYP2D6 and CYP3A metabolic phenotypes on the pharmacokinetics of tamoxifen and endoxifen

AIMS

Tamoxifen is considered a pro-drug of its active metabolite endoxifen. The major metabolic enzymes involved in endoxifen formation are CYP2D6 and CYP3A. There is considerable evidence that variability in activity of these enzymes influences endoxifen exposure and thereby may influence the clinical outcome of tamoxifen treatment. We aimed to quantify the impact of metabolic phenotype on the pharmacokinetics of tamoxifen and endoxifen.

METHODS

We assessed the CYP2D6 and CYP3A metabolic phenotypes in 40 breast cancer patients on tamoxifen treatment with a single dose of dextromethorphan as a dual phenotypic probe for CYP2D6 and CYP3A. The pharmacokinetics of dextromethorphan, tamoxifen and their relevant metabolites were analyzed using non-linear mixed effects modelling.

RESULTS

Population pharmacokinetic models were developed for dextromethorphan, tamoxifen and their metabolites. In the final model for tamoxifen, the dextromethorphan derived metabolic phenotypes for CYP2D6 as well as CYP3A significantly (P < 0.0001) explained 54% of the observed variability in endoxifen formation (inter-individual variability reduced from 55% to 25%).

CONCLUSIONS

We have shown that not only CYP2D6, but also CYP3A enzyme activity influences the tamoxifen to endoxifen conversion in breast cancer patients. Our developed model may be used to assess separately the impact of CYP2D6 and CYP3A mediated drug-drug interactions with tamoxifen without the necessity of administering this anti-oestrogenic drug and to support Bayesian guided therapeutic drug monitoring of tamoxifen in routine clinical practice.

Cytochrome P450 2D6*10 genotype affects the pharmacokinetics of dimemorfan in healthy Chinese subjects

This study aimed to investigate the effects of cytochrome P450 2D6*10 (100C > T, rs1065852) genotype on the pharmacokinetics of dimemorfan in healthy Chinese subjects. Data were evaluated from 24 subjects in two pharmacokinetic studies who received an oral dose of 40 mg of dimemorfan syrup (n = 12) or dimemorfan tablet (n = 12) after providing written informed consent and being divided into three groups: subjects with CYP2D6*10 CC (n = 5), CYP2D6*10 CT (n = 11) and CYP2D6*10 TT (n = 8). CC homozygotes and CT heterozygotes were defined to be C allele carriers. The CYP2D6*10 was genotyped by polymerase chain reaction-restriction fragment length polymorphism. Dimemorfan was measured by LC-MS/MS. There was significant difference in C max, AUC0-t , AUC0-inf, V z , and CL values of dimemorfan observed among the three CYP2D6*10 genotype groups (GLM, (a) P < 0.05, co-dominant model). CYP2D6*10 under the recessive model (CC + TC vs TT) was significantly associated with pharmacokinetics of dimemorfan ((c) P < 0.05). The C max values were significantly higher in subjects with CYP2D6*10 TT (8.06 ± 4.43 ng/mL) than CYP2D6*10 CC (3.41 ± 2.79 ng/mL), CYP2D6*10 CT (3.11 ± 2.47 ng/mL), so was AUC0-inf. V z /F and CL/F of subjects with CYP2D6*10 TT homozygotes were the lowest. We demonstrated that cytochrome P450 2D6*10 (100C > T, rs1065852) polymorphism can affect the pharmacokinetics of dimemorfan in humans, not dosage forms.

Common CYP2D6 polymorphisms affecting alternative splicing and transcription: long-range haplotypes with two regulatory variants modulate CYP2D6 activity

Cytochrome P450 2D6 (CYP2D6) is involved in the metabolism of 25% of clinically used drugs. Genetic polymorphisms cause substantial variation in CYP2D6 activity and serve as biomarkers guiding drug therapy. However, genotype-phenotype relationships remain ambiguous except for poor metabolizers carrying null alleles, suggesting the presence of yet unknown genetic variants. Searching for regulatory CYP2D6 polymorphisms, we find that a SNP defining the CYP2D6*2 allele, rs16947 [R296C, 17-60% minor allele frequency (MAF)], previously thought to convey normal activity, alters exon 6 splicing, thereby reducing CYP2D6 expression at least 2-fold. In addition, two completely linked SNPs (rs5758550/rs133333, MAF 13-42%) increase CYP2D6 transcription more than 2-fold, located in a distant downstream enhancer region (>100 kb) that interacts with the CYP2D6 promoter. In high linkage disequilibrium (LD) with each other, rs16947 and the enhancer SNPs form haplotypes that affect CYP2D6 enzyme activity in vivo. In a pediatric cohort of 164 individuals, rs16947 alone (minor haplotype frequency 28%) was associated with reduced CYP2D6 metabolic activity (measured as dextromethorphan/metabolite ratios), whereas rs5758550/rs133333 alone (frequency 3%) resulted in increased CYP2D6 activity, while haplotypes containing both rs16947 and rs5758550/rs133333 were similar to the wild-type. Other alleles used in biomarker panels carrying these variants such as CYP2D6*41 require re-evaluation of independent effects on CYP2D6 activity. The occurrence of two regulatory variants of high frequency and in high LD, residing on a long haplotype, highlights the importance of gene architecture, likely shaped by evolutionary selection pressures, in determining activity of encoded proteins.

Complexities of CYP2D6 gene analysis and interpretation

Cytochrome P450 2D6 (CYP2D6) plays an important role in the metabolism and bioactivation of about 25% of clinically used drugs including many antidepressants, antipsychotics and opioids. CYP2D6 activity is highly variably ranging from no activity in so-called poor metabolizers to ultrarapid metabolism at the other end of the extreme of the activity distribution. A large portion of this variability can be explained by the highly polymorphic nature of the CYP2D6 gene locus for which > 100 variants and subvariants identified to date. Allele frequencies vary markedly between ethnic groups; some have exclusively or predominantly only been observed in certain populations. Pharmacogenetic testing holds the promise of individualizing drug therapy by identifying patients with CYP2D6 diplotypes that puts them at an increased risk of experiencing dose-related adverse events or therapeutic failure. Inferring a patient's CYP2D6 metabolic capacity, or phenotype, however, is a challenging task due to the complexity of the CYP2D6 gene locus. Allelic variation includes SNPs, small insertions and deletions, gene copy number variation and rearrangements with CYP2D7, a highly related non-functional gene. This review provides a summary of the intricacies of CYP2D6 variation and genotype analysis, knowledge that is invaluable for the translation of genotype into clinically useful information.

Characterization of the CYP2D6 gene locus and metabolic activity in Indo- and Afro-Trinidadians: discovery of novel allelic variants

BACKGROUND

The highly polymorphic CYP2D6 gene has extensively been studied in many populations, but there is a void of knowledge regarding CYP2D6 pharmacogenetics and activity in populations with unique ancestries and admixture, such as those residing in Trinidad and Tobago.

MATERIALS & METHODS

167 healthy Indo- and 103 Afro-Trinidadians were phenotyped with dextromethorphan and extensively genotyped. Gene resequencing was performed to resolve cases with genotype/phenotype discordance.

RESULTS

CYP2D6 activity did not differ between the Indo-Trinidadians and Afro-Trinidadians. Poor metabolizers were, however, more frequent in the Indo-Trinidadians (4.19 vs 1.94%), and unique allele frequency patterns were observed. Two novel nonfunctional allelic variants were found among the Indo-Trinidadians in two discordant cases. CYP2D6*100 is characterized by a single nucleotide deletion and CYP2D6*101 by a 19-bp deletion; both cause frameshifts.

CONCLUSION

Our study underscores the importance of thoroughly characterizing the genetic make up of unique populations when considering pharmacogenetic testing for individualized therapy.

CYP2D6 genotype affects age-related decline in flecainide clearance: a population pharmacokinetic analysis

OBJECTIVE

To investigate the association between age-related decline in flecainide clearance and CYP2D6 genotype, we conducted a population pharmacokinetic analysis of flecainide using routine therapeutic drug monitoring data.

METHODS

Population pharmacokinetic analysis was performed on retrospective data from 163 genotyped patients treated with oral flecainide for supraventricular tachyarrhythmias. The CYP2D6 genotype was categorized as CYP2D6 homozygous extensive metabolizers (hom-EMs; n=57), heterozygous extensive metabolizers (het-EMs; n=79), and intermediate metabolizers and poor metabolizers (IMs/PMs; n=27).

RESULTS

Population pharmacokinetic analysis revealed that estimated glomerular filtration rate, body weight, female sex, and aging were important factors for estimating flecainide clearance. The metabolic clearance was decreased age dependently in a curvilinear fashion, where the lower clearance was observed in greater than 60 years for het-EMs and greater than 55 years for IMs/PMs. The reduction in metabolic clearance in elderly (70 years) patients compared with middle-aged (52 years) patients was different among the CYP2D6 genotype groups: 22.1 and 49.5% in CYP2D6 het-EMs and IMs/PMs, respectively, and no change in hom-EMs. A 11.4% reduction in estimated glomerular filtration rate in elderly patients compared with middle-aged patients corresponded to 6.1% decline in flecainide clearance. Overall, the age-related decline in flecainide clearance was 6.1% in hom-EMs, 16.3% in het-EMs, and 28.9% in IMs/PMs groups.

CONCLUSION

This study suggests that CYP2D6 genotype is a determinant factor of age-related decline in flecainide clearance.

A physiologically based pharmacokinetic model to predict disposition of CYP2D6 and CYP1A2 metabolized drugs in pregnant women

Conducting pharmacokinetic (PK) studies in pregnant women is challenging. Therefore, we asked if a physiologically based pharmacokinetic (PBPK) model could be used to evaluate different dosing regimens for pregnant women. We refined and verified our previously published pregnancy PBPK model by incorporating cytochrome P450 CYP1A2 suppression (based on caffeine PK) and CYP2D6 induction (based on metoprolol PK) into the model. This model accounts for gestational age-dependent changes in maternal physiology and hepatic CYP3A activity. For verification, the disposition of CYP1A2-metabolized drug theophylline (THEO) and CYP2D6-metabolized drugs paroxetine (PAR), dextromethorphan (DEX), and clonidine (CLO) during pregnancy was predicted. Our PBPK model successfully predicted THEO disposition during the third trimester (T3). Predicted mean postpartum to third trimester (PP:T3) ratios of THEO area under the curve (AUC), maximum plasma concentration, and minimum plasma concentration were 0.76, 0.95, and 0.66 versus observed values 0.75, 0.89, and 0.72, respectively. The predicted mean PAR steady-state plasma concentration (Css) ratio (PP:T3) was 7.1 versus the observed value 3.7. Predicted mean DEX urinary ratio (UR) (PP:T3) was 2.9 versus the observed value 1.9. Predicted mean CLO AUC ratio (PP:T3) was 2.2 versus the observed value 1.7. Sensitivity analysis suggested that a 100% induction of CYP2D6 during T3 was required to recover the observed PP:T3 ratios of PAR Css, DEX UR, and CLO AUC. Based on these data, it is prudent to conclude that the magnitude of hepatic CYP2D6 induction during T3 ranges from 100 to 200%. Our PBPK model can predict the disposition of CYP1A2, 2D6, and 3A drugs during pregnancy.

Detection and characterization of the CYP2D6*9x2 gene duplication in two Spanish populations: resolution of AmpliChip CYP450 test no-calls

BACKGROUND

CYP2D6 is a major drug-metabolizing enzyme. Polymorphic variation includes copy number variants such as gene deletions, duplications and multiplications of functional and nonfunctional gene units. In this article we describe the first systematic characterization of a CYP2D6*9x2 gene duplication. CYP2D6*9 is an allelic variant conferring reduced enzymatic activity. This novel gene duplication was discovered in two unrelated Spanish psychiatric patients. Both subjects were initially tested with the AmpliChip CYP450 test, which indicated the presence of a duplication and the CYP2D6*9 allele, but did not make a genotype call. The goal of the study was to resolve this issue by characterizing the CYP2D6 gene locus in these patients.

MATERIALS & METHODS

Both individuals and one offspring were regenotyped using our own CYP2D6 genotyping strategy employing long-range PCR and TaqMan-based SNP detection. In addition, gene resequencing and genotyping of duplication-specific long-range PCR products and quantitative gene copy number analysis was applied.

RESULTS

The duplication was mapped to the CYP2D6*9 allele and copy number analysis determined a CYP2D6*9x2 gene duplication in all three individuals. Because CYP2D6*9x2 is not recognized by the AmpliChip CYP450 test, this structural arrangement was responsible for the 'no call' on the AmpliChip CYP450 test report.

CONCLUSION

The full characterization of this allele will aid in the interpretation of AmpliChip CYP450 test results for clinical and research applications. Original submitted 8 June 2011; Revision submitted 18 July 2011.

Pharmacogenomics of CYP2D6: molecular genetics, interethnic differences and clinical importance

CYP2D6 has received intense attention since the beginning of the pharmacogenetic era in the 1970s. This is because of its involvement in the metabolism of more than 25% of the marketed drugs, the large geographical and inter-ethnic differences in the genetic polymorphism and possible drug-induced toxicity. Many interesting reviews have been published on CYP2D6 and this review aims to reinstate the importance of the genetic polymorphism of CYP2D6 in different populations as well as some clinical implications and important drug interactions.

Prediction of relative in vivo metabolite exposure from in vitro data using two model drugs: dextromethorphan and omeprazole

Metabolites can have pharmacological or toxicological effects, inhibit metabolic enzymes, and be used as probes of drug-drug interactions or specific cytochrome P450 (P450) phenotypes. Thus, better understanding and prediction methods are needed to characterize metabolite exposures in vivo. This study aimed to test whether in vitro data could be used to predict and rationalize in vivo metabolite exposures using two model drugs and P450 probes: dextromethorphan and omeprazole with their primary metabolites dextrorphan, 5-hydroxyomeprazole (5OH-omeprazole), and omeprazole sulfone. Relative metabolite exposures were predicted using metabolite formation and elimination clearances. For dextrorphan, the formation clearances of dextrorphan glucuronide and 3-hydroxymorphinan from dextrorphan in human liver microsomes were used to predict metabolite (dextrorphan) clearance. For 5OH-omeprazole and omeprazole sulfone, the depletion rates of the metabolites in human hepatocytes were used to predict metabolite clearance. Dextrorphan/dextromethorphan in vivo metabolite/parent area under the plasma concentration versus time curve ratio (AUC(m)/AUC(p)) was overpredicted by 2.1-fold, whereas 5OH-omeprazole/omeprazole and omeprazole sulfone/omeprazole were predicted within 0.75- and 1.1-fold, respectively. The effect of inhibition or induction of the metabolite's formation and elimination on the AUC(m)/AUC(p) ratio was simulated. The simulations showed that unless metabolite clearance pathways are characterized, interpretation of the metabolic ratios is exceedingly difficult. This study shows that relative in vivo metabolite exposure can be predicted from in vitro data and characterization of secondary metabolism of probe metabolites is critical for interpretation of phenotypic data.

Frequency of undetected CYP2D6 hybrid genes in clinical samples: impact on phenotype prediction

Cytochrome P450 2D6 (CYP2D6) is highly polymorphic. CYP2D6-2D7 hybrid genes can be present in samples containing CYP2D6*4 and CYP2D6*10 alleles. CYP2D7-2D6 hybrid genes can be present in samples with duplication signals and in samples with homozygous genotyping results. The frequency of hybrid genes in clinical samples is unknown. We evaluated 1390 samples for undetected hybrid genes by polymerase chain reaction (PCR) amplification, PCR fragment analysis, TaqMan copy number assays, DNA sequencing, and allele-specific primer extension assay. Of 508 CYP2D6*4-containing samples, 109 (21.5%) harbored CYP2D6*68 + *4-like, whereas 9 (1.8%) harbored CYP2D6*4N + *4-like. Of 209 CYP2D6*10-containing samples, 44 (21.1%) were found to have CYP2D6*36 + *10. Of 332 homozygous samples, 4 (1.2%) harbored a single CYP2D7-2D6 hybrid, and of 341 samples with duplication signals, 25 (7.3%) harbored an undetected CYP2D7-2D6 hybrid. Phenotype before and after accurate genotyping was predicted using a method in clinical use. The presence of hybrid genes had no effect on the phenotype prediction of CYP2D6*4- and CYP2D6*10-containing samples. Four of four (100%) homozygous samples containing a CYP2D7-2D6 gene had a change in predicted phenotype, and 23 of 25 (92%) samples with a duplication signal and a CYP2D7-2D6 gene had a change in predicted phenotype. Four novel genes were identified (CYP2D6*13A1 variants 1 and 2, CYP2D6*13G1, and CYP2D6*13G2), and two novel hybrid tandem structures consisting of CYP2D6*13B + *68×2 + *4-like and CYP2D6*13A1 variant 2 + *1×N were observed.

Carvedilol pharmacokinetics and pharmacodynamics in relation to CYP2D6 and ADRB pharmacogenetics

AIMS

Carvedilol is an effective treatment in hypertension and chronic heart failure. The medical impact of polymorphisms in CYP2D6 and in the β-adrenergic receptors ADRB1 and ADRB2 on the pharmacokinetics and pharmacodynamics of carvedilol is controversial.

METHODS

After carvedilol 25 mg was administered to 110 volunteers, concentrations were enantioselectively quantified and effects on resting and exercise-induced heart rate and blood pressure were analyzed using population pharmacokinetic, pharmacodynamic and pharmacogenetic modeling.

RESULTS

There were significant CYP2D6 allele-specific differences in carvedilol pharmacokinetics, but the CYP2D6 genotype had no effect on heart rate, blood pressure or adverse effects. ADRB1 Gly49 was associated with higher baseline heart rates and with greater carvedilol effects on exercise heart rates. Carriers of ADRB2 Gln27 had greater reduction in resting blood pressure by carvedilol compared with Glu27.

CONCLUSION

Carvedilol is a drug where CYP2D6-related pharmacokinetic variation is apparently not carried forward into pharmacodynamic variation. Although current knowledge does not allow utilizing ADRB1 and ADRB2 genotypes for clinical treatment decisions, our data should stimulate further research on the impact of these genotypes in health and disease.