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

Towards best use and regulatory acceptance of generic physiologically based kinetic (PBK) models for in vitro-to-in vivo extrapolation (IVIVE) in chemical risk assessment

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

Case study on the impact of the source of metabolism parameters in next generation physiologically based pharmacokinetic models: Implications for occupational exposures to trimethylbenzenes

Physiologically based pharmacokinetic (PBPK) models are a means of making important linkages between exposure assessment and in vitro toxicity. A key constraint on rapid application of PBPK models in risk assessment is traditional reliance on substance-specific in vivo toxicokinetic data to evaluate model quality. Bounding conditions, in silico, in vitro, and chemical read-across approaches have been proposed as alternative sources for metabolic clearance estimates. A case study to test consistency of predictive ability across these approaches was conducted using trimethylbenzenes (TMB) as prototype chemicals. Substantial concordance was found among TMB isomers with respect to accuracy (or inaccuracy) of approaches to estimating metabolism; for example, the bounding conditions never reproduced the human in vivo toxicokinetic data within two-fold. Using only approaches that gave acceptable prediction of in vivo toxicokinetics for the source compound (1,2,4-TMB) substantially narrowed the range of plausible internal doses for a given external dose for occupational, emergency response, and environmental/community health risk assessment scenarios for TMB isomers. Thus, risk assessments developed using the target compound models with a constrained subset of metabolism estimates (determined for source chemical models) can be used with greater confidence that internal dosimetry will be estimated with accuracy sufficient for the purpose at hand.

Development and evaluation of a high throughput inhalation model for organic chemicals

Currently it is difficult to prospectively estimate human toxicokinetics (particularly for novel chemicals) in a high-throughput manner. The R software package httk has been developed, in part, to address this deficiency, and the aim of this investigation was to develop a generalized inhalation model for httk. The structure of the inhalation model was developed from two previously published physiologically based models from Jongeneelen and Berge (Ann Occup Hyg 55:841-864, 2011) and Clewell et al. (Toxicol Sci 63:160-172, 2001), while calculated physicochemical data was obtained from EPA's CompTox Chemicals Dashboard. In total, 142 exposure scenarios across 41 volatile organic chemicals were modeled and compared to published data. The slope of the regression line of best fit between log-transformed simulated and observed blood and exhaled breath concentrations was 0.46 with an r2 = 0.45 and a root mean square error (RMSE) of direct comparison between the log-transformed simulated and observed values of 1.11. Approximately 5.1% (n = 108) of the data points analyzed were >2 orders of magnitude different than expected. The volatile organic chemicals examined in this investigation represent small, generally lipophilic molecules. Ultimately this paper details a generalized inhalation component that integrates with the httk physiologically based toxicokinetic model to provide high-throughput estimates of inhalation chemical exposures.

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.

Derivation of internal dose-based thresholds of toxicological concern for occupational inhalation exposure to systemically acting organic chemicals

This study aimed at deriving occupational thresholds of toxicological concern for inhalation exposure to systemically-acting organic chemicals using predicted internal doses. The latter were also used to evaluate the quantitative relationship between occupational exposure limit and internal dose. Three internal dose measures were identified for investigation: (i) the daily area under the venous blood concentration vs. time curve, (ii) the daily rate of the amount of parent chemical metabolized, and (iii) the maximum venous blood concentration at the end of an 8-hr work shift. A dataset of 276 organic chemicals with 8-hr threshold limit values-time-weighted average was compiled along with their molecular structure and Cramer classes (Class I: low toxicity, Class II: intermediate toxicity, Class III: suggestive of significant toxicity). Using a human physiologically-based pharmacokinetic model, the three identified dose metrics were predicted for an 8-hr occupational inhalation exposure to the threshold limit value for each chemical. Distributional analyses of the predicted dose metrics were performed to identify the percentile values corresponding to the occupational thresholds of toxicological concern. Also, simple linear regression analyses were performed to evaluate the relationship between the 8-hr threshold limit value and each of the predicted dose metrics, respectively. No threshold of toxicological concern could be derived for class II due to few chemicals. Based on the daily rate of the amount of parent chemical metabolized, the proposed internal dose-based occupational thresholds of toxicological concern were 5.61 × 10^-2 and 9 × 10^-4 mmol/d at the 10th percentile level for classes I and III, respectively, while they were 4.55 × 10^-1 and 8.5 × 10^-3 mmol/d at the 25th percentile level. Even though high and significant correlations were observed between the 8-hr threshold limit values and the predicted dose metrics, the one with the rate of the amount of chemical metabolized was remarkable regardless of the Cramer class (r² = 0.81; n = 276). The proposed internal dose-based occupational thresholds of toxicological concern are potentially useful for screening-level assessments as well as prioritization within an integrated occupational risk assessment framework.