A unified algorithm for predicting partition coefficients for PBPK modeling of drugs and environmental chemicals

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.

First-in-Human Predictions of Hepatic Clearance for Drugs With the Well-Stirred Model: Comparative Assessment Between Models of Fraction Unbound Based Either on the Free Drug Hypothesis, Albumin-Facilitated Hepatic Uptake or Dynamic Binding Kinetics

The well-stirred model (WSM) is commonly used to predict the hepatic clearance in vivo (CLH) of drugs. The necessary intrinsic clearance of the unbound drug (CLint-in vitro-unbound) is generated in the in vitro assays in the presence of microsomes or hepatocytes but in the absence of plasma proteins. The value of CLint-in vitro-unbound can be extrapolated with the fraction unbound determined in vitro in plasma (fup) only if the fraction unbound in vivo in liver is the same. However, this approach resulted to a systematic underprediction bias of CLH. With the goal of reducing this bias, two new models of fraction unbound were published in this journal. These models estimate the binding kinetics of the rates of association and dissociation of the drug-protein complex and propose that more dissociation in the liver compared to plasma will increase the fraction unbound available for the metabolism. Consequently, these two models generated higher values of fraction unbound, implying a lower underprediction bias of CLH with the WSM. The first model was developed by Poulin et al. and is referring to the value of fup that is adjusted (fu-adjusted) to quantify the effect of a full dissociation of the drug-protein complex at the hepatocyte membrane in accordance with the theory of the albumin-facilitated hepatic uptake. A second model was developed by Yan et al. who presented a dynamic fraction unbound (fu-dynamic) measuring the real dissociation kinetics of the drug-protein complex with a new in vitro assay in the presence and absence of a recombinant liver enzyme in plasma. Therefore, the objective of this study was to make the first comparative assessment between these two models. The results indicate that, in general, the WSM combined with the values of fu-adjusted was the most accurate approach for predicting CLH. The WSM combined with the values of fu-dynamic has underperformed particularly with the acidic and neutral drugs binding to the albumin and presenting a low metabolic turnover in vitro. Therefore, the new in vitro assay for fu-dynamic gave an underprediction bias of CLH for these drug properties. However, the values of fu-adjusted are significantly higher than those values of fu-dynamic, and, this resulted to no underprediction bias, which is reinforcing the theory of the ALB-facilitated hepatic uptake. For the other neutral and acidic drugs, the models of fu-dynamic and fu-adjusted are in closer agreement. Finally, for the basic drugs, the models of fu-adjusted and fu-dynamic as well as a third model only considering a pH gradient effect on fup are almost accurately equivalent.

Multimodal Model to Predict Tissue-to-Blood Partition Coefficients of Chemicals in Mammals and Fish

Tissue-to-blood partition coefficients (Ptb) are key parameters for assessing toxicokinetics of xenobiotics in organisms, yet their experimental data were lacking. Experimental methods for measuring Ptb values are inefficient, underscoring the urgent need for prediction models. However, most existing models failed to fully exploit Ptb data from diverse sources, and their applicability domain (AD) was limited. The current study developed a multimodal model capable of processing and integrating textual (categorical features) and numerical information (molecular descriptors/fingerprints) to simultaneously predict Ptb values across various species, tissues, blood matrices, and measurement methods. Artificial neural network algorithms with embedding layers were used for the multimodal modeling. The corresponding unimodal models were developed for comparison. Results showed that the multimodal model outperformed unimodal models. To enhance the reliability of the model, a method considering categorical features, weighted molecular similarity density, and weighted inconsistency in molecular activities of structure-activity landscapes was used to characterize the AD. The model constrained by the AD exhibited better prediction accuracy for the validation set, with the determination coefficient, root mean-square error, and mean absolute error being 0.843, 0.276, and 0.213 log units, respectively. The multimodal model coupled with the AD characterization can serve as an efficient tool for internal exposure assessment of chemicals.

Multi-task machine learning models for simultaneous prediction of tissue-to-blood partition coefficients of chemicals in mammals

Tissue-to-blood partition coefficients (Ptb) are crucial for assessing the distribution of chemicals in organisms. Given the lack of experimental data and laborious nature of experimental methods, there is an urgent need to develop efficient predictive models. With the help of machine learning algorithms, i,e., random forest (RF), and artificial neural network (ANN), this study developed multi-task (MT) models that can simultaneously predict Ptb values for various mammalian tissues, including liver, muscle, brain, lung, and adipose. Single-task (ST) models using partial least squares regression, RF, and ANN algorithms for each endpoint were established for comparison. Overall, the performances of MT models were superior to those of ST models. The MT model using ANN algorithms showed the highest prediction accuracy with determination coefficients ranging from 0.704 to 0.886, root mean square errors between 0.223 and 0.410, and mean absolute errors ranging from 0.178 to 0.285 log units. Results showed that lipophilicity and polarizability of molecules significantly influence their partition behavior in organisms. Applicability domains (ADs) of the models were characterized by weighted molecular similarity density, and weighted inconsistency in molecular activities of structure-activity landscapes. When constrained by ADs, the models displayed enhanced predictive accuracy, making them valuable tools for the risk assessment and management of chemicals.

A New Version of the Tissue Composition-Based Model for Improving the Mechanism-Based Prediction of Volume of Distribution at Steady-State for Neutral Drugs

In-vitro models are available in the literature for predicting the volume of distribution at steady-state (Vdss) of drugs. The mechanistic model refers to the tissue composition-based model (TCM), which includes important factors that govern Vdss such as drug physiochemistry and physiological data. The recognized TCM published by Rodgers and Rowland (TCM-RR) and a subsequent adjustment made by Simulations Plus Inc. (TCM-SP) have been shown to be generally less accurate with neutral compared to ionized drugs. Therefore, improving these models for neutral drugs becomes necessary. The objective of this study was to propose a new TCM for improving the prediction of Vdss for neutral drugs. The new TCM included two modifications of the published models (i) accentuate the effect of the blood-to-plasma ratio (BPR) that should cover permeated molecules across the biomembranes, which is lacking in these models for neutral compounds, and (ii) use a different approach to estimate the binding in tissues. The new TCM was validated with a large dataset of 202 commercial and proprietary compounds including preclinical and clinical data. All scenario datasets were predicted more accurately with the TCM-New, whereas all statistical parameters indicate that the TCM-New showed significant improvements in terms of accuracy over the TCM-RR and TCM-SP. Predictions of Vdss were frequently more accurate for the TCM-new with 83% within twofold error versus only 50% for the TCM-RR. And more than 95% of the predictions were within threefold error and patient interindividual differences can be predicted with the TCM-New, greatly exceeding the accuracy of the published models. Overall, the new TCM incorporating BPR significantly improved the Vdss predictions in animals and humans for neutral drugs, and, hence, has the potential to better support the drug discovery and facilitate the first-in-human predictions.

Co-PBK: a computational biomonitoring tool for assessing chronic internal exposure to chemicals and metabolites

Toxic chemicals are released into the environment through diverse human activities. An increasing number of chronic diseases are associated with ambient pollution, thus posing a threat to people. Given the high consumption of resources for human biomonitoring, this study proposed coupled physiologically-based kinetic (co-PBK) modeling matrices as a biomonitoring tool for simplifying chronic internal exposure estimates of environmental chemicals and their metabolites using naphthalene (NAP) and its metabolites (i.e., 1-OHN and 2-OHN) as simulation examples. According to the simulation of the steady-state mass among various organs/tissues via the co-PBK modeling matrices, fat had the highest potential bioaccumulation of NAP and its metabolites. With respect to body fluids, 1-OHN and 2-OHN tended to bioaccumulate more in the bile than in the urine. According to the sensitivity analysis, the calculated sensitivity factors for the first-order kinetics-based rate constants imply that due to the biotransformation process, target organs/tissues (e.g., liver and kidneys) would be continuously exposed to more NAP metabolites under chronic exposure. Meanwhile, 1-OHN may be more stably transported to the urine than 2-OHN for further human biomonitoring during long-term internal exposure. According to the case study of simulating population chronic exposure to NAP in Shenzhen, the co-PBK modeling estimated the population exposure to NAP with an intake rate of 8.77 × 10-2 mg d-1 and the aggregated urinary concentration of NAP metabolites of 2.60 μg L-1. Furthermore, the accuracy of the urinary levels between the real-world data and the values simulated by the co-PBK modeling was assessed and the root-mean-square error of c1-OHN,urine was found to be lower than that of c2-OHN,urine.

Determination of blood:air, urine:air and plasma:air partition coefficients of selected microbial volatile organic compounds

Partition coefficients (PCs) are essential parameters for understanding the toxicokinetics of chemicals in the human body since they are used in the description of different processes of absorption, distribution, and excretion in physiologically based pharmacokinetic (PBPK) models used in chemical exposure and risk assessment. The goal of this study was to determine urine:air, blood:air and plasma:air partition coefficients (PCs) of microbial volatile organic compounds (mVOCs) previously selected as having high potential as biomarkers of indoor mold exposure. To achieve this goal, the vial-equilibration technique was used, and quantification was performed using headspace gas chromatography tandem mass spectrometry (HS-GC-MS/MS) analysis. Matrix:air PCs of 19 different mVOCs have been successfully determined and their values ranged between 14 and 3586 for urine:air, 78 and 4721 for blood:air and 64 and 5604 for plasma:air PCs. Water:air PCs were also determined, and their values varied between 16 and 2210, showing a good correlation with urine:air PCs for 17 compounds of the selected mVOCs (R2 = 0.97, slope close to unity) indicating that water:air PCs below 103 may be a good surrogate for urine:air PCs. All studied mVOCs have high blood:air PCs (greater than 78) indicating strong pulmonary uptake. Due to their high blood:urine PCs, some mVOCs may be more easily measured in blood than in urine. This work is an important preliminary step toward the use of mVOCs as potential biomarkers of indoor mold exposure. The data obtained in this study will help to determine the most appropriate matrix to use in this biomonitoring approach and will eventually facilitate the development of PBPK models for these chemicals.

CETP inhibitor evacetrapib enters mouse brain tissue

High levels of plasma cholesterol, especially high levels of low-density lipoprotein cholesterol (LDL-C), have been associated with an increased risk of Alzheimer's disease. The cholesteryl ester transfer protein (CETP) in plasma distributes cholesteryl esters between lipoproteins and increases LDL-C in plasma. Epidemiologically, decreased CETP activity has been associated with sustained cognitive performance during aging, longevity, and a lower risk of Alzheimer's disease. Thus, pharmacological CETP inhibitors could be repurposed for the treatment of Alzheimer's disease as they are safe and effective at lowering CETP activity and LDL-C. Although CETP is mostly expressed by the liver and secreted into the bloodstream, it is also expressed by astrocytes in the brain. Therefore, it is important to determine whether CETP inhibitors can enter the brain. Here, we describe the pharmacokinetic parameters of the CETP inhibitor evacetrapib in the plasma, liver, and brain tissues of CETP transgenic mice. We show that evacetrapib crosses the blood-brain barrier and is detectable in brain tissue 0.5 h after a 40 mg/kg i.v. injection in a non-linear function. We conclude that evacetrapib may prove to be a good candidate to treat CETP-mediated cholesterol dysregulation in Alzheimer's disease.

Plasma Protein-Mediated Uptake and Contradictions to the Free Drug Hypothesis: A Critical Review

According to the free drug hypothesis (FDH), only free, unbound drug is available to interact with biological targets. This hypothesis is the fundamental principle that continues to explain the vast majority of all pharmacokinetic and pharmacodynamic processes. Under the FDH, the free drug concentration at the target site is considered the driver of pharmacodynamic activity and pharmacokinetic processes. However, deviations from the FDH are observed in hepatic uptake and clearance predictions, where observed unbound intrinsic hepatic clearance (CL_int,u) is larger than expected. Such deviations are commonly observed when plasma proteins are present and form the basis of the so-called plasma protein-mediated uptake effect (PMUE). This review will discuss the basis of plasma protein binding as it pertains to hepatic clearance based on the FDH, as well as several hypotheses that may explain the underlying mechanisms of PMUE. Notably, some, but not all, potential mechanisms remained aligned with the FDH. Finally, we will outline possible experimental strategies to elucidate PMUE mechanisms. Understanding the mechanisms of PMUE and its potential contribution to clearance underprediction is vital to improving the drug development process.

A "middle-out approach" for the prediction of human drug disposition from preclinical data using simplified physiologically based pharmacokinetic (PBPK) models

Simplified physiologically based pharmacokinetic (PBPK) models using estimated tissue-to-unbound plasma partition coefficients (Kpus) were previously investigated by fitting them to in vivo pharmacokinetic (PK) data. After optimization with preclinical data, the performance of these models for extrapolation of distribution kinetics to human were evaluated to determine the best approach for the prediction of human drug disposition and volume of distribution (Vss) using PBPK modeling. Three lipophilic bases were tested (diazepam, midazolam, and basmisanil) for which intravenous PK data were available in rat, monkey, and human. The models with Kpu scalars using k-means clustering were generally the best for fitting data in the preclinical species and gave plausible Kpu values. Extrapolations of plasma concentrations for diazepam and midazolam using these models and parameters obtained were consistent with the observed clinical data. For diazepam and midazolam, the human predictions of Vss after optimization in rats and monkeys were better compared with the Vss estimated from the traditional PBPK modeling approach (varying from 1.1 to 3.1 vs. 3.7-fold error). For basmisanil, the sparse preclinical data available could have affected the model performance for fitting and the subsequent extrapolation to human. Overall, this work provides a rational strategy to predict human drug distribution using preclinical PK data within the PBPK modeling strategy.

A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species

Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t_½) have been observed in some cases. Knowledge of chemical-specific t_½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t_½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t_½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t_½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t_½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.

Evaluation of a rapid, generic human gestational dose model

Chemical risk assessment considers potentially susceptible populations including pregnant women and developing fetuses. Humans encounter thousands of chemicals in their environments, few of which have been fully characterized. Toxicokinetic (TK) information is needed to relate chemical exposure to potentially bioactive tissue concentrations. Observational data describing human gestational exposures are unavailable for most chemicals, but physiologically based TK (PBTK) models estimate such exposures. Development of chemical-specific PBTK models requires considerable time and resources. As an alternative, generic PBTK approaches describe a standardized physiology and characterize chemicals with a set of standard physical and TK descriptors - primarily plasma protein binding and hepatic clearance. Here we report and evaluate a generic PBTK model of a human mother and developing fetus. We used a published set of formulas describing the major anatomical and physiological changes that occur during pregnancy to augment the High-Throughput Toxicokinetics (httk) software package. We simulated the ratio of concentrations in maternal and fetal plasma and compared to literature in vivo measurements. We evaluated the model with literature in vivo time-course measurements of maternal plasma concentrations in pregnant and non-pregnant women. Finally, we prioritized chemicals measured in maternal serum based on predicted fetal brain concentrations. This new model can be used for TK simulations of 859 chemicals with existing human-specific in vitro TK data as well as any new chemicals for which such data become available. This gestational model may allow for in vitro to in vivo extrapolation of point of departure doses relevant to reproductive and developmental toxicity.

A new headspace solid-phase microextraction coupled with gas chromatography-tandem mass spectrometry method for the simultaneous quantification of 21 microbial volatile organic compounds in urine and blood

Mold growth can cause the development of several metabolites including microbial volatile organic compounds (mVOCs). These latter may be considered as potential biomarkers of fungal presence and have been detected in human biological matrices such as urine and blood. Exposure to molds and their metabolites (e.g., mVOCs, mycotoxins) in occupational settings, is responsible for several health effects. Thus, this exposure cannot be neglected and must be evaluated. Herein, a method has been developed to quantify 21 mVOCs in urine and human blood by headspace solid phase micro-extraction (HS-SPME) coupled with gas chromatography-triple quadrupole mass spectrometry (GC-MS/MS). The parameters influencing the extraction process, such as the type of fiber, the incubation and extraction time and temperature and the desorption time, have been optimized to ensure better mVOCs extraction. The developed method showed good linearity over the concentration range of the compounds (R² ˃ 0.995) for all the mVOCs in all the matrices. The low limits of detection (LOD) ranging from 0.7 to 417 ng/L in urine and from 1 to 507 ng/L in blood, make the developed methods sensitive and effective for biomonitoring of exposure at low levels. Recoveries, at low and high concentrations, were between 87% and 120% in urine and between 83% and 118% in blood. The repeatability and the intermediate precision in terms of coefficients of variation (CV%) was lower than 13% and 8.58% respectively for all compounds in all matrices. These values show satisfactory accuracy and precision of the developed method. Thus, this practical, simple, and sensitive method is well suited for the simultaneous quantification of target mVOCs.

Opportunities and challenges of physiologically based pharmacokinetic modeling in drug delivery

Physiologically based pharmacokinetic (PBPK) modeling is an important in silico tool to bridge drug properties and in vivo PK behaviors during drug development. Over the recent decade, the PBPK method has been largely applied to drug delivery systems (DDS), including oral, inhaled, transdermal, ophthalmic, and complex injectable products. The related therapeutic agents have included small-molecule drugs, therapeutic proteins, nucleic acids, and even cells. Simulation results have provided important insights into PK behaviors of new dosage forms, which strongly support drug regulation. In this review, we comprehensively summarize recent progress in PBPK applications in drug delivery, which shows large opportunities for facilitating drug development. In addition, we discuss the challenges of applying this methodology from a practical viewpoint.

Use of Physiologically-Based Kinetics Modelling to Reliably Predict Internal Concentrations of the UV Filter, Homosalate, After Repeated Oral and Topical Application

Ethical and legal considerations have led to increased use of non-animal methods to evaluate the safety of chemicals for human use. We describe the development and qualification of a physiologically-based kinetics (PBK) model for the cosmetic UV filter ingredient, homosalate, to support its safety without the need of generating further animal data. The intravenous (IV) rat PBK model, using PK-Sim^®, was developed and validated using legacy in vivo data generated prior to the 2013 EU animal-testing ban. Input data included literature or predicted physicochemical and pharmacokinetic properties. The refined IV rat PBK model was subject to sensitivity analysis to identify homosalate-specific sensitive parameters impacting the prediction of C_max (more sensitive than AUC_(0-∞)). These were then considered, together with population modeling, to calculate the confidence interval (CI) 95% C_max and AUC_(0-∞). Final model parameters were established by visual inspection of the simulations and biological plausibility. The IV rat model was extrapolated to oral administration, and used to estimate internal exposures to doses tested in an oral repeated dose toxicity study. Next, a human PBK dermal model was developed using measured human in vitro ADME data and a module to represent the dermal route. Model performance was confirmed by comparing predicted and measured values from a US-FDA clinical trial (Identifier: NCT03582215, https://clinicaltrials.gov/). Final exposure estimations were obtained in a virtual population and considering the in vitro and input parameter uncertainty. This model was then used to estimate the C_max and AUC_(0–24 h) of homosalate according to consumer use in a sunscreen. The developed rat and human PBK models had a good biological basis and reproduced in vivo legacy rat and human clinical kinetics data. They also complied with the most recent WHO and OECD recommendations for assessing the confidence level. In conclusion, we have developed a PBK model which predicted reasonably well the internal exposure of homosalate according to different exposure scenarios with a medium to high level of confidence. In the absence of in vivo data, such human PBK models will be the heart of future completely non-animal risk assessments; therefore, valid approaches will be key in gaining their regulatory acceptance.

Clinical Trial Registration: https://clinicaltrials.gov/, identifier, NCT03582215

PBPK Modeling to Simulate the Fate of Compounds in Living Organisms

Pharmacokinetics study the fate of xenobiotics in a living organism. Physiologically based pharmacokinetic (PBPK) models provide realistic descriptions of xenobiotics' absorption, distribution, metabolism, and excretion processes. They model the body as a set of homogeneous compartments representing organs, and their parameters refer to anatomical, physiological, biochemical, and physicochemical entities. They offer a quantitative mechanistic framework to understand and simulate the time-course of the concentration of a substance in various organs and body fluids. These models are well suited for performing extrapolations inherent to toxicology and pharmacology (e.g., between species or doses) and for integrating data obtained from various sources (e.g., in vitro or in vivo experiments, structure-activity models). In this chapter, we describe the practical development and basic use of a PBPK model from model building to model simulations, through implementation with an easily accessible free software.

High-throughput PBTK models for in vitro to in vivo extrapolation

INTRODUCTION

Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals.

AREAS COVERED

This review covers generic physiologically based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on 'httk', a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data.

EXPERT OPINION

HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.

Evaluation and comparison of in vitro intrinsic clearance rates measured using cryopreserved hepatocytes from humans, rats, and rainbow trout

In vitro and in silico methods that can reduce the need for animal testing are being used with increasing frequency to assess chemical risks to human health and the environment. The rate of hepatic biotransformation is an important species-specific parameter for determining bioaccumulation potential and extrapolating in vitro bioactivity to in vivo effects. One approach to estimating hepatic biotransformation is to employ in vitro systems derived from liver tissue to measure chemical (substrate) depletion over time which can then be translated to a rate of intrinsic clearance (CL_int). In the present study, cryopreserved hepatocytes from humans, rats, and rainbow trout were used to measure CL_int values for 54 industrial and pesticidal chemicals at starting test concentrations of 0.1 and 1 μM. A data evaluation framework that emphasizes the behavior of Heat-Treated Controls (HTC) was developed to identify datasets suitable for rate reporting. Measured or estimated ("greater than" or "less than") CL_int values were determined for 124 of 226 (55 %) species-chemical-substrate concentration datasets with acceptable analytical chemistry. A large percentage of tested chemicals exhibited low HTC recovery values, indicating a substantial abiotic loss of test chemical over time. An evaluation of K_OW values for individual chemicals suggested that in vitro test performance declined with increasing chemical hydrophobicity, although differences in testing devices for mammals and fish also likely played a role. The current findings emphasize the value of negative controls as part of a rigorous approach to data quality assessment for in vitro substrate depletion studies. Changes in current testing protocols can be expected to result in the collection of higher quality data. However, poorly soluble chemicals are likely to remain a challenge for CL_int determination.

Olanzapine Mesoporous Nanostructured Lipid Carrier: Optimization, Characterization, In Vivo Assessment, and Physiologically Based Pharmacokinetic Modeling

A promising approach has been emerging to enhance dissolution of hydrophobicdrugsby encapsulation in mesoporous silica materials. Olanzapine is a practically insoluble antipsychotic drug which is subjected to excessive first pass effect and shows inadequate oral bioavailability. Therefore, mesoporous silica was used to improve bioavailability of olanzapine incorporated in nano-structured lipid carriers (NLCs). These systems were characterized for their particle size, polydispersity index (PDI), zeta potential, entrapment efficiency (EE) and differential scanning calorimetry (DSC) as well asits release profile. The optimized mesoporous NLC system displayed nano-spherical particles (120.56 nm), possessed high entrapment efficiency (88.46%) and the highest percentage of drug released after six hours (75.13%). The biological performance of the optimized system was assessed in comparison with the drug suspension in healthy albino rabbits. The optimized system showed significantly (P < 0.05) prolonged MRT (8.47 h), higher C_max (22.12± 0.40 ng/ml) and T_max (2.0 h) values compared to drug suspension. Physiologically based pharmacokinetic (PBPK) model was simulated and verified. All the predicted results were within 0.6 and 1-fold of the reported data. To set a conclusion, in vitro results as well as in vivo pharmacokinetic study and PBPK data showed an enhancement in bioavailability of the optimized NLCs system over the plain drug suspension. These results proved the potentiality of incorporating olanzapine in mesoporous NLC for a significant improvement in oral bioavailability of olanzapine.

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.

Systematic Analysis of the Whole-Body Tissue Distribution and Fatty Acid Compositions of Membrane Lipids in CD1 and NMRI Mice and Wistar Rats

Understanding the tissue distribution of phospholipids and glycerolipids in animal models enables promoting the pharmacokinetic study of drugs and related PK predictions. The measurement of lipid compositions in animal models, usually mice and rats, without a standardized approach hindered the accuracy of PBPK investigation. In this work, high resolution mass spectrometry was applied to profile the tissue distribution of phospholipids and glycerolipids in 12 organs/tissues of mice and rats. Using this method, not only the amounts of phospholipids and glycerolipids in each organ/tissue but also the fatty acid compositions were acquired. In order to explore the interspecies specificity of lipid distribution in different organs/tissues, three animal species including CD1 mice, NMRI mice, and Wister rats were used in this systematic study. Globally, more organ specificity was observed. It was found that the brain is the organ containing the most abundant phosphatidylserine lipids (PSs) in all three animal models, leading to brain tissues having the most concentrated acidic phospholipids. Diverse fatty acid compositions in each lipid class were clearly revealed. Certain tissues/organs also had a specific selection of unique fatty acid compositions, for example, unreferenced FA(18 : 2) in the brain. It turned out that the access of free fatty acids affects the incorporation of acyl chain in phospholipids and glycerolipids. In the analysis, ether lipids were also profiled with the observation of dominant ePEs in brain tissues. However, little interspecies difference was found for fatty acid constituents and tissues distribution of phospholipids and glycerolipids.

Current Approaches and Techniques in Physiologically Based Pharmacokinetic (PBPK) Modelling of Nanomaterials

There have been efforts to develop physiologically based pharmacokinetic (PBPK) models for nanomaterials (NMs). Since NMs have quite different kinetic behaviors, the applicability of the approaches and techniques that are utilized in current PBPK models for NMs is warranted. Most PBPK models simulate a size-independent endocytosis from tissues or blood. In the lungs, dosimetry and the air-liquid interface (ALI) models have sometimes been used to estimate NM deposition and translocation into the circulatory system. In the gastrointestinal (GI) tract, kinetics data are needed for mechanistic understanding of NM behavior as well as their absorption through GI mucus and their subsequent hepatobiliary excretion into feces. Following absorption, permeability (P_t) and partition coefficients (PCs) are needed to simulate partitioning from the circulatory system into various organs. Furthermore, mechanistic modelling of organ- and species-specific NM corona formation is in its infancy. More recently, some PBPK models have included the mononuclear phagocyte system (MPS). Most notably, dissolution, a key elimination process for NMs, is only empirically added in some PBPK models. Nevertheless, despite the many challenges still present, there have been great advances in the development and application of PBPK models for hazard assessment and risk assessment of NMs.

Global Sensitivity Analysis of the Rodgers and Rowland Model for Prediction of Tissue: Plasma Partitioning Coefficients: Assessment of the Key Physiological and Physicochemical Factors That Determine Small-Molecule Tissue Distribution

In physiologically based pharmacokinetic (PBPK) modelling, the large number of input parameters, limited amount of available data and the structural model complexity generally hinder simultaneous estimation of uncertain and/or unknown parameters. These parameters are generally subject to estimation. However, the approaches taken for parameter estimation vary widely. Global sensitivity analyses are proposed as a method to systematically determine the most influential parameters that can be subject to estimation. Herein, a global sensitivity analysis was conducted to identify the key drug and physiological parameters influencing drug disposition in PBPK models and to potentially reduce the PBPK model dimensionality. The impact of these parameters was evaluated on the tissue-to-unbound plasma partition coefficients (Kpus) predicted by the Rodgers and Rowland model using Latin hypercube sampling combined to partial rank correlation coefficients (PRCC). For most drug classes, PRCC showed that LogP and fraction unbound in plasma (fu_p) were generally the most influential parameters for Kpu predictions. For strong bases, blood:plasma partitioning was one of the most influential parameter. Uncertainty in tissue composition parameters had a large impact on Kpu and Vss predictions for all classes. Among tissue composition parameters, changes in Kpu outputs were especially attributed to changes in tissue acidic phospholipid concentrations and extracellular protein tissue:plasma ratio values. In conclusion, this work demonstrates that for parameter estimation involving PBPK models and dimensionality reduction purposes, less influential parameters might be assigned fixed values depending on the parameter space, while influential parameters could be subject to parameters estimation.

Quantitative Predictions for Molecular Initiating Events Using Three-Dimensional Quantitative Structure-Activity Relationships

The aim of human toxicity risk assessment is to determine a safe dose or exposure to a chemical for humans. This requires an understanding of the exposure of a person to a chemical and how much of the chemical is required to cause an adverse effect. To do this computationally, we need to understand how much of a chemical is required to perturb normal biological function in an adverse outcome pathway (AOP). The molecular initiating event (MIE) is the first step in an adverse outcome pathway and can be considered as a chemical interaction between a chemical toxicant and a biological molecule. Key chemical characteristics can be identified and used to model the chemistry of these MIEs. In this study, we do just this by using chemical substructures to categorize chemicals and 3D quantitative structure-activity relationships (QSARs) based on comparative molecular field analysis (CoMFA) to calculate molecular activity. Models have been constructed across a variety of human biological targets, the glucocorticoid receptor, mu opioid receptor, cyclooxygenase-2 enzyme, human ether-à-go-go related gene channel, and dopamine transporter. These models tend to provide molecular activity estimation well within one log unit and electronic and steric fields that can be visualized to better understand the MIE and biological target of interest. The outputs of these fields can be used to identify key aspects of a chemical's chemistry which can be changed to reduce its ability to activate a given MIE. With this methodology, the quantitative chemical activity can be predicted for a wide variety of MIEs, which can feed into AOP-based chemical risk assessments, and understanding of the chemistry behind the MIE can be gained.

An in vitro assay and artificial intelligence approach to determine rate constants of nanomaterial-cell interactions

In vitro assays and simulation technologies are powerful methodologies that can inform scientists of nanomaterial (NM) distribution and fate in humans or pre-clinical species. For small molecules, less animal data is often needed because there are a multitude of in vitro screening tools and simulation-based approaches to quantify uptake and deliver data that makes extrapolation to in vivo studies feasible. Small molecule simulations work because these materials often diffuse quickly and partition after reaching equilibrium shortly after dosing, but this cannot be applied to NMs. NMs interact with cells through energy dependent pathways, often taking hours or days to become fully internalized within the cellular environment. In vitro screening tools must capture these phenomena so that cell simulations built on mechanism-based models can deliver relationships between exposure dose and mechanistic biology, that is biology representative of fundamental processes involved in NM transport by cells (e.g. membrane adsorption and subsequent internalization). Here, we developed, validated, and applied the FORECAST method, a combination of a calibrated fluorescence assay (CF) with an artificial intelligence-based cell simulation to quantify rates descriptive of the time-dependent mechanistic biological interactions between NMs and individual cells. This work is expected to provide a means of extrapolation to pre-clinical or human biodistribution with cellular level resolution for NMs starting only from in vitro data.

A framework for application of quantitative property-property relationships (QPPRs) in physiologically based pharmacokinetic (PBPK) models for high-throughput prediction of internal dose of inhaled organic chemicals

New generation of toxicological tests and assessment strategies require validated toxicokinetic data or models that are lacking for most chemicals. This study aimed at developing a quantitative property-property relationship (QPPR)-based human physiologically based pharmacokinetic (PBPK) modeling framework for high-throughput predictions of inhalation toxicokinetics of organic chemicals. A PBPK model was parameterized with QPPR-derived values for hepatic clearance (CL_h) and partition coefficients (P) [blood:air (P_ba) and tissue:air (P_ta) and tissue:blood (P_tb)]. The model was initially applied to an evaluation dataset of 40 organic chemicals in the applicability domain, and then to an expanded dataset of 249 organic chemicals from diverse chemical classes. 'Batch' analyses were performed for rapid assessments of hundreds of chemicals. The simulations of inhalation toxicokinetics following an 8-h exposure to 1 ppm of each chemical were successful. The mean ratios of their predicted-to-experimental values were within a factor of 1.36-2.36 for P_tb and 1.18 for CL_h, for 80% of the chemicals in the evaluation dataset. The predicted 24-h area under the venous blood concentration-time curve (AUC₂₄) values were within the predicted envelopes obtained while using experimental values of P_ba and considering either no or maximal hepatic extraction. The reliability analysis (based on combined sensitivity and uncertainty analyses) indicated that AUC₂₄ predictions for 55% of the expanded dataset were moderately to highly reliable, with 46% exhibiting highly reliable values. Overall, the modeling framework suggests that molecular structure and chemical properties can together be effectively used to obtain first-cut estimates of the toxicokinetics of data-poor organic chemicals for screening and prioritization purposes.

Effect of dose and exposure protocol on the toxicokinetics and first-pass elimination of trichloroethylene and 1,1,1-trichloroethane

Trichloroethylene (TCE) and 1,1,1-trichloroethane (TRI) are frequent contaminants of drinking water and of groundwater at hazardous waste sites. There is relatively little information on the target organ deposition of TRI, despite its ingestion and common occurrence in humans. An important aim of the study was to delineate and contrast the toxicokinetics (TK) and bioavailability (F) of TRI and its well metabolized congener, TCE. Blood profiles were obtained from male Sprague-Dawley rats given aqueous emulsions of 6 or 48 mg TRI/kg and 10 or 50 mg TCE/kg as an oral bolus (po) or by gastric infusion (gi) over 2 h. TCE exhibited nonlinear TK, with a disproportionate increase in AUC and decrease in clearance and F with increase in dose. TRI exhibited linear TK. F did not vary significantly with TRI dose or dosage regimen. F values were substantially higher for TRI than for the respective TCE groups. TRI was distributed widely to tissues of rats gavaged with 6 mg TRI/kg, with accumulation in fat. This experiment yielded tissue uptake and elimination profiles and in vivo tissue:blood partition coefficients (PCs). Finally, additional rats were given 10 mg/kg of TCE and TRI po, ia and iv, so that first-pass hepatic (FP_h) and pulmonary (FP_p) elimination could be measured directly. Total and FP_h elimination of TCE exceeded that of TRI. TRI, with its higher air:blood PC, exhibited the higher FP_p. TCE and TRI, despite several common physical and chemical properties resulting in similar absorption and systemic distribution, displayed dissimilar dosage and dose rate effects on their TK.

Rapid experimental measurements of physicochemical properties to inform models and testing

The structures and physicochemical properties of chemicals are important for determining their potential toxicological effects, toxicokinetics, and route(s) of exposure. These data are needed to prioritize the risk for thousands of environmental chemicals, but experimental values are often lacking. In an attempt to efficiently fill data gaps in physicochemical property information, we generated new data for 200 structurally diverse compounds, which were rigorously selected from the USEPA ToxCast chemical library, and whose structures are available within the Distributed Structure-Searchable Toxicity Database (DSSTox). This pilot study evaluated rapid experimental methods to determine five physicochemical properties, including the log of the octanol:water partition coefficient (known as log(Kow) or logP), vapor pressure, water solubility, Henry's law constant, and the acid dissociation constant (pKa). For most compounds, experiments were successful for at least one property; log(Kow) yielded the largest return (176 values). It was determined that 77 ToxPrint structural features were enriched in chemicals with at least one measurement failure, indicating which features may have played a role in rapid method failures. To gauge consistency with traditional measurement methods, the new measurements were compared with previous measurements (where available). Since quantitative structure-activity/property relationship (QSAR/QSPR) models are used to fill gaps in physicochemical property information, 5 suites of QSPRs were evaluated for their predictive ability and chemical coverage or applicability domain of new experimental measurements. The ability to have accurate measurements of these properties will facilitate better exposure predictions in two ways: 1) direct input of these experimental measurements into exposure models; and 2) construction of QSPRs with a wider applicability domain, as their predicted physicochemical values can be used to parameterize exposure models in the absence of experimental data.

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.

High-throughput in-silico prediction of ionization equilibria for pharmacokinetic modeling

Chemical ionization plays an important role in many aspects of pharmacokinetic (PK) processes such as protein binding, tissue partitioning, and apparent volume of distribution at steady state (Vdss). Here, estimates of ionization equilibrium constants (i.e., pKa) were analyzed for 8132 pharmaceuticals and 24,281 other compounds to which humans might be exposed in the environment. Results revealed broad differences in the ionization of pharmaceutical chemicals and chemicals with either near-field (in the home) or far-field sources. The utility of these high-throughput ionization predictions was evaluated via a case-study of predicted PK Vdss for 22 compounds monitored in the blood and serum of the U.S. population by the U.S. Centers for Disease Control and Prevention National Health and Nutrition Examination Survey (NHANES). The chemical distribution ratio between water and tissue was estimated using predicted ionization states characterized by pKa. Probability distributions corresponding to ionizable atom types (IATs) were then used to analyze the sensitivity of predicted Vdss on predicted pKa using Monte Carlo methods. 8 of the 22 compounds were predicted to be ionizable. For 5 of the 8 the predictions based upon ionization are significantly different from what would be predicted for a neutral compound. For all but one (foramsulfuron), the probability distribution of predicted Vdss generated by IAT sensitivity analysis spans both the neutral prediction and the prediction using ionization. As new data sets of chemical-specific information on metabolism and excretion for hundreds of chemicals are being made available (e.g., Wetmore et al., 2015), high-throughput methods for calculating Vdss and tissue-specific PK distribution coefficients will allow the rapid construction of PK models to provide context for both biomonitoring data and high-throughput toxicity screening studies such as Tox21 and ToxCast.

Modeling Pharmacokinetics

Pharmacokinetics is the study of the fate of xenobiotics in a living organism. Physiologically based pharmacokinetic (PBPK) models provide realistic descriptions of xenobiotics' absorption, distribution, metabolism, and excretion processes. They model the body as a set of homogeneous compartments representing organs, and their parameters refer to anatomical, physiological, biochemical, and physicochemical entities. They offer a quantitative mechanistic framework to understand and simulate the time-course of the concentration of a substance in various organs and body fluids. These models are well suited for performing extrapolations inherent to toxicology and pharmacology (e.g., between species or doses) and for integrating data obtained from various sources (e.g., in vitro or in vivo experiments, structure-activity models). In this chapter, we describe the practical development and basic use of a PBPK model from model building to model simulations, through implementation with an easily accessible free software.

In vitro to in vivo extrapolation for high throughput prioritization and decision making

In vitro chemical safety testing methods offer the potential for efficient and economical tools to provide relevant assessments of human health risk. To realize this potential, methods are needed to relate in vitro effects to in vivo responses, i.e., in vitro to in vivo extrapolation (IVIVE). Currently available IVIVE approaches need to be refined before they can be utilized for regulatory decision-making. To explore the capabilities and limitations of IVIVE within this context, the U.S. Environmental Protection Agency Office of Research and Development and the National Toxicology Program Interagency Center for the Evaluation of Alternative Toxicological Methods co-organized a workshop and webinar series. Here, we integrate content from the webinars and workshop to discuss activities and resources that would promote inclusion of IVIVE in regulatory decision-making. We discuss properties of models that successfully generate predictions of in vivo doses from effective in vitro concentration, including the experimental systems that provide input parameters for these models, areas of success, and areas for improvement to reduce model uncertainty. Finally, we provide case studies on the uses of IVIVE in safety assessments, which highlight the respective differences, information requirements, and outcomes across various approaches when applied for decision-making.

Evaluation and calibration of high-throughput predictions of chemical distribution to tissues

Toxicokinetics (TK) provides critical information for integrating chemical toxicity and exposure assessments in order to determine potential chemical risk (i.e., the margin between toxic doses and plausible exposures). For thousands of chemicals that are present in our environment, in vivo TK data are lacking. The publicly available R package "httk" (version 1.8, named for "high throughput TK") draws from a database of in vitro data and physico-chemical properties in order to run physiologically-based TK (PBTK) models for 553 compounds. The PBTK model parameters include tissue:plasma partition coefficients (Kp) which the httk software predicts using the model of Schmitt (Toxicol In Vitro 22 (2):457-467, 2008). In this paper we evaluated and modified httk predictions, and quantified confidence using in vivo literature data. We used 964 rat Kp measured by in vivo experiments for 143 compounds. Initially, predicted Kp were significantly larger than measured Kp for many lipophilic compounds (log10 octanol:water partition coefficient > 3). Hence the approach for predicting Kp was revised to account for possible deficiencies in the in vitro protein binding assay, and the method for predicting membrane affinity was revised. These changes yielded improvements ranging from a factor of 10 to nearly a factor of 10,000 for 83 Kp across 23 compounds with only 3 Kp worsening by more than a factor of 10. The vast majority (92%) of Kp were predicted within a factor of 10 of the measured value (overall root mean squared error of 0.59 on log10-transformed scale). After applying the adjustments, regressions were performed to calibrate and evaluate the predictions for 12 tissues. Predictions for some tissues (e.g., spleen, bone, gut, lung) were observed to be better than predictions for other tissues (e.g., skin, brain, fat), indicating that confidence in the application of in silico tools to predict chemical partitioning varies depending upon the tissues involved. Our calibrated model was then evaluated using a second data set of human in vivo measurements of volume of distribution (Vss) for 498 compounds reviewed by Obach et al. (Drug Metab Dispos 36(7):1385-1405, 2008). We found that calibration of the model improved performance: a regression of the measured values as a function of the predictions has a slope of 1.03, intercept of - 0.04, and R2 of 0.43. Through careful evaluation of predictive methods for chemical partitioning into tissues, we have improved and calibrated these methods and quantified confidence for TK predictions in humans and rats.

European regulatory perspective on pediatric physiologically based pharmacokinetic models

Physiologically based pharmacokinetic (PBPK) models offer a mechanistic understanding of the disposition of the drug in the body. When well informed and conducted, this may lead to more efficient clinical studies in the vulnerable pediatric population. A review of pediatric submissions to European regulatory authorities has shown a limited number of recent examples. The use of PBPK models to inform pediatric drug development is encouraged. It is, however, important to consider the confidence in the predictions. A qualification of the PBPK platform for the intended purpose should be performed, and a learn-and-confirm approach is recommended. Further research in this area is encouraged to inform important parameters, and to increase availability of pediatric data for model qualification.

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.

Prediction of Deoxypodophyllotoxin Disposition in Mouse, Rat, Monkey, and Dog by Physiologically Based Pharmacokinetic Model and the Extrapolation to Human

Deoxypodophyllotoxin (DPT) is a potential anti-tumor candidate prior to its clinical phase. The aim of the study was to develop a physiologically based pharmacokinetic (PBPK) model consisting of 13 tissue compartments to predict DPT disposition in mouse, rat, monkey, and dog based on in vitro and in silico inputs. Since large interspecies difference was found in unbound fraction of DPT in plasma, we assumed that K_t:pl,u (unbound tissue-to-plasma concentration ratio) was identical across species. The predictions of our model were then validated by in vivo data of corresponding preclinical species, along with visual predictive checks. Reasonable matches were found between observed and predicted plasma concentrations and pharmacokinetic parameters in all four animal species. The prediction in the related seven tissues of mouse was also desirable. We also attempted to predict human pharmacokinetic profile by both the developed PBPK model and interspecies allometric scaling across mouse, rat and monkey, while dog was excluded from the scaling. The two approaches reached similar results. We hope the study will help in the efficacy and safety assessment of DPT in future clinical studies and provide a reference to the preclinical screening of similar compounds by PBPK model.

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.

A simple physiologically based pharmacokinetic model evaluating the effect of anti-nicotine antibodies on nicotine disposition in the brains of rats and humans

Physiologically based pharmacokinetic (PBPK) modeling was applied to investigate the effects of anti-nicotine antibodies on nicotine disposition in the brains of rats and humans. Successful construction of both rat and human models was achieved by fitting model outputs to published nicotine concentration time course data in the blood and in the brain. Key parameters presumed to have the most effect on the ability of these antibodies to prevent nicotine from entering the brain were selected for investigation using the human model. These parameters, which included antibody affinity for nicotine, antibody cross-reactivity with cotinine, and antibody concentration, were broken down into different, clinically-derived in silico treatment levels and fed into the human PBPK model. Model predictions suggested that all three parameters, in addition to smoking status, have a sizable impact on anti-nicotine antibodies' ability to prevent nicotine from entering the brain and that the antibodies elicited by current human vaccines do not have sufficient binding characteristics to reduce brain nicotine concentrations. If the antibody binding characteristics achieved in animal studies can similarly be achieved in human studies, however, nicotine vaccine efficacy in terms of brain nicotine concentration reduction is predicted to meet threshold values for alleviating nicotine dependence.

Equilibrium Sorption of Structurally Diverse Organic Ions to Bovine Serum Albumin

Reliable partitioning data are essential for assessing the bioaccumulation potential and the toxicity of chemicals. In contrast to neutral organic chemicals, the partitioning behavior of ionogenic organic chemicals (IOCs) is still a black box for environmental scientists. Partitioning to serum albumin, the major protein in blood plasma, strongly influences the freely dissolved concentration of many chemicals (including IOCs), which affects their transport and distribution in the body. Because consistent data sets for partitioning of IOCs are rarely available, bovine serum albumin-water partition coefficients (KBSA/w) were measured in this study for 45 anionic and 4 cationic organic chemicals, including various substituted benzoic and naphthoic acids, sulfonates and several pesticides and pharmaceuticals. The results of this study suggest that binding to BSA is substantially influenced by the three-dimensional structure of the chemicals and the position of substitutions on the sorbing molecules. For example, we found a difference of >1.5 log units between isomeric chemicals such as 3,4-dichlorobenzoic acid and 2,6-dichlorobenzoic acid, and 1-naphthoic acid and 2-naphthoic acid. Conventional modeling approaches (e.g., based on octanol-water partition coefficients) poorly predict log KBSA/w of organic ions (R(2) ≤ 0.5), partially because they do not capture the observed steric effects. Hence, alternative modeling strategies will be required for accurate prediction of serum albumin-water partition coefficients of organic ions.