Extrapolating In vitro Metabolic Interactions to Isolated Perfused Liver: Predictions of Metabolic Interactions between R-Bufuralol, Bunitrolol, and Debrisoquine

Estimation of toluene exposure in air from BMA (S-benzylmercapturic acid) urinary measures using a reverse dosimetry approach based on physiologically pharmacokinetic modeling

This study aimed to use a reverse dosimetry PBPK modeling approach to estimate toluene atmospheric exposure from urinary measurements of S-benzylmercapturic acid (BMA) in a small group of individuals and to evaluate the uncertainty associated to urinary spot-sampling compared to 24-h collected urine samples. Each exposure assessment technique was developed namely to estimate toluene air exposure from BMA measurements in 24-h urine samples (24-h-BMA) and from distributions of daily urinary BMA spot measurements (DUBSM). Model physiological parameters were described based upon age, weight, size and sex. Monte Carlo simulations with the PBPK model allowed converting DUBSM distribution (and 24-h-BMA) into toluene air levels. For the approach relying on DUBSM distribution, the ratio between the 95% probability of predicted toluene concentration and its 50% probability in each individual varied between 1.2 and 1.4, while that based on 24-h-BMA varied between 1.0 and 1.1. This suggests more variability in estimated exposure from spot measurements. Thus, estimating toluene exposure based on DUBSM distribution generated about 20% more uncertainty. Toluene levels estimated (0.0078-0.0138 ppm) are well below Health Canada's maximum chronic air guidelines. PBPK modeling and reverse dosimetry may be combined to interpret urinary metabolites data of VOCs and assess related uncertainties.

Role of Physiologically Based Kinetic modelling in addressing environmental chemical mixtures - A review

•

The availability and applicability of Physiologically Based Kinetic (PBK) models for mixtures is reviewed.

•

PBK models can support risk assessment of mixtures by incorporating the toxicokinetic processes.

•

Quantitative structure-activity relationship (QSAR) models can be used to fill data gaps in PBK modelling.

•

PBK models for mixtures can be improved by including various types of interactions.

The role of Physiologically Based Kinetic (PBK) modelling in assessing mixture toxicology has been growing for the last three decades. It has been widely used to investigate and address interactions in mixtures. This review describes the current state-of-the-art of PBK models for chemical mixtures and to evaluate the applications of PBK modelling for mixtures with emphasis on their role in chemical risk assessment. A total of 35 mixture PBK models were included after searching web resources (Scopus, PubMed, Web of Science, and Google Scholar), screening for duplicates, and excluding articles based on eligibility criteria. Binary mixtures and volatile organic compounds accounted for two-thirds of the chemical mixtures identified. The most common exposure route and modelled system were found to be inhalation and rats respectively. Twenty two (22) models were for binary mixtures, 5 for ternary mixtures, 3 for quaternary mixtures, and 5 for complex mixtures. Both bottom-up and top-down PBK modelling approaches are described. Whereas bottom-up approaches are based on a series of binary interactions, top-down approaches are based on the lumping of mixture components. Competitive inhibition is the most common type of interaction among the various types of mixtures, and usually becomes a concern at concentrations higher than environmental exposure levels. It leads to reduced biotransformation that either means a decrease in the amount of toxic metabolite formation or an increase in toxic parent chemical accumulation. The consequence is either lower or higher toxicity compared to that estimated for the mixture based on the additivity principle. Therefore, PBK modelling can play a central role in predicting interactions in chemical mixture risk assessment.

Reverse dosimetry modeling of toluene exposure concentrations based on biomonitoring levels from the Canadian health measures survey

Biomonitoring might provide useful estimates of population exposure to environmental chemicals. However, data uncertainties stemming from interindividual variability are common in large population biomonitoring surveys. Physiologically based pharmacokinetic (PBPK) models might be used to account for age- and gender-related variability in internal dose. The objective of this study was to reconstruct air concentrations consistent with blood toluene measures reported in the third Canadian Health Measures Survey using reverse dosimetry PBPK modeling techniques. Population distributions of model's physiological parameters were described based upon age, weight, and size for four subpopulations (12-19, 20-39, 40-59, and 60-79 years old). Monte Carlo simulations applied to PBPK modeling allowed converting the distributions of venous blood measures of toluene obtained from CHMS into related air levels. Based upon blood levels observed at the 50^th, 90^th and 95^th percentiles, corresponding air toluene concentrations were estimated for teenagers aged 12-19 years as being, respectively, 0.009, 0.04 and 0.06 ppm. Similarly, values were computed for adults aged 20-39 years (0.007, 0.036, and 0.06 ppm), 40-59 years (0.007, 0.036 and 0.06 ppm) and 60-79 years (0.006, 0.022 and 0.04 ppm). These estimations are well below Health Canada's maximum recommended chronic air guidelines for toluene. In conclusion, PBPK modeling and reverse dosimetry may be combined to help interpret biomonitoring data for chemical exposure in large population surveys and estimate the associated toxicological health risk.

Predictions of bisphenol A hepatic clearance in the isolated perfused rat liver (IPRL): impact of albumin binding and of co-administration with naproxen

1. This study aimed (i) to characterise hepatic clearance (CL) of bisphenol A (BPA) and naproxen (NAP) administered alone or in binary mixtures to highlight the influence of a binding to albumin (ALB) using an isolated perfused rat liver (IPRL) system; and (ii) to compare results of prediction algorithms with measured clearance rates. 2. The IPRL system and liver microsomes were used to determine the metabolic constants of BPA and NAP either in the presence or absence of ALB. In this study, the IPRL was used as proxy for the in vivo situation. Accordingly, diverse in vitro-to-in vivo and in vivo-to-in vivo extrapolations (IVIVEs) were made to predict CL of BPA determined in situ/in vivo with ALB from metabolic data determined without ALB by using different binding correction methods (i.e., direct and conventional scaling as well as a novel scaling considering an ALB-facilitated uptake mechanism). 3. The addition of ALB significantly influenced the liver kinetics of BPA and NAP either administered alone or in binary mixtures, which was reflected in the Michaelis-Menten constants. Analysis of concomitant exposures of BPA and NAP gave a fully competitive inhibition. Furthermore, the IVIVE method based on the ALB-facilitated uptake mechanism provided the most accurate predictions of CL_{in vivo} as compared with the other IVIVE methods when the impact of ALB is considered. 4. Our findings support the notion that high binding to ALB reduces the biotransformation of BPA and NAP when administered alone or in mixtures in the IPRL system. However, the free drug concentration in liver in vivo is probably higher than expected since the IVIVE method based on a potential ALB-facilitated uptake mechanism is the most robust prediction method. Overall, this study should improve the physiologically-based pharmacokinetic (PBPK) modelling of chemical-drug interactions.

The Role of Extracellular Binding Proteins in the Cellular Uptake of Drugs: Impact on Quantitative In Vitro-to-In Vivo Extrapolations of Toxicity and Efficacy in Physiologically Based Pharmacokinetic-Pharmacodynamic Research

A critical component in the development of physiologically based pharmacokinetic-pharmacodynamic (PBPK/PD) models for estimating target organ dosimetry in pharmacology and toxicology studies is the understanding of the uptake kinetics and accumulation of drugs and chemicals at the cellular level. Therefore, predicting free drug concentrations in intracellular fluid will contribute to our understanding of concentrations at the site of action in cells in PBPK/PD research. Some investigators believe that uptake of drugs in cells is solely driven by the unbound fraction; conversely, others argue that the protein-bound fraction contributes a significant portion of the total amount delivered to cells. Accordingly, the current literature suggests the existence of a so-called albumin-mediated uptake mechanism(s) for the protein-bound fraction (i.e., extracellular protein-facilitated uptake mechanisms) at least in hepatocytes and cardiac myocytes; however, such mechanism(s) and cells from other organs deserve further exploration. Therefore, the main objective of this present study was to discuss further the implication of potential protein-facilitated uptake mechanism(s) on drug distribution in cells under in vivo conditions. The interplay between the protein-facilitated uptake mechanism(s) and the effects of a pH gradient, metabolism, transport, and permeation limitation potentially occurring in cells was also discussed, as this should violate the basic assumption on similar free drug concentration in cells and plasma. This was made because the published equations used to calculate drug concentrations in cells in a PBPK/PD model did not consider potential protein-facilitated uptake mechanism(s). Consequently, we corrected some published equations for calculating the free drug concentrations in cells compared with plasma in PBPK/PD modeling studies, and we proposed a refined strategy for potentially performing more accurate quantitative in vitro-to-in vivo extrapolations (IVIVEs) of toxicity (efficacy) at the cellular level from data generated in cell assays. Overall, this present study may help to optimize the human dose prediction in preclinical and clinical studies, while prescribing drugs with narrow therapeutic windows that are highly bound to extracellular proteins and/or highly ionized at the physiological pH. This may facilitate building a more accurate safety (efficacy) profile for such drugs.

A paradigm shift in pharmacokinetic-pharmacodynamic (PKPD) modeling: rule of thumb for estimating free drug level in tissue compared with plasma to guide drug design

A basic assumption in pharmacokinetics-pharmacodynamics research is that the free drug concentration is similar in plasma and tissue, and, hence, in vitro plasma data can be used to estimate the in vivo condition in tissue. However, in a companion manuscript, it has been demonstrated that this assumption is violated for the ionized drugs. Nonetheless, these observations focus on in vitro static environments and do not challenge data with an in vivo dynamic system. Therefore, an extension from an in vitro to an in vivo system becomes the necessary next step. The objective of this study was to perform theoretical simulations of the free drug concentration in tissue and plasma by using a physiologically based pharmacokinetics (PBPK) model reproducing the in vivo conditions in human. Therefore, the effects of drug ionization, lipophilicity, and clearance have been taken into account in a dynamic system. This modeling exercise was performed as a proof of concept to demonstrate that free drug concentration in tissue and plasma may also differ in a dynamic system for passively permeable drugs that are ionized at the physiological pH. The PBPK model simulations indicated that free drug concentrations in tissue cells and plasma significantly differ for the ionized drugs because of the pH gradient effect between cells and interstitial space. Hence, a rule of thumb for potentially performing more accurate PBPK/PD modeling is suggested, which states that the free drug concentration in tissue and plasma will differ for the ionizable drugs in contrast to the neutral drugs. In addition to the pH gradient effect for the ionizable drugs, lipophilicity and clearance effects will increase or decrease the free drug concentration in tissue and plasma for each class of drugs; thus, higher will be the drug lipophilicity and clearance, lower would be the free drug concentration in plasma, and, hence, in tissue, in a dynamic in vivo system. Therefore, only considering the value of free fraction in plasma derived from a static in vitro environment might be biased to guide drug design (the old paradigm), and, hence, it is recommended to use a PBPK model to reproduce more accurately the in vivo condition in tissue (the new paradigm). This newly developed approach can be used to predict free drug concentration in diverse tissue compartments for small molecules in toxicology and pharmacology studies, which can be leveraged to optimize the pharmacokinetics drivers of tissue distribution based upon physicochemical and physiological input parameters in an attempt to optimize free drug level in tissue. Overall, this present study provides guidance on the application of plasma and tissue concentration information in PBPK/PD research in preclinical and clinical studies, which is in accordance with the recent literature.

Drug Distribution to Human Tissues: Prediction and Examination of the Basic Assumption in In Vivo Pharmacokinetics-Pharmacodynamics (PK/PD) Research

The tissue:plasma partition coefficients (Kp ) are good indicators of the extent of tissue distribution. Therefore, advanced tissue composition-based models were used to predict the Kp values of drugs under in vivo conditions on the basis of in vitro and physiological input data. These models, however, focus on animal tissues and do not challenge the predictions with human tissues for drugs. The first objective of this study was to predict the experimentally determined Kp values of seven human tissues for 26 drugs. In all, 95% of the predicted Kp values are within 2.5-fold error of the observed values in humans. Accordingly, these results suggest that the tissue composition-based model used in this study is able to provide accurate estimates of drug partitioning in the studied human tissues. Furthermore, as the Kp equals to the ratio of total concentration between tissue and plasma, or the ratio of unbound fraction between plasma (fup ) and tissue (fut ), this parameter Kp would deviate from the unity. Therefore, the second objective was to examine the corresponding relationships between fup and fut values experimentally determined in humans for several drugs. The results also indicate that fup may significantly deviate to fut ; the discrepancies are governed by the dissimilarities in the binding and ionization on both sides of the membrane, which were captured by the tissue composition-based model. Hence, this violated the basic assumption in in vivo pharmacokinetics-pharmacodynamics (PK/PD) research, since the free drug concentration in tissue and plasma was not equal particularly for the ionizable drugs due to the pH gradient effect on the fraction of unionized drug in plasma (fuip ) and tissue (fuit ) (i.e., fup × fuip × total plasma concentration = fut × fuit × total tissue concentration, and, hence, the free drug concentration in plasma and tissue differed by fuip/fuit). Therefore, this assumption should be adjusted for the ionized drugs, and, hence, a mathematical correction to the basic assumption of similar free drug concentration in plasma and tissues can be derived from the tissue composition-based model. Note that this assumption will be further challenged in a dynamic in vivo system in a companion manuscript. Overall, this study was a first attempt to predict the in vivo Kp values for specific human tissues by considering separately the effect of fup and fut , with the aim of facilitating the use of physiologically-based PK (PBPK) model in PK/PD studies.

Prediction of drug disposition on the basis of its chemical structure

The chemical structure of any drug determines its pharmacokinetics and pharmacodynamics. Detailed understanding of relationships between the drug chemical structure and individual disposition pathways (i.e., distribution and elimination) is required for efficient use of existing drugs and effective development of new drugs. Different approaches have been developed for this purpose, ranging from statistics-based quantitative structure-property (or structure-pharmacokinetic) relationships (QSPR) analysis to physiologically based pharmacokinetic (PBPK) models. This review critically analyzes currently available approaches for analysis and prediction of drug disposition on the basis of chemical structure. Models that can be used to predict different aspects of disposition are presented, including: (a) value of the individual pharmacokinetic parameter (e.g., clearance or volume of distribution), (b) efficiency of the specific disposition pathway (e.g., biliary drug excretion or cytochrome P450 3A4 metabolism), (c) accumulation in a specific organ or tissue (e.g., permeability of the placenta or accumulation in the brain), and (d) the whole-body disposition in the individual patients. Examples of presented pharmacological agents include "classical" low-molecular-weight compounds, biopharmaceuticals, and drugs encapsulated in specialized drug-delivery systems. The clinical efficiency of agents from all these groups can be suboptimal, because of inefficient permeability of the drug to the site of action and/or excessive accumulation in other organs and tissues. Therefore, robust and reliable approaches for chemical structure-based prediction of drug disposition are required to overcome these limitations. PBPK models are increasingly being used for prediction of drug disposition. These models can reflect the complex interplay of factors that determine drug disposition in a mechanistically correct fashion and can be combined with other approaches, for example QSPR-based prediction of drug permeability and metabolism, pharmacogenomic data and tools, pharmacokinetic-pharmacodynamic modeling approaches, etc. Moreover, the PBPK models enable detailed analysis of clinically relevant scenarios, for example the effect of the specific conditions on the time course of the analyzed drug in the individual organs and tissues, including the site of action. It is expected that further development of such combined approaches will increase their precision, enhance the effectiveness of drugs, and lead to individualized drug therapy for different patient populations (geriatric, pediatric, specific diseases, etc.).

Troubleshooting and improving the mouse and rat isolated perfused liver preparation

INTRODUCTION

Isolated perfused liver (IPL) model is not only widely performed in rats but is also used in mouse liver, although a detailed description of this procedure is absent. A comparison of the different techniques used on rats and mice will be discussed in this article association with a detailed description of the surgical and technical aspects needed to obtain and maintain the integrity of the livers during the organ isolation and perfusion.

METHODS

The surgical procedures, the IPL set-up, and the evaluation of hepatic function and damage will be described in relation to both rats and mice. In particular, the heparin dosage and administration, the portal vein cannulation avoiding portal leakage, the use of suprahepatic caval vein output, and the insertion of a cannula for bile collection will be reported. For the settings, the perfusion circuit, the perfusion solution, the temperature and the flow rate will be described, with particular regard to the balance between perfusion pressure and oxygen delivery. The monitoring of liver integrity by measuring oxygen concentration and calculating oxygen delivery rate and oxygen uptake rate, and recommendations for the collection of perfusate and bile samples will be considered. Accurate pH measurement with normalization, and the perfusion portal pressure assay by a calibrated water manometer will be also reported.

RESULTS AND DISCUSSION

This work analyzes the parameters crucial to performing a correct IPL both in rat and mouse, comparing our experience with the equivalent practice from other laboratories. An updated example of IPL applications in liver toxicology and pharmacology, physiology and pathophysiology, and liver graft preservation will be briefly presented, underlining how this technique provides essential information allowing a more accurate planning of the in vivo studies.

In vitro-in vivo extrapolation of clearance: modeling hepatic metabolic clearance of highly bound drugs and comparative assessment with existing calculation methods

In vitro-in vivo extrapolation (IVIVE) is an important method for estimating the hepatic metabolic clearance (CL) of drugs. This study highlights a problematic area observed when using microsomal data to predict in vivo CL of drugs that are highly bound to plasma proteins, and further explores mechanisms for human CL predictions by associating additional processes to IVIVE disconnect. Therefore, this study attempts to develop a novel IVIVE calculation method, which consists of adjusting the binding terms in a well-stirred liver model. A comparative assessment between the IVIVE method proposed here and previously published methods of Obach (1999. Drug Metab Dispos 27:1350-1359) and Berezhkovskiy (2010. J Pharm Sci 100:1167-1783) was also performed. The assessment was confined by the availability of measured in vitro and in vivo data in humans for 25 drugs highly bound to plasma proteins, for which it can be assumed that metabolism is the major route of elimination. Here, we argue that a difference in drug ionization and binding proteins such as albumin (AL) and alpha-1-acid glycoprotein (AAG) in plasma and liver also needs to be considered in IVIVE based on mechanistic studies. Therefore, converting unbound fraction in plasma to liver essentially increased the predicted CL values, which resulted in much more accurate estimates of in vivo CL as compared with the other IVIVE methods tested. The impact on CL estimate was more apparent for drugs binding to AL than to AAG. This is a mechanistic rational for explaining a considerable proportion of the divergence between previously estimated and observed CL values. Human CL was predicted within 1.5-fold, twofold, and threefold of the observed CL for 84%, 96%, and 100% of the compounds, respectively. Overall, this study demonstrates a significant improvement in the mechanism-based prediction of metabolic CL for these 25 highly bound drugs from in vitro data determined with microsomes, which should facilitate the application of physiologically based pharmacokinetic (PBPK) models in drug discovery and development.

Evaluating pharmacokinetic and pharmacodynamic interactions with computational models in supporting cumulative risk assessment

Simultaneous or sequential exposure to multiple chemicals may cause interactions in the pharmacokinetics (PK) and/or pharmacodynamics (PD) of the individual chemicals. Such interactions can cause modification of the internal or target dose/response of one chemical in the mixture by other chemical(s), resulting in a change in the toxicity from that predicted from the summation of the effects of the single chemicals using dose additivity. In such cases, conducting quantitative cumulative risk assessment for chemicals present as a mixture is difficult. The uncertainties that arise from PK interactions can be addressed by developing physiologically based pharmacokinetic (PBPK) models to describe the disposition of chemical mixtures. Further, PK models can be developed to describe mechanisms of action and tissue responses. In this article, PBPK/PD modeling efforts conducted to investigate chemical interactions at the PK and PD levels are reviewed to demonstrate the use of this predictive modeling framework in assessing health risks associated with exposures to complex chemical mixtures.