A semiphysiological population pharmacokinetic model for dynamic inhibition of liver and gut wall cytochrome P450 3A by voriconazole.

backgroundAccurate predictions of cytochrome P450 (CYP) 3A-mediated drug-drug interactions (DDIs) account for dynamic changes of CYP3A activity at both major expression sites (liver and gut wall) by considering the full pharmacokinetic profile of the perpetrator and the substrate. Physiological-based in vitro-in vivo extrapolation models have become of increasing interest. However, due to discrepancies between the predicted and observed magnitude of DDIs, the role of models fully based on in vivo data is still essential.

objectiveThe primary objective of this study was to develop a coupled dynamic model for the interaction of the CYP3A inhibitor voriconazole and the prototypical CYP3A substrate midazolam.

methodsRaw concentration data were obtained from a DDI study. Ten subjects were given either no pretreatment (control) or voriconazole twice daily orally. Midazolam was given either intravenously or orally after the last voriconazole dose and during control phases. Data analysis was performed by the population pharmacokinetic approach using non-linear mixed effects modelling (NONMEM 7.2.0). Model evaluation was performed using visual predictive checks and bootstrap analysis.

resultsA semiphysiological model was able to describe the pharmacokinetics of midazolam, its major metabolite and voriconazole simultaneously. By considering the temporal disposition of all three substances in the liver and gut wall, a time-varying CYP3A inhibition process was implemented. Only the incorporation of hypothetical enzyme site compartments resulted in an adequate fit, suggesting a sustained inhibitory effect through accumulation. Novel key features of this analysis are the identification of (1) an apparent sustained inhibitory effect by voriconazole due to a proposed quasi accumulation at the enzyme site, (2) a significantly reduced inhibitory potency of intravenous voriconazole for oral substrates, (3) voriconazole as a likely uridine diphosphate glucuronosyltransferase (UGT) 2B inhibitor and (4) considerable sources of interindividual variability.

conclusionThe proposed semiphysiological modelling approach generated a mechanistic description of the complex DDI occurring at major CYP3A expression sites and thus may serve as a powerful tool to maximise information acquired from clinical DDI studies. The model has been shown to draw precise and accurate predictions. Therefore, simulations based on this kind of models may be used for various clinical scenarios to improve pharmacotherapy.