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

Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans

Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration-time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration-time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.

A database of chemical absorption in human skin with mechanistic modeling applications

Whether from environmental and occupational hazards or from topical pharmaceuticals, the human skin comes into contact with various chemicals every day. In vivo experiments not only require large investments of both time and money, but in vivo experiments can also be unethical due to the need to intentionally or incidentally expose humans or animals to toxic chemicals. Comparatively, in vitro experiments offer ethical and financial advantages when combined with the opportunity to selectively choose chemicals for experimentation. With in vivo experimentation being so infeasible, many scientists have chosen to make their in vitro data available publicly. Using these data, a detailed database containing 73 chemicals was created with a robust set of descriptors to be used in connection with mathematical modeling to predict diffusion, permeability, and partition coefficients. This resulting database is tailored to be easily used in various coding languages.

Investigating open access new approach methods (NAM) to assess biological points of departure: A case study with 4 neurotoxic pesticides

Open access new approach methods (NAM) in the US EPA ToxCast program and NTP Integrated Chemical Environment (ICE) were used to investigate activities of four neurotoxic pesticides: endosulfan, fipronil, propyzamide and carbaryl. Concordance of in vivo regulatory points of departure (POD) adjusted for interspecies extrapolation (AdjPOD) to modelled human Administered Equivalent Dose (AEDHuman) was assessed using 3-compartment or Adult/Fetal PBTK in vitro to in vivo extrapolation. Model inputs were from Tier 1 (High throughput transcriptomics: HTTr, high throughput phenotypic profiling: HTPP) and Tier 2 (single target: ToxCast) assays. HTTr identified gene expression signatures associated with potential neurotoxicity for endosulfan, propyzamide and carbaryl in non-neuronal MCF-7 and HepaRG cells. The HTPP assay in U-2 OS cells detected potent effects on DNA endpoints for endosulfan and carbaryl, and mitochondria with fipronil (propyzamide was inactive). The most potent ToxCast assays were concordant with specific components of each chemical mode of action (MOA). Predictive adult IVIVE models produced fold differences (FD) < 10 between the AEDHuman and the measured in vivo AdjPOD. The 3-compartment model was concordant (i.e., smallest FD) for endosulfan, fipronil and carbaryl, and PBTK was concordant for propyzamide. The most potent AEDHuman predictions for each chemical showed HTTr, HTPP and ToxCast were mainly concordant with in vivo AdjPODs but assays were less concordant with MOAs. This was likely due to the cell types used for testing and/or lack of metabolic capabilities and pathways available in vivo. The Fetal PBTK model had larger FDs than adult models and was less predictive overall.

Human biomonitoring and toxicokinetics as key building blocks for next generation risk assessment

Human health risk assessment is historically built upon animal testing, often following Organisation for Economic Co-operation and Development (OECD) test guidelines and exposure assessments. Using combinations of human relevant in vitro models, chemical analysis and computational (in silico) approaches bring advantages compared to animal studies. These include a greater focus on the human species and on molecular mechanisms and kinetics, identification of Adverse Outcome Pathways and downstream Key Events as well as the possibility of addressing susceptible populations and additional endpoints. Much of the advancement and progress made in the Next Generation Risk Assessment (NGRA) have been primarily focused on new approach methodologies (NAMs) and physiologically based kinetic (PBK) modelling without incorporating human biomonitoring (HBM). The integration of toxicokinetics (TK) and PBK modelling is an essential component of NGRA. PBK models are essential for describing in quantitative terms the TK processes with a focus on the effective dose at the expected target site. Furthermore, the need for PBK models is amplified by the increasing scientific and regulatory interest in aggregate and cumulative exposure as well as interactions of chemicals in mixtures. Since incorporating HBM data strengthens approaches and reduces uncertainties in risk assessment, here we elaborate on the integrated use of TK, PBK modelling and HBM in chemical risk assessment highlighting opportunities as well as challenges and limitations. Examples are provided where HBM and TK/PBK modelling can be used in both exposure assessment and hazard characterization shifting from external exposure and animal dose/response assays to animal-free, internal exposure-based NGRA.

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.

A deep-learning approach for identifying prospective chemical hazards

With the aim of helping to set safe exposure limits for the general population, various techniques have been implemented to conduct risk assessments for chemicals and other environmental stressors; however, none of these tools facilitate the identification of completely new chemicals that are likely hazardous and elicit an adverse biological effect. Here, we detail a novel in silico, deep-learning framework that is designed to systematically generate structures for new chemical compounds that are predicted to be chemical hazards. To assess the utility of the framework, we applied the tool to four endpoints related to environmental toxicants and their impacts on human and animal health: (i) toxicity to honeybees, (ii) immunotoxicity, (iii) endocrine disruption via ER-α antagonism, and (iv) mutagenicity. In addition, we characterized the predicted potency of these compounds and examined their structural relationship to existing chemicals of concern. As part of the array of emerging new approach methodologies (NAMs), we anticipate that such a framework will be a significant asset to risk assessors and other environmental scientists when planning and forecasting. Though not in the scope of the present study, we expect that the methodology detailed here could also be useful in the de novo design of more environmentally-friendly industrial chemicals.

Characterizing Chemical Exposure Trends from NHANES Urinary Biomonitoring Data

BACKGROUND

Xenobiotic metabolites are widely present in human urine and can indicate recent exposure to environmental chemicals. Proper inference of which chemicals contribute to these metabolites can inform human exposure and risk. Furthermore, longitudinal biomonitoring studies provide insight into how chemical exposures change over time.

OBJECTIVES

We constructed an exposure landscape for as many human-exposure relevant chemicals over as large a time span as possible to characterize exposure trends across demographic groups and chemical types.

METHODS

We analyzed urine data of nine 2-y cohorts (1999-2016) from the National Health and Nutrition Examination Survey (NHANES). Chemical daily intake rates (in milligrams per kilogram bodyweight per day) were inferred, using the R package bayesmarker, from metabolite concentrations in each cohort individually to identify exposure trends. Trends for metabolites and parents were clustered to find chemicals with similar exposure patterns. Exposure variation by age, gender, and body mass index were also assessed.

RESULTS

Intake rates for 179 parent chemicals were inferred from 151 metabolites (96 measured in five or more cohorts). Seventeen metabolites and 44 parent chemicals exhibited fold-changes ≥ 10 between any two cohorts (deltamethrin, di-n -octyl phthalate, and di-isononyl phthalate had the greatest exposure increases). Di-2-ethylhexyl phthalate intake began decreasing in 2007, whereas both di-isobutyl and di-isononyl phthalate began increasing shortly before. Intake for four parabens was markedly higher in females, especially reproductive-age females, compared with males and children. Cadmium and arsenobetaine exhibited higher exposure for individuals > 65 years of age and lower for individuals < 20 years of age.

DISCUSSION

With appropriate analysis, NHANES indicates trends in chemical exposures over the past two decades. Decreases in exposure are observable as the result of regulatory action, with some being accompanied by increases in replacement chemicals. Age- and gender-specific variations in exposure were observed for multiple chemicals. Continued estimation of demographic-specific exposures is needed to both monitor and identify potential vulnerable populations. https://doi.org/10.1289/EHP12188.

Reproducibility of organ-level effects in repeat dose animal studies

This work estimates benchmarks for new approach method (NAM) performance in predicting organ-level effects in repeat dose studies of adult animals based on variability in replicate animal studies. Treatment-related effect values from the Toxicity Reference database (v2.1) for weight, gross, or histopathological changes in the adrenal gland, liver, kidney, spleen, stomach, and thyroid were used. Rates of chemical concordance among organ-level findings in replicate studies, defined by repeated chemical only, chemical and species, or chemical and study type, were calculated. Concordance was 39 - 88%, depending on organ, and was highest within species. Variance in treatment-related effect values, including lowest effect level (LEL) values and benchmark dose (BMD) values when available, was calculated by organ. Multilinear regression modeling, using study descriptors of organ-level effect values as covariates, was used to estimate total variance, mean square error (MSE), and root residual mean square error (RMSE). MSE values, interpreted as estimates of unexplained variance, suggest study descriptors accounted for 52-69% of total variance in organ-level LELs. RMSE ranged from 0.41 - 0.68 log10-mg/kg/day. Differences between organ-level effects from chronic (CHR) and subchronic (SUB) dosing regimens were also quantified. Odds ratios indicated CHR organ effects were unlikely if the SUB study was negative. Mean differences of CHR - SUB organ-level LELs ranged from -0.38 to -0.19 log10 mg/kg/day; the magnitudes of these mean differences were less than RMSE for replicate studies. Finally, in vitro to in vivo extrapolation (IVIVE) was employed to compare bioactive concentrations from in vitro NAMs for kidney and liver to LELs. The observed mean difference between LELs and mean IVIVE dose predictions approached 0.5 log10-mg/kg/day, but differences by chemical ranged widely. Overall, variability in repeat dose organ-level effects suggests expectations for quantitative accuracy of NAM prediction of LELs should be at least ± 1 log10-mg/kg/day, with qualitative accuracy not exceeding 70%.

Cutting-edge computational chemical exposure research at the U.S. Environmental Protection Agency

Exposure science is evolving from its traditional "after the fact" and "one chemical at a time" approach to forecasting chemical exposures rapidly enough to keep pace with the constantly expanding landscape of chemicals and exposures. In this article, we provide an overview of the approaches, accomplishments, and plans for advancing computational exposure science within the U.S. Environmental Protection Agency's Office of Research and Development (EPA/ORD). First, to characterize the universe of chemicals in commerce and the environment, a carefully curated, web-accessible chemical resource has been created. This DSSTox database unambiguously identifies >1.2 million unique substances reflecting potential environmental and human exposures and includes computationally accessible links to each compound's corresponding data resources. Next, EPA is developing, applying, and evaluating predictive exposure models. These models increasingly rely on data, computational tools like quantitative structure activity relationship (QSAR) models, and machine learning/artificial intelligence to provide timely and efficient prediction of chemical exposure (and associated uncertainty) for thousands of chemicals at a time. Integral to this modeling effort, EPA is developing data resources across the exposure continuum that includes application of high-resolution mass spectrometry (HRMS) non-targeted analysis (NTA) methods providing measurement capability at scale with the number of chemicals in commerce. These research efforts are integrated and well-tailored to support population exposure assessment to prioritize chemicals for exposure as a critical input to risk management. In addition, the exposure forecasts will allow a wide variety of stakeholders to explore sustainable initiatives like green chemistry to achieve economic, social, and environmental prosperity and protection of future generations.

A regression analysis using simple descriptors for multiple dermal datasets: Going from individual membranes to the full skin

In silico methods to estimate and/or quantify skin absorption of chemicals as a function of chemistry are needed to realistically predict pharmacological, occupational, and environmental exposures. The Potts-Guy equation is a well-established approach, using multi-linear regression analysis describing skin permeability (Kp) in terms of the octanol/water partition coefficient (logP) and molecular weight (MW). In this work, we obtained regression equations for different human datasets relevant to environmental and cosmetic chemicals. Since the Potts-Guy equation was published in 1992, we explored recent datasets that include different skin layers, such as dermatomed (including dermis to a defined thickness) and full skin. Our work was consistent with others who have observed that fits to the Potts-Guy equation are stronger for experiments focused on the epidermis. Permeability estimates for dermatomed skin and full skin resulted in low regression coefficients when compared to epidermis datasets. An updated regression equation uses a combination of fitted permeability values obtained with a published 2D compartmental model previously evaluated. The resulting regression equation was: logKp = -2.55 + 0.65logP - 0.0085MW, R2  = 0.91 (applicability domain for all datasets: MW ranges from 18 to >584 g/mol and -4 to >5 for logP). This approach demonstrates the advantage of combining mechanistic with structural activity relationships in a single modeling approach. This combination approach results in an improved regression fit when compared to permeability estimates obtained using the Potts-Guy approach alone. The analysis presented in this work assumes a one-compartment skin absorption route; future modeling work will consider adding multiple compartments.

Application of cell painting for chemical hazard evaluation in support of screening-level chemical assessments

'Cell Painting' is an imaging-based high-throughput phenotypic profiling (HTPP) method in which cultured cells are fluorescently labeled to visualize subcellular structures (i.e., nucleus, nucleoli, endoplasmic reticulum, cytoskeleton, Golgi apparatus / plasma membrane and mitochondria) and to quantify morphological changes in response to chemicals or other perturbagens. HTPP is a high-throughput and cost-effective bioactivity screening method that detects effects associated with many different molecular mechanisms in an untargeted manner, enabling rapid in vitro hazard assessment for thousands of chemicals. Here, 1201 chemicals from the ToxCast library were screened in concentration-response up to ~100 μM in human U-2 OS cells using HTPP. A phenotype altering concentration (PAC) was estimated for chemicals active in the tested range. PACs tended to be higher than lower bound potency values estimated from a broad collection of targeted high-throughput assays, but lower than the threshold for cytotoxicity. In vitro to in vivo extrapolation (IVIVE) was used to estimate administered equivalent doses (AEDs) based on PACs for comparison to human exposure predictions. AEDs for 18/412 chemicals overlapped with predicted human exposures. Phenotypic profile information was also leveraged to identify putative mechanisms of action and group chemicals. Of 58 known nuclear receptor modulators, only glucocorticoids and retinoids produced characteristic profiles; and both receptor types are expressed in U-2 OS cells. Thirteen chemicals with profile similarity to glucocorticoids were tested in a secondary screen and one chemical, pyrene, was confirmed by an orthogonal gene expression assay as a novel putative GR modulating chemical. Most active chemicals demonstrated profiles not associated with a known mechanism-of-action. However, many structurally related chemicals produced similar profiles, with exceptions such as diniconazole, whose profile differed from other active conazoles. Overall, the present study demonstrates how HTPP can be applied in screening-level chemical assessments through a series of examples and brief case studies.

Use of computational toxicology models to predict toxicological points of departure: A case study with triazine herbicides

BACKGROUND

Atrazine simazine and propazine, widely used triazine herbicides on food crops and in residential areas, disrupt the neuroendocrine system raising human health concerns. USEPA developed a PBPK model based on triazine common Mode of Action (MOA)-suppression of luteinizing hormone surge in female rats-to generate human regulatory points of departure (POD: mg/kg/day). We compared triazine Human Administered Equivalent Dose (AED_Human mg/kg/day) predictions from open access computational tools to the PBPK PODs to assess concordance.

METHODS

Computational tools were the following: ToxCast/Tox21 in vitro assays; Toxicogenomic databases to assess concordance with ToxCast/Tox21 targets; integrated chemical environment (ICE) models with ToxCast/Tox21 inputs to predict AED_Human PODs and population-based age-refined high throughput toxicokinetics (HTTK-Pop) to compare to age-related PBPK PODs.

RESULTS

ToxCast/Tox21 assays identified critical targets in the triazine common MOA and gene databases; ICE AED_Human predictions were mainly concordant with the USEPA PBPK PODs quantitatively. Low fold-differences between PBPK POD and ICE AED_Human predictions indicated that the ICE models are health-protective. HTTK-Pop age-refinements were within 10-fold of the USEPA PBPK PODs.

CONCLUSIONS

CompTox tools were used to identify assay targets in the MOA and identify potential molecular initiating targets in the adverse outcome pathway for potential use in risk assessment.

A Pragmatic Framework for the Application of New Approach Methodologies in One Health Toxicological Risk Assessment

Globally, industries and regulatory authorities are faced with an urgent need to assess the potential adverse effects of chemicals more efficiently by embracing new approach methodologies (NAMs). NAMs include cell and tissue methods (in vitro), structure-based/toxicokinetic models (in silico), methods that assess toxicant interactions with biological macromolecules (in chemico), and alternative models. Increasing knowledge on chemical toxicokinetics (what the body does with chemicals) and toxicodynamics (what the chemicals do with the body) obtained from in silico and in vitro systems continues to provide opportunities for modernizing chemical risk assessments. However, directly leveraging in vitro and in silico data for derivation of human health-based reference values has not received regulatory acceptance due to uncertainties in extrapolating NAM results to human populations, including metabolism, complex biological pathways, multiple exposures, interindividual susceptibility and vulnerable populations. The objective of this article is to provide a standardized pragmatic framework that applies integrated approaches with a focus on quantitative in vitro to in vivo extrapolation (QIVIVE) to extrapolate in vitro cellular exposures to human equivalent doses from which human reference values can be derived. The proposed framework intends to systematically account for the complexities in extrapolation and data interpretation to support sound human health safety decisions in diverse industrial sectors (food systems, cosmetics, industrial chemicals, pharmaceuticals etc.). Case studies of chemical entities, using new and existing data, are presented to demonstrate the utility of the proposed framework while highlighting potential sources of human population bias and uncertainty, and the importance of Good Method and Reporting Practices.

Identifying xenobiotic metabolites with in silico prediction tools and LCMS suspect screening analysis

Understanding the metabolic fate of a xenobiotic substance can help inform its potential health risks and allow for the identification of signature metabolites associated with exposure. The need to characterize metabolites of poorly studied or novel substances has shifted exposure studies towards non-targeted analysis (NTA), which often aims to profile many compounds within a sample using high-resolution liquid-chromatography mass-spectrometry (LCMS). Here we evaluate the suitability of suspect screening analysis (SSA) liquid-chromatography mass-spectrometry to inform xenobiotic chemical metabolism. Given a lack of knowledge of true metabolites for most chemicals, predictive tools were used to generate potential metabolites as suspect screening lists to guide the identification of selected xenobiotic substances and their associated metabolites. Thirty-three substances were selected to represent a diverse array of pharmaceutical, agrochemical, and industrial chemicals from Environmental Protection Agency's ToxCast chemical library. The compounds were incubated in a metabolically-active in vitro assay using primary hepatocytes and the resulting supernatant and lysate fractions were analyzed with high-resolution LCMS. Metabolites were simulated for each compound structure using software and then combined to serve as the suspect screening list. The exact masses of the predicted metabolites were then used to select LCMS features for fragmentation via tandem mass spectrometry (MS/MS). Of the starting chemicals, 12 were measured in at least one sample in either positive or negative ion mode and a subset of these were used to develop the analysis workflow. We implemented a screening level workflow for background subtraction and the incorporation of time-varying kinetics into the identification of likely metabolites. We used haloperidol as a case study to perform an in-depth analysis, which resulted in identifying five known metabolites and five molecular features that represent potential novel metabolites, two of which were assigned discrete structures based on in silico predictions. This workflow was applied to five additional test chemicals, and 15 molecular features were selected as either reported metabolites, predicted metabolites, or potential metabolites without a structural assignment. This study demonstrates that in some-but not all-cases, suspect screening analysis methods provide a means to rapidly identify and characterize metabolites of xenobiotic chemicals.

A geospatial modeling approach to quantifying the risk of exposure to environmental chemical mixtures via a common molecular target

In the real world, individuals are exposed to chemicals from sources that vary over space and time. However, traditional risk assessments based on in vivo animal studies typically use a chemical-by-chemical approach and apical disease endpoints. New approach methodologies (NAMs) in toxicology, such as in vitro high-throughput (HTS) assays generated in Tox21 and ToxCast, can more readily provide mechanistic chemical hazard information for chemicals with no existing data than in vivo methods. In this paper, we establish a workflow to assess the joint action of 41 modeled ambient chemical exposures in the air from the USA-wide National Air Toxics Assessment by integrating human exposures with hazard data from curated HTS (cHTS) assays to identify counties where exposure to the local chemical mixture may perturb a common biological target. We exemplify this proof-of-concept using CYP1A1 mRNA up-regulation. We first estimate internal exposure and then convert the inhaled concentration to a steady state plasma concentration using physiologically based toxicokinetic modeling parameterized with county-specific information on ages and body weights. We then use the estimated blood plasma concentration and the concentration-response curve from the in vitro cHTS assay to determine the chemical-specific effects of the mixture components. Three mixture modeling methods were used to estimate the joint effect from exposure to the chemical mixture on the activity levels, which were geospatially mapped. Finally, a Monte Carlo uncertainty analysis was performed to quantify the influence of each parameter on the combined effects. This workflow demonstrates how NAMs can be used to predict early-stage biological perturbations that can lead to adverse health outcomes that result from exposure to chemical mixtures. As a result, this work will advance mixture risk assessment and other early events in the effects of chemicals.

The impact of environmental contaminants on extracellular vesicles and their key molecular regulators: A literature and database-driven review

Exposure to environmental chemicals is now well recognized as a significant factor contributing to the global burden of disease; however, there remain critical gaps in understanding the types of biological mechanisms that link environmental chemicals to adverse health outcomes. One type of mechanism that remains understudied involves extracellular vesicles (EVs), representing small cell-derived particles capable of carrying molecular signals such as RNAs, miRNAs, proteins, lipids, and chemicals through biological fluids and imparting beneficial, neutral, or negative effects on target cells. In fact, evidence is just now starting to grow that supports the role of EVs in various disease etiologies. This review aims to (1) Provide a landscape of the current understanding of the functional relationship between EVs and environmental chemicals; (2) Summarize current knowledge of EV regulatory processes including production, packaging, and release; and (3) Conduct a database-driven analysis of known chemical-gene interactions to predict and prioritize environmentally relevant chemicals that may impact EV regulatory genes and thus EV regulatory processes. This approach to predicting environmentally relevant chemicals that may alter EVs provides a novel method for evidence-based hypothesis generation for future studies evaluating the link between environmental exposures and EVs.

Three perspectives on the prediction of chemical effects in ecosystems

The increasing production, use and emission of synthetic chemicals into the environment represents a major driver of global change. The large number of synthetic chemicals, limited knowledge on exposure patterns and effects in organisms and their interaction with other global change drivers hamper the prediction of effects in ecosystems. However, recent advances in biomolecular and computational methods are promising to improve our capacity for prediction. We delineate three idealised perspectives for the prediction of chemical effects: the suborganismal, organismal and ecological perspective, which are currently largely separated. Each of the outlined perspectives includes essential and complementary theories and tools for prediction but captures only part of the phenomenon of chemical effects. Links between the perspectives may foster predictive modelling of chemical effects in ecosystems and extrapolation between species. A major challenge for the linkage is the lack of data sets simultaneously covering different levels of biological organisation (here referred to as biological levels) as well as varying temporal and spatial scales. Synthesising the three perspectives, some central aspects and associated types of data seem particularly necessary to improve prediction. First, suborganism- and organism-level responses to chemicals need to be recorded and tested for relationships with chemical groups and organism traits. Second, metrics that are measurable at many biological levels, such as energy, need to be scrutinised for their potential to integrate across levels. Third, experimental data on the simultaneous response over multiple biological levels and spatiotemporal scales are required. These could be collected in nested and interconnected micro- and mesocosm experiments. Lastly, prioritisation of processes involved in the prediction framework needs to find a balance between simplification and capturing the essential complexity of a system. For example, in some cases, eco-evolutionary dynamics and interactions may need stronger consideration. Prediction needs to move from a static to a real-world eco-evolutionary view.

Simulating toxicokinetic variability to identify susceptible and highly exposed populations

BACKGROUND

Toxicokinetic (TK) data needed for chemical risk assessment are not available for most chemicals. To support a greater number of chemicals, the U.S. Environmental Protection Agency (EPA) created the open-source R package "httk" (High Throughput ToxicoKinetics). The "httk" package provides functions and data tables for simulation and statistical analysis of chemical TK, including a population variability simulator that uses biometrics data from the National Health and Nutrition Examination Survey (NHANES).

OBJECTIVE

Here we modernize the "HTTK-Pop" population variability simulator based on the currently available data and literature. We provide explanations of the algorithms used by "httk" for variability simulation and uncertainty propagation.

METHODS

We updated and revised the population variability simulator in the "httk" package with the most recent NHANES biometrics (up to the 2017-18 NHANES cohort). Model equations describing glomerular filtration rate (GFR) were revised to more accurately represent physiology and population variability. The model output from the updated "httk" package was compared with the current version.

RESULTS

The revised population variability simulator in the "httk" package now provides refined, more relevant, and better justified estimations.

SIGNIFICANCE

Fulfilling the U.S. EPA's mission to provide open-source data and models for evaluations and applications by the broader scientific community, and continuously improving the accuracy of the "httk" package based on the currently available data and literature.

Literature review and evaluation of biomarkers, matrices and analytical methods for chemicals selected in the research program Human Biomonitoring for the European Union (HBM4EU)

Humans are potentially exposed to a large amount of chemicals present in the environment and in the workplace. In the European Human Biomonitoring initiative (Human Biomonitoring for the European Union = HBM4EU), acrylamide, mycotoxins (aflatoxin B1, deoxynivalenol, fumonisin B1), diisocyanates (4,4'-methylenediphenyl diisocyanate, 2,4- and 2,6-toluene diisocyanate), and pyrethroids were included among the prioritized chemicals of concern for human health. For the present literature review, the analytical methods used in worldwide biomonitoring studies for these compounds were collected and presented in comprehensive tables, including the following parameter: determined biomarker, matrix, sample amount, work-up procedure, available laboratory quality assurance and quality assessment information, analytical techniques, and limit of detection. Based on the data presented in these tables, the most suitable methods were recommended. According to the paradigm of biomonitoring, the information about two different biomarkers of exposure was evaluated: a) internal dose = parent compounds and metabolites in urine and blood; and b) the biologically effective = dose measured as blood protein adducts. Urine was the preferred matrix used for deoxynivalenol, fumonisin B1, and pyrethroids (biomarkers of internal dose). Markers of the biological effective dose were determined as hemoglobin adducts for diisocyanates and acrylamide, and as serum-albumin-adducts of aflatoxin B1 and diisocyanates. The analyses and quantitation of the protein adducts in blood or the metabolites in urine were mostly performed with LC-MS/MS or GC-MS in the presence of isotope-labeled internal standards. This review also addresses the critical aspects of the application, use and selection of biomarkers. For future biomonitoring studies, a more comprehensive approach is discussed to broaden the selection of compounds.

Exposure forecasting - ExpoCast - for data-poor chemicals in commerce and the environment

Estimates of exposure are critical to prioritize and assess chemicals based on risk posed to public health and the environment. The U.S. Environmental Protection Agency (EPA) is responsible for regulating thousands of chemicals in commerce and the environment for which exposure data are limited. Since 2009 the EPA's ExpoCast ("Exposure Forecasting") project has sought to develop the data, tools, and evaluation approaches required to generate rapid and scientifically defensible exposure predictions for the full universe of existing and proposed commercial chemicals. This review article aims to summarize issues in exposure science that have been addressed through initiatives affiliated with ExpoCast. ExpoCast research has generally focused on chemical exposure as a statistical systems problem intended to inform thousands of chemicals. The project exists as a companion to EPA's ToxCast ("Toxicity Forecasting") project which has used in vitro high-throughput screening technologies to characterize potential hazard posed by thousands of chemicals for which there are limited toxicity data. Rapid prediction of chemical exposures and in vitro-in vivo extrapolation (IVIVE) of ToxCast data allow for prioritization based upon risk of adverse outcomes due to environmental chemical exposure. ExpoCast has developed (1) integrated modeling approaches to reliably predict exposure and IVIVE dose, (2) highly efficient screening tools for chemical prioritization, (3) efficient and affordable tools for generating new exposure and dose data, and (4) easily accessible exposure databases. The development of new exposure models and databases along with the application of technologies like non-targeted analysis and machine learning have transformed exposure science for data-poor chemicals. By developing high-throughput tools for chemical exposure analytics and translating those tools into public health decisions ExpoCast research has served as a crucible for identifying and addressing exposure science knowledge gaps.

Evaluation of a rapid, generic human gestational dose model

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

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.

Advances in computational methods along the exposure to toxicological response paradigm

Human health risk assessment is a function of chemical toxicity, bioavailability to reach target biological tissues, and potential environmental exposure. These factors are complicated by many physiological, biochemical, physical and lifestyle factors. Furthermore, chemical health risk assessment is challenging in view of the large, and continually increasing, number of chemicals found in the environment. These challenges highlight the need to prioritize resources for the efficient and timely assessment of those environmental chemicals that pose greatest health risks. Computational methods, either predictive or investigative, are designed to assist in this prioritization in view of the lack of cost prohibitive in vivo experimental data. Computational methods provide specific and focused toxicity information using in vitro high throughput screening (HTS) assays. Information from the HTS assays can be converted to in vivo estimates of chemical levels in blood or target tissue, which in turn are converted to in vivo dose estimates that can be compared to exposure levels of the screened chemicals. This manuscript provides a review for the landscape of computational methods developed and used at the U.S. Environmental Protection Agency (EPA) highlighting their potentials and challenges.

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 of the InTelligence And Machine LEarning (TAME) Toolkit for Introductory Data Science, Chemical-Biological Analyses, Predictive Modeling, and Database Mining for Environmental Health Research

Research in environmental health is becoming increasingly reliant upon data science and computational methods that can more efficiently extract information from complex datasets. Data science and computational methods can be leveraged to better identify relationships between exposures to stressors in the environment and human disease outcomes, representing critical information needed to protect and improve global public health. Still, there remains a critical gap surrounding the training of researchers on these in silico methods. We aimed to address this gap by developing the inTelligence And Machine lEarning (TAME) Toolkit, promoting trainee-driven data generation, management, and analysis methods to “TAME” data in environmental health studies. Training modules were developed to provide applications-driven examples of data organization and analysis methods that can be used to address environmental health questions. Target audiences for these modules include students, post-baccalaureate and post-doctorate trainees, and professionals that are interested in expanding their skillset to include recent advances in data analysis methods relevant to environmental health, toxicology, exposure science, epidemiology, and bioinformatics/cheminformatics. Modules were developed by study coauthors using annotated script and were organized into three chapters within a GitHub Bookdown site. The first chapter of modules focuses on introductory data science, which includes the following topics: setting up R/RStudio and coding in the R environment; data organization basics; finding and visualizing data trends; high-dimensional data visualizations; and Findability, Accessibility, Interoperability, and Reusability (FAIR) data management practices. The second chapter of modules incorporates chemical-biological analyses and predictive modeling, spanning the following methods: dose-response modeling; machine learning and predictive modeling; mixtures analyses; -omics analyses; toxicokinetic modeling; and read-across toxicity predictions. The last chapter of modules was organized to provide examples on environmental health database mining and integration, including chemical exposure, health outcome, and environmental justice indicators. Training modules and associated data are publicly available online (https://uncsrp.github.io/Data-Analysis-Training-Modules/). Together, this resource provides unique opportunities to obtain introductory-level training on current data analysis methods applicable to 21st century science and environmental health.

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.

Application of an Accessible Interface for Pharmacokinetic Modeling and In Vitro to In Vivo Extrapolation

Regulatory toxicology testing has traditionally relied on in vivo methods to inform decision-making. However, scientific, practical, and ethical considerations have led to an increased interest in the use of in vitro and in silico methods to fill data gaps. While in vitro experiments have the advantage of rapid application across large chemical sets, interpretation of data coming from these non-animal methods can be challenging due to the mechanistic nature of many assays. In vitro to in vivo extrapolation (IVIVE) has emerged as a computational tool to help facilitate this task. Specifically, IVIVE uses physiologically based pharmacokinetic (PBPK) models to estimate tissue-level chemical concentrations based on various dosing parameters. This approach is used to estimate the administered dose needed to achieve in vitro bioactivity concentrations within the body. IVIVE results can be useful to inform on metrics such as margin of exposure or to prioritize potential chemicals of concern, but the PBPK models used in this approach have extensive data requirements. Thus, access to input parameters, as well as the technical requirements of applying and interpreting models, has limited the use of IVIVE as a routine part of in vitro testing. As interest in using non-animal methods for regulatory and research contexts continues to grow, our perspective is that access to computational support tools for PBPK modeling and IVIVE will be essential for facilitating broader application and acceptance of these techniques, as well as for encouraging the most scientifically sound interpretation of in vitro results. We highlight recent developments in two open-access computational support tools for PBPK modeling and IVIVE accessible via the Integrated Chemical Environment (https://ice.ntp.niehs.nih.gov/), demonstrate the types of insights these tools can provide, and discuss how these analyses may inform in vitro-based decision making.

Quo vadis blood protein adductomics?

Chemicals are measured regularly in air, food, the environment, and the workplace. Biomonitoring of chemicals in biological fluids is a tool to determine the individual exposure. Blood protein adducts of xenobiotics are a marker of both exposure and the biologically effective dose. Urinary metabolites and blood metabolites are short term exposure markers. Stable hemoglobin adducts are exposure markers of up to 120 days. Blood protein adducts are formed with many xenobiotics at different sites of the blood proteins. Newer methods apply the techniques developed in the field of proteomics. Larger adducted peptides with 20 amino acids are used for quantitation. Unfortunately, at present the methods do not reach the limits of detection obtained with the methods looking at single amino acid adducts or at chemically cleaved adducts. Therefore, to progress in the field new approaches are needed.