ABSTRACT:Conventional methods to forecast CYP3A-mediated drug-drug interactions have not employed stochastic approaches that integrate pharmacokinetic (PK) variability and relevant covariates to predict inhibition in terms of probability and uncertainty. Empirical approaches to predict the extent of inhibition may not account for nonlinear or non-steady-state conditions, such as first-pass effects or accumulation of inhibitor concentration with multiple dosing. A physiologically based PK model was developed to predict the inhibition of CYP3A by ketoconazole (KTZ), using midazolam (MDZ) as the substrate. The model integrated PK models of MDZ and KTZ, in vitro inhibition kinetics of KTZ, and the variability and uncertainty associated with these parameters. This model predicted the time-and dose-dependent inhibitory effect of KTZ on MDZ oral clearance. The predictive performance of the model was validated using the results of five published KTZ-MDZ studies. The model improves the accuracy of predicting the inhibitory effect of increasing KTZ dosing on MDZ PK by incorporating a saturable KTZ efflux from the site of enzyme inhibition in the liver. The results of simulations using the model supported the KTZ dose of 400 mg once daily as the optimal regimen to achieve maximum inhibition by KTZ. Sensitivity analyses revealed that the most influential variable on the prediction of inhibition was the fractional clearance of MDZ mediated by CYP3A. The model may be used prospectively to improve the quantitative prediction of CYP3A inhibition and aid the optimization of study designs for CYP3A-mediated drug-drug interaction studies in drug development.Metabolism-based pharmacokinetic (PK) interactions are well recognized as a source of clinically significant adverse drug reactions Huang and Lesko, 2004). The early forecast of clinically significant metabolism-based pharmacokinetic drug-drug interactions is an increasingly important aspect of drug development. The assessment of drug-drug interaction potential before the conduct of a clinical trial often involves simple algebraic calculations. One such calculation is the ratio of the expected clinical exposure to the in vitro inhibition constant (K i ) of a new drug entity for a specific cytochrome P450 (Sahajwalla et al., 1999;Bjornsson et al., 2003). The primary sources of quantitative errors associated with this prediction approach are attributed to the experimental procedures used to determine the in vitro parameters, the variability in the intrinsic factors associated with PK properties of the inhibitor (e.g., absorption kinetics, plasma protein binding, and tissue partition coefficients), or the substrate (i.e., the fraction of total clearance attributed to the elimination pathway of interest), the extrinsic factors associated with the clinical study designs (e.g., dosing scheme, sampling design, and population demographics), and, most notably, the uncertainty in the effective concentration of the inhibitor at the enzyme site. In addition, complex drug disposition proper...