Hepatocyte Composition-Based Model as a Mechanistic Tool for Predicting the Cell Suspension: Aqueous Phase Partition Coefficient of Drugs in In Vitro Metabolic Studies

Application of partition coefficient methods to predict tissue:plasma affinities in common farm animals: Influence of ionisation state

Tissue affinities are conventionally determined from in vivo steady-state tissue and plasma or plasma-water chemical concentration data. In silico approaches were initially developed for preclinical species but standardly applied and tested in human physiologically-based kinetic (PBK) models. Recently, generic PBK models for farm animals have been made available and require partition coefficients as input parameters. In the current investigation, data for species-specific tissue compositions have been collected, and prediction of chemical distribution in various tissues of livestock species for cattle, chicken, sheep and swine have been performed. Overall, tissue composition was very similar across the four farm animal species. However, small differences were observed in moisture, fat and protein content in the various organs within each species. Such differences could be attributed to factors such as variations in age, breed, and weight of the animals and general conditions of the animal itself. With regards to the predictions of tissue:plasma partition coefficients, 80 %, 71 %, 77 % of the model predictions were within a factor 10 using the methods of Berezhkovskiy (2004), Rodgers and Rowland (2006) and Schmitt (2008). The method of Berezhkovskiy (2004) was often providing the most reliable predictions except for swine, where the method of Schmitt (2008) performed best. In addition, investigation of the impact of chemical classes on prediction performance, all methods had very similar reliability. Notwithstanding, no clear pattern regarding specific chemicals or tissues could be detected for the values predicted outside a 10-fold change in certain chemicals or specific tissues. This manuscript concludes with the need for future research, particularly focusing on lipophilicity and species differences in protein binding.

Non-specific binding of compounds in in vitro metabolism assays: a comparison of microsomal and hepatocyte binding in different species and an assessment of the accuracy of prediction models

Non-specific binding in in vitro metabolism systems leads to an underestimation of the true intrinsic metabolic clearance of compounds being studied. Therefore in vitro binding needs to be accounted for when extrapolating in vitro data to predict the in vivo metabolic clearance of a compound. While techniques exist for experimentally determining the fraction of a compound unbound in in vitro metabolism systems, early in drug discovery programmes computational approaches are often used to estimate the binding in the in vitro system.Experimental fraction unbound data (n = 60) were generated in liver microsomes (fu_mic) from five commonly used pre-clinical species (rat, mouse, dog, minipig, monkey) and humans. Unbound fraction in incubations with mouse, rat or human hepatocytes was determined for the same 60 compounds. These data were analysed to determine the relationship between experimentally determined binding in the different matrices and across different species. In hepatocytes there was a good correlation between fraction unbound in human and rat (r²=0.86) or mouse (r²=0.82) hepatocytes. Similar correlations were observed between binding in human liver microsomes and microsomes from rat, mouse, dog, Göttingen minipig or monkey liver microsomes (r² of >0.89, n = 51 - 52 measurements in different species). Physicochemical parameters (logP, pKa and logD) were predicted for all evaluated compounds. In addition, logP and/or logD were measured for a subset of compounds.Binding to human hepatocytes predicted using 5 different methods was compared to the measured data for a set of 59 compounds. The best methods evaluated used measured microsomal binding in human liver microsomes to predict hepatocyte binding. The collated physicochemical data were used to predict the human fu_mic using four different in silico models for a set of 53-60 compounds. The correlation (r²) and root mean square error between predicted and observed microsomal binding was 0.69 & 0.20, 0.47 & 0.23, 0.56 & 0.21 and 0.54 & 0.26 for the Turner-Simcyp, Austin, Hallifax-Houston and Poulin models, respectively. These analyses were extended to include measured literature values for binding in human liver microsomes for a larger set of compounds (n=697). For the larger dataset of compounds, microsomal binding was well predicted for neutral compounds (r²=0.67 - 0.70) using the Poulin, Austin, or Turner-Simcyp methods but not for acidic or basic compounds (r²<0.5) using any of the models. While the lipophilicity-based models can be used, the in vitro binding should be measured for compounds where more certainty is needed, using appropriately calibrated assays and possibly established weak, moderate, and strong binders as reference compounds to allow comparison across databases.

A Novel Experimental and Theoretical Method for Estimating Albumin-Mediated Hepatic Uptake Based on the Albumin Binding Fraction in Plasma and Human PK Prediction Using a Physiologically-Based Pharmacokinetic Approach

Recently, protein-facilitated uptake has been suggested to be an important factor in the precise prediction of the pharmacokinetic (PK) profiles of drugs. In our previous study, a physiologically-based pharmacokinetic (PBPK) approach considering the mechanism of albumin-mediated hepatic uptake was developed for predicting human PK profiles. It was assumed that drugs affected by albumin-mediated hepatic uptake would bind only to albumin, which means that there would be over-estimation of the contribution of protein-facilitated uptake for a drug that could bind to multiple proteins. In this study, we developed a method that can evaluate the albumin binding fraction in plasma considering the affinity for other proteins. Based on the albumin binding fraction, the contribution of albumin-mediated hepatic uptake was theoretically estimated, and then the human PK profiles were predicted by our proposed PBPK approach incorporating this mechanism. As a result, the predicted human PK profiles agreed well with the observed ones, and the absolute average fold error of PK parameters was almost within a 1.5-fold error on average. These findings show the importance of considering protein-facilitated uptake and also suggest that our proposed PBPK approach can be useful in scientific discussions with regulatory authorities.

Atropselective Partitioning of Polychlorinated Biphenyls in a HepG2 Cell Culture System: Experimental and Modeling Results

Cell culture models are used to study the toxicity of polychlorinated biphenyls (PCBs); however, it is typically unknown how much PCB enters the cells and, for chiral PCBs, if the partitioning is atropselective. We investigated the partitioning of racemic PCB 91, PCB 95, PCB 132, and PCB 136 in HepG2 cells following a 72 h incubation. PCBs were present in the cell culture medium (60.7-88.8%), cells (8.0-14.6%), and dishes (2.3-7.8%) and displayed atropisomeric enrichment in the cells (enantiomeric fraction [EF] = 0.55-0.77) and dishes (EF = 0.53-0.68). Polyparameter linear free energy relationships coupled with a composition-based model provided a good estimate of the PCB levels in the cells and cell culture medium. The free concentration was subsequently used to extrapolate from the nominal cell culture concentration to PCB tissue levels and vice versa. This approach can be used for in vitro-in vivo extrapolations for all 209 PCB congeners. However, this model (and modified models based on descriptors incorporating atropselective interactions, i.e., relative retention times on chiral columns) did not predict the atropselective partitioning in the cell culture system. Improved chemical descriptors that account for the atropselective binding of PCBs to biological macromolecules are, therefore, needed to predict the atropselective partitioning of PCBs in biological systems.

Application of the Tissue Composition-Based Model to Minipig for Predicting the Volume of Distribution at Steady State and Dermis-to-Plasma Partition Coefficients of Drugs Used in the Physiologically Based Pharmacokinetics Model in Dermatology

The minipig continues to build a reputation as a viable alternative large animal model to predict humans in dermatology and toxicology studies. Therefore, it is essential to describe and predict the pharmacokinetics in that species to speed up the clinical candidate selection. Essential input parameters in whole-body physiologically based pharmacokinetic models are the tissue-to-plasma partition coefficients and the resulting volume of distribution at steady-state (V_ss). Mechanistic in vitro- and in silico-based models used for predicting these parameters of tissue distribution of drugs refer to the tissue composition-based model (TCM). Robust TCMs were initially developed for some preclinical species (e.g., rat and dog) and human; however, there is currently no model available for the minipig. Therefore, the objective of this present study was to develop a TCM for the minipig and to estimate the corresponding tissue composition data. Drug partitioning into the tissues was predominantly governed by lipid and protein binding effects in addition to drug solubilization and pH gradient effects in the aqueous phase on both sides of the biological membranes; however, some more complex tissue distribution processes such as drug binding to the collagen-laminin material in dermis and a restricted drug partitioning into membranes of tissues for compounds that are amphiphilic and contain sulfur atom(s) were also challenged. The model was validated by predicting V_ss and the dermis-to-plasma partition coefficients (K_p-dermis) of 68 drugs. The prediction of K_p-dermis was extended to humans for comparison with the minipig. The results indicate that the extended TCM provided generally good agreements with observations in the minipig showing that it is also applicable to this preclinical species. In general, up to 86% and 100% of the predicted V_ss values are respectively within 2-fold and 3-fold errors compared to the experimentally determined values, whereas these numbers are 78% and 94% for K_p-dermis when the anticipated outlier compounds are not included. Binding data to dermis are comparable between minipigs and humans. Overall, this study is a first step toward developing a mechanistic TCM for the minipig, with the aim of increasing the use of physiologically based pharmacokinetic models of drugs for that species in addition to rats, dogs, and humans because such models are used in preclinical and clinical transdermal studies.

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.

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 V_ss can be predicted for most drugs from plasma protein binding and microsomal partitioning.

Combining Chimeric Mice with Humanized Liver, Mass Spectrometry, and Physiologically-Based Pharmacokinetic Modeling in Toxicology

Species differences exist in terms of drug oxidation activities, which are mediated mainly by cytochrome P450 (P450) enzymes. To overcome the problem of species extrapolation, transchromosomic mice containing a human P450 3A cluster or chimeric mice transplanted with human hepatocytes have been introduced into the human toxicology research area. In this review, drug metabolism and disposition mediated by humanized livers in chimeric mice are summarized in terms of biliary/urinary excretions of phthalate and bisphenol A and plasma clearances of the human cocktail probe drugs caffeine, warfarin, omeprazole, metoprolol, and midazolam. Simulation of human plasma concentrations of the teratogen thalidomide and its human metabolites is possible with a simplified physiologically based pharmacokinetic model based on data obtained in chimeric mice, in accordance with reported clinical thalidomide concentrations. In addition, in vivo nonspecific hepatic protein binding parameters of metabolically activated ¹⁴C-drug candidate and hepatotoxic medicines in humanized liver mice can be analyzed by accelerator mass spectrometry and are useful for predictions in humans.

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.

Human biofluid concentrations of mono(2-ethylhexyl)phthalate extrapolated from pharmacokinetics in chimeric mice with humanized liver administered with di(2-ethylhexyl)phthalate and physiologically based pharmacokinetic modeling

Di(2-ethylhexyl)phthalate (DEHP) is a reproductive toxicant in male rodents. The aim of the current study was to extrapolate the pharmacokinetics and toxicokinetics of mono(2-ethylhexyl)phthalate (MEHP, a primary metabolite of DEHP) in humans by using data from oral administration of DEHP to chimeric mice transplanted with human hepatocytes. MEHP and its glucuronide were detected in plasma from control mice and chimeric mice after single oral doses of 250mg DEHP/kg body weight. Biphasic plasma concentration-time curves of MEHP and its glucuronide were seen only in control mice. MEHP and its glucuronide were extensively excreted in urine within 24h in mice with humanized liver. In contrast, fecal excretion levels of MEHP glucuronide were high in control mice compared with those with humanized liver. Adjusted animal biomonitoring equivalents from chimeric mice studies were scaled to human biomonitoring equivalents using known species allometric scaling factors and in vitro metabolic clearance data with a simple physiologically based pharmacokinetic (PBPK) model. Estimated urine MEHP concentrations in humans were consistent with reported concentrations. This research illustrates how chimeric mice transplanted with human hepatocytes in combination with a simple PBPK model can assist evaluations of pharmacokinetics or toxicokinetics of the primary or secondary metabolites of DEHP.

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.

Application of mass balance models and the chemical activity concept to facilitate the use of in vitro toxicity data for risk assessment

Practical, financial, and ethical considerations related to conducting extensive animal testing have resulted in various initiatives to promote and expand the use of in vitro testing data for chemical evaluations. Nominal concentrations in the aqueous phase corresponding to an effect (or biological activity) are commonly reported and used to characterize toxicity (or biological response). However, the true concentration in the aqueous phase can be substantially different from the nominal. To support in vitro test design and aid the interpretation of in vitro toxicity data, we developed a mass balance model that can be parametrized and applied to represent typical in vitro test systems. The model calculates the mass distribution, freely dissolved concentrations, and cell/tissue concentrations corresponding to the initial nominal concentration and experimental conditions specified by the user. Chemical activity, a metric which can be used to assess the potential for baseline toxicity to occur, is also calculated. The model is first applied to a set of hypothetical chemicals to illustrate the degree to which test conditions (e.g., presence or absence of serum) influence the distribution of the chemical in the test system. The model is then applied to set of 1194 real substances (predominantly from the ToxCast chemical database) to calculate the potential range of concentrations and chemical activities under assumed test conditions. The model demonstrates how both concentrations and chemical activities can vary by orders of magnitude for the same nominal concentration.

The role of drug metabolizing enzymes in clearance

INTRODUCTION

Metabolism is one of the most important clearance pathways representing the major clearance route of 75% drugs. The four most common drug metabolizing enzymes (DME) that contribute significantly to elimination pathways of new chemical entities are cytochrome P450s, UDP-glucuronosyltransferases, aldehyde oxidase and sulfotransferases. Accurate prediction of human in vivo clearance by these enzymes, using both in vitro and in vivo tools, is critical for the success of drug candidates in human translation.

AREAS COVERED

Important recent advances of key DME are reviewed and highlighted in the following areas: major isoforms, tissue distribution, generic polymorphism, substrate specificity, species differences, mechanism of catalysis, in vitro-in vivo extrapolation and the importance of using optimal assay conditions and relevant animal models.

EXPERT OPINION

Understanding the clearance mechanism of a compound is the first step toward successful prediction of human clearance. It is critical to apply appropriate in vitro and in vivo methodologies and physiologically based models in human translation. While high-confidence prediction for P450-mediated clearance has been achieved, the accuracy of human clearance prediction is significantly lower for other enzyme classes. More accurate predictive methods and models are being developed to address these challenges.