An integrated QSAR-PBPK modelling approach for predicting the inhalation toxicokinetics of mixtures of volatile organic chemicals in the rat

Novel Models for Accurate Estimation of Air-Blood Partitioning: Applications to Individual Compounds and Complex Mixtures of Neutral Organic Compounds

The air-blood partition coefficient (Kab) is extensively employed in human health risk assessment for chemical exposure. However, current Kab estimation approaches either require an extensive number of parameters or lack precision. In this study, we present two novel and parsimonious models to accurately estimate Kab values for individual neutral organic compounds, as well as their complex mixtures. The first model, termed the GC×GC model, was developed based on the retention times of nonpolar chemical analytes on comprehensive two-dimensional gas chromatography (GC×GC). This model is unique in its ability to estimate the Kab values for complex mixtures of nonpolar organic chemicals. The GC×GC model successfully accounted for the Kab variance (R2 = 0.97) and demonstrated strong prediction power (RMSE = 0.31 log unit) for an independent set of nonpolar chemical analytes. Overall, the GC×GC model can be used to estimate Kab values for complex mixtures of neutral organic compounds. The second model, termed the partition model (PM), is based on two types of partition coefficients: octanol to water (Kow) and air to water (Kaw). The PM was able to effectively account for the variability in Kab data (n = 344), yielding an R2 value of 0.93 and root-mean-square error (RMSE) of 0.34 log unit. The predictive power and explanatory performance of the PM were found to be comparable to those of the parameter-intensive Abraham solvation models (ASMs). Additionally, the PM can be integrated into the software EPI Suite, which is widely used in chemical risk assessment for initial screening. The PM provides quick and reliable estimation of Kab compared to ASMs, while the GC×GC model is uniquely suited for estimating Kab values for complex mixtures of neutral organic compounds. In summary, our study introduces two novel and parsimonious models for the accurate estimation of Kab values for both individual compounds and complex mixtures.

The Application of a Physiologically Based Toxicokinetic Model in Health Risk Assessment

Toxicokinetics plays a crucial role in the health risk assessments of xenobiotics. Classical compartmental models are limited in their ability to determine chemical concentrations in specific organs or tissues, particularly target organs or tissues, and their limited interspecific and exposure route extrapolation hinders satisfactory health risk assessment. In contrast, physiologically based toxicokinetic (PBTK) models quantitatively describe the absorption, distribution, metabolism, and excretion of chemicals across various exposure routes and doses in organisms, establishing correlations with toxic effects. Consequently, PBTK models serve as potent tools for extrapolation and provide a theoretical foundation for health risk assessment and management. This review outlines the construction and application of PBTK models in health risk assessment while analyzing their limitations and future perspectives.

Machine learning and artificial intelligence in physiologically based pharmacokinetic modeling

Physiologically based pharmacokinetic (PBPK) models are useful tools in drug development and risk assessment of environmental chemicals. PBPK model development requires the collection of species-specific physiological, and chemical-specific absorption, distribution, metabolism, and excretion (ADME) parameters, which can be a time-consuming and expensive process. This raises a need to create computational models capable of predicting input parameter values for PBPK models, especially for new compounds. In this review, we summarize an emerging paradigm for integrating PBPK modeling with machine learning (ML) or artificial intelligence (AI)-based computational methods. This paradigm includes 3 steps (1) obtain time-concentration PK data and/or ADME parameters from publicly available databases, (2) develop ML/AI-based approaches to predict ADME parameters, and (3) incorporate the ML/AI models into PBPK models to predict PK summary statistics (eg, area under the curve and maximum plasma concentration). We also discuss a neural network architecture "neural ordinary differential equation (Neural-ODE)" that is capable of providing better predictive capabilities than other ML methods when used to directly predict time-series PK profiles. In order to support applications of ML/AI methods for PBPK model development, several challenges should be addressed (1) as more data become available, it is important to expand the training set by including the structural diversity of compounds to improve the prediction accuracy of ML/AI models; (2) due to the black box nature of many ML models, lack of sufficient interpretability is a limitation; (3) Neural-ODE has great potential to be used to generate time-series PK profiles for new compounds with limited ADME information, but its application remains to be explored. Despite existing challenges, ML/AI approaches will continue to facilitate the efficient development of robust PBPK models for a large number of chemicals.

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.

Development of a generic lifelong physiologically based biokinetic model for exposome studies

The current study within the frame of the HEALS project aims at the development of a lifelong physiologically based biokinetic (PBBK) model for exposome studies. The aim was to deliver a comprehensive modelling framework for addressing a large chemical space. Towards this aim, the delivered model can easily adapt parameters from existing ad-hoc models or complete the missing compound specific parameters using advanced quantitative structure activity relationship (QSAR). All major human organs are included, as well as arterial, venous, and portal blood compartments. Xenobiotics and their metabolites are linked through the metabolizing tissues. This is mainly the liver, but also other sites of metabolism might be considered (intestine, brain, skin, placenta) based on the presence or not of the enzymes involved in the metabolism of the compound of interest. Each tissue is described by three mass balance equations for (a) red blood cells, (b) plasma and interstitial tissue and (c) cells respectively. The anthropometric parameters of the models are time dependent, so as to provide a lifetime internal dose assessment, as well as to describe the continuously changing physiology of the mother and the developing fetus. An additional component of flexibility is that the biokinetic processes that relate to metabolism are related with either Michaelis-Menten kinetics, as well as intrinsic clearance kinetics. The capability of the model is demonstrated in the assessment of internal exposure and the prediction of expected biomonitored levels in urine for three major compounds within the HEALS project, namely bisphenol A (BPA), Bis(2-ethylhexyl) phthalate (DEHP) and cadmium (Cd). The results indicated that the predicted urinary levels fit very well with the ones from human biomonitoring (HBM) studies; internal exposure to plasticizers is very low (in the range of ng/L), while internal exposure to Cd is in the range of μg/L.

Novel and Predictive QSAR Model for Steroidal and Nonsteroidal 5α- Reductase Type II Inhibitors

AIMS AND OBJECTIVE

In this study, a novel quantitative structure activity relationship (QSAR) model has been developed for inhibitors of human 5-alpha reductase type II, which are used to treat benign prostate hypertrophy (BPH).

METHODS

The dataset consisted of 113 compounds-mainly nonsteroidal-with known inhibitory concentration. Then 3D structures of compounds were optimized and molecular structure descriptors were calculated. The stepwise multiple linear regression was used to select descriptors encoding the inhibitory activity of the compounds. Multiple linear regression (MLR) was used to build up the linear QSAR model.

RESULTS

The results obtained revealed that the descriptors which best describe the activity were atom type electropological state, carbon type, radial distribution function (RDF), barysz matrix and molecular linear free energy relation. The suggested model could achieve satisfied square correlation coefficient of R² = 0.72, higher than of many previous studies, indicating its superiority. Rigid validation criteria were met using external data with Q² ˃ 0.5 and R² = 0.75, reflecting the predictive power of the model.

CONCLUSION

The QSAR model was applied for screening botanical components of herbal preparations used to treat BPH, and could predict the activity of some, among others, making reasonable attribution to the proposed effect of these preparations. Gamma tocopherol was found to be an active inhibitor, in consistence with many previous studies, anticipating the power of this model in the prediction of new candidate molecules and suggesting further investigations.

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

•

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

•

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

•

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

•

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

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

Physiology-based toxicokinetic modelling in the frame of the European Human Biomonitoring Initiative

Given the opportunities provided by internal dosimetry modelling in the interpretation of human biomonitoring (HBM) data, the assessment of the links between exposure to chemicals and observed HBM data can be effectively supported by PBTK modelling. This paper gives a comprehensive review of available human PBTK models for compounds selected as a priority by the European Human Biomonitoring Initiative (HBM4EU). We highlight their advantages and deficiencies and suggest steps for advanced internal dose modelling. The review of the available PBTK models highlighted the conceptual differences between older models compared to the ones developed recently, reflecting commensurate differences in research questions. Due to the lack of coordinated strategies for deriving useful biomonitoring data for toxicokinetic properties, significant problems in model parameterisation still remain; these are further increased by the lack of human toxicokinetic data due to ethics issues. Finally, questions arise as well as to the extent they are really representative of interindividual variability. QSARs for toxicokinetic properties is a complementary approach for PBTK model parameterisation, especially for data poor chemicals. This approach could be expanded to model chemico-biological interactions such as intestinal absorption and renal clearance; this could serve the development of more complex generic PBTK models that could be applied to newly derived chemicals. Another gap identified is the framework for mixture interaction terms among compounds that could eventually interact in metabolism. From the review it was concluded that efforts should be shifted toward the development of generic multi-compartmental and multi-route models, supported by targeted biomonitoring coupled with parameterisation by both QSAR approach and experimental (in-vivo and in-vitro) data for newly developed and data poor compounds.

Health risk assessment of trihalomethanes mixtures from daily water-related activities via multi-pathway exposure based on PBPK model

In this study, the concentrations of trihalomethanes (THMs) in tap water and direct drinking water were analyzed, and based on the human behavior patterns and building parameters, the concentrations of THMs in indoor air were simulated with the water-air concentration conversion model. In addition, concentrations of THMs in human tissues were predicted based on physiologically based pharmacokinetic (PBPK) model, and the health risk of THMs for participants were estimated. Furthermore, the carcinogenic risk of mixtures according to the method proposed by USEPA and PBPK model based method were calculated and compared. The concentrations of chloroform, bromodichloromethane, dibromochloromethane and bromoform in tap water were 11.28-16.21, 4.83-6.28, 0.81-1.32 and 0.08-0.21 μg/L, and those in direct drinking water were 3.29-6.88, 0.35-0.47, 0.03-0.08 and 0.04-0.08 μg/L, respectively. The results of water-air concentration conversion model demonstrated that pollutants in air had a strong correlation with water-related activities. Multi-pathway PBPK model showed that THMs concentrations in liver, kidney and richly perfused tissue were higher than those in other tissues. The results of risk assessment showed that the mean risk levels of mixtures were 1.69 × 10^-4 and 1.72 × 10^-4 calculated by the USEPA recommended method and PBPK based method, which seriously exceeded the acceptable level. TCM and BDCM were the major risk factors, and inhalation was the main exposure route of THMs.

Target and Tissue Selectivity Prediction by Integrated Mechanistic Pharmacokinetic-Target Binding and Quantitative Structure Activity Modeling

Selectivity is an important attribute of effective and safe drugs, and prediction of in vivo target and tissue selectivity would likely improve drug development success rates. However, a lack of understanding of the underlying (pharmacological) mechanisms and availability of directly applicable predictive methods complicates the prediction of selectivity. We explore the value of combining physiologically based pharmacokinetic (PBPK) modeling with quantitative structure-activity relationship (QSAR) modeling to predict the influence of the target dissociation constant (K _D) and the target dissociation rate constant on target and tissue selectivity. The K _D values of CB1 ligands in the ChEMBL database are predicted by QSAR random forest (RF) modeling for the CB1 receptor and known off-targets (TRPV1, mGlu5, 5-HT1a). Of these CB1 ligands, rimonabant, CP-55940, and Δ⁸-tetrahydrocanabinol, one of the active ingredients of cannabis, were selected for simulations of target occupancy for CB1, TRPV1, mGlu5, and 5-HT1a in three brain regions, to illustrate the principles of the combined PBPK-QSAR modeling. Our combined PBPK and target binding modeling demonstrated that the optimal values of the K _D and k _off for target and tissue selectivity were dependent on target concentration and tissue distribution kinetics. Interestingly, if the target concentration is high and the perfusion of the target site is low, the optimal K _D value is often not the lowest K _D value, suggesting that optimization towards high drug-target affinity can decrease the benefit-risk ratio. The presented integrative structure-pharmacokinetic-pharmacodynamic modeling provides an improved understanding of tissue and target selectivity.

Toluene and methylethylketone: effect of combined exposure on their metabolism in rat

1. Multiple exposures are ubiquitous in industrial environments. In this article, we highlight the risks faced by workers and complete the data available on the metabolic impact of a common mixture: toluene (TOL) and methylethylketone (MEK). 2. Rats were exposed by inhalation under controlled conditions either to each solvent individually, or to mixtures of the two. How the interaction between the two solvents affected their fate in the blood and brain, their main relevant urinary metabolites (o-cresol, benzylmercapturic acid for TOL and 2,3-butanediols for MEK) and their hepatic metabolism were investigated. 3. Although the cytochrome P450 concentration was unchanged, and the activities of CYP1A2 and CYP2E1 isoforms were not additively or synergistically induced by co-exposure, TOL metabolism was inhibited by the presence of MEK (and vice versa). Depending on the relative proportions of each compound in the mixture, this sometimes resulted in a large increase in blood and brain concentrations. Apart from extreme cases (unbalanced mixtures), the amount of o-cresol and benzylmercapturic acid (and to a lesser extent 2,3-butanediols) excreted were proportional to the blood solvent concentrations. 4. In a co-exposure context, ortho-cresol and benzylmercapturic acid can be used as urinary biomarkers in biomonitoring for employees to relatively accurately assess TOL exposure.

NanoEHS beyond Toxicity - Focusing on Biocorona

The first phase of environmental health and safety of nanomaterials (nanoEHS) studies has been mainly focused on evidence-based investigations that probe the impact of nanoparticles, nanomaterials and nano-enabled products on biological and ecological systems. The integration of multiple disciplines, including colloidal science, nanomaterial science, chemistry, toxicology/immunology and environmental science, is necessary to understand the implications of nanotechnology for both human health and the environment. While strides have been made in connecting the physicochemical properties of nanomaterials with their hazard potential in tiered models, fundamental understanding of nano-biomolecular interactions and their implications for nanoEHS is largely absent from the literature. Research on nano-biomolecular interactions within the context of natural systems not only provides important clues for deciphering nanotoxicity and nanoparticle-induced pathology, but also presents vast new opportunities for screening beneficial material properties and designing greener products from bottom up. This review highlights new opportunities concerning nano-biomolecular interactions beyond the scope of toxicity.

Informatics and Data Analytics to Support Exposome-Based Discovery for Public Health

The complexity of the human exposome-the totality of environmental exposures encountered from birth to death-motivates systematic, high-throughput approaches to discover new environmental determinants of disease. In this review, we describe the state of science in analyzing the human exposome and provide recommendations for the public health community to consider in dealing with analytic challenges of exposome-based biomedical research. We describe extant and novel analytic methods needed to associate the exposome with critical health outcomes and contextualize the data-centered challenges by drawing parallels to other research endeavors such as human genomics research. We discuss efforts for training scientists who can bridge public health, genomics, and biomedicine in informatics and statistics. If an exposome data ecosystem is brought to fruition, it will likely play a role as central as genomic science has had in molding the current and new generations of biomedical researchers, computational scientists, and public health research programs.

Application of Biologically Based Lumping To Investigate the Toxicokinetic Interactions of a Complex Gasoline Mixture

People are often exposed to complex mixtures of environmental chemicals such as gasoline, tobacco smoke, water contaminants, or food additives. We developed an approach that applies chemical lumping methods to complex mixtures, in this case gasoline, based on biologically relevant parameters used in physiologically based pharmacokinetic (PBPK) modeling. Inhalation exposures were performed with rats to evaluate the performance of our PBPK model and chemical lumping method. There were 109 chemicals identified and quantified in the vapor in the chamber. The time-course toxicokinetic profiles of 10 target chemicals were also determined from blood samples collected during and following the in vivo experiments. A general PBPK model was used to compare the experimental data to the simulated values of blood concentration for 10 target chemicals with various numbers of lumps, iteratively increasing from 0 to 99. Large reductions in simulation error were gained by incorporating enzymatic chemical interactions, in comparison to simulating the individual chemicals separately. The error was further reduced by lumping the 99 nontarget chemicals. The same biologically based lumping approach can be used to simplify any complex mixture with tens, hundreds, or thousands of constituents.

Using physiologically based pharmacokinetic models to estimate the health risk of mixtures of trihalomethanes from reclaimed water

To estimate the health risk of mixture of trihalomethanes (THMs) from reclaimed water during toilet flushing, the interaction-based Hazard Index (HI(interaction-based)) and the mixture carcinogenic risk (CRM) according to tissue dose were conducted through the integrated use of both the exposure concentrations model and the physiologically based pharmacokinetic (PBPK) model of THMs. Monte Carlo simulations were employed to implement the probabilistic risk analysis and sensitivity analysis. Nine samples were analyzed, which were collected from J Water Reclamation Plant (JWRP) in Tianjin of China. The results indicated that the mean HI(interaction-based) (=0.85) was lower than the acceptable risk level (=1). The probability that the HI(interaction-based) exceeded the acceptable risk level is 22.97%. For carcinogenic risk, the CRM ranges from 9.41×10(-7) to 3.54×10(-5), with a mean of 5.49×10(-6). Moreover, the probability of exceeding the acceptable risk level (1×10(-6)) is near 100%. And the values of HI(interaction-based) from sample no. 1, 5, and 7 exceeded 1, while the values of CRM for all samples exceeded 1×10(-6). Consequently, the reclaimed water used for flushing toilets should be paid more attention, though non-carcinogenic effect is relatively small. Furthermore, the concentrations of DBCM had greater impact on both the carcinogenic and non-carcinogenic risk based on sensitivity analysis.

Assessing human variability in kinetics for exposures to multiple environmental chemicals: a physiologically based pharmacokinetic modeling case study with dichloromethane, benzene, toluene, ethylbenzene, and m-xylene

The objective of this study was to compare the magnitude of interindividual variability in internal dose for inhalation exposure to single versus multiple chemicals. Physiologically based pharmacokinetic models for adults (AD), neonates (NEO), toddlers (TODD), and pregnant women (PW) were used to simulate inhalation exposure to "low" (RfC-like) or "high" (AEGL-like) air concentrations of benzene (Bz) or dichloromethane (DCM), along with various levels of toluene alone or toluene with ethylbenzene and xylene. Monte Carlo simulations were performed and distributions of relevant internal dose metrics of either Bz or DCM were computed. Area under the blood concentration of parent compound versus time curve (AUC)-based variability in AD, TODD, and PW rose for Bz when concomitant "low" exposure to mixtures of increasing complexities occurred (coefficient of variation (CV) = 16-24%, vs. 12-15% for Bz alone), but remained unchanged considering DCM. Conversely, AUC-based CV in NEO fell (15 to 5% for Bz; 12 to 6% for DCM). Comparable trends were observed considering production of metabolites (AMET), except for NEO's CYP2E1-mediated metabolites of Bz, where an increased CV was observed (20 to 71%). For "high" exposure scenarios, Cmax-based variability of Bz and DCM remained unchanged in AD and PW, but decreased in NEO (CV= 11-16% to 2-6%) and TODD (CV= 12-13% to 7-9%). Conversely, AMET-based variability for both substrates rose in every subpopulation. This study analyzed for the first time the impact of multiple exposures on interindividual variability in toxicokinetics. Evidence indicates that this impact depends upon chemical concentrations and biochemical properties, as well as the subpopulation and internal dose metrics considered.

PBTK modelling platforms and parameter estimation tools to enable animal-free risk assessment: recommendations from a joint EPAA--EURL ECVAM ADME workshop

Information on toxicokinetics is critical for animal-free human risk assessment. Human external exposure must be translated into human tissue doses and compared with in vitro actual cell exposure associated to effects (in vitro-in vivo comparison). Data on absorption, distribution, metabolism and excretion in humans (ADME) could be generated using in vitro and QSAR tools. Physiologically-based toxicokinetic (PBTK) computer modelling could serve to integrate disparate in vitro and in silico findings. However, there are only few freely-available PBTK platforms currently available. And although some ADME parameters can be reasonably estimated in vitro or in silico, important gaps exist. Examples include unknown or limited applicability domains and lack of (high-throughput) tools to measure penetration of barriers, partitioning between blood and tissues and metabolic clearance. This paper is based on a joint EPAA--EURL ECVAM expert meeting. It provides a state-of-the-art overview of the availability of PBTK platforms as well as the in vitro and in silico methods to parameterise basic (Tier 1) PBTK models. Five high-priority issues are presented that provide the prerequisites for wider use of non-animal based PBTK modelling for animal-free chemical risk assessment.

Perspectives for integrating human and environmental risk assessment and synergies with socio-economic analysis

For more than a decade, the integration of human and environmental risk assessment (RA) has become an attractive vision. At the same time, existing European regulations of chemical substances such as REACH (EC Regulation No. 1907/2006), the Plant Protection Products Regulation (EC regulation 1107/2009) and Biocide Regulation (EC Regulation 528/2012) continue to ask for sector-specific RAs, each of which have their individual information requirements regarding exposure and hazard data, and also use different methodologies for the ultimate risk quantification. In response to this difference between the vision for integration and the current scientific and regulatory practice, the present paper outlines five medium-term opportunities for integrating human and environmental RA, followed by detailed discussions of the associated major components and their state of the art. Current hazard assessment approaches are analyzed in terms of data availability and quality, and covering non-test tools, the integrated testing strategy (ITS) approach, the adverse outcome pathway (AOP) concept, methods for assessing uncertainty, and the issue of explicitly treating mixture toxicity. With respect to exposure, opportunities for integrating exposure assessment are discussed, taking into account the uncertainty, standardization and validation of exposure modeling as well as the availability of exposure data. A further focus is on ways to complement RA by a socio-economic assessment (SEA) in order to better inform about risk management options. In this way, the present analysis, developed as part of the EU FP7 project HEROIC, may contribute to paving the way for integrating, where useful and possible, human and environmental RA in a manner suitable for its coupling with SEA.

Quantitative Property-Property Relationship for Screening-Level Prediction of Intrinsic Clearance of Volatile Organic Chemicals in Rats and Its Integration within PBPK Models to Predict Inhalation Pharmacokinetics in Humans

The objectives of this study were (i) to develop a screening-level Quantitative property-property relationship (QPPR) for intrinsic clearance (CL_int) obtained from in vivo animal studies and (ii) to incorporate it with human physiology in a PBPK model for predicting the inhalation pharmacokinetics of VOCs. CL_int, calculated as the ratio of the in vivo V_max (μmol/h/kg bw rat) to the K_m (μM), was obtained for 26 VOCs from the literature. The QPPR model resulting from stepwise linear regression analysis passed the validation step (R² = 0.8; leave-one-out cross-validation Q² = 0.75) for CL_int normalized to the phospholipid (PL) affinity of the VOCs. The QPPR facilitated the calculation of CL_int (L PL/h/kg bw rat) from the input data on log P_ow, log blood: water PC and ionization potential. The predictions of the QPPR as lower and upper bounds of the 95% mean confidence intervals (LMCI and UMCI, resp.) were then integrated within a human PBPK model. The ratio of the maximum (using LMCI for CL_int) to minimum (using UMCI for CL_int) AUC predicted by the QPPR-PBPK model was 1.36 ± 0.4 and ranged from 1.06 (1,1-dichloroethylene) to 2.8 (isoprene). Overall, the integrated QPPR-PBPK modeling method developed in this study is a pragmatic way of characterizing the impact of the lack of knowledge of CL_int in predicting human pharmacokinetics of VOCs, as well as the impact of prediction uncertainty of CL_int on human pharmacokinetics of VOCs.

Bayesian Analysis of a Lipid-Based Physiologically Based Toxicokinetic Model for a Mixture of PCBs in Rats

A lipid-based physiologically based toxicokinetic (PBTK) model has been developed for a mixture of six polychlorinated biphenyls (PCBs) in rats. The aim of this study was to apply population Bayesian analysis to a lipid PBTK model, while incorporating an internal exposure-response model linking enzyme induction and metabolic rate. Lipid-based physiologically based toxicokinetic models are a subset of PBTK models that can simulate concentrations of highly lipophilic compounds in tissue lipids, without the need for partition coefficients. A hierarchical treatment of population metabolic parameters and a CYP450 induction model were incorporated into the lipid-based PBTK framework, and Markov-Chain Monte Carlo was applied to in vivo data. A mass balance of CYP1A and CYP2B in the liver was necessary to model PCB metabolism at high doses. The linked PBTK/induction model remained on a lipid basis and was capable of modeling PCB concentrations in multiple tissues for all dose levels and dose profiles.

Translational research to develop a human PBPK models tool kit-volatile organic compounds (VOCs)

Toxicity and exposure evaluations remain the two of the key components of human health assessment. While improvement in exposure assessment relies on a better understanding of human behavior patterns, toxicity assessment still relies to a great extent on animal toxicity testing and human epidemiological studies. Recent advances in computer modeling of the dose-response relationship and distribution of xenobiotics in humans to important target tissues have advanced our abilities to assess toxicity. In particular, physiologically based pharmacokinetic (PBPK) models are among the tools than can enhance toxicity assessment accuracy. Many PBPK models are available to the health assessor, but most are so difficult to use that health assessors rarely use them. To encourage their use these models need to have transparent and user-friendly formats. To this end the Agency for Toxic Substances and Disease Registry (ATSDR) is using translational research to increase PBPK model accessibility, understandability, and use in the site-specific health assessment arena. The agency has initiated development of a human PBPK tool-kit for certain high priority pollutants. The tool kit comprises a series of suitable models. The models are recoded in a single computer simulation language and evaluated for use by health assessors. While not necessarily being state-of-the-art code for each chemical, the models will be sufficiently accurate to use for screening purposes. This article presents a generic, seven-compartment PBPK model for six priority volatile organic compounds (VOCs): benzene (BEN), carbon tetrachloride (CCl(4)), dichloromethane (DCM), perchloroethylene (PCE), trichloroethylene (TCE), and vinyl chloride (VC). Limited comparisons of the generic and original model predictions to published kinetic data were conducted. A goodness of fit was determined by calculating the means of the sum of the squared differences (MSSDs) for simulation vs. experimental kinetic data using the generic and original models. Using simplified solvent exposure assumptions for oral ingestion and inhalation, steady-state blood concentrations of each solvent were simulated for exposures equivalent to the ATSDR Minimal Risk Levels (MRLs). The predicted blood levels were then compared to those reported in the National Health and Nutrition Examination Survey (NHANES). With the notable exception of BEN, simulations of combined oral and inhalation MRLs using our generic VOC model yielded blood concentrations well above those reported for the 95th percentile blood concentrations for the U.S. populations, suggesting no health concerns. When the PBPK tool kit is fully developed, risk assessors will have a readily accessible tool for evaluating human exposure to a variety of environmental pollutants.

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

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