How precise is quantitative prediction of pharmacokinetic effects due to drug-drug interactions and genotype from in vitro data? A comprehensive analysis on the example CYP2D6 and CYP2C19 substrates

https://doi.org/10.1016/j.pharmthera.2020.107629Get rights and content

Abstract

Drug-drug interactions (DDI) and genomic variation (PG) can lead to dangerously high blood and tissue concentrations with some drugs but may be negligible with other drugs. Using a quantitative metaanalysis, we analyzed on the example of CYP2D6 and CYP2C19 substrates, how well the effects of DDI and PG can be predicted by in vitro methods. In addition, we analyzed the quantitative effect of prototypic inhibitors of the two enzymes in relation to their genetic deficiency. More than 600 published studies were screened which compared either human pharmacokinetics with and without comedication, or which compared human pharmacokinetics of deficient with extensive metabolizers, or which assessed metabolism by in vitro approaches. With human liver microsomes, the in vitro to in vivo agreement of fractional clearances was reasonably high if loss of substrate was quantified in the in vitro assays performed with and without enzyme specific inhibitors. Also a generally very high correlation between the clinical pharmacokinetic effects of inherited deficiency and inhibition by drug-drug interactions could be demonstrated. Most cases of poor correlation were explained by the lack of CYP2D6 versus CYP2C19 specificity of fluoxetine or by a poor knowledge of the quantitative contribution of the metabolic pathways if metabolite formation was quantified in the in vitro assays. The good correspondence of the in vitro data with clinical DDI and clinical PG studies may be a good basis for future application of these methods in drug development and drug therapy.

Introduction

Pharmacokinetic drug-drug interactions (DDI) may lead to serious adverse reactions and have caused numerous market withdrawals in the past (Mullins et al., 1998; Wienkers & Heath, 2005). Similarly, pharmacogenetic variation (PG) may compromise a safe and effective therapy with drugs. Therefore, consideration of DDI and PG in drug development and in daily medical care of patients is essential. But this is complicated by the fact, that the quantitative relevance of DDI or PG on systemic drug exposure widely differs between the different drugs. There are drugs for which DDI and PG effects are negligibly small, for other drugs DDI and PG effects are moderate and for some drugs they are very significant (Stingl, Brockmöller, & Viviani, 2013). Thus, a quantitative understanding and prediction of the magnitude of these effects is important, and this is particularly relevant to provide concise recommendations for an individually adjusted drug therapy (Kirchheiner, Fuhr, & Brockmoller, 2005).

Cytochrome P450 enzymes 2D6 and 2C19, together, are involved in biotransformation of about 25% of all small molecule drugs. About 8% and 3% of Europeans are so called poor metabolizers (PM) of CYP2D6 and CYP2C19 substrates, respectively (Ingelman-Sundberg, Oscarson, & McLellan, 1999). These PM are “natural knock-outs” providing unique opportunities to study the biological effects of the two enzymes in human beings. In the analysis presented here we utilize the unique opportunity of having a zero-activity comparison group (the PM) to assess the validity of frequently used in vivo methods for prediction of the fractional clearance via a specific enzyme.

In vitro methods, so-called reaction phenotyping, aim to identify the enzymes involved in the metabolism of a specific drug and their quantitative contribution to biotransformation best defined by fractional clearance (Houston & Galetin, 2003; Rodrigues, 1999; Venkatakrishnan et al., 2001a, Venkatakrishnan et al., 2001b; Wood, Houston, & Hallifax, 2017). Thus, those assays aim to predict the enzyme-specific contribution to in vivo pharmacokinetics in humans. Typically used reaction phenotyping methods (Table 1) are either performed as substrate depletion assay, quantifying loss of the substrate, or as metabolite formation assay quantifying formation of one or several metabolites. Both approaches have their advantages and disadvantages (Yang, Atkinson, & Di, 2016).

In vivo, the difference in total clearance between poor and extensive metabolizers (EM) corresponds to the partial clearance via the respective polymorphic enzyme (Kirchheiner et al., 2004; Kirchheiner, Meineke, Müller, Roots, & Brockmoller, 2002). In vitro data has been used to predict that difference in systemic exposure between extensive and poor metabolizers of CYP2D6 substrates (Gibbs, Hyland, & Youdim, 2006). Vice versa, data from in vitro approaches may then be validated using clinical pharmacokinetic data. Comedication with strong inhibitors can have effects similar to genetic deficiency (Rost, Brockmoller, Esdorn, & Roots, 1995). In vitro predictions of fractional clearance via a specific enzyme may therefore also be validated with in vivo data from drug-drug interaction studies (DDI).

Although the potential of in vitro methods to predict the in vivo pharmacokinetics has been shown for numerous single drugs, there are few comprehensive analyses on many drugs. And there are even less analyses comparing both, pharmacogenetics and drug-drug interactions, and how well both can be predicted by the in vitro methods. In vitro, there are several methods, but it is controversial which one is the best. Thus, here we want to extend research aiming to validate in vitro reaction phenotyping methods with in vivo data from PG studies and/or with in vivo data from DDI studies (Parmentier et al., 2019).

Thus, here we compared data from reaction phenotyping assays with in vivo PG and in vivo DDI data for substrates of CYP2D6 and CYP2C19. Our search strategy aimed to identify all available in vitro and in vivo data on CYP2D6 and CYP2C19 substrates suitable for our analysis. First, we compared reaction phenotyping results with clinical PG study data, second, we compared reaction phenotyping results with clinical drug-drug interaction data and, finally, we compared clinical PG data with drug-drug interaction data. First, this should allow to compare validity of different reaction phenotyping methods. Second, it should be possible to compare the in vitro assays' ability to predict PG as compared to DDI effects. Third, it will also be possible to compare in vivo PG effects with in vivo DDI effects. Finally, joint consideration of in vitro data, clinical PG and clinical DDI data may result in the best fractional clearance data.

Section snippets

Search criteria

We searched for all drugs known to be substrates of CYP2D6 and/or CYP2C19 (Fig. 1). Inclusion into the analyses was possible only if at least one in vitro study and one in vivo study (PG or DDI) was available. But the supplemental data as well as the in vivo only analyses also include in vivo data without corresponding in vitro data. For racemic drugs we included the enantiomers if reported in any publication. Our search was explicitly not restricted to major substrates of the two enzymes but

Requirements for data inclusion

We included all in vitro data on reaction phenotyping with human liver microsomes or single recombinant enzymes as long as fractional clearance was reported or could be extracted.

In metabolite formation assays a fractional clearance for every combination of CYP enzyme and every primary metabolite has to be taken into consideration. To calculate the overall fractional clearance of the parent compound here, it was necessary to weight fractional clearances by the fraction of parent drug eliminated

Processing of in vitro data

If data from more than one laboratory were available, we calculated means and standard deviations of fractional clearance for every combination of drug and enzyme.

Fractional clearances via CYP2D6 and/or CYP2C19 were multiplied with the extrarenally excreted fraction of the drug because in vitro reaction phenotyping does not take into account unchanged renal excretion. Unchanged renal elimination was taken from various sources (Brunton, Chabner, & Knollman, 2011; Haefeli, 2019). For enantiomers

Processing of in vivo data

For vivo studies, geometric means (if not available, alternatively arithmetic means or medians) of the parameter with standard deviations were gathered. We calculated the PG PK ratio AUCPM/AUCEM (eq. 1) reflecting the ratio of systemic exposure to the drug in poor metabolizers (PM) over extensive metabolizers (EM). And we calculated the DDI PK ratio AUCinhibited/AUCnot inhibited (eq. 2) reflecting the ratio of systemic exposure with enzyme-specific inhibition versus no inhibition.

Processing of data for in vitro in vivo comparison

The fractional clearance as estimated in vitro was compared against the fractional clearances (yi^) as derived from the in vivo PG results according to yi^≈1-(AUCEM/AUCPM), eq. 6, and from the in vivo DDI results according to yi^≈1-(AUCnot inhibited/AUCinhibited), eq. 7 (Gibbs et al., 2006; VandenBrink, Foti, Rock, Wienkers, & Wahlstrom, 2012). With therapeutic drug monitoring data (TDM) calculations were performed accordingly, for studies with clearance information reciprocal calculation was

Statistical power calculations

To perform comparable statistical power calculations over all clinical PG and DDI studies, we utilized Cohen's d which was estimated for the PG studies according to the following equation (eq. 8):d=x¯1x¯2n11s12+n21s22/n1+n22

With x¯1and x¯2 representing the mean AUCs (or trough levels in TDM) in PM and EM (for clearance the other way around) and n and s the corresponding sample sizes and estimations for standard deviations. As shown, calculations were performed with different variances and

Software

Data was prepared using basic R functions and the sqldf and dplyr packages. Figures were generated with ggplot2 and by the use of ggrepel and gridExtra. Additional steps in data acquisition were performed with Python.

Results: Data included

Screening for possibly eligible studies (Fig. 1) identified 173 in vitro, 306 PG and 144 DDI studies. Based on the predefined criteria, we could include 41 in vitro studies. The reaction phenotyping methods (Table 1) used most frequently were chemical inhibition (50%), the abundance factor approach (31%) and the relative activity factor approach (18%). All specific data from these studies is available in supplementary file A (https://data.mendeley.com/datasets/c7jhvd7g9p/1). Concerning the

Correlation between in vitro and clinical PG and DDI data

Fractional clearance via CYP2D6 calculated from the PM/EM ratios in clinical PG studies correlated very well with the predictions from chemical inhibition/substrate depletion in vitro assays (Fig. 2 A). Only one drug, sertraline, was far out of the 95% confidence interval for the correlation coefficient, possibly explained by an unspecific inhibition of other enzymes. Fractional clearance calculated from clinical DDI studies correlated equally well with the predictions from chemical

Correlation of clinical pharmacokinetic PG and DDI data

Drug metabolizing enzymes like CYP2D6 and CYP2C19 studied here may lack any activity due to genetic deficiency or due to enzyme-inhibitory comedication. The relationship between CYP2D6 genetic deficiency and enzyme inhibition could be analyzed for three inhibitors, fluoxetine, paroxetine and quinidine, which were used in several studies. As illustrated (Fig. 4), the effect of the 3 inhibitors tended to be less than the genotype effect, indicating that there was no complete CYP2D6 inhibition in

Discussion

An individually optimally adjusted drug dosing requires a quantitative understanding of the impact of single enzymes or membrane transporters on the pharmacokinetics in human beings. Utilizing the fact that many humans have complete inherited deficiency of CYP2D6 and CYP2C19 activities, clinical pharmacogenetic studies in poor metabolizers versus extensive metabolizers provide a unique and valid gold standard on the relative contribution of these two enzymes (i.e. the fractional clearance) for

Declaration of competing interest

The authors declare that there are no conflicts of interest. No external funding was obtained for this analysis.

References (26)

  • M. Ingelman-Sundberg et al.

    Polymorphic human cytochrome P450 enzymes: An opportunity for individualized drug treatment

    Trends in Pharmacological Sciences

    (1999)
  • K.L. Rost et al.

    Phenocopies of poor metabolizers of omeprazole caused by liver disease and drug treatment

    Journal of Hepatology

    (1995)
  • J. Yadav et al.

    Time-dependent enzyme inactivation: Numerical analyses of in vitro data and prediction of drug-drug interactions

    Pharmacology & Therapeutics

    (2020)
  • V. Bergk et al.

    Information deficits in the summary of product characteristics preclude an optimal management of drug interactions: A comparison with evidence from the literature

    European Journal of Clinical Pharmacology

    (2005)
  • L. Brunton et al.

    Goodman and Gilman's the pharmacological basis of therapeutics

    (2011)
  • FDA

    In vitro metabolism- and transporter- mediated drug-drug interaction studies: guidance for industry

    (2017)
  • A.M. Filppula et al.

    Improved predictions of time-dependent drug-drug interactions by determination of cytosolic drug concentrations

    Scientific Reports

    (2019)
  • J.P. Gibbs et al.

    Minimizing polymorphic metabolism in drug discovery: Evaluation of the utility of in vitro methods for predicting pharmacokinetic consequences associated with CYP2D6 metabolism

    Drug Metabolism and Disposition

    (2006)
  • Haefeli, W. E. (2019). DOSING tool for drug utilization. Retrieved from...
  • J.B. Houston et al.

    Progress towards prediction of human pharmacokinetic parameters from in vitro technologies

    Drug Metabolism Reviews

    (2003)
  • J. Kirchheiner et al.

    Pharmacogenetics-based therapeutic recommendations—Ready for clinical practice?

    Nature Reviews. Drug Discovery

    (2005)
  • J. Kirchheiner et al.

    Impact of the ultrarapid metabolizer genotype of cytochrome P450 2D6 on metoprolol pharmacokinetics and pharmacodynamics

    Clinical Pharmacology and Therapeutics

    (2004)
  • J. Kirchheiner et al.

    Contributions of CYP2D6, CYP2C9 and CYP2C19 to the biotransformation of E- and Z-doxepin in healthy volunteers

    Pharmacogenetics

    (2002)
  • Cited by (2)

    View full text