Predicting blood:air partition coefficients using basic physicochemical properties

A systems toxicology approach for identification of disruptions in cholesterol homeostasis after aggregated exposure to mixtures of perfluorinated compounds in humans

Per- and polyfluoroalkyl substances (PFAS) are used in various household and industrial products. In humans, positive associations were reported between PFAS, including perfluorsulfonic acid and perfluorooctanoic acid, and cholesterol, a cardiometabolic risk factor. Animal studies show the opposite. Human-centered approaches are needed to better understand the effects of PFAS mixtures on cholesterol. Here, a systems toxicology approach is described, using a gene-centered cholesterol biokinetic model. PFAS exposure-gene expression relations from published data were introduced into the model. An existing PFAS physiologically based kinetic model was augmented with lung and dermal compartments and integrated with the cholesterol model to enable exposure-effect modeling. The final model was populated with data reflecting lifetime mixture exposure from: tolerable weekly intake values; the environment; high occupational exposures (ski waxing, PFAS industry). Results indicate that low level exposures (tolerable weekly intake, environmental) did not change cholesterol. In contrast, occupational exposures clearly resulted in internal PFAS exposure and disruption of cholesterol homeostasis, largely in line with epidemiological observations. Despite model limitations (eg, dynamic range, directionality), changes in cholesterol homeostasis were predicted for ski waxers, hitherto unknown from epidemiological studies. Here, future studies involving lipid metabolism could improve risk assessment.

A revision of the multiple-path particle dosimetry model focusing on tobacco product aerosol dynamics

To assess the health impact of inhaled aerosols, it is necessary to understand aerosol dynamics and the associated dosimetry in the human respiratory tract. Although several studies have measured or simulated the dosimetry of aerosol constituents, the respiratory tract focus areas have been limited. In particular, the aerosols generated from tobacco products are complex composites and simulating their dynamics in the respiratory tract is challenging. To assess the dosimetry of the aerosol constituents of tobacco products, we developed a revised version of the Multiple-Path Particle Dosimetry (MPPD) model, which employs (1) new geometry based on CT-scanned human respiratory tract data, (2) convective mixing in the oral cavity and deep lung, and (3) constituent partitioning between the tissue and air, and clearance. The sensitivity analysis was conducted using aerosols composed of four major constituents of electronic cigarette (EC) aerosols to investigate the parameters that have a significant impact on the results. In addition, the revised model was run with 4 and 10 constituents in ECs and conventional cigarettes (CCs), respectively. Sensitivity analysis revealed that the new modeling and the physicochemical properties of constituents had a considerable impact on the simulated aerosol concentration and dosimetry. The simulations could be carried out within 3 min even when 10 constituents of CC aerosols were analyzed simultaneously. The revised model based on MPPD is an efficient and easy-to-use tool for understanding the aerosol dynamics of CC and EC constituents and their effect on the human body.

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.

Exposure Prioritization (Ex Priori): A Screening-Level High-Throughput Chemical Prioritization Tool

To estimate potential chemical risk, tools are needed to prioritize potential exposures for chemicals with minimal data. Consumer product exposures are a key pathway, and variability in consumer use patterns is an important factor. We designed Ex Priori, a flexible dashboard-type screening-level exposure model, to rapidly visualize exposure rankings from consumer product use. Ex Priori is Excel-based. Currently, it is parameterized for seven routes of exposure for 1108 chemicals present in 228 consumer product types. It includes toxicokinetics considerations to estimate body burden. It includes a simple framework for rapid modeling of broad changes in consumer use patterns by product category. Ex Priori rapidly models changes in consumer user patterns during the COVID-19 pandemic and instantly shows resulting changes in chemical exposure rankings by body burden. Sensitivity analysis indicates that the model is sensitive to the air emissions rate of chemicals from products. Ex Priori's simple dashboard facilitates dynamic exploration of the effects of varying consumer product use patterns on prioritization of chemicals based on potential exposures. Ex Priori can be a useful modeling and visualization tool to both novice and experienced exposure modelers and complement more computationally intensive population-based exposure models.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Comparison of in vivo derived and scaled in vitro metabolic rate constants for several volatile organic compounds (VOCs)

Metabolic rate parameters estimation using in vitro data is necessary due to numbers of chemicals for which data are needed, trend towards minimizing laboratory animal use, and limited opportunity to collect data in human subjects. We evaluated how well metabolic rate parameters derived from in vitro data predict overall in vivo metabolism for a set of environmental chemicals for which well validated and established methods exist. We compared values of VmaxC derived from in vivo vapor uptake studies with estimates of VmaxC scaled up from in vitro hepatic microsomal metabolism studies for VOCs for which data were available in male F344 rats. For 6 of 7 VOCs, differences between the in vivo and scaled up in vitro VmaxC estimates were less than 2.6-fold. For bromodichloromethane (BDCM), the in vivo derived VmaxC was approximately 4.4-fold higher than the in vitro derived and scaled up VmaxC. The more rapid rate of BDCM metabolism estimated based in vivo studies suggests other factors such as extrahepatic metabolism, binding or other non-specific losses making a significant contribution to overall clearance. Systematic and reliable utilization of scaled up in vitro biotransformation rate parameters in PBPK models will require development of methods to predict cases in which extrahepatic metabolism and binding as well as other factors are likely to be significant contributors.

Topical review on monitoring tetrahydrocannabinol in breath

Legalization of cannabis for recreational use has compelled governments to seek new tools to accurately monitor Δ9-tetrahydrocannabinol (Δ9-THC) and understand its effect on impairment. Various methods have been employed to measure Δ9-THC, and its respective metabolites, in different biological matrices. Recently, breath analysis has gained interest as a non-invasive method for the detection of chemicals that are either produced as part of biological processes or are absorbed from the environment. Existing breath analyzers function by analyzing previously collected samples or by direct real-time analysis. Portable hand-held devices are of particular interest for law enforcement and personal use. This paper reviews and compares both commercially available and prototype devices that proclaim Δ9-THC detection in exhaled breath using methods such as Field Asymmetric Ion Mobility Spectrometry, Semiconductor-Enriched Single-Walled Carbon Nanotube chemiresistors, Liquid Chromatography Tandem-mass Spectrometry, microfluidic-based artificial olfaction, and optical-based gas sensing.

Predicting estrogen receptor binding of chemicals using a suite of in silico methods - Complementary approaches of (Q)SAR, molecular docking and molecular dynamics

With the aim of obtaining reliable estimates of Estrogen Receptor (ER) binding for diverse classes of compounds, a weight of evidence approach using estimates from a suite of in silico models was assessed. The predictivity of a simple Majority Consensus of (Q)SAR models was assessed using a test set of compounds with experimental Relative Binding Affinity (RBA) data. Molecular docking was also carried out and the binding energies of these compounds to the ERα receptor were determined. For a few selected compounds, including a known full agonist and antagonist, the intrinsic activity was determined using low-mode molecular dynamics methods. Individual (Q)SAR model predictivity varied, as expected, with some models showing high sensitivity, others higher specificity. However, the Majority Consensus (Q)SAR prediction showed a high accuracy and reasonably balanced sensitivity and specificity. Molecular docking provided quantitative information on strength of binding to the ERα receptor. For the 50 highest binding affinity compounds with positive RBA experimental values, just 5 of them were predicted to be non-binders by the Majority QSAR Consensus. Furthermore, agonist-specific assay experimental values for these 5 compounds were negative, which indicates that they may be ER antagonists. We also showed different scenarios of combining (Q)SAR results with Molecular docking classification of ER binding based on cut-off values of binding energies, providing a rational combined strategy to maximize terms of toxicological interest.

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.

Critical Review of Volatile Organic Compound Analysis in Breath and In Vitro Cell Culture for Detection of Lung Cancer

Breath analysis is a promising technique for lung cancer screening. Despite the rapid development of breathomics in the last four decades, no consistent, robust, and validated volatile organic compound (VOC) signature for lung cancer has been identified. This review summarizes the identified VOC biomarkers from both exhaled breath analysis and in vitro cultured lung cell lines. Both clinical and in vitro studies have produced inconsistent, and even contradictory, results. Methodological issues that lead to these inconsistencies are reviewed and discussed in detail. Recommendations on addressing specific issues for more accurate biomarker studies have also been made.

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.

Derivation of Occupational Thresholds of Toxicological Concern for Systemically Acting Noncarcinogenic Organic Chemicals

Many substances in workplace do not have occupational exposure limits. The threshold of toxicological concern (TTC) principle is part of the hierarchy of approaches useful in occupational health risk assessment. The aim of this study was to derive occupational TTCs (OTTCs) reflecting the airborne concentrations below which no significant risk to workers would be anticipated. A reference dataset consisting of the 8-h threshold limit values-Time-Weighted Average for 280 organic substances was compiled. Each substance was classified into low (class I), intermediate (class II), or high (class III) hazard categories as per Cramer rules. For each chemical, n-octanol:water partition coefficient and vapor pressure along with the molecular weight were used to predict the blood:air partition coefficient. The blood:air partition coefficient along with data on water solubility and ventilation rate allowed the prediction of pulmonary retention factor and absorbed dose in workers. For each Cramer class, the distribution of the predicted doses was analyzed to identify the various percentile values corresponding to the OTTC. Accordingly, for Cramer classes I-III, the OTTCs derived in this study correspond to 0.15, 0.0085, and 0.006 mmol/d, respectively, at the 10th percentile level, while these values were 1.5, 0.09 and 0.03 mmol/d at the 25th percentile level. The proposed OTTCs are not meant to replace the traditional occupational exposure limits, but can be used in data-poor situations along with exposure estimates to support screening level risk assessment and prioritization.

Modeling Human Exposure to Indoor Contaminants: External Source to Body Tissues

Information on human indoor exposure is necessary to assess the potential risk to individuals from many chemicals of interest. Dynamic indoor and human physicologically based pharmacokinetic (PBPK) models of the distribution of nonionizing, organic chemical concentrations in indoor environments resulting in delivered tissue doses are developed, described and tested. The Indoor model successfully reproduced independently measured, reported time-dependent air concentrations of chloroform released during showering and of 2-butyoxyethanol following use of a volatile surface cleaner. The Indoor model predictions were also comparable to those from a higher tier consumer model (ConsExpo 4.1) for the surface cleaner scenario. The PBPK model successful reproduced observed chloroform exhaled air concentrations resulting from an inhalation exposure. Fugacity based modeling provided a seamless description of the partitioning, fluxes, accumulation and release of the chemical in indoor media and tissues of the exposed subject. This has the potential to assist in health risk assessments, provided that appropriate physical/chemical property, usage characteristics, and toxicological information are available.

The acute effects of daily nicotine intake on heart rate--a toxicokinetic and toxicodynamic modelling study

Joint physiologically-based toxicokinetic and toxicodynamic (PBTK/TD) modelling was applied to simulate concentration-time profiles of nicotine, a well-known stimulant, in the human body following single and repeated dosing. Both kinetic and dynamic models were first calibrated by using in vivo literature data for the Caucasian population. The models were then used to estimate the blood and liver concentrations of nicotine in terms of the Area Under Curve (AUC) and the peak concentration (Cmax) for selected exposure scenarios based on inhalation (cigarette smoking), oral intake (nicotine lozenges) and dermal absorption (nicotine patches). The model simulations indicated that whereas frequent cigarette smoking gives rise to high AUC and Cmax in blood, the use of nicotine-rich dermal patches leads to high AUC and Cmax in the liver. Venous blood concentrations were used to estimate one of the most common acute effects, mean heart rate, both at rest and during exercise. These estimations showed that cigarette smoking causes a high peak heart rate, whereas dermal absorption causes a high mean heart rate over 48h. This study illustrates the potential of using PBTK/TD modelling in the safety assessment of nicotine-containing products.

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.