Designing QSARs for Parameters of High-Throughput Toxicokinetic Models Using Open-Source Descriptors

Interpretation of thyroid-relevant bioactivity data for comparison to in vivo exposures: A prioritization approach for putative chemical inhibitors of in vitro deiodinase activity

Many ToxCast assay endpoints can be mapped to molecular initiating events (MIEs) within the thyroid adverse outcome pathway (AOP) network. Herein, we provide a framework for interpretation of thyroid-relevant bioactivity data across MIEs. As a proof-of-concept, we used ToxCast data on the inhibition of deiodinase (DIO) enzymes, which convert thyroid hormones between active and inactive forms, and identified substances most likely to inhibit DIO enzymes. Data from 4 relevant cell-free in vitro assays are available for > 2000 chemicals in single concentration screening and 327 chemicals in multi-concentration screening. We filtered to identify chemicals that demonstrated inhibition for each DIO enzyme less likely to be confounded by assay interference, refining the list of putatively active chemicals from 523 to 135. In vitro bioactivity data were then used to estimate administered equivalent doses (AEDs) using a novel high-throughput toxicokinetic (HTTK) model for in vitro to in vivo extrapolation (IVIVE) of dose. To consider potential thyroid-disrupting activity in an appropriate life-stage and dose context, we extended an existing human maternal-fetal HTTK model to allow for simulations involving the first trimester of pregnancy. For many chemicals, using modeled fetal tissue concentrations produced lower AED estimates than using modeled maternal plasma concentrations alone, at least partially due to conservative assumptions in our HTTK model of complete gestation. This extensible approach for MIE groups of thyroid-related bioactivity data from ToxCast may inform further screening or analyses for potential adverse outcomes during pregnancy and development.

Combined In Vitro and In Silico Workflow to Deliver Robust, Transparent, and Contextually Rigorous Models of Bioactivity

New approach methodologies (NAMs) are an increasing priority in the field of toxicology to fill data gaps and reduce time and resources in chemical safety assessment. We describe an NAMs workflow that integrates an in vitro high-throughput bioassay with an in silico computational model. In defining this workflow, we propose, as a crucial step of in silico development, the identification of explicit "purpose contexts": a priori definitions of the scope and intent of an in silico solution, which provide natural targets for the mechanistic interpretation, validation, and output design of the model. By inspecting data from an in vitro assay measuring the displacement of fluorescent probe 8-anilino-1-naphthalenesulfonic acid (ANSA) from the serum transport protein transthyretin (TTR) as a proxy for potential disruption of thyroxine (T4) binding, in collaboration with the experimenters, we developed three relevant purpose contexts for this in silico modeling effort: (1) examination and confirmation of the in vitro assay principle via orthogonal information, (2) immediate integration with the in vitro experimental cycle to reduce costs and enhance hit rates, and (3) ultimate replacement of the use of single-concentration screening as a prioritization strategy for bioactivity testing of bulk chemical libraries. From these purpose contexts, we derived the foundations of a robust and transparent quantitative structure-activity relationship (QSAR) model that is constructively fit for purpose, characterized by first-principles mechanistic analysis, strict data quality evaluation, contextually rigorous performance testing and, finally, delivery of a quantitative recommendation schedule to simultaneously improve in vitro hit rates and in silico model learning potential.

Integration of new approach methods for the assessment of data-poor chemicals

The use of new approach methods (NAMs), including high-throughput, in vitro bioactivity data, in setting a point-of-departure (POD) will accelerate the pace of human health hazard assessments. Combining hazard and exposure predictions into a bioactivity:exposure ratio (BER) for use in risk-based prioritization and utilizing NAM-based bioactivity flags to indicate potential hazards of interest for further prediction or mechanism-based screening together comprise a prospective approach for management of substances with limited traditional toxicity testing data. In this work, we demonstrate a NAM-based assessment case study conducted via the Accelerating the Pace of Chemical Risk Assessment initiative, a consortium of international research and regulatory scientists. The primary objective was to develop a reusable and adaptable approach for addressing chemicals with limited traditional toxicity data using a NAM-based POD, BER, and bioactivity-based flags for indication of putative endocrine, developmental, neurological, and immunosuppressive effects via data generation and interpretation for 200 substances. Multiple data streams, including in silico and in vitro NAMs, were used. High-throughput transcriptomics and phenotypic profiling data, as well as targeted biochemical and cell-based assays, were combined with generic high-throughput toxicokinetic models parameterized with chemical-specific data to estimate dose for comparison to exposure predictions. This case study further enables regulatory scientists from different international purviews to utilize efficient approaches for prospective chemical management, addressing hazard and risk-based data needs, while reducing the need for animal studies. This work demonstrates the feasibility of using a battery of toxicodynamic and toxicokinetic NAMs to provide a NAM-based POD for screening-level assessment.

Enabling transparent toxicokinetic modeling for public health risk assessment

Toxicokinetic modeling describes the absorption, distribution, metabolism, and elimination of chemicals by the body. Chemical-specific in vivo toxicokinetic data is often unavailable for the thousands of chemicals in commerce. However, predictions from generalized toxicokinetic models allow for extrapolation from in vitro toxicological data, obtained via new approach methods (NAMs), to predict in vivo human health outcomes and provide key information on chemicals for public health risk assessment. The httk R package provides an open-source software tool containing a suite of generalized toxicokinetic models covering various exposure scenarios, a library of chemical-specific data from peer-reviewed high-throughput toxicokinetic (HTTK) studies, and other utility functions to parameterize and evaluate toxicokinetic models. Generalized HTTK models in httk use the open-source language MCSim to describe the compartmental and physiologically based toxicokinetics (PBTK). New HTTK models may be integrated into httk with a model description code file (C script generated via MCSim) and a model documentation file (R script). httk provides a series of functionalities such as model parameterization, in vivo-derived data for evaluating model predictions, unit conversion, Monte Carlo simulations for uncertainty propagation and biological variability, and other model utilities. Here, we describe in detail how to add new HTTK models into the httk package to leverage its pre-existing data and functionality. As a demonstration, we describe the integration of a gas inhalation PBTK model. The intention of httk is to provide a transparent, open-source tool for toxicokinetics, bioinformatics, and public health risk assessment that makes use of publicly available data on more than one thousand chemicals.

The environmental neuroactive chemicals list of prioritized substances for human biomonitoring and neurotoxicity testing: A database and high-throughput toxicokinetics approach

There is a diversity of chemicals to which humans are potentially exposed. Few of these chemicals have linked human biomonitoring data, and most have very limited neurotoxicity testing. Of particular concern are environmental exposures impacting children, who constitute a population of heightened susceptibility due to rapid neural growth and plasticity, yet lack biomonitoring data compared to other age/population subgroups. This study set out to develop a prioritized list of neuroactive substances, titled the Environmental NeuRoactIve CHemicals (ENRICH) list, to be used as a defined screening library in the evaluation of human biological samples, with emphasis on early childhood-relevant environmental exposures. In silico database mining approaches were used to prioritize chemicals based upon likelihood of neuroactivity, human exposure, and feasible detection in biological samples. Evidence of neuroactivity was compiled across in vitro high-throughput screening, animal testing, and/or human epidemiological findings. Chemicals were considered for their likelihood of human exposure and detection presence in biological samples (including metabolites), with additional evidence indicating presence within child-relevant products. The resulting list of 1827 chemicals were ranked using a Chemical Prioritization Index. Manual inclusion/exclusion criteria were employed for the top-ranking chemical candidates to ensure that chemicals were within the study's scope (i.e., environmentally relevant) and, for the purposes of biomonitoring, had properties amenable to mass spectrometry methods. These elements were combined to produce the ENRICH list of 250 top-ranking chemicals, spanning pesticides and those used in home maintenance, personal care, cleaning products, vehicles, arts and crafts, and consumer electronics, among other sources. Chemicals were additionally evaluated for high-throughput toxicokinetics to predict how much of a chemical and/or its metabolite(s) may reach urine, as an example biological matrix for practical use in biomonitoring efforts. This novel study couples databases and in silico-based predictions to prioritize chemicals in the environment with potential neurological impacts for future study.

Incorporating new approach methods (NAMs) data in dose-response assessments: The future is now!

Regulatory dose-response assessments traditionally rely on in vivo data and default assumptions. New Approach Methods (NAMs) present considerable opportunities to both augment traditional dose-response assessments and accelerate the evaluation of new/data-poor chemicals. This review aimed to determine the potential utilization of NAMs through a unified conceptual framework that compartmentalizes derivation of toxicity values into five sequential Key Dose-response Modules (KDMs): (1) point-of-departure (POD) determination, (2) test system-to-human (e.g. inter-species) toxicokinetics and (3) toxicodynamics, (4) human population (intra-species) variability in toxicodynamics, and (5) toxicokinetics. After using several "traditional" dose-response assessments to illustrate this framework, a review is presented where existing NAMs, including in silico, in vitro, and in vivo approaches, might be applied across KDMs. Further, the false dichotomy between "traditional" and NAMs-derived data sources is broken down by organizing dose-response assessments into a matrix where each KDM has Tiers of increasing precision and confidence: Tier 0: Default/generic values, Tier 1: Computational predictions, Tier 2: Surrogate measurements, and Tier 3: Direct measurements. These findings demonstrated that although many publications promote the use of NAMs in KDMs (1) for POD determination and (5) for human population toxicokinetics, the proposed matrix of KDMs and Tiers reveals additional immediate opportunities for NAMs to be integrated across other KDMs. Further, critical needs were identified for developing NAMs to improve in vitro dosimetry and quantify test system and human population toxicodynamics. Overall, broadening the integration of NAMs across the steps of dose-response assessment promises to yield higher throughput, less animal-dependent, and more science-based toxicity values for protecting human health.

Potential impact of organophosphate esters on thyroid eye disease based on machine learning and molecular docking

Organophosphate esters (OPEs) are widely used as flame retardants and plasticizers in daily commodities and building materials. Some OPEs, acting as agonists of the thyroid-stimulating hormone receptor (TSHR), may contribute to the development of thyroid eye disease (TED). This study analyzes the serum and urine of patients and control groups, using machine learning and molecular docking to investigate the potential impact of OPEs on TED. Results indicate significantly higher concentrations of OPEs and di-OPEs of TED patients compared to controls (Mann-Whitney U test, p < 0.05). Aryl OPEs exhibit the strongest binding affinity with TSHR. We developed a predictive model for OPE-TSHR affinity to explore the impact of OPE structural features on TSHR activity and effectively capture the complex relationships between changes in OPE side chains and their effects on TSHR. Predictions from the USEPA's database indicate that 28 % of 1011 OPEs have a tendency to bind with TSHR. Furthermore, a high-accuracy classification model successfully identified key substructures associated with high affinity for TSHR. This study not only enhances our understanding of the complex relationship between the structural diversity of OPEs and their thyroid impact but also offers molecular design insights to prevent releasing OPEs with high thyroid harm potential into the environment.

Chemical space exploration with Molpher: Generating and assessing a glucocorticoid receptor ligand library

Computational exploration of chemical space is crucial in modern cheminformatics research for accelerating the discovery of new biologically active compounds. In this study, we present a detailed analysis of the chemical library of potential glucocorticoid receptor (GR) ligands generated by the molecular generator, Molpher. To generate the targeted GR library and construct the classification models, structures from the ChEMBL database as well as from the internal IMG library, which was experimentally screened for biological activity in the primary luciferase reporter cell assay, were utilized. The composition of the targeted GR ligand library was compared with a reference library that randomly samples chemical space. A random forest model was used to determine the biological activity of ligands, incorporating its applicability domain using conformal prediction. It was demonstrated that the GR library is significantly enriched with GR ligands compared to the random library. Furthermore, a prospective analysis demonstrated that Molpher successfully designed compounds, which were subsequently experimentally confirmed to be active on the GR. A collection of 34 potential new GR ligands was also identified. Moreover, an important contribution of this study is the establishment of a comprehensive workflow for evaluating computationally generated ligands, particularly those with potential activity against targets that are challenging to dock.

Discovery Phase Agrochemical Predictive Safety Assessment Using High Content In Vitro Data to Estimate an In Vivo Toxicity Point of Departure

Utilization of in vitro (cellular) techniques, like Cell Painting and transcriptomics, could provide powerful tools for agrochemical candidate sorting and selection in the discovery process. However, using these models generates challenges translating in vitro concentrations to the corresponding in vivo exposures. Physiologically based pharmacokinetic (PBPK) modeling provides a framework for quantitative in vitro to in vivo extrapolation (IVIVE). We tested whether in vivo (rat liver) transcriptomic and apical points of departure (PODs) could be accurately predicted from in vitro (rat hepatocyte or human HepaRG) transcriptomic PODs or HepaRG Cell Painting PODs using PBPK modeling. We compared two PBPK models, the ADMET predictor and the httk R package, and found httk to predict the in vivo PODs more accurately. Our findings suggest that a rat liver apical and transcriptomic POD can be estimated utilizing a combination of in vitro transcriptome-based PODs coupled with PBPK modeling for IVIVE. Thus, high content in vitro data can be translated with modest accuracy to in vivo models of ultimate regulatory importance to help select agrochemical analogs in early stage discovery program.

Investigation of the Binding Fraction of PFAS in Human Plasma and Underlying Mechanisms Based on Machine Learning and Molecular Dynamics Simulation

More than 7000 per- and polyfluorinated alkyl substances (PFAS) have been documented in the U.S. Environmental Protection Agency's CompTox Chemicals database. These PFAS can be used in a broad range of industrial and consumer applications but may pose potential environmental issues and health risks. However, little is known about emerging PFAS bioaccumulation to assess their chemical safety. This study focuses specifically on the large and high-quality data set of fluorochemicals from the related environmental and pharmaceutical chemicals databases, and machine learning (ML) models were developed for the classification prediction of the unbound fraction of compounds in plasma. A comprehensive evaluation of the ML models shows that the best blending model yields an accuracy of 0.901 for the test set. The predictions suggest that most PFAS (∼92%) have a high binding fraction in plasma. Introduction of alkaline amino groups is likely to reduce the binding affinities of PFAS with plasma proteins. Molecular dynamics simulations indicate a clear distinction between the high and low binding fractions of PFAS. These computational workflows can be used to predict the bioaccumulation of emerging PFAS and are also helpful for the molecular design of PFAS to prevent the release of high-bioaccumulation compounds into the environment.

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.

Category-Based Toxicokinetic Evaluations of Data-Poor Per- and Polyfluoroalkyl Substances (PFAS) using Gas Chromatography Coupled with Mass Spectrometry

Concern over per- and polyfluoroalkyl substances (PFAS) has increased as more is learned about their environmental presence, persistence, and bioaccumulative potential. The limited monitoring, toxicokinetic (TK), and toxicologic data available are inadequate to inform risk across this diverse domain. Here, 73 PFAS were selected for in vitro TK evaluation to expand knowledge across lesser-studied PFAS alcohols, amides, and acrylates. Targeted methods developed using gas chromatography-tandem mass spectrometry (GC-MS/MS) were used to measure human plasma protein binding and hepatocyte clearance. Forty-three PFAS were successfully evaluated in plasma, with fraction unbound (f_up) values ranging from 0.004 to 1. With a median f_up of 0.09 (i.e., 91% bound), these PFAS are highly bound but exhibit 10-fold lower binding than legacy perfluoroalkyl acids recently evaluated. Thirty PFAS evaluated in the hepatocyte clearance assay showed abiotic loss, with many exceeding 60% loss within 60 min. Metabolic clearance was noted for 11 of the 13 that were successfully evaluated, with rates up to 49.9 μL/(min × million cells). The chemical transformation simulator revealed potential (bio)transformation products to consider. This effort provides critical information to evaluate PFAS for which volatility, metabolism, and other routes of transformation are likely to modulate their environmental fates.

Occupational exposure to veterinary antibiotics: Pharmacokinetics of enrofloxacin in humans after dermal, inhalation and oral uptake - A clinical study

Occupational exposure to veterinary antibiotics in hen houses at poultry feeding farms was demonstrated by biomonitoring campaigns in the past. The objective of this study was to investigate pharmacokinetics of three uptake routes: dermal, oral and inhaled. In an open-label cross-over study, six healthy volunteers were exposed to single occupational relevant doses of enrofloxacin. Plasma and urine samples were analysed for enrofloxacin and ciprofloxacin. Physiologically based pharmacokinetic (PBPK) modelling based on bioanalysis data showed underestimation for the elimination rate in comparison to experimental data pointing towards a lack of sufficient ADME information and limitations of available physico-chemical properties of the parent drug. The data obtained in this study indicate that oral uptake with its various sources, e.g. airborne enrofloxacin, direct hand-mouth contact, is the major source for occupational exposure to enrofloxacin in hen houses. Dermal exposure was considered negligible.

From Protein Sequence to Structure: The Next Frontier in Cross-Species Extrapolation for Chemical Safety Evaluations

Computational screening for potentially bioactive molecules using advanced molecular modeling approaches including molecular docking and molecular dynamic simulation is mainstream in certain fields like drug discovery. Significant advances in computationally predicting protein structures from sequence information have also expanded the availability of structures for nonmodel species. Therefore, the objective of the present study was to develop an analysis pipeline to harness the power of these bioinformatics approaches for cross-species extrapolation for evaluating chemical safety. The Sequence Alignment to Predict Across Species Susceptibility (SeqAPASS) tool compares protein-sequence similarity across species for conservation of known chemical targets, providing an initial line of evidence for extrapolation of toxicity knowledge. However, with the development of structural models from tools like the Iterative Threading ASSEmbly Refinement (ITASSER), analyses of protein structural conservation can be included to add further lines of evidence and generate protein models across species. Models generated through such a pipeline could then be used for advanced molecular modeling approaches in the context of species extrapolation. Two case examples illustrating this pipeline from SeqAPASS sequences to I-TASSER-generated protein structures were created for human liver fatty acid-binding protein (LFABP) and androgen receptor (AR). Ninety-nine LFABP and 268 AR protein models representing diverse species were generated and analyzed for conservation using template modeling (TM)-align. The results from the structural comparisons were in line with the sequence-based SeqAPASS workflow, adding further evidence of LFABL and AR conservation across vertebrate species. The present study lays the foundation for expanding the capabilities of the web-based SeqAPASS tool to include structural comparisons for species extrapolation, facilitating more rapid and efficient toxicological assessments among species with limited or no existing toxicity data. Environ Toxicol Chem 2023;42:463-474. © 2022 SETAC. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Machine learning and artificial intelligence in physiologically based pharmacokinetic modeling

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

A Machine Learning Model to Estimate Toxicokinetic Half-Lives of Per- and Polyfluoro-Alkyl Substances (PFAS) in Multiple Species

Per- and polyfluoroalkyl substances (PFAS) are a diverse group of man-made chemicals that are commonly found in body tissues. The toxicokinetics of most PFAS are currently uncharacterized, but long half-lives (t½) have been observed in some cases. Knowledge of chemical-specific t½ is necessary for exposure reconstruction and extrapolation from toxicological studies. We used an ensemble machine learning method, random forest, to model the existing in vivo measured t½ across four species (human, monkey, rat, mouse) and eleven PFAS. Mechanistically motivated descriptors were examined, including two types of surrogates for renal transporters: (1) physiological descriptors, including kidney geometry, for renal transporter expression and (2) structural similarity of defluorinated PFAS to endogenous chemicals for transporter affinity. We developed a classification model for t½ (Bin 1: <12 h; Bin 2: <1 week; Bin 3: <2 months; Bin 4: >2 months). The model had an accuracy of 86.1% in contrast to 32.2% for a y-randomized null model. A total of 3890 compounds were within domain of the model, and t½ was predicted using the bin medians: 4.9 h, 2.2 days, 33 days, and 3.3 years. For human t½, 56% of PFAS were classified in Bin 4, 7% were classified in Bin 3, and 37% were classified in Bin 2. This model synthesizes the limited available data to allow tentative extrapolation and prioritization.

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.

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.

The chemical landscape of high-throughput new approach methodologies for exposure

The rapid characterization of risk to humans and ecosystems from exogenous chemicals requires information on both hazard and exposure. The U.S. Environmental Protection Agency's ToxCast program and the interagency Tox21 initiative have screened thousands of chemicals in various high-throughput (HT) assay systems for in vitro bioactivity. EPA's ExpoCast program is developing complementary HT methods for characterizing the human and ecological exposures necessary to interpret HT hazard data in a real-world risk context. These new approach methodologies (NAMs) for exposure include computational and analytical tools for characterizing multiple components of the complex pathways chemicals take from their source to human and ecological receptors. Here, we analyze the landscape of exposure NAMs developed in ExpoCast in the context of various chemical lists of scientific and regulatory interest, including the ToxCast and Tox21 libraries and the Toxic Substances Control Act (TSCA) inventory. We examine the landscape of traditional and exposure NAM data covering chemical use, emission, environmental fate, toxicokinetics, and ultimately external and internal exposure. We consider new chemical descriptors, machine learning models that draw inferences from existing data, high-throughput exposure models, statistical frameworks that integrate multiple model predictions, and non-targeted analytical screening methods that generate new HT monitoring information. We demonstrate that exposure NAMs drastically improve the coverage of the chemical landscape compared to traditional approaches and recommend a set of research activities to further expand the development of HT exposure data for application to risk characterization. Continuing to develop exposure NAMs to fill priority data gaps identified here will improve the availability and defensibility of risk-based metrics for use in chemical prioritization and screening. IMPACT: This analysis describes the current state of exposure assessment-based new approach methodologies across varied chemical landscapes and provides recommendations for filling key data gaps.

Estimating provisional margins of exposure for data-poor chemicals using high-throughput computational methods

Current computational technologies hold promise for prioritizing the testing of the thousands of chemicals in commerce. Here, a case study is presented demonstrating comparative risk-prioritization approaches based on the ratio of surrogate hazard and exposure data, called margins of exposure (MoEs). Exposures were estimated using a U.S. EPA’s ExpoCast predictive model (SEEM3) results and estimates of bioactivity were predicted using: 1) Oral equivalent doses (OEDs) derived from U.S. EPA’s ToxCast high-throughput screening program, together with in vitro to in vivo extrapolation and 2) thresholds of toxicological concern (TTCs) determined using a structure-based decision-tree using the Toxtree open source software. To ground-truth these computational approaches, we compared the MoEs based on predicted noncancer TTC and OED values to those derived using the traditional method of deriving points of departure from no-observed adverse effect levels (NOAELs) from in vivo oral exposures in rodents. TTC-based MoEs were lower than NOAEL-based MoEs for 520 out of 522 (99.6%) compounds in this smaller overlapping dataset, but were relatively well correlated with the same (r² = 0.59). TTC-based MoEs were also lower than OED-based MoEs for 590 (83.2%) of the 709 evaluated chemicals, indicating that TTCs may serve as a conservative surrogate in the absence of chemical-specific experimental data. The TTC-based MoE prioritization process was then applied to over 45,000 curated environmental chemical structures as a proof-of-concept for high-throughput prioritization using TTC-based MoEs. This study demonstrates the utility of exploiting existing computational methods at the pre-assessment phase of a tiered risk-based approach to quickly, and conservatively, prioritize thousands of untested chemicals for further study.

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.

Artificial Intelligence-Based Toxicity Prediction of Environmental Chemicals: Future Directions for Chemical Management Applications

Recently, research on the development of artificial intelligence (AI)-based computational toxicology models that predict toxicity without the use of animal testing has emerged because of the rapid development of computer technology. Various computational toxicology techniques that predict toxicity based on the structure of chemical substances are gaining attention, including the quantitative structure-activity relationship. To understand the recent development of these models, we analyzed the databases, molecular descriptors, fingerprints, and algorithms considered in recent studies. Based on a selection of 96 papers published since 2014, we found that AI models have been developed to predict approximately 30 different toxicity end points using more than 20 toxicity databases. For model development, molecular access system and extended-connectivity fingerprints are the most commonly used molecular descriptors. The most used algorithm among the machine learning techniques is the random forest, while the most used algorithm among the deep learning techniques is a deep neural network. The use of AI technology in the development of toxicity prediction models is a new concept that will aid in achieving a scientific accord and meet regulatory applications. The comprehensive overview provided in this study will provide a useful guide for the further development and application of toxicity prediction 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.

Considerations for Improving Metabolism Predictions for In Vitro to In Vivo Extrapolation

High-throughput (HT) in vitro to in vivo extrapolation (IVIVE) is an integral component in new approach method (NAM)-based risk assessment paradigms, for rapidly translating in vitro toxicity assay results into the context of in vivo exposure. When coupled with rapid exposure predictions, HT-IVIVE supports the use of HT in vitro assays for risk-based chemical prioritization. However, the reliability of prioritization based on HT bioactivity data and HT-IVIVE can be limited as the domain of applicability of current HT-IVIVE is generally restricted to intrinsic clearance measured primarily in pharmaceutical compounds. Further, current approaches only consider parent chemical toxicity. These limitations occur because current state-of-the-art HT prediction tools for clearance and metabolite kinetics do not provide reliable data to support HT-IVIVE. This paper discusses current challenges in implementation of IVIVE for prioritization and risk assessment and recommends a path forward for addressing the most pressing needs and expanding the utility of IVIVE.

EFSA Panel on Plant Protection Products and their Residues (EFSA PPR Panel), Hernandez‐Jerez AF, Adriaanse P, Aldrich A, Berny P, Coja T, Duquesne S, Focks A, Marinovich M, Millet M, Pelkonen O, Pieper S, Tiktak A, Topping CJ, Widenfalk A, Wilks M, Wolterink G, Gundert‐Remy U, Louisse J, Rudaz S, Testai E, Lostia A, Dorne J, Parra Morte JM. Scientific Opinion of the Scientific Panel on Plant Protection Products and their Residues (PPR Panel) on testing and interpretation of comparative in vitro metabolism studies

EFSA asked the Panel on Plant Protection Products and their residues to deliver a Scientific Opinion on testing and interpretation of comparative in vitro metabolism studies for both new active substances and existing ones. The main aim of comparative in vitro metabolism studies of pesticide active substances is to evaluate whether all significant metabolites formed in the human in vitro test system, as a surrogate of the in vivo situation, are also present at comparable level in animal species tested in toxicological studies and, therefore, if their potential toxicity has been appropriately covered by animal studies. The studies may also help to decide which animal model, with regard to a particular compound, is the most relevant for humans. In the experimental strategy, primary hepatocytes in suspension or culture are recommended since hepatocytes are considered the most representative in vitro system for prediction of in vivo metabolites. The experimental design of 3 × 3 × 3 (concentrations, time points, technical replicates, on pooled hepatocytes) will maximise the chance to identify unique (UHM) and disproportionate (DHM) human metabolites. When DHM and UHM are being assessed, test item-related radioactivity recovery and metabolite profile are the most important parameters. Subsequently, structural characterisation of the assigned metabolites is performed with appropriate analytical techniques. In toxicological assessment of metabolites, the uncertainty factor approach is the first alternative to testing option, followed by new approach methodologies (QSAR, read-across, in vitro methods), and only if these fail, in vivo animal toxicity studies may be performed. Knowledge of in vitro metabolites in human and animal hepatocytes would enable toxicological evaluation of all metabolites of concern, and, furthermore, add useful pieces of information for detection and evaluation of metabolites in different matrices (crops, livestock, environment), improve biomonitoring efforts via better toxicokinetic understanding, and ultimately, develop regulatory schemes employing physiologically based or physiology-mimicking in silico and/or in vitro test systems to anticipate the exposure of humans to potentially hazardous substances in plant protection products.

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.

Opportunities and challenges related to saturation of toxicokinetic processes: Implications for risk assessment

Top dose selection for repeated dose animal studies has generally focused on identification of apical endpoints, use of the limit dose, or determination of a maximum tolerated dose (MTD). The intent is to optimize the ability of toxicity tests performed in a small number of animals to detect effects for hazard identification. An alternative approach, the kinetically derived maximum dose (KMD), has been proposed as a mechanism to integrate toxicokinetic (TK) data into the dose selection process. The approach refers to the dose above which the systemic exposures depart from being proportional to external doses. This non-linear external-internal dose relationship arises from saturation or limitation of TK process(es), such as absorption or metabolism. The importance of TK information is widely acknowledged when assessing human health risks arising from exposures to environmental chemicals, as TK determines the amount of chemical at potential sites of toxicological responses. However, there have been differing opinions and interpretations within the scientific and regulatory communities related to the validity and application of the KMD concept. A multi-stakeholder working group, led by the Health and Environmental Sciences Institute (HESI), was formed to provide an opportunity for impacted stakeholders to address commonly raised scientific and technical issues related to this topic and, more specifically, a weight of evidence approach is recommended to inform design and dose selection for repeated dose animal studies. Commonly raised challenges related to the use of TK data for dose selection are discussed, recommendations are provided, and illustrative case examples are provided to address these challenges or refute misconceptions.

Prediction Models for Agonists and Antagonists of Molecular Initiation Events for Toxicity Pathways Using an Improved Deep-Learning-Based Quantitative Structure-Activity Relationship System

In silico approaches have been studied intensively to assess the toxicological risk of various chemical compounds as alternatives to traditional in vivo animal tests. Among these approaches, quantitative structure-activity relationship (QSAR) analysis has the advantages that it is able to construct models to predict the biological properties of chemicals based on structural information. Previously, we reported a deep learning (DL) algorithm-based QSAR approach called DeepSnap-DL for high-performance prediction modeling of the agonist and antagonist activity of key molecules in molecular initiating events in toxicological pathways using optimized hyperparameters. In the present study, to achieve high throughput in the DeepSnap-DL system-which consists of the preparation of three-dimensional molecular structures of chemical compounds, the generation of snapshot images from the three-dimensional chemical structures, DL, and statistical calculations-we propose an improved DeepSnap-DL approach. Using this improved system, we constructed 59 prediction models for the agonist and antagonist activity of key molecules in the Tox21 10K library. The results indicate that modeling of the agonist and antagonist activity with high prediction performance and high throughput can be achieved by optimizing suitable parameters in the improved DeepSnap-DL system.

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

INTRODUCTION

Toxicity data are unavailable for many thousands of chemicals in commerce and the environment. Therefore, risk assessors need to rapidly screen these chemicals for potential risk to public health. High-throughput screening (HTS) for in vitro bioactivity, when used with high-throughput toxicokinetic (HTTK) data and models, allows characterization of these thousands of chemicals.

AREAS COVERED

This review covers generic physiologically based toxicokinetic (PBTK) models and high-throughput PBTK modeling for in vitro-in vivo extrapolation (IVIVE) of HTS data. We focus on 'httk', a public, open-source set of computational modeling tools and in vitro toxicokinetic (TK) data.

EXPERT OPINION

HTTK benefits chemical risk assessors with its ability to support rapid chemical screening/prioritization, perform IVIVE, and provide provisional TK modeling for large numbers of chemicals using only limited chemical-specific data. Although generic TK model design can increase prediction uncertainty, these models provide offsetting benefits by increasing model implementation accuracy. Also, public distribution of the models and data enhances reproducibility. For the httk package, the modular and open-source design can enable the tool to be used and continuously improved by a broad user community in support of the critical need for high-throughput chemical prioritization and rapid dose estimation to facilitate rapid hazard assessments.