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

Upcyte Human Hepatocytes: a Potent In Vitro Tool for the Prediction of Hepatic Clearance of Metabolically Stable Compounds

In vitro models based on primary human hepatocytes (PHH) have been advanced for clearance (CL) prediction of metabolically stable compounds, representing state-of-the-art assay systems for drug discovery and development. Yet, limited cell availability and large interindividual variability of metabolic profiles remain shortcomings of PHH. Upcyte human hepatocytes (UHH) represent a novel hepatic cell system derived from PHH, exhibiting proliferative capacity for approximately 35 population doublings. UHH from three donors were evaluated during culture for up to 18 days, investigating relative mRNA expression and in situ enzyme activity of cytochrome P450s (P450s), UDP-glucuronosyltransferases, and sulfotransferases. Furthermore, UHH were used for predicting hepatic CL of 21 marketed low to intermediate CL drugs. In a typical experiment, expansion from 3.9 × 10(6) up to 8.5 × 10(7) cells was achieved during subculture. When maintained at confluence, transcripts of major P450s were expressed at donor-specific levels with sustained activities for the majority of isoforms, showing generally low CYP1A2 and high CYP2B6 activity levels. For donor 151-03, CL prediction based on depletion experiments resulted in an average fold error of 2.0, and 80% of compounds being predicted within twofold to in vivo CL for a subset of 10 low CL drugs. UHH showed sustained and consistent activity of drug-metabolizing enzymes (DME), resulting in highly reproducible CL prediction performance. In conclusion, UHH show promising potential as alternative to PHH for standardized in vitro applications in discovery research based on their stable, hepatocyte-like DME phenotype and virtually unlimited cell availability.

Content and activity of human liver microsomal protein and prediction of individual hepatic clearance in vivo

The lack of information concerning individual variation in content and activity of human liver microsomal protein is one of the most important obstacles for designing personalized medicines. We demonstrated that the mean value of microsomal protein per gram of liver (MPPGL) was 39.46 mg/g in 128 human livers and up to 19-fold individual variations existed. Meanwhile, the metabolic activities of 10 cytochrome P450 (CYPs) were detected in microsomes and liver tissues, respectively, which showed huge individual variations (200-fold). Compared with microsomes, the activities of liver tissues were much suitable to express the individual variations of CYP activities. Furthermore, individual variations in the in vivo clearance of tolbutamide were successfully predicted with the individual parameter values. In conclusion, we offer the values for MPPGL contents in normal liver tissues and build a new method to assess the in vitro CYP activities. In addition, large individual variations exist in predicted hepatic clearance of tolbutamide. These findings provide important physiological parameters for physiologically-based pharmacokinetics models and thus, establish a solid foundation for future development of personalized medicines.

In Silico Modeling for the Prediction of Dose and Pathway-Related Adverse Effects in Humans From In Vitro Repeated-Dose Studies

Long-term repeated-dose toxicity is mainly assessed in animals despite poor concordance of animal data with human toxicity. Nowadays advanced human in vitro systems, eg, metabolically competent HepaRG cells, are used for toxicity screening. Extrapolation of in vitro toxicity to in vivo effects is possible by reverse dosimetry using pharmacokinetic modeling. We assessed long-term repeated-dose toxicity of bosentan and valproic acid (VPA) in HepaRG cells under serum-free conditions. Upon 28-day exposure, the EC50 values for bosentan and VPA decreased by 21- and 33-fold, respectively. Using EC(10) as lowest threshold of toxicity in vitro, we estimated the oral equivalent doses for both test compounds using a simplified pharmacokinetic model for the extrapolation of in vitro toxicity to in vivo effect. The model predicts that bosentan is safe at the considered dose under the assumed conditions upon 4 weeks exposure. For VPA, hepatotoxicity is predicted for 4% and 47% of the virtual population at the maximum recommended daily dose after 3 and 4 weeks of exposure, respectively. We also investigated the changes in the central carbon metabolism of HepaRG cells exposed to orally bioavailable concentrations of both drugs. These concentrations are below the 28-day EC(10) and induce significant changes especially in glucose metabolism and urea production. These metabolic changes may have a pronounced impact in susceptible patients such as those with compromised liver function and urea cycle deficiency leading to idiosyncratic toxicity. We show that the combination of modeling based on in vitro repeated-dose data and metabolic changes allows the prediction of human relevant in vivo toxicity with mechanistic insights.

Microengineered cell and tissue systems for drug screening and toxicology applications: Evolution of in-vitro liver technologies

The liver performs many key functions, the most prominent of which is serving as the metabolic hub of the body. For this reason, the liver is the focal point of many investigations aimed at understanding an organism's toxicological response to endogenous and exogenous challenges. Because so many drug failures have involved direct liver toxicity or other organ toxicity from liver generated metabolites, the pharmaceutical industry has constantly sought superior, predictive in-vitro models that can more quickly and efficiently identify problematic drug candidates before they incur major development costs, and certainly before they are released to the public. In this broad review, we present a survey and critical comparison of in-vitro liver technologies along a broad spectrum, but focus on the current renewed push to develop "organs-on-a-chip". One prominent set of conclusions from this review is that while a large body of recent work has steered the field towards an ever more comprehensive understanding of what is needed, the field remains in great need of several key advances, including establishment of standard characterization methods, enhanced technologies that mimic the in-vivo cellular environment, and better computational approaches to bridge the gap between the in-vitro and in-vivo results.

In vitro metabolism of the lignan (-)-grandisin, an anticancer drug candidate, by human liver microsomes

(-)-grandisin is a tetrahydrofuran lignan that displays important biological properties, such as trypanocidal, anti-inflammatory, cytotoxic, and antitumor activities, suggesting its utility as a potential drug candidate. One important step in drug development is metabolic characterization and metabolite identification. To perform a biotransformation study of (-)-grandisin and to determine its kinetic properties in humans, a high performance liquid chromatography (HPLC) method was developed and validated. After HPLC method validation, the kinetic properties of (-)-grandisin were determined. (-)-grandisin metabolism obeyed Michaelis-Menten kinetics. The maximal reaction rate (Vmax ) was 3.96 ± 0.18 µmol/mg protein/h, and the Michaelis-Menten constant (Km ) was 8.23 ± 0.99 μM. In addition, the structures of the metabolites derived from (-)-grandisin were characterized via gas chromatography-mass spectrometry (GC-MS) and liquid chromatography-mass spectrometry (LC-MS) analysis. Four metabolites, 4-O-demethylgrandisin, 3-O-demethylgrandisin, 4,4'-di-O-demethylgrandisin, and a metabolite that may correspond to either 3,4-di-O-demethylgrandisin or 3,5-di-O-demethylgrandisin, were detected. CYP2C9 isoform was the main responsible for the formation of the metabolites. These metabolites have not been previously described, demonstrating the necessity of assessing (-)-grandisin metabolism using human-derived materials.

Analysis of the impact of controlled release formulations on oral drug absorption, gut wall metabolism and relative bioavailability of CYP3A substrates using a physiologically-based pharmacokinetic model

Controlled release (CR) formulations are usually designed to achieve similar exposure (AUC) levels as the marketed immediate release (IR) formulation. However, the AUC is often lower following CR compared to IR formulations. There are a few exceptions when the CR formulations have shown higher AUC. This study investigated the impact of CR formulations on oral drug absorption and CYP3A4-mediated gut wall metabolism. A review of the current literature on relative bioavailability (Frel) between CR and IR formulations of CYP3A substrates was conducted. This was followed by a systematic analysis to assess the impact of the release characteristics and the drug-specific factors (including metabolism and permeability) on oral bioavailability employing a physiologically-based pharmacokinetic (PBPK) modelling and simulation approach. From the literature review, only three CYP3A4 substrates showed higher Frel when formulated as CR. Several scenarios were investigated using the PBPK approach; in most of them, the oral absorption of CR formulations was lower as compared to the IR formulations. However, for highly permeable compounds that were CYP3A4 substrates the reduction in absorption was compensated by an increase in the fraction that escapes from first pass metabolism in the gut wall (FG), where the magnitude was dependent on CYP3A4 affinity. The systematic simulations of various interplays between different parameters demonstrated that BCS class 1 highly-cleared CYP3A4 substrates can display up to 220% higher relative bioavailability when formulated as CR compared to IR, in agreement with the observed data collected from the literature. The results and methodology of this study can be employed during the formulation development process in order to optimize drug absorption, especially for CYP3A4 substrates.

In vitro biotransformation rates in fish liver S9: effect of dosing techniques

In vitro biotransformation assays are currently being explored to improve estimates of bioconcentration factors of potentially bioaccumulative organic chemicals in fish. The present study compares thin-film and solvent-delivery dosing techniques as well as single versus multiple chemical dosing for measuring biotransformation rates of selected polycyclic aromatic hydrocarbons in rainbow trout (Oncorhynchus mykiss) liver S9. The findings show that biotransformation rates of very hydrophobic substances can be accurately measured in thin-film sorbent-dosing assays from concentration-time profiles in the incubation medium but not from those in the sorbent phase because of low chemical film-to-incubation-medium mass-transfer rates at the incubation temperature of 13.5 °C required for trout liver assays. Biotransformation rates determined by thin-film dosing were greater than those determined by solvent-delivery dosing for chrysene (octanol-water partition coefficient [KOW ] =10(5.60) ) and benzo[a]pyrene (KOW  =10(6.04) ), whereas there were no statistical differences in pyrene (KOW  =10(5.18) ) biotransformation rates between the 2 methods. In sorbent delivery-based assays, simultaneous multiple-chemical dosing produced biotransformation rates that were not statistically different from those measured in single-chemical dosing experiments for pyrene and benzo[a]pyrene but not for chrysene. In solvent-delivery experiments, multiple-chemical dosing produced biotransformation rates that were much smaller than those in single-chemical dosing experiments for all test chemicals. While thin-film sorbent-phase and solvent delivery-based dosing methods are both suitable methods for measuring biotransformation rates of substances of intermediate hydrophobicity, thin-film sorbent-phase dosing may be more suitable for superhydrophobic chemicals.

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.).

Enzyme-transporter interplay in the formation and clearance of abundant metabolites of faldaprevir found in excreta but not in circulation

Faldaprevir is a hepatitis C virus protease inhibitor that effectively reduces viral load in patients. Since faldaprevir exhibits slow metabolism in vitro and low clearance in vivo, metabolism was expected to be a minor clearance pathway. The human [(14)C] absorption, distribution, metabolism, and excretion study revealed that two monohydroxylated metabolites (M2a and M2b) were the most abundant excretory metabolites in feces, constituting 41% of the total administered dose. To deconvolute the formation and disposition of M2a and M2b in humans and determine why the minor change in structure [the addition of 16 atomic mass units (amu)] produced chemical entities that were excreted and were not present in the circulation, multiple in vitro test systems were used. The results from these in vitro studies clarified the formation and clearance of M2a and M2b. Faldaprevir is metabolized primarily in the liver by CYP3A4/5 to form M2a and M2b, which are also substrates of efflux transporters (P-glycoprotein and breast cancer resistance protein). The role of transporters is considered important for M2a and M2b as they demonstrate low permeability. It is proposed that both metabolites are efficiently excreted via bile into feces and do not enter the systemic circulation to an appreciable extent. If these metabolites permeate to blood, they can be readily taken up into hepatocytes from the circulation by uptake transporters (likely organic anion transporting polypeptides). These results highlight the critical role of drug-metabolizing enzymes and multiple transporters in the process of the formation and clearance of faldaprevir metabolites. Faldaprevir metabolism also provides an interesting case study for metabolites that are exclusively excreted in feces but are of clinical relevance.

Meeting the challenge of predicting hepatic clearance of compounds slowly metabolized by cytochrome P450 using a novel hepatocyte model, HepatoPac

Generating accurate in vitro intrinsic clearance data is an important aspect of predicting in vivo human clearance. Primary hepatocytes in suspension are routinely used to predict in vivo clearance; however, incubation times have typically been limited to 4-6 hours, which is not long enough to accurately evaluate the metabolic stability of slowly metabolized compounds. HepatoPac is a micropatterened hepatocyte-fibroblast coculture system that can be used for continuous incubations of up to 7 days. This study evaluated the ability of human HepatoPac to predict the in vivo clearance (CL) of 17 commercially available compounds with low to intermediate clearance (<12 ml/min per kg). In vitro half-life for disappearance of each compound was converted to hepatic clearance using the well stirred model, with and without correction for plasma protein binding. Hepatic CL, using three individual donors, was accurately predicted for 10 of 17 compounds (59%; predicted clearance within 2-fold of observed human in vivo clearance values). The accuracy of prediction increased to 76% (13 of 17 compounds) with an acceptance criterion defined as within 3-fold. When considering only low clearance compounds (<5 ml/min per kg), which represented 10 of the 17 compounds, the accuracy of prediction was 60% within 2-fold and 90% within 3-fold. In addition, the turnover of three slowly metabolized compounds (alprazolam, meloxicam, and tolbutamide) in HepatoPac was directly compared with turnover in suspended hepatocytes. The turnover of alprazolam and tolbutamide was approximately 2-fold greater using HepatoPac compared with suspended hepatocytes, which was roughly in line with the extrapolated values (correcting for the longer incubation time and lower cell number with HepatoPac). HepatoPac, but not suspended hepatocytes, demonstrated significant turnover of meloxicam. These results demonstrate the utility of HepatoPac for prediction of in vivo hepatic clearance, particularly with low clearance compounds.

Recent advances in 2D and 3D in vitro systems using primary hepatocytes, alternative hepatocyte sources and non-parenchymal liver cells and their use in investigating mechanisms of hepatotoxicity, cell signaling and ADME

This review encompasses the most important advances in liver functions and hepatotoxicity and analyzes which mechanisms can be studied in vitro. In a complex architecture of nested, zonated lobules, the liver consists of approximately 80 % hepatocytes and 20 % non-parenchymal cells, the latter being involved in a secondary phase that may dramatically aggravate the initial damage. Hepatotoxicity, as well as hepatic metabolism, is controlled by a set of nuclear receptors (including PXR, CAR, HNF-4α, FXR, LXR, SHP, VDR and PPAR) and signaling pathways. When isolating liver cells, some pathways are activated, e.g., the RAS/MEK/ERK pathway, whereas others are silenced (e.g. HNF-4α), resulting in up- and downregulation of hundreds of genes. An understanding of these changes is crucial for a correct interpretation of in vitro data. The possibilities and limitations of the most useful liver in vitro systems are summarized, including three-dimensional culture techniques, co-cultures with non-parenchymal cells, hepatospheres, precision cut liver slices and the isolated perfused liver. Also discussed is how closely hepatoma, stem cell and iPS cell-derived hepatocyte-like-cells resemble real hepatocytes. Finally, a summary is given of the state of the art of liver in vitro and mathematical modeling systems that are currently used in the pharmaceutical industry with an emphasis on drug metabolism, prediction of clearance, drug interaction, transporter studies and hepatotoxicity. One key message is that despite our enthusiasm for in vitro systems, we must never lose sight of the in vivo situation. Although hepatocytes have been isolated for decades, the hunt for relevant alternative systems has only just begun.

Lysosomal sequestration (trapping) of lipophilic amine (cationic amphiphilic) drugs in immortalized human hepatocytes (Fa2N-4 cells)

Lipophilic (logP > 1) and amphiphilic drugs (also known as cationic amphiphilic drugs) with ionizable amines (pKa > 6) can accumulate in lysosomes, a process known as lysosomal trapping. This process contributes to presystemic extraction by lysosome-rich organs (such as liver and lung), which, together with the binding of lipophilic amines to phospholipids, contributes to the large volume of distribution characteristic of numerous cardiovascular and central nervous system drugs. Accumulation of lipophilic amines in lysosomes has been implicated as a cause of phospholipidosis. Furthermore, elevated levels of lipophilic amines in lysosomes can lead to high organ-to-blood ratios of drugs that can be mistaken for active drug transport. In the present study, we describe an in vitro fluorescence-based method (using the lysosome-specific probe LysoTracker Red) to identify lysosomotropic agents in immortalized hepatocytes (Fa2N-4 cells). A diverse set of compounds with various physicochemical properties were tested, such as acids, bases, and zwitterions. In addition, the partitioning of the nonlysosomotropic atorvastatin (an anion) and the lysosomotropics propranolol and imipramine (cations) were quantified in Fa2N-4 cells in the presence or absence of various lysosomotropic or nonlysosomotropic agents and inhibitors of lysosomal sequestration (NH4Cl, nigericin, and monensin). Cellular partitioning of propranolol and imipramine was markedly reduced (by at least 40%) by NH4Cl, nigericin, or monensin. Lysosomotropic drugs also inhibited the partitioning of propranolol by at least 50%, with imipramine partitioning affected to a lesser degree. This study demonstrates the usefulness of immortalized hepatocytes (Fa2N-4 cells) for determining the lysosomal sequestration of lipophilic amines.

Prediction of hepatic plasma clearance in vivo from gel-entrapped rat and human hepatocytes

This paper aimed to evaluate the applicability of gel-entrapped rat and human hepatocytes in the prediction of hepatic plasma clearance (CLh,plasma) in vivo. The in vitro intrinsic clearances (CLint,in vitro) for the selected compounds were determined from the substrate disappearance rate, and further used to predict CLh,plasma using 3 classical mathematical models (well-stirred, parallel-tube, and dispersion) and disregarding drug binding. As a result, the predicted values from gel-entrapped rat hepatocytes were mostly within 2 SE of the literature data with a high correlation coefficient (R2) of 0.88–0.91. The predicted data with human hepatocytes also fitted well with the clinical data, indicating a high accuracy in prediction of in-vivo clearance. With respect to the mathematical model for predicting CLh,plasma, the parallel-tube and dispersion models produced a better prediction than the well-stirred model, and we suggest using the parallel-tube model because it is less complex mathematically. In conclusion, gel-entrapped hepatocytes predicted the drug clearance well and seemed to be a useful tool in the process of drug discovery.

Applications of minimal physiologically-based pharmacokinetic models

Conventional mammillary models are frequently used for pharmacokinetic (PK) analysis when only blood or plasma data are available. Such models depend on the quality of the drug disposition data and have vague biological features. An alternative minimal-physiologically-based PK (minimal-PBPK) modeling approach is proposed which inherits and lumps major physiologic attributes from whole-body PBPK models. The body and model are represented as actual blood and tissue (usually total body weight) volumes, fractions (f ( d )) of cardiac output with Fick's Law of Perfusion, tissue/blood partitioning (K ( p )), and systemic or intrinsic clearance. Analyzing only blood or plasma concentrations versus time, the minimal-PBPK models parsimoniously generate physiologically-relevant PK parameters which are more easily interpreted than those from mammillary models. The minimal-PBPK models were applied to four types of therapeutic agents and conditions. The models well captured the human PK profiles of 22 selected beta-lactam antibiotics allowing comparison of fitted and calculated K ( p ) values. Adding a classical hepatic compartment with hepatic blood flow allowed joint fitting of oral and intravenous (IV) data for four hepatic elimination drugs (dihydrocodeine, verapamil, repaglinide, midazolam) providing separate estimates of hepatic intrinsic clearance, non-hepatic clearance, and pre-hepatic bioavailability. The basic model was integrated with allometric scaling principles to simultaneously describe moxifloxacin PK in five species with common K ( p ) and f ( d ) values. A basic model assigning clearance to the tissue compartment well characterized plasma concentrations of six monoclonal antibodies in human subjects, providing good concordance of predictions with expected tissue kinetics. The proposed minimal-PBPK modeling approach offers an alternative and more rational basis for assessing PK than compartmental models.

The influence of old age and poloxamer-407 on the hepatic disposition of diazepam in the isolated perfused rat liver

BACKGROUND AND PURPOSE

The normal liver sinusoidal endothelium is thin and punctuated with fenestrations 50-200 nm in diameter that filter endobiotics and xenobiotics. Defenestration of the liver sinusoidal endothelium in old age and after pre-treatment with poloxamer-407 (P407) has been shown to prevent the transfer of small chylomicrons across the liver sinusoidal endothelium. The aim of this study was to investigate the impact of liver sinusoidal endothelium fenestrations on the hepatic uptake of the highly protein-bound drug diazepam. We hypothesized that defenestration will reduce the hepatic extraction of drugs which are highly bound to albumin.

METHODOLOGY

The isolated perfused rat liver (IPRL) model and multiple indicator dilution technique were used to investigate the effect of fenestrations in the liver sinusoidal endothelium on the hepatic disposition of diazepam in old and young rats, and in young rats treated with P407 or vehicle. A bolus dose of (14)C-diazpeam and non-extracted tracers ((3)H-sucrose and Evans blue) was injected into the portal vein. The single-pass recovery of diazepam and markers and the apparent volume of distribution were determined.

RESULTS

Scanning electron microscopy confirmed reduced porosity of the liver sinusoidal endothelial cells in P407-treated rats and old rats compared to young and control rats. The fractional recovery of diazepam was significantly increased in P407-treated rats compared to controls (0.20 ± 0.16, n = 12, P407; 0.08 ± 0.05, n = 8, controls; p = 0.0029), and in old rats compared to young rats (0.15 ± 0.03, n = 11, old; 0.10 ± 0.02, n = 11, young; p = 0.0004) following a single pass.

CONCLUSION

Defenestration due to age-related pseudocapillarization and treatment with P407 resulted in reduced hepatic extraction of diazepam after a single pass through the IPRL. These results highlight the importance of the liver sinusoidal endothelium in the ultrafiltration of highly protein-bound drugs, and may also provide an additional mechanism for reduced hepatic clearance of diazepam in conditions associated with defenestration.

Comparative assessment of In Vitro-In Vivo extrapolation methods used for predicting hepatic metabolic clearance of drugs

The purpose of this study was to perform a comparative analysis of various in vitro--in vivo extrapolation (IVIVE) methods used for predicting hepatic metabolic clearance (CL) of drugs on the basis of intrinsic CL data determined in microsomes. Five IVIVE methods were evaluated: the "conventional and conventional bias-corrected methods" using the unbound fraction in plasma (fu(p) ), the "Berezhkovskiy method" in which the fu(p) is adjusted for drug ionization, the "Poulin et al. method" using the unbound fraction in liver (fu(liver) ), and the "direct scaling method," which does not consider any binding corrections. We investigated the effects of the following scenarios on the prediction of CL: the use of preclinical or human datasets, the extent of plasma protein binding, the magnitude of CL in vivo, and the extent of drug disposition based on biopharmaceutics drug disposition classification system (BDDCS) categorization. A large and diverse dataset of 139 compounds was collected, including those from the literature and in house from Genentech. The results of this study confirm that the Poulin et al. method is robust and showed the greatest accuracy as compared with the other IVIVE methods in the majority of prediction scenarios studied here. The difference across the prediction methods is most pronounced for (a) albumin-bound drugs, (b) highly bound drugs, and (c) low CL drugs. Predictions of CL showed relevant interspecies differences for BDDCS class 2 compounds; the direct scaling method showed the greatest predictivity for these compounds, particularly for a reduced dataset in rat that have unexpectedly high CL in vivo. This result is a reflection of the direct scaling method's natural tendency to overpredict the true metabolic CL. Overall, this study should facilitate the use of IVIVE correlation methods in physiologically based pharmacokinetics (PBPK) model.

Cutting Edge PBPK Models and Analyses: Providing the Basis for Future Modeling Efforts and Bridges to Emerging Toxicology Paradigms

Physiologically based Pharmacokinetic (PBPK) models are used for predictions of internal or target dose from environmental and pharmacologic chemical exposures. Their use in human risk assessment is dependent on the nature of databases (animal or human) used to develop and test them, and includes extrapolations across species, experimental paradigms, and determination of variability of response within human populations. Integration of state-of-the science PBPK modeling with emerging computational toxicology models is critical for extrapolation between in vitro exposures, in vivo physiologic exposure, whole organism responses, and long-term health outcomes. This special issue contains papers that can provide the basis for future modeling efforts and provide bridges to emerging toxicology paradigms. In this overview paper, we present an overview of the field and introduction for these papers that includes discussions of model development, best practices, risk-assessment applications of PBPK models, and limitations and bridges of modeling approaches for future applications. Specifically, issues addressed include: (a) increased understanding of human variability of pharmacokinetics and pharmacodynamics in the population, (b) exploration of mode of action hypotheses (MOA), (c) application of biological modeling in the risk assessment of individual chemicals and chemical mixtures, and (d) identification and discussion of uncertainties in the modeling process.

Evaluation of hepatic clearance prediction using in vitro data: emphasis on fraction unbound in plasma and drug ionisation using a database of 107 drugs

Underprediction of in vivo intrinsic clearance (CL(int)) of unbound drug from human hepatic in vitro systems using physiological extrapolation methodology is accepted as a common outcome. Poulin et al. (2012. J Pharm Sci 101:838-851) recently proposed an approach involving determination of effective fraction unbound in plasma (fu(p)) based on albumin-facilitated hepatic uptake of acidic/neutral drugs which improved prediction accuracy and precision for 25 drugs highly bound to plasma proteins. This approach includes correction of unbound drug according to the ionisation fraction either side of the plasma membrane based on pH difference. Here, we assessed the proposed method using a larger database of predictions of CL(int) for 107 drugs involving hepatocytes (89 drugs) and microsomes (64 drugs). The proposed method was similarly effective in minimising average prediction bias (to within twofold), unlike the conventional fu(p) correction method. However, precision was similar between methods and there was no evidence in the larger database that prediction bias was associated with fu(p). Prediction bias for hepatocytes was clearance dependent by either method, indicating important sources of bias from in vitro methodology. Therefore, to progress beyond empirical correction of bias, there is further need of mechanistic elucidation to improve prediction methodology.