PROTOCOLS article

Front. Pharmacol., 27 December 2016

Sec. Drug Metabolism and Transport

Volume 7 - 2016 | https://doi.org/10.3389/fphar.2016.00517

Optimized Cocktail Phenotyping Study Protocol Using Physiological Based Pharmacokinetic Modeling and In silico Assessment of Metabolic Drug–Drug Interactions Involving Modafinil

  • 1. Department of Clinical Pharmacology, School of Medicine, Flinders University Adelaide, SA, Australia

  • 2. Precision Medicine Group, Flinders Center for Innovation in Cancer, School of Medicine, Flinders University AdelaideSA, Australia

Abstract

In vivo cocktail pathway phenotyping (ICPP) is routinely used to assess the metabolic drug–drug interaction (mDDI) potential of new drug candidates (NDC) during drug development. However, there are a number of potential limitations to this approach and the use of validated drug cocktails and study protocols is essential. Typically ICPP mDDI studies assess only the impact of interactions following multiple postulated perpetrator doses and hence the emphasis in terms of validation of these studies has been ensuring that there are no interactions between probe substrates. Studies assessing the comparative impact of single and multiple doses of the postulated perpetrator have the potential to provide richer information regarding both the clinical impact and mechanism of mDDIs. Using modafinil as a model compound, we sought to develop an optimized ICPP mDDI study protocol to evaluate the potential magnitude and clinical relevance of mDDIs using a physiologically based pharmacokinetic modeling approach.

Introduction

In vivo cocktail pathway phenotyping (ICPP) is a relatively new approach that facilitates the simultaneous, but independent, assessment of the activity of multiple metabolic pathways. ICPP is ideally suited for the assessment of the magnitude of in vivo metabolic drug–drug interactions (mDDIs) involving cytochrome P450 (CYP) enzymes (Snyder et al., 2014). Indeed, ICPP is routinely used to assess the mDDI potential of new drug candidates (NDC) during drug development. However, there are a number of potential limitations to this approach, such as the capacity for interactions between probe drugs, increased risk of adverse effects caused by probe drugs, and increased analytical complexity (Tanaka et al., 2003). As such, the use of validated drug cocktails and study protocols is essential.

Frye et al. (1997) the “Pittsburgh cocktail” was reported as the first validated in vivo cocktail for the assessment of CYP enzyme activities. Today, validated drug cocktails include the “Cooperstown” (Streetman et al., 2000a,b), “Karolinska” (Christensen et al., 2003), and “Inje” (Ryu et al., 2007) cocktails. Each cocktail has advantages and disadvantages, and in some cases they have been modified to expand their capacity (Chainuvati et al., 2003; Snyder et al., 2014) or address limitations. The four probe drugs in the Cooperstown cocktail; caffeine (CYP1A2), omeprazole (CYP2C19), dextromethorphan (CYP2D6), and intravenous midazolam (CYP3A), are generally considered the core ‘probes’ of the CYP cocktail approach. These drugs do not interact with each other and are not associated with significant adverse effects. However, the intravenous administration of midazolam (probe for CYP3A) precludes the assessment of gastrointestinal CYP3A enzyme activity. Furthermore, this cocktail does not contain a probe drug for CYP2C9, which accounts for approximately 18% of the CYP protein in human liver (Miners and Birkett, 1998). These limitations are partially ameliorated by using warfarin as a probe for CYP2C9 in the “Cooperstown 5+1” cocktail (Chainuvati et al., 2003). The Karolinska 5-drug cocktail uses caffeine (CYP1A2), losartan (CYP2C9), omeprazole (CYP2C19), debrisoquine (CYP2D6), and quinine (CYP3A) as probe drugs. However, the simultaneous administration of the cocktail drugs causes a significant increase of metabolic ratio of debrisoquine. Therefore, this approach requires the separate oral intake of debrisoquine, with additional logistical complexity. The “Inje” cocktail comprises caffeine (CYP1A2), losartan (CYP2C9), omeprazole (CYP2C19), dextromethorphan (CYP2D6) and oral midazolam (CYP3A). These probe drugs do not interact with each other and are not associated with significant adverse effects. Furthermore, the use of these drugs as probes facilitates the assessment of the effects of mDDIs on all major drug metabolizing CYP enzymes including both hepatic and gastrointestinal CYP3A.

Historically ICPP based mDDI studies have only assessed the impact of interactions following multiple (3 to 5) doses of the postulated ‘perpetrator,’ where an interaction perpetrator is a compound that has the potential to alter the metabolic clearance of other drugs. As such, the emphasis in terms of validation of these studies has been ensuring that there are no interactions between probe substrates. Increasingly, however, it is recognized that the mechanism of mDDIs both in terms of inhibition (i.e., mechanism based versus reversible) and activation (i.e., induction of enzyme expression versus enhanced substrate binding) of CYP activity can have important clinical implications. While multiple dose mDDI studies are useful in describing the effect of steady state dosing of the postulated perpetrator on CYP activity, they provide limited insights regarding the mechanism of mDDIs. Studies assessing the comparative impact of single and multiple doses of the postulated perpetrator, providing information both on the clinical impact and mechanism of mDDIs, have not been routinely performed and appropriate validated dosing protocols to facilitate such studies remain to be defined.

Modafinil, a drug able to both induce and inhibit CYP activities in vitro, particularly CYP3A (Robertson et al., 2000), is an ideal candidate to assess the comparative impact of single dose versus steady state dosing on CYP activity using an ICPP mDDI study design. We sought to develop an optimized ICPP mDDI study protocol to evaluate the potential magnitude and clinical relevance of mDDIs perpetrated by modafinil using a physiologically based pharmacokinetic (PBPK) modeling approach.

Methods

Structural Model

Simulations were performed using the Simcyp® Simulator (version 15.1; Jamei et al., 2009) with a ‘minimal PBPK model’ comprising a liver compartment and a merged compartment representing all other organs (Howgate et al., 2006; Polasek et al., 2009; Wattanachai et al., 2011). The differential equations used by the simulator describing enzyme kinetics and the impact of co-variates have been described previously (Rowland Yeo et al., 2010).

Population Model

Unless specified otherwise, simulations were performed using the Simcyp Healthy Volunteer population profile with virtual study cohorts comprising an equal distribution of healthy males and females aged 21–40 years old. Simulation data are presented as the geometric mean and 95% confidence interval (CI) for 120 estimations (10 studies of 12 individuals).

Development of Modafinil Model

Compound (Substrate and Inhibitor) Profile

The physicochemical, blood binding, absorption, distribution, elimination, and interaction data used to construct the modafinil compound (substrate and inhibitor) profiles are summarized in Table 1. A schematic depicting the minimal PBPK model used to simulate and predict the pharmacokinetics of modafinil is shown in Figure 1. All parameters were based on published literature values or were model predicted based on the physicochemical characteristics of the drug. Modafinil hepatic microsomal intrinsic clearance (CLint) was back-calculated from clinically observed oral clearances (CLPO; Wong et al., 1998a) using the retrograde model function in Simcyp (Yamazaki et al., 2015). The interaction (‘inhibitor’) profile was created based on published in vitro microsomal inhibition and hepatocyte induction data for CYP 1A2, 2C9, 2C19, 2D6 and 3A4 (Robertson et al., 2000). The presence of endogenous fatty acids in microsomal preparations is known to effect the in vitro determination of kinetic parameters for multiple CYP enzymes including CYP 1A2, 2C9 and 2C19 (Rowland et al., 2008; Wattanachai et al., 2012, 2015). As such, published Ki values for CYP 1A2, 2C9, and 2C19 were scaled based on the known effects of endogenous fatty acids on these enzymes.

Table 1

Physiochemical properties
Molecular weight273.4
log Po:w1.53
pKa (monoprotic base)8.84
B/P (SimCYP predicted)0.887
fup0.4 (Wong et al., 1998a)

In vivo pharmacokinetic properties

fa (Simcyp predicted)0.99
ka (1/h; Simcyp predicted)3.93
Qgut (L/h; Simcyp predicted)15.56
Vss (L/kg; Simcyp predicted)0.72
CLpo (L/h)3.23 (Wong et al., 1998a)

In vitro inhibition parameters (Ki;μM; Robertson et al., 2000)

CYP1A2750
CYP2C9750
CYP2C197.8
CYP2D61,500
CYP3A4632

In vitro induction parameters (IndC50; μM/IndMax; Robertson et al., 2000)

CYP1A275/1.9
CYP2C9N/A
CYP2C19N/A
CYP2D6N/A
CYP3A410/2.5

Substrate and inhibitor parameter values used for modafinil substrate profile.

Po:w, neutral species octanol: buffer partition coefficient; B/P, blood-to-plasma partition ratio; fup, fraction unbound in plasma; fa, fraction available from dosage form; ka, first-order absorption rate constant; Vss, volume of distribution at steady state; CLpo, oral clearance; Ki, inhibition constant; IndC50, inducer concentration that causes half maximal induction; IndMax, maximal fold induction.

FIGURE 1

Validation of Substrate Profile

The modafinil substrate profile (simulated exposure) was validated by comparison of the pharmacokinetics parameters maximal plasma concentration (CMax), time to maximal plasma concentration (tMax), area under the plasma concentration time curve (AUC) and elimination half-life (t1/2) from simulations with data from age and gender matched participants from three clinical studies (n = 6–12) that were not used in the development of the modafinil compound profile (Wong et al., 1998a, 1999b; Robertson and Hellriegel, 2003).

Stepwise Procedure

Compound Profiles

Simulations performed to assess CYP 1A2, 2D6, 3A4 and 2C19 activities used validated Simcyp substrate profiles for caffeine, dextromethorphan, midazolam and omeprazole, respectively (Jamei et al., 2009). Simulations performed to assess CYP2C9 activity used a losartan substrate profile based on published parameters (Brantley et al., 2014). Probe doses were 100 mg orally (PO) of caffeine, 25 mg PO of losartan, 20 mg PO of enteric coated omeprazole, 30 mg PO of dextromethorphan and 1 mg PO of midazolam.

Validation of Study Dosing Protocol

The study dosing protocol was validated using PBPK modeling in terms of the minimal washout period between probe doses to facilitate assessment of CYP activity under three conditions (baseline, single modafinil dose and steady state modafinil dosing) and the minimal number of modafinil doses to achieve steady state. Following a series of preliminary simulations, two dosing protocols were assessed in full; a “Day 0, 1, 7 protocol,” where probes were dosed at 0, 24, and 168 h, with assessment of exposure at 200 uniformly distributed sampling time points over 192 h, and a “Day 0, 2, 8 protocol,” where probes were dosed at 0, 48, and 192 h, with assessment of exposure at 200 uniformly distributed sampling time points over 216 h (Figure 2). Residual concentration at 0.5 h prior to the second dose, AUC and CMax for each probe in the absence of modafinil dosing were compared between the protocols. Trough modafinil concentrations (0.5 h prior to the subsequent dose) were evaluated over 7 days to determine the minimal dosing duration required to achieve repeat trough concentrations within 5%.

FIGURE 2

Simulated Modafinil ICPP mDDI Study

The simulated impact of a single dose of modafinil (200 mg; PO) administered 1 h prior to the second probe dose and steady state dosing of modafinil [200 mg PO once daily (QD) for 7 days] on caffeine, losartan, omeprazole, dextromethorphan and midazolam exposure was assessed using the Day 0, 2, 8 protocol. The capacity for modafinil to perpetrate mDDIs was predicted on the basis of the simulated probe AUC ratios, with and without modafinil (single dose or steady state). The geometric mean of the AUC ratio for each probe was estimated using a mixed effects model of logarithmically transformed data. Time period was included as a fixed effect and participant as a random effect. Back transformation was utilized to provide a point estimate and 95% CI for the AUC ratio. Based on an equivalence approach a lack of effect of modafinil on CYP activity was demonstrated if the 95% CI for the estimated AUC ratio for the probe was contained within the range 0.85 to 1.2.

Anticipated (Simulation) Results

Optimisation and Validation of ICPP mDDI Model

Modafinil Substrate Profile

Comparison of simulated and observed pharmacokinetic parameters reported in Table 2 demonstrate that the modafinil substrate profile accurately estimates exposure to this drug following a single dose (Wong et al., 1999a; Robertson and Hellriegel, 2003) and steady state dosing (Wong et al., 1999b). Simulated and observed (Wong et al., 1999b) plasma concentration-time profiles for modafinil (200 mg PO QD) over 7 days are shown in Figure 3. Estimated mean modafinil trough concentrations were 1.00, 1.30, 1.45, 1.47, 1.52, and 1.53 mg/L on days 2, 3, 4, 5, 6, and 7, respectively. Repeatable trough concentrations within 5% of the highest trough concentration were obtained from day 5 onwards. Simulations confirm that steady state exposure to modafinil is achieved within the 7 days perpetrator dosing period recommended by the Food and Drug Administration (FDA) for the assessment of induction and mechanism based inhibition mDDIs (FDA, 2012).

Table 2

Population (age range)DoseStudyCMax (mg/L)tMax (h)AUC (mg/L h)t1/2 (h)
Males (22–37 y/o; Wong et al., 1999a)200 mgObserved4.21.55712
Simulation3.0 (2.8–3.1)1.4 (1.3–1.4)59 (55–62)16 (14–17)
Females (19–40 y/o; Wong et al., 1999a)200 mgObserved5.21.56110
Simulation3.3 (1.3–1.4)1.4 (1.3–1.4)60 (57–64)14 (13–15)
Males (20–39 y/o; Wong et al., 1999b)200 mg QD for7 daysObserved6.42.57917
Simulation4.4 (4.2–4.7)1.3 (1.2–1.4)77 (72–83)17 (16 –19)

Comparison of modafinil pharmacokinetics in age and gender matched simulations (mean and 95 % CI) and in vivo clinical studies.

AUC, area under the plasma concentration-time curve; CMax, maximal plasma concentration; tMax, time taken to achieve maximal plasma concentration; t1/2, elimination half-life.

FIGURE 3

Study Dosing Protocol

Baseline and second dose mean plasma concentration-time profiles for each probe obtained using the Day 0, 1, 7 and Day 0, 2, 8 protocols are shown in Figure 4. Corresponding residual concentration, AUC and CMax values are described in Table 3. Simulations demonstrated significant residual caffeine and dextromethorphan concentrations prior to the second dose using the Day 0, 1, 7 protocols of 186 and 2.3 μg/L, respectively. Consistent with this observation, comparison of exposure (AUC and CMax) for probes dosed on Day 0 and Day 1 in the absence of modafinil, demonstrate that the second dose (Day 1) mean AUC and CMax for caffeine were increased by 9.9 and 9.8 %, respectively, while the second dose (Day 1) mean AUC and CMax for dextromethorphan were increased by 18.1 and 17%, respectively. Simulations performed using the Day 0, 2, 8 protocol demonstrated that estimated second dose (Day 2) mean AUC and CMax values for all probes were within 5% of the corresponding first dose (Day 0) value. Mean AUC and CMax values for probes administered 7 days after the second dose (i.e., Day 7 or 8) were also invariably within 5% of the first dose (Day 0) value with residual baseline probe concentrations of 0 μg/L (not shown).

Table 3

ProbeProtocol[Probe]Residual (μg/L)
AUC (μg/L h)
CMax (μg/L)
23.5 h47.5 hFirst dose24–72 h48–96 hFirst doseDay 2Day 3
CaffeineDay 0, 1, 7186 (143–229)16,97018,653#(16,518–21,083)2,128 (1,986–2,282)2,337#(2,192–2,491)
Day 0, 2, 834 (26–45)(14,868–18,961)16,773 (14,858–18,936)2,145
DextromethorphanDay 0, 1, 72.3 (1.8–2.9)72 (59–88)85# (72–101)9.9 (8.7–10.6)11.6# (10.7–12.9)(2,003–2,299)
Day 0, 2, 81.2 (0.8–1.6)74 (61–89)10.4 (9.3–12.4)
LosartanDay 0, 1, 70.4 (0.0–0.9)575 (529–624)576 (530–625)165 (152–178)165 (152–178)
Day 0, 2, 80.0 (0.0–0.2)575 (529–624)167 (154–180)
MidazolamDay 0, 1, 70.1 (0.0–0.2)11.4 (10.1–13.0)11.6 (10.2–13.3)4.2 (3.8–4.7)4.5 (4.1–4.9)
Day 0, 2, 80.0 (0.0–0.0)11.4 (10.1–13.0)4.5 (4.1–4.9)
OmeprazoleDay 0, 1, 70.0 (0.0–0.0)391 (341–449)391 (341–449)420 (376–465)420 (376–465)
Day 0, 2, 80.0 (0.0–0.0)391 (341–449)420 (376–465)

Comparison of mean (and 95% CI) pharmacokinetic parameters defining probe exposure between Day 0, 1 and 7 and Day 0, 2 and 8 dosing protocols.

# Mean parameter estimate 10% greater than corresponding first dose value.

FIGURE 4

Assessment of Modafinil as an mDDI Perpetrator

Simulated mean probe AUC values and ratios with 95% CIs in the presence and absence of a single dose of modafinil and steady state dosing of modafinil are shown in Tables 4 and 5, respectively. Following a single oral dose of modafinil, the mean AUC omeprazole increased by 56% (391 to 610 μg/L h). However, the magnitude of this predicted mDDI was partially attenuated following dosing of modafinil to steady state; 38% increase in omeprazole AUC from 391 to 455 μg/L h (Figure 5). These data indicate that modafinil may perpetrate clinically relevant inhibitory mDDIs when co-administered with drugs primarily metabolized by CYP2C19. The partially attenuated interaction simulated here following dosing of modafinil to steady state is likely due to an increase in CYP3A4 catalyzed omeprazole metabolism (see Discussion). Dosing of modafinil to steady state (7 days) resulted in a 50% decrease in midazolam mean AUC (11.2 to 5.6 μg/L h; Figure 6). These data indicate that modafinil causes clinically relevant induction of CYP3A4 with steady state dosing. No change in midazolam AUC was observed following a single dose of modafinil. Similarly, no significant effect of modafinil on CYP1A2 (caffeine), CYP2C9 (losartan) or CYP2D6 (dextromethorphan) was observed following either a single dose or dosing to steady state (probe drug Δ AUC < 5 %).

Table 4

EnzymeSubstrate (dose)AUC (μg/L h)
Without modafinilWith modafinilRatio
CYP1A2Caffeine (100 mg)16,787 (14,867–18,956)16,713 (14,806–18,865)1.00 (0.99–1.00)
CYP2C9Losartan (25 mg)575 (529–624)571 (526–620)0.99 (0.99–1.00)
CYP2C19Omeprazole (20 mg)391 (341–449)610 (537–691)1.56 (1.52–1.59)
CYP2D6Dextromethorphan (30 mg)71.6 (58.8–87.2)71.5 (58.9–86.8)1.00 (0.99–1.00)
CYP3A4Midazolam (1 mg)11.4 (10.5–13.0)10.9 (9.6–12.4)0.96 (0.95–0.96)

Area under the plasma-concentration time curve (mean and 95 % CI) for probes in the absence and presence of a single dose of modafinil (200 mg PO).

Table 5

EnzymeSubstrate (dose)AUC (μg/L h)
Without modafinilWith modafinilRatio
CYP1A2Caffeine (100 mg)16,556 (14710–18633)15,941 (14166–17939)0.96 (0.96–0.97)
CYP2C9Losartan (25 mg)573 (528–623)526 (485–571)0.92 (0.91–0.93)
CYP2C19Omeprazole (20 mg)391 (341–449)455 (413–501)1.38 (1.32–1.44)
CYP2D6Dextromethorphan (30 mg)70.2 (58.2–84.7)68.5 (57.2–82.1)0.98 (0.97–0.98)
CYP3A4Midazolam (5 mg)11.2 (9.9–12.8)5.6 (4.9–6.4)0.50 (0.47–0.53)

Area under the plasma-concentration time curve (mean and 95% CI) for probes in the absence and presence of modafinil (200 mg PO QD) dosed to steady state (7 days).

FIGURE 5

FIGURE 6

Discussion

This study reports a validated dosing protocol for the administration of Inje cocktail probes that facilitates the assessment of the magnitude of potential mDDIs following a single perpetrator dose and steady state perpetrator dosing. Consistent with the estimated half-life of modafinil in healthy volunteers (∼15 h; Robertson and Hellriegel, 2003), simulated steady-state concentrations were obtained with 5 days of dosing. Repeatable simulated mean trough concentrations (defined as within 5%) were observed from Day 5 onwards. The minimal duration of perpetrator dosing considered sufficient to facilitate assessment of potential induction and mechanism-based inactivation of CYP enzymes e.g., time-dependent changes in CYP activity is 7 days (FDA, 2012). Simulations demonstrate that this dosing duration was also sufficient to ensure steady state modafinil concentrations.

In order to assess changes in CYP activity caused by the administration of the postulated mDDI perpetrator (single or steady state dosing), ICPP based mDDI studies require administration of probes prior to perpetrator administration to establish baseline CYP activities (Day 0). Simulations demonstrate that when assessing the impact of a single dose of the postulated mDDI perpetrator, repeat probe dosing on Day 1 (i.e., the Day 0, 1 and 7 protocol) is suitable for evaluation of CYP 2C9 (losartan), 2C19 (omeprazole) and 3A4 (midazolam) activities since no residual probe is present and AUC and CMax values are equivalent to the first dose. However, this protocol is not suitable for assessment of CYP 1A2 (caffeine) or 2D6 (dextromethorphan) activities due to the presence of significant residual probe in the systemic circulation at the time of the second dose (24 h), leading to an overestimation of AUC and CMax for these probes (Table 3). Use of this protocol will bias activity assessment toward inhibitory mDDIs involving CYP 1A2 and 2D6 increasing the risk of either detecting an inhibitory interaction (false positive) or failing to detect induction (false negative).

Simulations demonstrate that when assessing the single dose effect of the postulated mDDI perpetrator, administration of probes 48 h after the initial dose (Day 0, 2 and 8 protocol) is suitable for evaluation of all CYP activities since essentially no residual probe is present and AUC and CMax values are comparable between the first and second dose (Table 3). On this basis, the Day 0, 2 and 8 protocol was used for the assessment of modafinil as a perpetrator of mDDIs involving CYP and is recommended for ICPP mDDI studies performed to assess single dose and steady state dosing effects of the a postulated mDDI perpetrator.

The capacity of modafinil to inhibit and induce CYP activity has been demonstrated in vitro. Classical (i.e., [I]/Ki) extrapolation of in vitro data suggests that modafinil may induce expression of CYP 1A2 and 3A4 and inhibit CYP 2C19 and 3A4 activities in a reversible, competitive manner (Robertson et al., 2000). On the basis of these classically extrapolated in vitro data, modafinil is classified as a ‘moderate’ inducer of CYP 3A4 by the US FDA. Increasingly, however, it is accepted that there are several limitations (such as the appropriate selection of [I] (Rowland et al., 2006)) associated with the use of classical extrapolation of in vitro data when assessing mDDIs that preclude sensible assessment of the clinical relevance of novel interactions. Indeed, the substantial benefits of mechanistic (PBPK) extrapolation of mDDI assessments have been extensively reported and adopted (Sheiner and Steimer, 2000; Danhof et al., 2008; Rowland Yeo et al., 2010; Rowland et al., 2011) and this approach is recommended for the preclinical assessment of mDDIs for regulatory review (Zhao et al., 2011; FDA, 2012).

On the basis of the in vitro interaction profile (warfarin; CYP2C9; Robertson et al., 2002a), or potential for concomitant use with other vigilance (dexamphetamine, methylphenidate) or sedative (triazolam) agents (Wong et al., 1998b; Hellriegel et al., 2001, 2002; Robertson et al., 2002b), a limited number of in vivo studies assessing the mDDI potential of modafinil have been reported. Notably, these study ‘victim’ drugs, including warfarin, which was administered as a racemic mix or r- and s-enantiomers (metabolised by CYP 3A4 and 2C9, respectively) have been substrates for multiple metabolic pathways, and as such these studies do not readily facilitate the direct assessment of the effects of modafinil on individual CYP enzymes. Only one in vivo interaction study has reported the capacity of modafinil to induce drug metabolizing CYP (Moachon et al., 1996). This study used antipyrine, a pan-CYP substrate, with 7 days of modafinil dosing at 100 to 500 mg BD. The results of the study suggest that modafinil may be a weak general inducer of CYP at doses ≥400 mg/day (double the recommended daily dose), but no definitive conclusions were possible. The study design did not facilitate assessment of effects on individual CYP enzymes or account for concurrent inhibition of these enzymes.

In the current study, a clinically relevant increase in omeprazole exposure was predicted following both a single modafinil dose and dosing of modafinil to steady state. Notably, the apparent magnitude of the inhibitory interaction was attenuated with dosing of modafinil to steady state; omeprazole AUC ratios following a single modafinil dose and steady state modafinil dosing were 1.56 (1.52–1.59) and 1.38 (1.32–1.44), respectively. While omeprazole is primarily (88%) cleared by CYP2C19 catalyzed demethylation and hydroxylation, the remainder of omeprazole clearance occurs via CYP3A4 catalyzed sulfone formation. Following repeated dosing of modafinil, which induces CYP3A4 while inhibiting CYP2C19, the inhibitory effect on CYP2C19 is partially offset by induction of CYP3A4. The resulting change in fraction metabolized (fm) for omeprazole is shown in Figure 7. These data highlight both the advantages of mechanistic extrapolation of mDDI data when assessing complex interactions involving victim drugs that are metabolized by multiple enzymes (Rowland Yeo et al., 2010), and the need to consider even minor metabolic pathways when considering the selection of probe substrates as induction of a minor metabolic pathway can result in a significant contribution to clearance when induced. Where multiple enzyme involvement results in the formation of unique metabolites for each enzyme, this potential confounder may be overcome by the consideration of an enzyme specific metabolic ratio. Simulations also predict that steady state dosing of modafinil may cause clinically relevant mDDIs resulting in reduced exposure to drugs metabolized by CYP3A4. This observation supports the current FDA classification of modafinil as a ‘moderate inducer’ of this enzyme (FDA, 2012).

FIGURE 7

Conclusion

These data support consideration of the risk of clinically relevant mDDIs when co-administered modafinil with drugs that are primarily cleared by CYP 2C19 and 3A4 catalyzed metabolic pathways. In particular, given the major role of CYP3A4 in the clearance and metabolic activation of a myriad of drugs from various therapeutic classes including the oral contraceptive pill and orally administered non-cytotoxic anticancer drugs (e.g., erlotinib, pazopanib and sunitinib; Rowland et al., 2016) potential induction of CYP3A4 by modafinil should be considered when evaluating new indications for this drug such as the management of chemotherapy induced fatigue and cognitive impairment in order to avoid a potential decrease in therapeutic response for the co-prescribed drug that may result in therapeutic failure.

Statements

Author contributions

Participated in research design: AgR, AM, MS, and AdR. Conducted experiments: AgR and AdR. Performed data analysis: AgR, AM and AdR. Wrote or contributed to the writing of the manuscript: AgR, AM, AH, MS, and AdR.

Acknowledgments

This work was supported by seed funding from the Flinders Foundation (Grant ID 3000007902).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  • 1

    BrantleyS. J.GuffordB. T.DuaR.FediukD. J.GrafT. N.ScarlettY. V.et al (2014). Physiologically based pharmacokinetic modeling framework for quantitative prediction of an herb-drug interaction.CPT Pharmacometrics Syst. Pharmacol.3:e107. 10.1038/psp.2013.69

  • 2

    ChainuvatiS.NafzigerA. N.LeederJ. S.GaedigkA.KearnsG. L.SellersE.et al (2003). Combined phenotypic assessment of cytochrome P450 1A2, 2C9, 2C19, 2D6, and 3A. N-acetyltransferase-2, and xanthine oxidase activities with the cooperstown 5+1 cocktail.Clin. Pharmacol. Ther.74437447.

  • 3

    ChristensenM.AnderssonK.DalénP.MirghaniR. A.MuirheadG. J.NordmarkA.et al (2003). The Karolinska cocktail for phenotyping of five human cytochrome P450 enzymes.Clin. Pharmacol. Ther.73517528. 10.1016/S0009-9236(03)00050-X

  • 4

    DanhofM.de LangeE. C.Della PasquaO. E.PloegerB. A.VoskuylR. A. (2008). Mechanism-based pharmacokinetic-pharmacodynamic (PK-PD) modeling in translational drug research.Trends Pharmacol. Sci.29186191. 10.1016/j.tips.2008.01.007

  • 5

    FDA (2012). Draft Guidance/Guidance for Industry in Drug Interaction Studies–Study Design, Data Analysis, Implications for Dosing, and Labeling Recommendations.Silver Spring, MA: FDA.

  • 6

    FryeR. F.MatzkeG. R.AdedoyinA.PorterJ. A.BranchR. A. (1997). Validation of the five-drug “Pittsburgh cocktail” approach for assessment of selective regulation of drug-metabolizing enzymes.Clin. Pharmacol. Ther.62365376. 10.1016/S0009-9236(97)90114-4

  • 7

    HellriegelE. T.AroraS.NelsonM.RobertsonP.Jr. (2001). Steady-state pharmacokinetics and tolerability of modafinil given alone or in combination with methylphenidate in healthy volunteers.J. Clin. Pharmacol.41895904. 10.1177/00912700122010690

  • 8

    HellriegelE. T.AroraS.NelsonM.RobertsonP.Jr. (2002). Steady-state pharmacokinetics and tolerability of modafinil administered alone or in combination with dextroamphetamine in healthy volunteers.J. Clin. Pharmacol.42450460. 10.1177/00912700222011382

  • 9

    HowgateE. M.Rowland YeoK.ProctorN. J.TuckerG. T.Rostami-HodjeganA. (2006). Prediction of in vivo drug clearance from in vitro data. I: impact of inter-individual variability.Xenobiotica36473497.

  • 10

    JameiM.MarciniakS.FengK.BarnettmA.TuckerG.Rostami-HodjeganA. (2009). The simcyp® population-based ADME simulator.Exp. Opin. Drug Metab. Toxicol.5211223. 10.1517/17425250802691074

  • 11

    MinersJ. O.BirkettD. J. (1998). Cytochrome P4502C9: an enzyme of major importance in human drug metabolism.Br. J. Clin. Pharmacol.45525538. 10.1046/j.1365-2125.1998.00721.x

  • 12

    MoachonG.KanmacherI.ClenetM. (1996). Pharmacokinetic profile of modafinil.Drugs Today32327337.

  • 13

    PolasekT. M.PolakS.DoogueM. P.Rostami-HodjeganA.MinersJ. O. (2009). Assessment of inter-individual variability in predicted phenytoin clearance.Eur. J. Clin. Pharmacol.6512031210. 10.1007/s00228-009-0703-y

  • 14

    RobertsonP.Jr.HellriegelE. T. (2003). Clinical pharmacokinetic profile of modafinil.Clin. Pharmacokinet.42123137. 10.2165/00003088-200342020-00002

  • 15

    RobertsonP.DeCoryH. H.MadanA.ParkinsonA. (2000). In vitro inhibition and induction of human hepatic cytochrome P450 enzymes by modafinil.Drug Metab. Dispos.28664671.

  • 16

    RobertsonP.Jr.HellriegelE. T.AroraS.NelsonM. (2002a). Effect of modafinil at steady state on the single-dose pharmacokinetic profile of warfarin in healthy volunteers.J. Clin. Pharmacol.42205214. 10.1177/00912700222011120

  • 17

    RobertsonP.Jr.HellriegelE. T.AroraS.NelsonM. (2002b). Effect of modafinil on the pharmacokinetics of ethinyl estradiol and triazolam in healthy volunteers.Clin. Pharmacol. Ther.714656. 10.1067/mcp.2002.121217

  • 18

    RowlandA.ElliotD.WilliamsJ.MackenzieP. I.DickinsonR. G.MinersJ. O. (2006). In vitro characterization of lamotrigine N2-glucuronidation and the lamotrigine - valproic acid interaction.Drug Metab. Dispos.3410551062.

  • 19

    RowlandA.ElliotD. J.KnightsK. M.MackenzieP. I.MinersJ. O. (2008). The “albumin effect” and in vitro-in vivo extrapolation: sequestration of long-chain unsaturated fatty acids enhances phenytoin hydroxylation by human liver microsomal and recombinant cytochrome P450 2C9.Drug Metab. Dispos.36870877. 10.1124/dmd.107.019885

  • 20

    RowlandA.van DykM.MangoniA. A.MinersJ. O.McKinnonR. A.WieseM. D.et al (2016). Kinase inhibitor pharmacokinetics: comprehensive summary and roadmap for addressing inter-individual variability in exposure.Exp. Opin. Drug Metab. Toxicol.133149. 10.1080/17425255.2016.1229303

  • 21

    RowlandM.PeckC.TuckerG. (2011). Physiologically-based pharmacokinetics in drug development and regulatory science.Ann. Rev. Pharmacol. Toxicol.514573. 10.1146/annurev-pharmtox-010510-100540

  • 22

    Rowland YeoK.JameiM.YangJ.TuckerG. T.Rostami-HodjeganA. (2010). Physiologically based mechanistic modelling to predict complex drug-drug interactions involving simultaneous competitive and time-dependent enzyme inhibition by parent compound and its metabolite in both liver and gut - the effect of diltiazem on the time-course of exposure to triazolam.Eur. J. Pharm. Sci.39298309. 10.1016/j.ejps.2009.12.002

  • 23

    RyuJ. Y.SongI. S.SunwooY. E.ShonJ. H.LiuK. H.ChaI. J.et al (2007). Development of the “Inje cocktail” for high-throughput evaluation of five human cytochrome P450 isoforms in vivo.Clin. Pharmacol. Ther.82531540. 10.1038/sj.clpt.6100187

  • 24

    SheinerL. B.SteimerJ.-L. (2000). Pharmacokinetic/pharmacodynamic modeling in drug development.Ann. Rev. Pharmacol. Toxicol.406795. 10.1146/annurev.pharmtox.40.1.67

  • 25

    SnyderB. D.RowlandA.PolasekT. M.MinersJ. O.DoogueM. P. (2014). Evaluation of felodipine as a potential perpetrator of pharmacokinetic drug-drug interactions.Eur. J. Clin. Pharmacol.7011151122. 10.1007/s00228-014-1716-8

  • 26

    StreetmanD. S.BertinoJ. S.Jr.NafzigerA. N. (2000a). Phenotyping of drug-metabolizing enzymes in adults: a review of in-vivo cytochrome P450 phenotyping probes.Pharmacogenetics10187216. 10.1097/00008571-200004000-00001

  • 27

    StreetmanD. S.BleakleyJ. F.KimJ. S.NafzigerA. N.LeederJ. S.GaedigkA.et al (2000b). Combined phenotypic assessment of CYP1A2, CYP2C19, CYP2D6, CYP3A, N-acetyltransferase-2, and xanthine oxidase with the “Cooperstown cocktail”.Clin. Pharmacol. Ther.68375383.

  • 28

    TanakaE.KurataN.YasuharaH. (2003). How useful is the “cocktail approach” for evaluating human hepatic drug metabolizing capacity using cytochrome P450 phenotyping probes in vivo?J. Clin. Pharm. Ther.28157165. 10.1046/j.1365-2710.2003.00486.x

  • 29

    WattanachaiN.ElliotD. J.UchaipichatV.TassaneeyakulW.MinersJ. O. (2015). Effect of albumin on human cytochromes P450 kinetics: extrapolation of in vivo clearance from in vitro data.Thai J. Pharmacol.32:4.

  • 30

    WattanachaiN.PolasekT. M.HeathT. M.UchaipichatV.TassaneeyakulW.TassaneeyakulW.et al (2011). In vitro-in vivo extrapolation of CYP2C8-catalyzed paclitaxel 6alpha-hydroxylation: effects of albumin on in vitro kinetic parameters and assessment of interindividual variability in predicted clearance.Eur. J. Clin. Pharmacol.67815824. 10.1007/s00228-011-1001-z

  • 31

    WattanachaiN.TassaneeyakulW.RowlandA.ElliotD. J.BowalgahaK.KnightsK. M.et al (2012). Effect of albumin on human liver microsomal and recombinant CYP1A2 activities: impact on in vitro-in vivo extrapolation of drug clearance.Drug Metab. Dispos.40982989. 10.1124/dmd.111.044057

  • 32

    WongY. N.KingS. P.LaughtonW. B.McCormickG. C.GrebowP. E. (1998a). Single-dose pharmacokinetics of modafinil and methylphenidate given alone or in combination in healthy male volunteers.J. Clin. Pharmacol.38276282. 10.1002/j.1552-4604.1998.tb04395.x

  • 33

    WongY. N.KingS. P.SimcoeD.GormanS.LaughtonW.McCormickG. C.et al (1999a). Open-label, single-dose pharmacokinetic study of modafinil tablets: influence of age and gender in normal subjects.J. Clin. Pharmacol.39281288.

  • 34

    WongY. N.SimcoeD.HartmanL. N.LaughtonW. B.KingS. P.McCormickG. C.et al (1999b). A double-blind, placebo-controlled, ascending-dose evaluation of the pharmacokinetics and tolerability of modafinil tablets in healthy male volunteers.J. Clin. Pharmacol.393040. 10.1177/00912709922007534

  • 35

    WongY. N.WangL.HartmanL.SimcoeD.ChenY.LaughtonW.et al (1998b). Comparison of the single-dose pharmacokinetics and tolerability of modafinil and dextroamphetamine administered alone or in combination in healthy male volunteers.J. Clin. Pharmacol.38971978. 10.1002/j.1552-4604.1998.tb04395.x

  • 36

    YamazakiS.JohnsonT. R.SmithB. J. (2015). Prediction of drug-drug interactions with crizotinib as the CYP3A substrate using a physiologically based pharmacokinetic model.Drug Metab. Dispos.4314171429. 10.1124/dmd.115.064618

  • 37

    ZhaoP.ZhangL.GrilloJ. A.LiuQ.BullockJ. M.MoonY. J.et al (2011). Applications of physiologically based pharmacokinetic (PBPK) modeling and simulation during regulatory review.Clin. Pharmacol. Ther.89259267. 10.1038/clpt.2010.298

Summary

Keywords

physiological based pharmacokinetic modeling, cocktail phenotyping, metabolic drug–drug interactions, modafinil, study protocol

Citation

Rowland A, Mangoni AA, Hopkins A, Sorich MJ and Rowland A (2016) Optimized Cocktail Phenotyping Study Protocol Using Physiological Based Pharmacokinetic Modeling and In silico Assessment of Metabolic Drug–Drug Interactions Involving Modafinil. Front. Pharmacol. 7:517. doi: 10.3389/fphar.2016.00517

Received

11 November 2016

Accepted

14 December 2016

Published

27 December 2016

Volume

7 - 2016

Edited by

David Stec, University of Mississippi Medical Center, USA

Reviewed by

Takato Hiranita, National Center for Toxicological Research (FDA), USA; Wayne Louis Backes, LSU Health Sciences Center New Orleans, USA

Updates

Copyright

*Correspondence: Andrew Rowland,

This article was submitted to Drug Metabolism and Transport, a section of the journal Frontiers in Pharmacology

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics