Stereoselective pharmacokinetics of stable isotope (+/-)-[13C]-pantoprazole: Implications for a rapid screening phenotype test of CYP2C19 activity

AIMS: We have previously shown that the (±)-[(13) C]-pantoprazole breath test is a promising noninvasive probe of CYP2C19 activity. As part of that trial, plasma, breath test indices and CYP2C19 (*2, *3, and *17) genotype were collected. Here, we examined whether [(13) C]-pantoprazole exhibits enantioselective pharmacokinetics and whether this enantioselectivity is correlated with indices of breath test. METHODS: Plasma (-)- and (+)-[(13) C]-pantoprazole that were measured using a chiral HPLC were compared between CYP2C19 genotypes and correlated with breath test indices. RESULTS: The AUC( 0-∞) of (+)-[(13) C]-pantoprazole in PM (*2/*2, n = 4) was 10.1- and 5.6-fold higher that EM (*1/*1or *17, n = 10) and IM (*1/*2or *3, n = 10) of CYP2C19, respectively (P < 0.001). The AUC( 0-∞) of (-)-[(13) C]-pantoprazole only significantly differed between PMs and EMs (1.98-fold; P = 0.05). The AUC( 0-∞) ratio of (+)-/(-)-[(13) C]-pantoprazole was 3.45, 0.77, and 0.67 in PM, IM, and EM genotypes, respectively. Breath test index, delta over baseline show significant correlation with AUC( 0-∞) of (+)-[(13) C]-pantoprazole (Pearson's r = 0.62; P < 0.001). CONCLUSIONS: [(13) C]-pantoprazole exhibits enantioselective elimination. (+)-[(13) C]-pantoprazole is more dependent on CYP2C19 metabolic status and may serve as a more attractive probe of CYP2C19 activity than (-)-[(13) C]-pantoprazole or the racemic mixture

Rifampin modulation of xeno- and endobiotic conjugating enzyme mRNA expression and associated microRNAs in human hepatocytes

Rifampin is a pleiotropic inducer of multiple drug metabolizing enzymes and transporters. This work utilized a global approach to evaluate rifampin effects on conjugating enzyme gene expression with relevance to human xeno- and endo-biotic metabolism. Primary human hepatocytes from 7 subjects were treated with rifampin (10 μmol/L, 24 hours). Standard methods for RNA-seq library construction, EZBead preparation, and NextGen sequencing were used to measure UDP-glucuronosyl transferase UGT, sulfonyltransferase SULT, N acetyltransferase NAT, and glutathione-S-transferase GST mRNA expression compared to vehicle control (0.01% MeOH). Rifampin-induced (>1.25-fold) mRNA expression of 13 clinically important phase II drug metabolizing genes and repressed (>1.25-fold) the expression of 3 genes (P < .05). Rifampin-induced miRNA expression changes correlated with mRNA changes and miRNAs were identified that may modulate conjugating enzyme expression. NAT2 gene expression was most strongly repressed (1.3-fold) by rifampin while UGT1A4 and UGT1A1 genes were most strongly induced (7.9- and 4.8-fold, respectively). Physiologically based pharmacokinetic modeling (PBPK) was used to simulate the clinical consequences of rifampin induction of CYP3A4- and UGT1A4-mediated midazolam metabolism. Simulations evaluating isolated UGT1A4 induction predicted increased midazolam N-glucuronide exposure (~4-fold) with minimal reductions in parent midazolam exposure (~10%). Simulations accounting for simultaneous induction of both CYP3A4 and UGT1A4 predicted a ~10-fold decrease in parent midazolam exposure with only a ~2-fold decrease in midazolam N-glucuronide metabolite exposure. These data reveal differential effects of rifampin on the human conjugating enzyme transcriptome and potential associations with miRNAs that form the basis for future mechanistic studies to elucidate the interplay of conjugating enzyme regulatory elements

Population Pharmacokinetic Modeling To Estimate the Contributions of Genetic and Nongenetic Factors to Efavirenz Disposition

Efavirenz pharmacokinetics is characterized by large between-subject variability, which determines both therapeutic response and adverse effects. Some of the variability in efavirenz pharmacokinetics has been attributed to genetic variability in cytochrome P450 genes that alter efavirenz metabolism, such as CYP2B6 and CYP2A6. While the effects of additional patient factors have been studied, such as sex, weight, and body mass index, the extent to which they contribute to variability in efavirenz exposure is inconsistently reported. The aim of this analysis was to develop a pharmacometric model to quantify the contribution of genetic and nongenetic factors to efavirenz pharmacokinetics. A population-based pharmacokinetic model was developed using 1,132 plasma efavirenz concentrations obtained from 73 HIV-seronegative volunteers administered a single oral dose of 600 mg efavirenz. A two-compartment structural model with absorption occurring by zero- and first-order processes described the data. Allometric scaling adequately described the relationship between fat-free mass and apparent oral clearance, as well as fat mass and apparent peripheral volume of distribution. Inclusion of fat-free mass and fat mass in the model mechanistically accounted for correlation between these disposition parameters and sex, weight, and body mass index. Apparent oral clearance of efavirenz was reduced by 25% and 51% in subjects predicted to have intermediate and slow CYP2B6 metabolizer status, respectively. The final pharmacokinetic model accounting for fat-free mass, fat mass, and CYP2B6 metabolizer status was consistent with known mechanisms of efavirenz disposition, efavirenz physiochemical properties, and pharmacokinetic theory. (This study has been registered at ClinicalTrials.gov under identifier NCT00668395.

Genome-Wide Discovery of Drug-Dependent Human Liver Regulatory Elements

Inter-individual variation in gene regulatory elements is hypothesized to play a causative role in adverse drug reactions and reduced drug activity. However, relatively little is known about the location and function of drug-dependent elements. To uncover drug-associated elements in a genome-wide manner, we performed RNA-seq and ChIP-seq using antibodies against the pregnane X receptor (PXR) and three active regulatory marks (p300, H3K4me1, H3K27ac) on primary human hepatocytes treated with rifampin or vehicle control. Rifampin and PXR were chosen since they are part of the CYP3A4 pathway, which is known to account for the metabolism of more than 50% of all prescribed drugs. We selected 227 proximal promoters for genes with rifampin-dependent expression or nearby PXR/p300 occupancy sites and assayed their ability to induce luciferase in rifampin-treated HepG2 cells, finding only 10 (4.4%) that exhibited drug-dependent activity. As this result suggested a role for distal enhancer modules, we searched more broadly to identify 1,297 genomic regions bearing a conditional PXR occupancy as well as all three active regulatory marks. These regions are enriched near genes that function in the metabolism of xenobiotics, specifically members of the cytochrome P450 family. We performed enhancer assays in rifampin-treated HepG2 cells for 42 of these sequences as well as 7 sequences that overlap linkage-disequilibrium blocks defined by lead SNPs from pharmacogenomic GWAS studies, revealing 15/42 and 4/7 to be functional enhancers, respectively. A common African haplotype in one of these enhancers in the GSTA locus was found to exhibit potential rifampin hypersensitivity. Combined, our results further suggest that enhancers are the predominant targets of rifampin-induced PXR activation, provide a genome-wide catalog of PXR targets and serve as a model for the identification of drug-responsive regulatory elements

A Penalized Mixture Model Approach in Genotype/Phenotype Association Analysis for Quantitative Phenotypes

A mixture normal model has been developed to partition genotypes in predicting quantitative phenotypes. Its estimation and inference are performed through an EM algorithm. This approach can conduct simultaneous genotype clustering and hypothesis testing. It is a valuable method for predicting the distribution of quantitative phenotypes among multi-locus genotypes across genes or within a gene. This mixture model’s performance is evaluated in data analyses for two pharmacogenetics studies. In one example, thirty five CYP2D6 genotypes were partitioned into three groups to predict pharmacokinetics of a breast cancer drug, Tamoxifen, a CYP2D6 substrate (p-value = 0.04). In a second example, seventeen CYP2B6 genotypes were categorized into three clusters to predict CYP2B6 protein expression (p-value = 0.002). The biological validities of both partitions are examined using established function of CYP2D6 and CYP2B6 alleles. In both examples, we observed genotypes clustered in the same group to have high functional similarities. The power and recovery rate of the true partition for the mixture model approach are investigated in statistical simulation studies, where it outperforms another published method

Efavirenz inhibits the human ether-a-go-go related current (hERG) and induces QT interval prolongation in CYP2B6*6*6 allele carriers

Background Efavirenz (EFV) has been associated with torsade de pointes despite marginal QT interval lengthening. Since EFV is metabolized by the cytochrome P450 (CYP) 2B6 enzyme, we hypothesized that EFV would lengthen the rate-corrected QT (QTcF) interval in carriers of the CYP2B6*6 decreased functional allele. Objective The primary objective of this study was to evaluate EFV-associated QT interval changes with regard to CYP2B6 genotype and to explore mechanisms of QT interval lengthening. Methods EFV was administered to healthy volunteers (n=57) as a single 600 mg dose followed by multiple doses to steady-state. Subjects were genotyped for known CYP2B6 alleles and ECGs and EFV plasma concentrations were obtained serially. Whole-cell, voltage-clamp experiments were performed on cells stably expressing hERG and exposed to EFV in the presence and absence of CYP2B6 expression. Results EFV demonstrated a gene-dose effect and exceeded the FDA criteria for QTcF interval prolongation in CYP2B6*6/*6 carriers. The largest mean time-matched differences ΔΔQTcF were observed at 6 hrs (14 ms; 95% CI [1; 27]), 12 hrs (18 ms; 95% CI [−4; 40] and 18 hrs (6 ms; 95% CI [−1; 14]) in the CYP2B6*6/*6 genotype. EFV concentrations exceeding 0.4 µg/mL significantly inhibited outward hERG tail currents (P<0.05). Conclusions This study demonstrates that homozygous carriers of CYP2B6*6 allele may be at increased risk for EFV-induced QTcF interval prolongation via inhibition of hERG

Rapid Identification of the Hepatic Cytochrome P450 2C19 Activity Using a Novel and Noninvasive [ 13 C]Pantoprazole Breath Test

ABSTRACT We tested the hypothesis that the stable isotope [ 13 C]pantoprazole is O-demethylated by cytochrome P450 CYP2C19 and that the 13 CO 2 produced and exhaled in breath as a result can serve as a safe, rapid, and noninvasive phenotyping marker of CYP2C19 activity in vivo. Healthy volunteers who had been genotyped for the CYP2C19*2, CYP2C19*3, and CYP2C19*17 alleles were administered a single oral dose of [ 13 C]pantoprazole sodium-sesquihydrate (100 mg) with 2.1 g of sodium bicarbonate. Exhaled 13 CO 2 and 12 CO 2 were measured by IR spectroscopy before (baseline) and 2.5 to 120 min after dosing. Ratios of 13 CO 2 / 12 CO 2 after [ 13 C]pantoprazole relative to 13 CO 2 / 12 CO 2 at baseline were expressed as change over baseline (DOB). Maximal DOB, DOB 15 to DOB 120 , and area under the DOB versus time curve (AUC 0 -120 and AUC 0 -ϱ ) were significantly different among three genotype groups (CYP2C19*1/ *1, n ϭ 10; CYP2C19*1/*2 or CYP2C19*1/*3, n ϭ 10; and CYP2C19*2/*2, n ϭ 5) with predicted extensive metabolizers (EMs), intermediate metabolizers (IMs), and poor metabolizers (PMs) of CYP2C19, respectively (Kruskal-Wallis test, p Ͻ 0.01); linear regression analysis indicated a gene-dose effect relationship (r 2 ranged between 0.236 and 0.522; all p Ͻ 0.05). These breath test indices were significantly lower in PMs than IMs (p Ͻ 0.05) or EMs (p Ͻ 0.01) of CYP2C19. [ 13 C]Pantoprazole plasma exposure showed significant inverse correlation with breath test indices in the respective subjects (Pearson r ϭ Ϫ0.74; p ϭ 0.038). These feasibility data suggest that the [ 13 C]pantoprazole breath test is a reliable, rapid, and noninvasive probe of CYP2C19 and seems to be a useful tool to optimize drug therapy metabolized by CYP2C19

Drug–gene and drug–drug interactions associated with tramadol and codeine therapy in the INGENIOUS trial

Background: Tramadol and codeine are metabolized by CYP2D6 and are subject to drug-gene and drug-drug interactions. Methods: This interim analysis examined prescribing behavior and efficacy in 102 individuals prescribed tramadol or codeine while receiving pharmaco-genotyping as part of the INGENIOUS trial (NCT02297126). Results: Within 60 days of receiving tramadol or codeine, clinicians more frequently prescribed an alternative opioid in ultrarapid and poor metabolizers (odds ratio: 19.0; 95% CI: 2.8-160.4) as compared with normal or indeterminate metabolizers (p = 0.01). After adjusting the CYP2D6 activity score for drug-drug interactions, uncontrolled pain was reported more frequently in individuals with reduced CYP2D6 activity (odds ratio: 0.50; 95% CI: 0.25-0.94). Conclusion: Phenoconversion for drug-drug and drug-gene interactions is an important consideration in pharmacogenomic implementation; drug-drug interactions may obscure the potential benefits of genotyping

17-Hydroexemestane: A Potent Inhibitor of CYP19 (Aromatase) and Substrate of CYP3A

17-hydroexemestane is the major metabolite of exemestane in vivo. Previous studies have shown that 17-hydroexemestane is androgenic and bone protective. Due to structure similarities, we hypothesized that, like exemestane, 17-hydroexemestane is an inhibitor of aromatase (CYP19). Our aim was to assess the potency (IC50) of 17-hydroexemestane toward CYP19 inhibition, and to determine the specific CYPsresponsible for 17-hydroexemestane metabolism. Using recombinant human CYP19, we investigated the ability of exemestane and 17-hydroexemestane to block the formation of estradiol from testosterone. We found that 17-hydroexemestane potently inhibited aromatase. IC 50 values for the inhibition of CYP19 by exemestane and 17-hydroexemestane were 1.5 μM and 3 μM, respectively. Furthermore, using recombinant human P450s, human liver microsomes, and HPLC analytical techniques, we identified one major metabolite (MIII) of 17-hydroexemestane in the human liver microsomal incubate. In a bank of 15 well-characterized HLMs, MIII formation rate was significantly correlated with the activity of CYP3A (rs= 0.78, p=0.001).In a panel of baculovirus-expressed CYP enzymes, only CYP3A4 and CYP3A5 catalyzed MIII formation at the highest rate. In sum, these in vitro data suggest that 17-hydroexemestane is a potent inhibitor of CYP19 and that CYP3A plays a major role in its metabolism. Whether genetic polymorphisms and drug interactions involving these enzymes may contribute to the disposition and action of 17-hydroexemestane in breast cancer patients remains to be elucidated