Simulation of urinary excretion of 1-hydroxypyrene in various scenarios of exposure to polycyclic aromatic hydrocarbons with a generic, cross-chemical predictive PBTK-model

Profiling population-wide exposure to environmental chemicals: A case study of naphthalene

Long-term exposure to environmental chemicals can detrimentally impact human health, and understanding the relationship between age distribution and levels of external and internal exposure is crucial. Nonetheless, existing methods for assessing population-wide exposure across age groups are limited. To bridge this research gap, we introduced a modeling approach designed to assess both chronic external and internal exposure to chemicals at the population level. The external and internal exposure assessments were quantified in terms of the average daily dose (ADD) and steady-state blood concentration of the environmental chemical, respectively, which were categorized by age and gender groups. The modeling process was presented within a spreadsheet framework, affording users the capability to execute population-wide exposure analyses across a spectrum of chemicals. Our simulation outcomes underscored a salient trend: younger age groups, particularly infants and children, exhibited markedly higher ADD values and blood concentrations of environmental chemicals compared to their older counterparts. This observation is due to the elevated basal metabolic rate per unit of body weight characteristic of younger individuals, coupled with their diminished biotransformation kinetics of xenobiotics within their livers. These factors collectively contribute to increased intake rates of environmental chemicals per unit of body weight through air and food consumption, along with heightened bioaccumulation of these chemicals within their bodies (e.g., blood). Furthermore, we augmented the precision of the external and internal exposure assessment by incorporating the age distribution across the population. The simulation outcomes unveiled that, to estimate the central tendency of the population's exposure levels, employing the baseline value group (age group 21-30) or the surrogate age of 25 serves as a simple yet dependable approach. However, for comprehensive population protection, our recommendation aligns with conducting exposure assessments for the younger age groups (age group 0-11). Future studies should integrate individual-level exposure assessment, analyze vulnerable population groups, and refine population structures within our developed model.

Novel Strategy to Assess the Neurotoxicity of Organic Solvents Such as Glycol Ethers: Protocol for Combining In Vitro and In Silico Methods With Human-Controlled Exposure Experiments

BACKGROUND

Chemicals are not required to be tested systematically for their neurotoxic potency, although they may contribute to the development of several neurological diseases. The absence of systematic testing may be partially explained by the current Organisation for Economic Co-operation and Development (OECD) Test Guidelines, which rely on animal experiments that are expensive, laborious, and ethically debatable. Therefore, it is important to understand the risks to exposed workers and the general population exposed to domestic products. In this study, we propose a strategy to test the neurotoxicity of solvents using the commonly used glycol ethers as a case study.

OBJECTIVE

This study aims to provide a strategy that can be used by regulatory agencies and industries to rank solvents according to their neurotoxicity and demonstrate the use of toxicokinetic modeling to predict air concentrations of solvents that are below the no observed adverse effect concentrations (NOAECs) for human neurotoxicity determined in in vitro assays.

METHODS

The proposed strategy focuses on a complex 3D in vitro brain model (BrainSpheres) derived from human-induced pluripotent stem cells (hiPSCs). This model is accompanied by in vivo, in vitro, and in silico models for the blood-brain barrier (BBB) and in vitro models for liver metabolism. The data are integrated into a toxicokinetic model. Internal concentrations predicted using this toxicokinetic model are compared with the results from in vivo human-controlled exposure experiments for model validation. The toxicokinetic model is then used in reverse dosimetry to predict air concentrations, leading to brain concentrations lower than the NOAECs determined in the hiPSC-derived 3D brain model. These predictions will contribute to the protection of exposed workers and the general population with domestic exposures.

RESULTS

The Swiss Centre for Applied Human Toxicology funded the project, commencing in January 2021. The Human Ethics Committee approval was obtained on November 16, 2022. Zebrafish experiments and in vitro methods started in February 2021, whereas recruitment of human volunteers started in 2022 after the COVID-19 pandemic-related restrictions were lifted. We anticipate that we will be able to provide a neurotoxicity testing strategy by 2026 and predicted air concentrations for 6 commonly used propylene glycol ethers based on toxicokinetic models incorporating liver metabolism, BBB leakage parameters, and brain toxicity.

CONCLUSIONS

This study will be of great interest to regulatory agencies and chemical industries needing and seeking novel solutions to develop human chemical risk assessments. It will contribute to protecting human health from the deleterious effects of environmental chemicals.

INTERNATIONAL REGISTERED REPORT IDENTIFIER (IRRID)

DERR1-10.2196/50300.

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.

Development of a Generic PBK Model for Human Biomonitoring with an Application to Deoxynivalenol

Toxicokinetic modelling provides a powerful tool in relating internal human exposure (i.e., assessed through urinary biomarker levels) to external exposure. Chemical specific toxicokinetic models are available; however, this specificity prevents their application to similar contaminants or to other routes of exposure. For this reason, we investigated whether a generic physiological-based kinetic (PBK) model might be a suitable alternative for a biokinetic model of deoxynivalenol (DON). IndusChemFate (ICF) was selected as a generic PBK model, which could be fit for purpose. Being suited for simulating multiple routes of exposure, ICF has particularly been used to relate the inhalation and dermal exposure of industrial chemicals to their urinary excretion. For the first time, the ICF model was adapted as a generic model for the human biomonitoring of mycotoxins, thereby extending its applicability domain. For this purpose, chemical-specific data for DON and its metabolites were collected directly from the literature (distribution and metabolism) or indirectly (absorption and excretion) by fitting the ICF model to previously described urinary excretion data. The obtained results indicate that this generic model can be used to model the urinary excretion of DON and its glucuronidated metabolites following dietary exposure to DON. Additionally, the present study establishes the basis for further development of the model to include an inhalation exposure route alongside the oral exposure route.

A Model Template Approach for Rapid Evaluation and Application of Physiologically Based Pharmacokinetic Models: Extension to Volatile Organic Compounds

Chemical risk assessors use physiologically based pharmacokinetic (PBPK) models to perform dosimetric calculations, including extrapolations between exposure scenarios, species, and populations of interest. Assessors should complete a thorough quality assurance (QA) review to ensure biological accuracy and correct implementation prior to using these models. This process can be time-consuming, and we developed a PBPK model template that allows for faster, more efficient QA review. The model template consists of a single model "superstructure" with equations and logic commonly found in PBPK models, allowing users to implement a wide variety of chemical-specific PBPK models. QA review can be completed more quickly than for conventional PBPK model implementations because the general model equations have already been reviewed and only parameters describing chemical-specific model and exposure scenarios need review for any given model implementation. We have expanded a previous version of the PBPK model template by adding features commonly included in PBPK models for volatile organic compounds (VOCs). We included multiple options for representing concentrations in blood, describing metabolism, and modeling gas exchange processes to allow for inhalation exposures. We created PBPK model template implementations of published models for seven VOCs: dichloromethane, methanol, chloroform, styrene, vinyl chloride, trichloroethylene, and carbon tetrachloride. Simulations performed using our template implementations matched published simulation results to a high degree of accuracy (maximum observed percent error: 1%). Thus, the model template approach can now be applied to a broader class of chemical-specific PBPK models while continuing to bolster efficiency of QA processes that should be conducted prior to using models for risk assessment applications.

Human skin responses to environmental pollutants: A review of current scientific models

Whatever the exposure route, chemical, physical and biological pollutants modify the whole organism response, leading to nerve, cardiac, respiratory, reproductive, and skin system pathologies. Skin acts as a barrier for preventing pollutant modifications. This review aims to present the available scientific models, which help investigate the impact of pollution on the skin. The research question was "Which experimental models illustrate the impact of pollution on the skin in humans?" The review covered a period of 10 years following a PECO statement on in vitro, ex vivo, in vivo and in silico models. Of 582 retrieved articles, 118 articles were eligible. In oral and inhalation routes, dermal exposure had an important impact at both local and systemic levels. Healthy skin models included primary cells, cell lines, co-cultures, reconstructed human epidermis, and skin explants. In silico models estimated skin exposure and permeability. All pollutants affected the skin by altering elasticity, thickness, the structure of epidermal barrier strength, and dermal extracellular integrity. Some specific models concerned wound healing or the skin aging process. Underlying mechanisms were an exacerbated inflammatory skin reaction with the modulation of several cytokines and oxidative stress responses, ending with apoptosis. Pathological skin models revealed the consequences of environmental pollutants on psoriasis, atopic dermatitis, and tumour development. Finally, scientific models were used for evaluating the safety and efficacy of potential skin formulations in preventing the skin aging process or skin irritation after repeated contact. The review gives an overview of scientific skin models used to assess the effects of pollutants. Chemical and physical pollutants were mainly represented while biological contaminants were little studied. In future developments, cell hypoxia and microbiota models may be considered as more representative of clinical situations. Models considering humidity and temperature variations may reflect the impact of these changes.

A Model Template Approach for Rapid Evaluation and Application of Physiologically Based Pharmacokinetic Models for Use in Human Health Risk Assessments: A Case Study on Per- and Polyfluoroalkyl Substances

Physiologically based pharmacokinetic (PBPK) models are commonly used in risk assessments to perform inter- and intraspecies extrapolations as well as to extrapolate between different dosing scenarios; however, they must first undergo quality assurance review, which can be a time-consuming process, especially when model code is not readily available. We developed and implemented (using R and MCSim) a PBPK model template capable of replicating published model results for several chemical-specific PBPK models. This model template allows for faster quality assurance review because the general model equations only need to be reviewed once, and application to a specific chemical then only requires reviewing input parameters. The model template can implement PBPK models with oral and intravenous exposure routes, varying numbers of tissue compartments, renal reabsorption, and multiple elimination pathways, including fecal, urinary, and biliary. Using the model template, we reproduced published model simulation results for perfluorohexanesulfonic acid, perfluorononanoic acid, perfluorodecanoic acid, perfluorooctanoate, and perflouorooctane sulfonate. We also show that the template can be a useful tool for identifying potential model errors. Thus, the model template allows for faster evaluation and review of published PBPK models and provides a proof of concept for using this approach with broader classes of chemical-specific PBPK models.

A Scoping Review of Technologies and Their Applicability for Exposome-Based Risk Assessment in the Oil and Gas Industry

INTRODUCTION

Oil and gas workers have been shown to be at increased risk of chronic diseases including cancer, asthma, chronic obstructive pulmonary disease, and hearing loss, among others. Technological advances may be used to assess the external (e.g. personal sensors, smartphone apps and online platforms, exposure models) and internal exposome (e.g. physiologically based kinetic modeling (PBK), biomonitoring, omics), offering numerous possibilities for chronic disease prevention strategies and risk management measures. The objective of this study was to review the literature on these technologies, by focusing on: (i) evaluating their applicability for exposome research in the oil and gas industry, and (ii) identifying key challenges that may hamper the successful application of such technologies in the oil and gas industry.

METHOD

A scoping review was conducted by identifying peer-reviewed literature with searches in MEDLINE/PubMed and SciVerse Scopus. Two assessors trained on the search strategy screened retrieved articles on title and abstract. The inclusion criteria used for this review were: application of the aforementioned technologies at a workplace in the oil and gas industry or, application of these technologies for an exposure relevant to the oil and gas industry but in another occupational sector, English language and publication period 2005-end of 2019.

RESULTS

In total, 72 articles were included in this scoping review with most articles focused on omics and bioinformatics (N = 22), followed by biomonitoring and biomarkers (N = 20), external exposure modeling (N = 11), PBK modeling (N = 10), and personal sensors (N = 9). Several studies were identified in the oil and gas industry on the application of PBK models and biomarkers, mainly focusing on workers exposed to benzene. The application of personal sensors, new types of exposure models, and omics technology are still in their infancy with respect to the oil and gas industry. Nevertheless, applications of these technologies in other occupational sectors showed the potential for application in this sector.

DISCUSSION AND CONCLUSION

New exposome technologies offer great promise for personal monitoring of workers in the oil and gas industry, but more applied research is needed in collaboration with the industry. Current challenges hindering a successful application of such technologies include (i) the technological readiness of sensors, (ii) the availability of data, (iii) the absence of standardized and validated methods, and (iv) the need for new study designs to study the development of disease during working life.

In vitro evaluation of the metabolic stability of nine fragrance chemicals in trout and human hepatocytes

In vitro metabolic stability of nine fragrance chemicals: p-tolyl acetate, cashmeran, ethylene brassylate, celestolide, galaxolide, traseolide, ambretone, tonalide and pentadecanolide, was evaluated in trout and human hepatocytes. The compounds were incubated with trout hepatocytes at 12°C and human hepatocytes at 37°C. Quantification of compound disappearance with time was performed using gas chromatography/mass spectrometry. in vivo hepatic intrinsic clearance values were calculated from the in vitro data. Significant metabolism was observed with trout hepatocytes for five of the nine fragrance chemicals, while all nine were metabolized significantly with human hepatocytes. Previously published models were used to examine expected bioaccumulation and persistence in whole organisms. Calculated half-lives due to metabolism of the nine chemicals are significantly shorter for humans than trout: <1 hour and <1 day, respectively. For all chemicals with demonstrated hepatic metabolism, the models indicate a lack of accumulation. For those where metabolism was demonstrated in trout, calculated bioconcentration factors would not be classified as bioaccumulative under prevailing regulatory systems.

Newborn Exposure to Polycyclic Aromatic Hydrocarbons Through Parenteral Nutrition

BACKGROUND AND OBJECTIVES

Parenteral nutrition (PN) administered to newborns (NB) may be contaminated with polycyclic aromatic hydrocarbons (PAHs) and may therefore increase the contact with these toxicants in very early life stages. The aim of the study is to determine to what extent, if any, commercial products for PN are contaminated with PAHs and to determine whether these contaminants, when present in the bag content, are delivered to NB and whether 1-hydroxypyrene (1-HP), the pyrene metabolite, can be detected in the urine of exposed NB.

METHODS

Commercial products and the bags administered to 10 NB during their period in the NICU were analyzed for the 16 priority US Environmental Protection Agency PAHs. Urine samples were collected and analyzed for their 1-HP content. Urine samples of a control group composed of 8 breastfed NB were also analyzed for the determination of 1-HP.

RESULTS

From 9 different commercial products used to compound PN bags, 6 were contaminated with PAHs, with total concentrations varying from 0.02 to 10.56 mg/L. In the bags administered to the NB, this sum varied from 0.01 to 6.30 mg/L with a mean of 2.62 mg/L. Therefore, for each 100 mL PN, an average load of 0.26 mg PAHs was observed. The majority of the urine samples taken from NB in the study group (80%) contained 1-HP, but it was not detected in the urine of any baby in the control group.

CONCLUSIONS

The contamination of PN with PAHs poses a critical toxicological risk. The elevated contaminant concentrations and the parenteral way of administration make this source of PAHs considerably worse than any other, including maternal exposure to environmental pollution or tobacco.

A review of the current state of the art of physiologically-based tests for measuring human dermal in vitro bioavailability of polycyclic aromatic hydrocarbons (PAH) in soil

Polycyclic Aromatic Hydrocarbons are classed as Persistent Organic Pollutants, a large group of compounds that share similar characteristics. They are lipophilic, resistant to degradation in the environment and harmful to human and environmental health. Soil has been identified as the primary reservoir for Polycyclic Aromatic Hydrocarbons in the United Kingdom. This study reviews the literature associated with, or is relevant to, the measurement and modelling of dermal absorption of Polycyclic Aromatic Hydrocarbons from soils. The literature illustrates the use of in vivo, in vitro and in silico methods from a wide variety of scientific disciplines including occupational and environmental exposure, medical, pharmaceutical and cosmetic research and associated mathematical modelling. The review identifies a number of practical shortcomings which must be addressed if dermal bioavailability tests are to be applied to laboratory analysis of contaminated soils for human health risk assessment.

Cancer risk of petrochemical workers exposed to airborne PAHs in industrial Lanzhou City, China

This paper reports the connections between red blood cells abnormality risk of petrochemical workers and their exposure to airborne polycyclic aromatic hydrocarbons (PAHs). Urinary 1-hydroxypyrene (1-OHP), as the biomarker of PAHs exposure, was adopted to assess the exposure risk of the petrochemical workers to PAHs in Xigu, the west suburb of Lanzhou where petrochemical industries are located. Fifty-three workers, sub-grouped to 36 petrochemical workers and 17 office workers, participated in this investigation. Logistic regression model and spearman correlation analysis were performed to estimate the associations between PAHs exposure levels and red blood cells abnormality risk of petrochemical workers. Strong associations between some red cell indices (MCH, MCHC, RDW) and 1-OHP concentration were found. Results also show that the red blood cells abnormality risk increased with increasing PAHs exposure level. Compared with office workers, risk level of red blood cells abnormality in petrochemical workers was higher by 41.7 % (OR, 1.417; 95 % CI: 0.368-5.456) than that in office workers. This result was verified by the tissue-to-human blood partition coefficient for pyrene and 1-OHP. The quantitative assessments of the potential health risk through inhalation exposure to PAHs were conducted using the Incremental Lifetime Cancer Risk (ILCR) model. It was found the ILCR from inhalation exposure to PAHs for the petrochemical workers ranged from 10(-5) to 10(-4) with 95 % probability, indicating that petrochemical plant workers were under a high potential cancer risk level.

Use of physiologically-based pharmacokinetic modeling to simulate the profiles of 3-hydroxybenzo(a)pyrene in workers exposed to polycyclic aromatic hydrocarbons

Biomathematical modeling has become an important tool to assess xenobiotic exposure in humans. In the present study, we have used a human physiologically-based pharmacokinetic (PBPK) model and an simple compartmental toxicokinetic model of benzo(a)pyrene (BaP) kinetics and its 3-hydroxybenzo(a)pyrene (3-OHBaP) metabolite to reproduce the time-course of this biomarker of exposure in the urine of industrially exposed workers and in turn predict the most plausible exposure scenarios. The models were constructed from in vivo experimental data in rats and then extrapolated from animals to humans after assessing and adjusting the most sensitive model parameters as well as species specific physiological parameters. Repeated urinary voids from workers exposed to polycyclic aromatic hydrocarbons (PAHs) have been collected over the course of a typical workweek and during subsequent days off work; urinary concentrations of 3-OHBaP were then determined. Based on the information obtained for each worker (BaP air concentration, daily shift hours, tasks, protective equipment), the time courses of 3-OHBaP in the urine of the different workers have been simulated using the PBPK and toxicokinetic models, considering the various possible exposure routes, oral, dermal and inhalation. Both models were equally able to closely reproduce the observed time course of 3-OHBaP in the urine of workers and predicted similar exposure scenarios. Simulations of various scenarios suggest that the workers under study were exposed mainly by the dermal route. Comparison of measured air concentration levels of BaP with simulated values needed to obtain a good approximation of observed time course further pointed out that inhalation was not the main route of exposure for most of the studied workers. Both kinetic models appear as a useful tool to interpret biomonitoring data of PAH exposure on the basis of 3-OHBaP levels.

Understanding the linked kinetics of benzo(a)pyrene and 3-hydroxybenzo(a)pyrene biomarker of exposure using physiologically-based pharmacokinetic modelling in rats

3-hydroxybenzo(a)pyrene (3-OHBaP) in urine has been proposed as a biomarker of occupational exposure to polycyclic aromatic hydrocarbons. However, to reconstruct exposure doses in workers from biomarker measurements, a thorough knowledge of the kinetics of the benzo(a)pyrene (BaP) and 3-OHBaP given different routes of exposure is needed. A rat physiologically-based pharmacokinetic model of BaP and 3-OHBaP was built. Organs (tissues) represented as compartments were based on in vivo experimental data in rats. Tissue: blood partition coefficients, permeability coefficients, metabolism rates, excretion parameters, and absorption fractions and rates for different routes-of-entry were obtained directly from published in vivo time courses of BaP and 3-OHBaP in blood, various tissues and excreta of rats. The latter parameter values were best-fitted by least square procedures and Monte Carlo simulations. Sensitivity analyses were then carried out to ensure the stability of the model and the key parameters driving the overall modeled kinetics. This modeling pointed out critical determinants of the kinetics: (1) hepatic metabolism of BaP and 3-OHBaP elimination rate as the most sensitive parameters; (2) the strong partition of BaP in lungs compared to other tissues, followed by adipose tissues and liver; (3) the strong partition of 3-OHBaP in kidneys; (4) diffusion-limited tissue transfers of BaP in lungs and 3-OHBaP in lungs, adipose tissues and kidneys; (5) significant entero-hepatic recycling of 3-OHBaP. Very good fits to various sets of experimental data in rats from four different routes-of-entry (intravenous, oral, dermal and inhalation) were obtained with the model.