A Continuous Intestinal Absorption Model to Predict Drug Enterohepatic Recirculation in Healthy Humans: Nalbuphine as a Model Substrate.

Nalbuphine (NAL) is a κ-agonist/µ-antagonist opioid being developed as an oral extended formulation (ER) for the treatment of chronic cough in idiopathic pulmonary fibrosis and itch in prurigo nodularis. NAL is extensively glucuronidated and likely undergoes enterohepatic recirculation (EHR). The purpose of this work is to develop pharmacokinetic models for NAL absorption and enterohepatic recirculation (EHR). Clinical pharmacokinetic (PK) data sets in healthy subjects from three trials that included IV, oral solution, and ER tablets in fed and fasted state and two published trials were used to parametrize a novel partial differential equation (PDE)-based model, termed "PDE-EHR" model. Experimental inputs included in vitro dissolution and permeability data. The model incorporates a continuous intestinal absorption framework, explicit liver and gall bladder compartments, and compartments for systemic drug disposition. The model was fully PDE-based with well-stirred compartments achieved by rapid diffusion. The PDE-EHR model accurately reproduces NAL concentration-time profiles for all clinical data sets. NAL disposition simulations required inclusion of both parent and glucuronide recirculation. Inclusion of intestinal P-glycoprotein efflux in the simulations suggests that NAL is not expected to be a victim or perpetrator of P-glycoprotein-mediated drug interactions. The PDE-EHR model is a novel tool to predict EHR and food/formulation effects on drug PK. The results strongly suggest that even intravenous dosing studies be conducted in fasted subjects when EHR is suspected. The modeling effort is expected to aid in improved prediction of dosing regimens and drug disposition in patient populations.

Process and System Clearances in Pharmacokinetic Models: Our Basic Clearance Concepts Are Correct.

Clearances are important parameters in pharmacokinetic (PK) models. All clearances in PK models are either process clearances that include diffusion, transport, and metabolism clearances or system clearances that include organ and systemic clearance. Clearance and volume of distribution are two independent parameters that characterize drug disposition in both individual compartments and systems of compartments. In this minireview, we show that systemic and organ clearances are net clearances that can be easily derived by partition analysis. When drugs are eliminated from the central compartment by first-order processes, systemic clearance is constant. When drugs are eliminated from a peripheral compartment, instantaneous systemic clearance will vary with time. However, average clearance and clearance at steady state will be constant and will equal dose divided by area under the curve. We show that peripheral elimination will not have a large impact on most pharmacokinetic analyses and that standard models of organ and systemic clearance are useful and appropriate. SIGNIFICANCE STATEMENT: There are two basic kinds of clearances used in pharmacokinetic models, process and system clearances. We show that organ and systemic clearances are net clearances with blood or plasma as the driving concentration. For linear pharmacokinetics, clearance is constant for elimination from the central compartment but varies with time for peripheral elimination. Despite the different kinds of clearance parameters and models, standard clearance models and concepts remain valid.

The Rat Continuous Intestine Model Predicts the Impact of Particle Size and Transporters on the Oral Absorption of Glyburide.

Oral drug absorption is known to be impacted by the physicochemical properties of drugs, properties of oral formulations, and physiological characteristics of the intestine. The goal of the present study was to develop a mathematical model to predict the impact of particle size, feeding time, and intestinal transporter activity on oral absorption. A previously published rat continuous intestine absorption model was extended for solid drug absorption. The impact of active pharmaceutical ingredient particle size was evaluated with glyburide (GLY) as a model drug. Two particle size suspensions of glyburide were prepared with average particle sizes of 42.7 and 4.1 µm. Each suspension was dosed as a single oral gavage to male Sprague Dawley rats, and concentration-time (C-t) profiles of glyburide were measured with liquid chromatography coupled with tandem mass spectrometry. A continuous rat intestine absorption model was extended to include drug dissolution and was used to predict the absorption kinetics of GLY depending on particle size. Additional literature datasets of rat GLY formulations with particle sizes ranging from 0.25 to 4.0 µm were used for model predictions. The model predicted reasonably well the absorption profiles of GLY based on varying particle size and varying feeding time. The model predicted inhibition of intestinal uptake or efflux transporters depending on the datasets. The three datasets used formulations with different excipients, which may impact the transporter activity. Model simulations indicated that the model provides a facile framework to predict the impact of transporter inhibition on drug C-t profiles. Model simulations can also be conducted to evaluate the impact of an altered intestinal lumen environment. In conclusion, the rat continuous intestine absorption model may provide a useful tool to predict the impact of varying drug formulations on rat oral absorption profiles.

Predicting Impact of Food and Feeding Time on Oral Absorption of Drugs with a Novel Rat Continuous Intestinal Absorption Model.

Intricacies in intestinal physiology, drug properties, and food effects should be incorporated into models to predict complex oral drug absorption. A previously published human continuous intestinal absorption model based on the convection-diffusion equation was modified specifically for the male Sprague-Dawley rat in this report. Species-specific physiologic conditions along intestinal length - experimental velocity and pH under fasted and fed conditions, were measured and incorporated into the intestinal absorption model. Concentration-time (C-t) profiles were measured upon a single intravenous and peroral (PO) dose for three drugs: amlodipine (AML), digoxin (DIG), and glyburide (GLY). Absorption profiles were predicted and compared with experimentally collected data under three feeding conditions: 12-hour fasted rats were provided food at two specific times after oral drug dose (1 hour and 2 hours for AML and GLY; 0.5 hours and 1 hour for DIG), or they were provided food for the entire study. Intravenous versus PO C-t profiles suggested absorption even at later times and informed design of appropriate mathematical input functions based on experimental feeding times. With this model, AML, DIG, and GLY oral C-t profiles for all feeding groups were generally well predicted, with exposure overlap coefficients in the range of 0.80-0.97. Efflux transport for DIG and uptake and efflux transport for GLY were included, modeling uptake transporter inhibition in the presence of food. Results indicate that the continuous intestinal rat model incorporates complex physiologic processes and feeding times relative to drug dose into a simple framework to provide accurate prediction of oral absorption. SIGNIFICANCE STATEMENT: A novel rat continuous intestinal model predicts drug absorption with respect to time and intestinal length. Feeding time relative to dose was modeled as a key effect. Experimental fasted/fed intestinal pH and velocity, efflux and uptake transporter expression along intestinal length, and uptake transporter inhibition in the presence of food were modeled. The model uses the pharmacokinetic profiles of three model drugs and provides a novel framework to study food effects on absorption.

Complex Cytochrome P450 Kinetics Due to Multisubstrate Binding and Sequential Metabolism. Part 1. Theoretical Considerations.

Complexities in P450-mediated metabolism kinetics include multisubstrate binding, multiple-product formation, and sequential metabolism. Saturation curves and intrinsic clearances were simulated for single-substrate and multisubstrate models using derived velocity equations and numerical solutions of ordinary differential equations (ODEs). Multisubstrate models focused on sigmoidal kinetics because of their dramatic impact on clearance predictions. These models were combined with multiple-product formation and sequential metabolism, and simulations were performed with random error. Use of single-substrate models to characterize multisubstrate data can result in inaccurate kinetic parameters and poor clearance predictions. Comparing results for use of standard velocity equations with ODEs clearly shows that ODEs are more versatile and provide better parameter estimates. It would be difficult to derive concentration-velocity relationships for complex models, but these relationships can be easily modeled using numerical methods and ODEs. SIGNIFICANCE STATEMENT: The impact of multisubstrate binding, multiple-product formation, and sequential metabolism on the P450 kinetics was investigated. Numerical methods are capable of characterizing complicated P450 kinetics.

Complex Cytochrome P450 Kinetics Due to Multisubstrate Binding and Sequential Metabolism. Part 2. Modeling of Experimental Data.

Three CYP3A4 substrates, midazolam, ticlopidine, and diazepam, display non-Michaelis-Menten kinetics, form multiple primary metabolites, and are sequentially metabolized to secondary metabolites. We generated saturation curves for these compounds and analyzed the resulting datasets using a number of single-substrate and multisubstrate binding models. These models were parameterized using rate equations and numerical solutions of the ordinary differential equations. Multisubstrate binding models provided results superior to single-substrate models, and simultaneous modeling of multiple metabolites provided better results than fitting the individual datasets independently. Although midazolam datasets could be represented using standard two-substrate models, more complex models that include explicit enzyme-product complexes were needed to model the datasets for ticlopidine and diazepam. In vivo clearance predictions improved markedly with the use of in vitro parameters from the complex models versus the Michaelis-Menten equation. The results highlight the need to use sufficiently complex kinetic schemes instead of the Michaelis-Menten equation to generate accurate kinetic parameters. SIGNIFICANCE STATEMENT: The metabolism of midazolam, ticlopidine, and diazepam by CYP3A4 results in multiple metabolites and sequential metabolism. This study evaluates the use of rate equations and numerical methods to characterize the in vitro enzyme kinetics. Use of complex cytochrome P450 kinetic models is necessary to obtain accurate parameter estimates for predicting in vivo disposition.

A hybrid model to evaluate the impact of active uptake transport on hepatic distribution of atorvastatin in rats.

1. Mathematical modeling remains a useful tool to study the impact of transporters on overall and intracellular drug disposition. The impact of organic anion transporter protein mediated uptake on atorvastatin systemic and intracellular pharmacokinetics, specifically distribution volume, was studied in rats with mathematical modeling and conducting in vivo pharmacokinetic studies for atorvastatin in presence and absence of rifampicin. A previously developed 5-compartment explicit membrane model for the liver was combined with a compartmental model to develop a semi-physiological hybrid model for atorvastatin disposition. 2. Rifampicin treatment resulted in a decrease in systemic clearance and steady-state distribution volume, and an increase in half-life of atorvastatin. The hybrid model predicted higher unbound intracellular liver atorvastatin concentrations than unbound plasma concentrations in both rifampicin treated and untreated groups, indicating involvement of uptake transporters. The intracellular unbound concentrations during the distributive phase were unaffected by rifampicin. The dependence of clearance on blood flow as well as hepatic uptake for atorvastatin (a moderate-to-high extraction ratio drug) can explain this lack of change in intracellular concentrations if there is incomplete inhibition of transport at the tested rifampicin dose. 3. The hybrid model successfully allowed the evaluation of effect of active uptake on intracellular and plasma atorvastatin disposition.

Numerical analysis of time-dependent inhibition kinetics: comparison between rat liver microsomes and rat hepatocyte data for mechanistic model fitting.

Time-dependent inhibition (TDI) may confound drug interaction predictions. Recently, models were generated for an array of TDI kinetic schemes using numerical analysis of microsomal assays. Additionally, a distinct terminal inactivation step was identified for certain mechanism based inhibitors (MBI) following reversible metabolite intermediate complex (MIC) formation. Longer hepatocyte incubations potentially allow analysis of slow TDI and terminal inactivation. In the experiments presented here, we compared the quality of TDI parameterization by numerical analysis between hepatocyte and microsomal data. Rat liver microsomes (RLM), suspended rat hepatocytes (SRH) and sandwich-cultured rat hepatocytes (SCRH) were incubated with the prototypical CYP3A MBI troleandomycin and the substrate midazolam. Data from RLM provided a better model fit as compared to SRH. Increased CYP3A expression after dexamethasone (DEX) induction improved the fit for RLM and SRH. A novel sequential kinetic scheme, defining inhibitor metabolite production prior to MIC formation, improved the fit compared to direct MIC formation. Furthermore, terminal inactivation rate constants were parameterized for RLM and SRH samples with DEX-induced CYP3A. The low expression of CYP3A and experimental error in SCRH resulted in poor data for model fitting. Overall, RLM generated data better suited for elucidation of TDI mechanisms by numerical analysis.

Impact of Lipid Partitioning on the Design, Analysis, and Interpretation of Microsomal Time-Dependent Inactivation.

Nonspecific drug partitioning into microsomal membranes must be considered for in vitro-in vivo correlations. This work evaluated the effect of including lipid partitioning in the analysis of complex TDI kinetics with numerical methods. The covariance between lipid partitioning and multiple inhibitor binding was evaluated. Simulations were performed to test the impact of lipid partitioning on the interpretation of TDI kinetics, and experimental TDI datasets for paroxetine (PAR) and itraconazole (ITZ) were modeled. For most kinetic schemes, modeling lipid partitioning results in statistically better fits. For MM-IL simulations (KI,u = 0.1 µM, kinact = 0.1 minute-1), concurrent modeling of lipid partitioning for an fumic range (0.01, 0.1, and 0.5) resulted in better fits compared with post hoc correction (AICc: -526 vs. -496, -579 vs. -499, and -636 vs. -579, respectively). Similar results were obtained with EII-IL. Lipid partitioning may be misinterpreted as double binding, leading to incorrect parameter estimates. For the MM-IL datasets, when fumic = 0.02, MM-IL, and EII model fits were indistinguishable (δAICc = 3). For less partitioned datasets (fumic = 0.1 or 0.5), the inclusion of partitioning resulted in better models. The inclusion of lipid partitioning can lead to markedly different estimates of KI,u and kinact A reasonable alternate experimental design is nondilution TDI assays, with post hoc fumic incorporation. The best fit models for PAR (MIC-M-IL) and ITZ (MIC-EII-M-IL and MIC-EII-M-Seq-IL) were consistent with their reported mechanism and kinetics. Overall, experimental fumic values should be concurrently incorporated into TDI models with complex kinetics, when dilution protocols are used.

Prediction of Tissue-Plasma Partition Coefficients Using Microsomal Partitioning: Incorporation into Physiologically based Pharmacokinetic Models and Steady-State Volume of Distribution Predictions.

Drug distribution is a necessary component of models to predict human pharmacokinetics. A new membrane-based tissue-plasma partition coefficient (K p) method (K p,mem) to predict unbound tissue to plasma partition coefficients (K pu) was developed using in vitro membrane partitioning [fraction unbound in microsomes (f um)], plasma protein binding, and log P The resulting K p values were used in a physiologically based pharmacokinetic (PBPK) model to predict the steady-state volume of distribution (V ss) and concentration-time (C-t) profiles for 19 drugs. These results were compared with K p predictions using a standard method [the differential phospholipid K p prediction method (K p,dPL)], which differentiates between acidic and neutral phospholipids. The K p,mem method was parameterized using published rat K pu data and tissue lipid composition. The K pu values were well predicted with R 2 = 0.8. When used in a PBPK model, the V ss predictions were within 2-fold error for 12 of 19 drugs for K p,mem versus 11 of 19 for Kp,dPL With one outlier removed for K p,mem and two for K p,dPL, the V ss predictions for R 2 were 0.80 and 0.79 for the K p,mem and K p,dPL methods, respectively. The C-t profiles were also predicted and compared. Overall, the K p,mem method predicted the V ss and C-t profiles equally or better than the K p,dPL method. An advantage of using f um to parameterize membrane partitioning is that f um data are used for clearance prediction and are, therefore, generated early in the discovery/development process. Also, the method provides a mechanistically sound basis for membrane partitioning and permeability for further improving PBPK models. SIGNIFICANCE STATEMENT: A new method to predict tissue-plasma partition coefficients was developed. The method provides a more mechanistic basis to model membrane partitioning.

Improved Predictions of Drug-Drug Interactions Mediated by Time-Dependent Inhibition of CYP3A.

Time-dependent inactivation (TDI) of cytochrome P450s (CYPs) is a leading cause of clinical drug-drug interactions (DDIs). Current methods tend to overpredict DDIs. In this study, a numerical approach was used to model complex CYP3A TDI in human-liver microsomes. The inhibitors evaluated included troleandomycin (TAO), erythromycin (ERY), verapamil (VER), and diltiazem (DTZ) along with the primary metabolites N-demethyl erythromycin (NDE), norverapamil (NV), and N-desmethyl diltiazem (NDD). The complexities incorporated into the models included multiple-binding kinetics, quasi-irreversible inactivation, sequential metabolism, inhibitor depletion, and membrane partitioning. The resulting inactivation parameters were incorporated into static in vitro-in vivo correlation (IVIVC) models to predict clinical DDIs. For 77 clinically observed DDIs, with a hepatic-CYP3A-synthesis-rate constant of 0.000 146 min-1, the average fold difference between the observed and predicted DDIs was 3.17 for the standard replot method and 1.45 for the numerical method. Similar results were obtained using a synthesis-rate constant of 0.000 32 min-1. These results suggest that numerical methods can successfully model complex in vitro TDI kinetics and that the resulting DDI predictions are more accurate than those obtained with the standard replot approach.

Continuous Intestinal Absorption Model Based on the Convection-Diffusion Equation.

Prediction of the rate and extent of drug absorption upon oral dosing needs models that capture the complexities of both the drug molecule and intestinal physiology. We report here the development of a continuous intestinal absorption model based on the convection-diffusion equation. The model includes explicit enterocyte apical membrane and intracellular lipid radial compartments along the length of the intestine. Physiologic functions along length x are built into the model and include velocity, diffusion, surface areas, and pH of the intestine. Also included are expression levels of the intestinal active uptake transporter OATP2B1 and efflux transporter P-gp. Oral dosing of solution as well as solid (with a dissolution function) was modeled for several drugs. The fraction absorbed (FA) and concentration-time (C-t) profiles were predicted and compared with clinical data. Overall, FA was well predicted upon oral (n = 21) or colonic dosing (n = 11), with four outliers. The overall accuracy (prediction of the correct bin) was 81% with outliers and 90% without outliers. Of the nine solution dosing data sets, six drugs were very well predicted with an exposure overlap coefficient (EOC) > 0.9 and predicted Cmax and Tmax values similar to those observed. Of the six solid dose formulations evaluated, the EOC values were > 0.9 for all drugs except budesonide. The observed precipitation of nifedipine at high doses was predicted by the model. Most of the poor predictions were for drugs that are known to be transporter substrates. As proof of concept, incorporating OATP2B1 and P-gp markedly improved the EOC and predicted Cmax and Tmax for fexofenadine. Finally, the continuous intestinal model accurately recapitulated the known relationships between drug absorption and permeability, solubility, and particle size. Together, these results indicate that this preliminary intestinal absorption model offers a simple and straightforward framework to build in complexities such as drug permeability, lipid partitioning, solubility, metabolism, and transport for improved prediction of the rate and extent of drug absorption.

Drug Distribution Part 2. Predicting Volume of Distribution from Plasma Protein Binding and Membrane Partitioning.

PURPOSE: Volume of distribution is an important pharmacokinetic parameter in the distribution and half-life of a drug. Protein binding and lipid partitioning together determine drug distribution. METHODS: Here we present a simple relationship that estimates the volume of distribution with the fraction of drug unbound in both plasma and microsomes. Model equations are based upon a two-compartment system and the experimental fractions unbound in plasma and microsomes represent binding to plasma proteins and cellular lipids, respectively. RESULTS: The protein and lipid binding components were parameterized using a dataset containing human in vitro and in vivo parameters for 63 drugs. The resulting equation explains ~84% of the variance in the log of the volume of distribution with an average fold-error of 1.6, with 3 outliers. CONCLUSIONS: These results suggest that Vss can be predicted for most drugs from plasma protein binding and microsomal partitioning.

Drug Distribution. Part 1. Models to Predict Membrane Partitioning.

PURPOSE: Tissue partitioning is an important component of drug distribution and half-life. Protein binding and lipid partitioning together determine drug distribution. METHODS: Two structure-based models to predict partitioning into microsomal membranes are presented. An orientation-based model was developed using a membrane template and atom-based relative free energy functions to select drug conformations and orientations for neutral and basic drugs. RESULTS: The resulting model predicts the correct membrane positions for nine compounds tested, and predicts the membrane partitioning for n = 67 drugs with an average fold-error of 2.4. Next, a more facile descriptor-based model was developed for acids, neutrals and bases. This model considers the partitioning of neutral and ionized species at equilibrium, and can predict membrane partitioning with an average fold-error of 2.0 (n = 92 drugs). CONCLUSIONS: Together these models suggest that drug orientation is important for membrane partitioning and that membrane partitioning can be well predicted from physicochemical properties.

On the Nature of Physiologically-Based Pharmacokinetic Models -A Priori or A Posteriori? Mechanistic or Empirical?

Physiologically-based pharmacokinetic (PBPK) models explicitly incorporate tissue-specific blood flows, partition coefficients, and metabolic processes. Since PBPK models are derived using physiologic parameters and interactions of the compound with tissue components, these models are considered to be "bottom up" as opposed to "top down". Modeling approaches can be characterized as either a posteriori (observational) or a priori (based solely on theory). Furthermore, approaches can be mechanistic (structure and components based on mechanisms) or empirical (based on observations alone). Both "bottom up" and "top down" approaches can incorporate either empirical or mechanistic components. In this perspective, we discuss some of the methods and assumptions of current PBPK modeling approaches. Specifically, we discuss drug partitioning into phospholipids and neutral lipids, use of blood-plasma ratios to estimate basic drug tissue partitioning, and clearance of neutral and acidic drugs. Based on these discussions, we believe that current PBPK models are mechanistic but a posteriori and semi-empirical.

Intracellular Unbound Atorvastatin Concentrations in the Presence of Metabolism and Transport.

Accurate prediction of drug target activity and rational dosing regimen design require knowledge of drug concentrations at the target. It is important to understand the impact of processes such as membrane permeability, partitioning, and active transport on intracellular drug concentrations. The present study aimed to predict intracellular unbound atorvastatin concentrations and characterize the effect of enzyme-transporter interplay on these concentrations. Single-pass liver perfusion studies were conducted in rats using atorvastatin (ATV, 1 µM) alone at 4°C and at 37°C in presence of rifampin (RIF, 20 µM) and 1-aminobenzotriazole (ABT, 1 mM), separately and in combination. The unbound intracellular ATV concentration was predicted with a five-compartment explicit membrane model using the parameterized diffusional influx clearance, active basolateral uptake clearance, and metabolic clearance. Chemical inhibition of uptake and metabolism at 37°C proved to be better controls relative to studies at 4°C. The predicted unbound intracellular concentration at the end of the 50-minute perfusion in the +ABT , +ABT+RIF, and the ATV-only groups was 6.5 µM, 0.58 µM, and 5.14 µM, respectively. The predicted total liver concentrations and amount recovered in bile were within 0.94-1.3 fold of the observed value in all groups. The fold difference in total liver concentration did not always extrapolate to the fold difference in predicted unbound concentration across groups. Together, these results support the use of compartmental modeling to predict intracellular concentrations in dynamic organ-based systems. These predictions can provide insight into the role of uptake transporters and metabolizing enzymes in determining drug tissue concentrations.

Mechanism-Based Inhibition of CYP3A4 by Podophyllotoxin: Aging of an Intermediate Is Important for in Vitro/in Vivo Correlations.

An in vitro observation of time-dependent inhibition (TDI) of metabolic enzymes often results in removing a potential drug from the drug pipeline. However, the accepted method for predicting TDIs of the important drug metabolizing cytochrome P450 enzymes often overestimates the drug interaction potential. Better models that take into account the complexities of the cytochrome P450 enzyme system will lead to better predictions. Herein we report the use of our previously described models for complex kinetics of podophyllotoxin. Spectral characterization of the kinetics indicates that an intermediate MI complex is formed, which slowly progresses to an essentially irreversible MI complex. The intermediate MI complex can release free enzyme during the time course of a typical 30 min TDI experiment. This slow rate of MI complex conversion results in an overprediction of the kinact value if this process is not included in the analysis of the activity versus time profile. In vitro kinetic experiments in rat liver microsomes predicted a lack of drug interaction between podophyllotoxin and midazolam. In vivo rat pharmacokinetic studies confirmed this lack of drug interaction.

A physiologically based pharmacokinetic model to predict the pharmacokinetics of highly protein-bound drugs and the impact of errors in plasma protein binding.

Predicting the pharmacokinetics of highly protein-bound drugs is difficult. Also, since historical plasma protein binding data were often collected using unbuffered plasma, the resulting inaccurate binding data could contribute to incorrect predictions. This study uses a generic physiologically based pharmacokinetic (PBPK) model to predict human plasma concentration-time profiles for 22 highly protein-bound drugs. Tissue distribution was estimated from in vitro drug lipophilicity data, plasma protein binding and the blood: plasma ratio. Clearance was predicted with a well-stirred liver model. Underestimated hepatic clearance for acidic and neutral compounds was corrected by an empirical scaling factor. Predicted values (pharmacokinetic parameters, plasma concentration-time profile) were compared with observed data to evaluate the model accuracy. Of the 22 drugs, less than a 2-fold error was obtained for the terminal elimination half-life (t1/2 , 100% of drugs), peak plasma concentration (Cmax , 100%), area under the plasma concentration-time curve (AUC0-t , 95.4%), clearance (CLh , 95.4%), mean residence time (MRT, 95.4%) and steady state volume (Vss , 90.9%). The impact of fup errors on CLh and Vss prediction was evaluated. Errors in fup resulted in proportional errors in clearance prediction for low-clearance compounds, and in Vss prediction for high-volume neutral drugs. For high-volume basic drugs, errors in fup did not propagate to errors in Vss prediction. This is due to the cancellation of errors in the calculations for tissue partitioning of basic drugs. Overall, plasma profiles were well simulated with the present PBPK model. Copyright © 2016 John Wiley & Sons, Ltd.

"Target-Site" Drug Metabolism and Transport.

The recent symposium on "Target-Site" Drug Metabolism and Transport that was sponsored by the American Society for Pharmacology and Experimental Therapeutics at the 2014 Experimental Biology meeting in San Diego is summarized in this report. Emerging evidence has demonstrated that drug-metabolizing enzyme and transporter activity at the site of therapeutic action can affect the efficacy, safety, and metabolic properties of a given drug, with potential outcomes including altered dosing regimens, stricter exclusion criteria, or even the failure of a new chemical entity in clinical trials. Drug metabolism within the brain, for example, can contribute to metabolic activation of therapeutic drugs such as codeine as well as the elimination of potential neurotoxins in the brain. Similarly, the activity of oxidative and conjugative drug-metabolizing enzymes in the lung can have an effect on the efficacy of compounds such as resveratrol. In addition to metabolism, the active transport of compounds into or away from the site of action can also influence the outcome of a given therapeutic regimen or disease progression. For example, organic anion transporter 3 is involved in the initiation of pancreatic ß-cell dysfunction and may have a role in how uremic toxins enter pancreatic ß-cells and ultimately contribute to the pathogenesis of gestational diabetes. Finally, it is likely that a combination of target-specific metabolism and cellular internalization may have a significant role in determining the pharmacokinetics and efficacy of antibody-drug conjugates, a finding which has resulted in the development of a host of new analytical methods that are now used for characterizing the metabolism and disposition of antibody-drug conjugates. Taken together, the research summarized herein can provide for an increased understanding of potential barriers to drug efficacy and allow for a more rational approach for developing safe and effective therapeutics.

A numerical method for analysis of in vitro time-dependent inhibition data. Part 1. Theoretical considerations.

Inhibition of cytochromes P450 by time-dependent inhibitors (TDI) is a major cause of clinical drug-drug interactions. It is often difficult to predict in vivo drug interactions based on in vitro TDI data. In part 1 of these manuscripts, we describe a numerical method that can directly estimate TDI parameters for a number of kinetic schemes. Datasets were simulated for Michaelis-Menten (MM) and several atypical kinetic schemes. Ordinary differential equations were solved directly to parameterize kinetic constants. For MM kinetics, much better estimates of KI can be obtained with the numerical method, and even IC50 shift data can provide meaningful estimates of TDI kinetic parameters. The standard replot method can be modified to fit non-MM data, but normal experimental error precludes this approach. Non-MM kinetic schemes can be easily incorporated into the numerical method, and the numerical method consistently predicts the correct model at errors of 10% or less. Quasi-irreversible inactivation and partial inactivation can be modeled easily with the numerical method. The utility of the numerical method for the analyses of experimental TDI data is provided in our companion manuscript in this issue of Drug Metabolism and Disposition (Korzekwa et al., 2014b).