Advancing internal exposure and physiologically-based toxicokinetic modeling for 21st-century risk assessments

Human risk assessment through development and application of a physiologically based toxicokinetic model for 4-tert-octylphenol

4-tert-octylphenol (4-tert-OP) is an ecologically hazardous substance, and exposure to it in the environment has been consistently reported in the past. Despite the hazards and widespread exposure to 4-tert-OP, tools for scientific assessment of 4-tert-OP exposure risk level in humans are lacking. The main purpose of this study was to develop a physiologically-based-toxicokinetic (PBTK) model for 4-tert-OP and to perform quantitative risk assessment of 4-tert-OP in various population groups using the established model. Based on the results of toxicokinetic experiments on male rats, the PBTK model for 4-tert-OP was established and verified, and this was converted to a model for humans through interspecies extrapolation. Based on the previously reported no-observed-adverse-effect-levels for rats, it was possible to estimate the 4-tert-OP reference dose in humans through reverse dosimetry using the model. Biomonitoring data derived from various population groups were applied to the human PBTK model to calculate external exposures and margin of safety for 4-tert-OP for each population group. The PBTK model established in this study adequately explained the toxicokinetic experimental values at acceptable levels and was able to quantitatively predict the 4-tert-OP exposure level in the testes related to male reproductive toxicity. In addition, the degree of external exposure to 4-tert-OP could be scientifically estimated based on biomonitoring values derived from various biological matrices. The reference doses for systemic and reproductive toxicity caused by 4-tert-OP in male humans were calculated to be 0.16 and 1.12 mg/kg/day, respectively. The mean external exposure to 4-tert-OP in each population group estimated based on plasma and urine biomonitoring data was 0.04-66.24 mg/kg/day, showing very large exposure diversity between groups. Exposure risks to 4-tert-OP in populations ranged from safe to risky, suggesting the need for continued monitoring and risk management of 4-tert-OP worldwide. This study provides valuable scientific insight regarding the 4-tert-OP human risk assessment.

A Predictive Toxicokinetic Model for Nickel Leaching from Vascular Stents

In vitro testing methods offer valuable insights into the corrosion vulnerability of metal implants and enable prompt comparison between devices. However, they fall short in predicting the extent of leaching and the biodistribution of implant byproducts under in vivo conditions. Physiologically based toxicokinetic (PBTK) models are capable of quantitatively establishing such correlations and therefore provide a powerful tool in advancing nonclinical methods to test medical implants and assess patient exposure to implant debris. In this study, we present a multicompartment PBTK model and a simulation engine for toxicological risk assessment of vascular stents. The mathematical model consists of a detailed set of constitutive equations that describe the transfer of nickel ions from the device to peri-implant tissue and circulation and the nickel mass exchange between blood and the various tissues/organs and excreta. Model parameterization was performed using (1) in-house-produced data from immersion testing to compute the device-specific diffusion parameters and (2) full-scale animal in situ implantation studies to extract the mammalian-specific biokinetic functions that characterize the time-dependent biodistribution of the released ions. The PBTK model was put to the test using a simulation engine to estimate the concentration-time profiles, along with confidence intervals through probabilistic Monte Carlo, of nickel ions leaching from the implanted devices and determine if permissible exposure limits are exceeded. The model-derived output demonstrated prognostic conformity with reported experimental data, indicating that it may provide the basis for the broader use of modeling and simulation tools to guide the optimal design of implantable devices in compliance with exposure limits and other regulatory requirements.

Big Question to Developing Solutions: A Decade of Progress in the Development of Aquatic New Approach Methodologies from 2012 to 2022

In 2012, 20 key questions related to hazard and exposure assessment and environmental and health risks of pharmaceuticals and personal care products in the natural environment were identified. A decade later, this article examines the current level of knowledge around one of the lowest-ranking questions at that time, number 19: "Can nonanimal testing methods be developed that will provide equivalent or better hazard data compared with current in vivo methods?" The inclusion of alternative methods that replace, reduce, or refine animal testing within the regulatory context of risk and hazard assessment of chemicals generally faces many hurdles, although this varies both by organism (human-centric vs. other), sector, and geographical region or country. Focusing on the past 10 years, only works that might reasonably be considered to contribute to advancements in the field of aquatic environmental risk assessment are highlighted. Particular attention is paid to methods of contemporary interest and importance, representing progress in (1) the development of methods which provide equivalent or better data compared with current in vivo methods such as bioaccumulation, (2) weight of evidence, or (3) -omic-based applications. Evolution and convergence of these risk assessment areas offer the basis for fundamental frameshifts in how data are collated and used for the protection of taxa across the breadth of the aquatic environment. Looking to the future, we are at a tipping point, with a need for a global and inclusive approach to establish consensus. Bringing together these methods (both new and old) for regulatory assessment and decision-making will require a concerted effort and orchestration. Environ Toxicol Chem 2024;43:559-574. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Alternatives to animal testing in toxicity testing: Current status and future perspectives in food safety assessments

The development of alternative methods to animal testing has gained great momentum since Russel and Burch introduced the "3Rs" concept of Reduction, Refinement, and Replacement of animals in safety testing in 1959. Several alternatives to animal testing have since been introduced, including but not limited to in vitro and in chemico test systems, in silico models, and computational models (e.g., [quantitative] structural activity relationship models, high-throughput screens, organ-on-chip models, and genomics or bioinformatics) to predict chemical toxicity. Furthermore, several agencies have developed robust integrated testing strategies to determine chemical toxicity. The cosmetics sector is pioneering the adoption of alternative methodologies for safety evaluations, and other sectors are aiming to completely abandon animal testing by 2035. However, beyond the use of in vitro genetic testing, agencies regulating the food industry have been slow to implement alternative methodologies into safety evaluations compared with other sectors; setting health-based guidance values for food ingredients requires data from systemic toxicity, and to date, no standalone validated alternative models to assess systemic toxicity exist. The abovementioned models show promise for assessing systemic toxicity with further research. In this paper, we review the current alternatives and their applicability and limitations in food safety evaluations.

A physiologically based pharmacokinetic model to optimize the dosage regimen and withdrawal time of cefquinome in pigs

Cefquinome is widely used to treat respiratory tract diseases of swine. While extra-label dosages of cefquinome could improve clinical efficacy, they might lead to excessively high residues in animal-derived food. In this study, a physiologically based pharmacokinetic (PBPK) model was calibrated based on the published data and a microdialysis experiment to assess the dosage efficiency and food safety. For the microdialysis experiment, in vitro/in vivo relative recovery and concentration-time curves of cefquinome in the lung interstitium were investigated. This PBPK model is available to predict the drug concentrations in the muscle, kidney, liver, plasma, and lung interstitial fluid. Concentration-time curves of 1000 virtual animals in different tissues were simulated by applying sensitivity and Monte Carlo analyses. By integrating pharmacokinetic/pharmacodynamic target parameters, cefquinome delivered at 3-5 mg/kg twice daily is advised for the effective control of respiratory tract infections of nursery pig, which the bodyweight is around 25 kg. Based on the predicted cefquinome concentrations in edible tissues, the withdrawal interval is 2 and 3 days for label and the extra-label doses, respectively. This study provides a useful tool to optimize the dosage regimen of cefquinome against respiratory tract infections and predicts the concentration of cefquinome residues in edible tissues. This information would be helpful to improve the food safety and guide rational drug usage.

Evaluating the food safety and risk assessment evidence-base of polyethylene terephthalate oligomers: A systematic evidence map

BACKGROUND

The presence of polyethylene terephthalate (PET) oligomers in food contact materials (FCMs) is well-documented. Consumers are exposed through their migration into foods and beverages; however, there is no specific guidance for their safety evaluation.

OBJECTIVES

This systematic evidence map (SEM) aims to identify and organize existing knowledge and associated gaps in hazard and exposure information on 34 PET oligomers to support regulatory decision-making.

METHODS

The methodology for this SEM was recently registered. A systematic search in bibliographic and gray literature sources was conducted and studies evaluated for inclusion according to the Populations, Exposures, Comparators, Outcomes, and Study type (PECOS) framework. Inclusion criteria were designed to record hazard and exposure information for all 34 PET oligomers and coded into the following evidence streams: human, animal, organism (non-animal), ex vivo, in vitro, in silico, migration, hydrolysis, and absorption, distribution, metabolism, excretion/toxicokinetics/pharmacokinetics (ADME/TK/PK) studies. Relevant information was extracted from eligible studies and synthesized according to the protocol.

RESULTS

Literature searches yielded 7445 unique records, of which 96 were included. Data comprised migration (560 entries), ADME/TK/PK-related (253 entries), health/bioactivity (98 entries) and very few hydrolysis studies (7 entries). Cyclic oligomers were studied more frequently than linear PET oligomers. In vitro results indicated that hydrolysis of cyclic oligomers generated a mixture of linear oligomers, but not monomers, potentially allowing their absorption in the gastrointestinal tract. Cyclic dimers, linear trimers and the respective smaller oligomers exhibit physico-chemical properties making oral absorption more likely. Information on health/bioactivity effects of oligomers was almost non-existent, except for limited data on mutagenicity.

CONCLUSIONS

This SEM revealed substantial deficiencies in the available evidence on ADME/TK/PK, hydrolysis, and health/bioactivity effects of PET oligomers, currently preventing appropriate risk assessment. It is essential to develop more systematic and tiered approaches to address the identified research needs and assess the risks of PET oligomers.

New Approach Methodologies for the Endocrine Activity Toolbox: Environmental Assessment for Fish and Amphibians

Multiple in vivo test guidelines focusing on the estrogen, androgen, thyroid, and steroidogenesis pathways have been developed and validated for mammals, amphibians, or fish. However, these tests are resource-intensive and often use a large number of laboratory animals. Developing alternatives for in vivo tests is consistent with the replacement, reduction, and refinement principles for animal welfare considerations, which are supported by increasing mandates to move toward an "animal-free" testing paradigm worldwide. New approach methodologies (NAMs) hold great promise to identify molecular, cellular, and tissue changes that can be used to predict effects reliably and more efficiently at the individual level (and potentially on populations) while reducing the number of animals used in (eco)toxicological testing for endocrine disruption. In a collaborative effort, experts from government, academia, and industry met in 2020 to discuss the current challenges of testing for endocrine activity assessment for fish and amphibians. Continuing this cross-sector initiative, our review focuses on the current state of the science regarding the use of NAMs to identify chemical-induced endocrine effects. The present study highlights the challenges of using NAMs for safety assessment and what work is needed to reduce their uncertainties and increase their acceptance in regulatory processes. We have reviewed the current NAMs available for endocrine activity assessment including in silico, in vitro, and eleutheroembryo models. New approach methodologies can be integrated as part of a weight-of-evidence approach for hazard or risk assessment using the adverse outcome pathway framework. The development and utilization of NAMs not only allows for replacement, reduction, and refinement of animal testing but can also provide robust and fit-for-purpose methods to identify chemicals acting via endocrine mechanisms. Environ Toxicol Chem 2023;42:757-777. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Quantitative in vitro to in vivo extrapolation of genotoxicity data provides protective estimates of in vivo dose

Genotoxicity assessment is a critical component in the development and evaluation of chemicals. Traditional genotoxicity assays (i.e., mutagenicity, clastogenicity, and aneugenicity) have been limited to dichotomous hazard classification, while other toxicity endpoints are assessed through quantitative determination of points-of-departures (PODs) for setting exposure limits. The more recent higher-throughput in vitro genotoxicity assays, many of which also provide mechanistic information, offer a powerful approach for determining defined PODs for potency ranking and risk assessment. In order to obtain relevant human dose context from the in vitro assays, in vitro to in vivo extrapolation (IVIVE) models are required to determine what dose would elicit a concentration in the body demonstrated to be genotoxic using in vitro assays. Previous work has demonstrated that application of IVIVE models to in vitro bioactivity data can provide PODs that are protective of human health, but there has been no evaluation of how these models perform with in vitro genotoxicity data. Thus, the Genetic Toxicology Technical Committee, under the Health and Environmental Sciences Institute, conducted a case study on 31 reference chemicals to evaluate the performance of IVIVE application to genotoxicity data. The results demonstrate that for most chemicals considered here (20/31), the PODs derived from in vitro data and IVIVE are health protective relative to in vivo PODs from animal studies. PODs were also protective by assay target: mutations (8/13 chemicals), micronuclei (9/12), and aneugenicity markers (4/4). It is envisioned that this novel testing strategy could enhance prioritization, rapid screening, and risk assessment of genotoxic chemicals.

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.

Quantitative in vitro to in vivo extrapolation for developmental toxicity potency of valproic acid analogues

BACKGROUND

The developmental toxicity potential (dTP) concentration from the devTOX quickPredict (devTOX^qP ) assay, a metabolomics-based human induced pluripotent stem cell assay, predicts a chemical's developmental toxicity potency. Here, in vitro to in vivo extrapolation (IVIVE) approaches were applied to address whether the devTOX^qP assay could quantitatively predict in vivo developmental toxicity lowest effect levels (LELs) for the prototypical teratogen valproic acid (VPA) and a group of structural analogues.

METHODS

VPA and a series of structural analogues were tested with the devTOX^qP assay to determine dTP concentration and we estimated the equivalent administered doses (EADs) that would lead to plasma concentrations equivalent to the in vitro dTP concentrations. The EADs were compared to the LELs in rat developmental toxicity studies, human clinical doses, and EADs reported using other in vitro assays. To evaluate the impact of different pharmacokinetic (PK) models on IVIVE outcomes, we compared EADs predicted using various open-source and commercially available PK and physiologically based PK (PBPK) models. To evaluate the effect of in vitro kinetics, an equilibrium distribution model was applied to translate dTP concentrations to free medium concentrations before subsequent IVIVE analyses.

RESULTS

The EAD estimates for the VPA analogues based on different PK/PBPK models were quantitatively similar to in vivo data from both rats and humans, where available, and the derived rank order of the chemicals was consistent with observed in vivo developmental toxicity. Different models were identified that provided accurate predictions for rat prenatal LELs and conservative estimates of human safe exposure. The impact of in vitro kinetics on EAD estimates is chemical-dependent. EADs from this study were within range of predicted doses from other in vitro and model organism data.

CONCLUSIONS

This study highlights the importance of pharmacokinetic considerations when using in vitro assays and demonstrates the utility of the devTOX^qP human stem cell-based platform to quantitatively assess a chemical's developmental toxicity potency.

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.

Estimation of the bio-accessibility of methylmercury from aquatic foods using a PBTK model with an approximate Bayesian computation method in Chinese pregnant women

BACKGROUND

Fish are the primary source of protein and docosahexaenoic acid (DHA) for pregnant women and children, but methylmercury (MeHg) pollution is the potential hazard of fish consumption. In risk assessments, the bio-accessibility of MeHg is usually assumed to be 100%, which could lead to overestimation of dietary exposure.

METHOD

An existing PBTK model was adapted to estimate parameters of the bio-accessibility based on MeHg exposure data from a cohort of 397 Chinese pregnant women. The posterior distributions of parameters were determined by using the ABC - MCMC. RMSEP and Spearman's rank correlation coefficients (Rho) were calculated to determine the goodness of model fitting. The Monte Carlo analysis was performed for the parameter distributions to estimate the model variability.

RESULT

The median of daily MeHg intake and maternal MeHg levels were 0.018 μg/kg bw and 3.01 μg/kg in the early and middle terms of pregnancy. The estimated bio-accessibility of freshwater fish, marine fish and others were 46.1, 17.3 and 58.2%, separately. The RMSEP improved from 11.18 to 2.54 and the Rho improved from 0.19 to 0.22 after bio-accessibility optimization. The model variability was estimated to be 2.6.

CONCLUSION

The bio-accessibility estimated in this study was comparable to that determined in previous in vitro studies. The optimized model could improve the prediction performance on the MeHg body burden by dietary exposure.

The Evolution of Regulatory Toxicology: Where is the Gardener?

There is a need for paradigm change in the methodology employed for toxicological testing and assessment. It could be said that this change is well on its way, through an evolutionary progress analogous to that of natural selection. Darwin's Theory of Evolution has defined the idea of evolution and descendancy since the last third of the 19th century. Increasingly, this concept of 'evolution' is being applied beyond the field of biology. This Comment article discusses the progress of toxicological testing in the context of 'evolutionary pressure' and deliberates how this process can help foster the development, implementation and acceptance of mechanistic and human-relevant methods in this field. By comparing the current regulatory landscape in toxicity testing and assessment to specific elements in Charles Darwin's evolutionary theory, we aim to better understand the needs and requirements for the future.

Nanosafety: An Evolving Concept to Bring the Safest Possible Nanomaterials to Society and Environment

The use of nanomaterials has been increasing in recent times, and they are widely used in industries such as cosmetics, drugs, food, water treatment, and agriculture. The rapid development of new nanomaterials demands a set of approaches to evaluate the potential toxicity and risks related to them. In this regard, nanosafety has been using and adapting already existing methods (toxicological approach), but the unique characteristics of nanomaterials demand new approaches (nanotoxicology) to fully understand the potential toxicity, immunotoxicity, and (epi)genotoxicity. In addition, new technologies, such as organs-on-chips and sophisticated sensors, are under development and/or adaptation. All the information generated is used to develop new in silico approaches trying to predict the potential effects of newly developed materials. The overall evaluation of nanomaterials from their production to their final disposal chain is completed using the life cycle assessment (LCA), which is becoming an important element of nanosafety considering sustainability and environmental impact. In this review, we give an overview of all these elements of nanosafety.

Understanding the combined effects of multiple stressors: A new perspective on a longstanding challenge

Wildlife populations and their habitats are exposed to an expanding diversity and intensity of stressors caused by human activities, within the broader context of natural processes and increasing pressure from climate change. Estimating how these multiple stressors affect individuals, populations, and ecosystems is thus of growing importance. However, their combined effects often cannot be predicted reliably from the individual effects of each stressor, and we lack the mechanistic understanding and analytical tools to predict their joint outcomes. We review the science of multiple stressors and present a conceptual framework that captures and reconciles the variety of existing approaches for assessing combined effects. Specifically, we show that all approaches lie along a spectrum, reflecting increasing assumptions about the mechanisms that regulate the action of single stressors and their combined effects. An emphasis on mechanisms improves analytical precision and predictive power but could introduce bias if the underlying assumptions are incorrect. A purely empirical approach has less risk of bias but requires adequate data on the effects of the full range of anticipated combinations of stressor types and magnitudes. We illustrate how this spectrum can be formalised into specific analytical methods, using an example of North Atlantic right whales feeding on limited prey resources while simultaneously being affected by entanglement in fishing gear. In practice, case-specific management needs and data availability will guide the exploration of the stressor combinations of interest and the selection of a suitable trade-off between precision and bias. We argue that the primary goal for adaptive management should be to identify the most practical and effective ways to remove or reduce specific combinations of stressors, bringing the risk of adverse impacts on populations and ecosystems below acceptable thresholds.

A modular mechanistic framework for estimating exposure to SVOCs: Next steps for modeling emission and partitioning of plasticizers and PFAS

Estimates of human exposure to semi-volatile organic compounds (SVOCs) such as phthalates, phthalate alternatives, and some per- and polyfluoroalkyl substances (PFAS) are required for the risk-based evaluation of chemicals. Recently, a modular mechanistic modeling framework to rapidly predict SVOC emission and partitioning in indoor environments has been presented, in which several mechanistically consistent source emission categories (SECs) were identified. However, not all SECs have well-developed emission models. In addition, data on model parameters are missing even for frequently studied SVOCs. These knowledge gaps impede the comprehensive prediction of the fate of SVOCs indoors. In this paper, sets of high-priority phthalates, phthalate alternatives, and PFAS were identified based on chemical occurrence indoors and additional selection criteria. These high-priority chemicals served as the basis for exploring model parameter availability for existing indoor SVOC emission and partitioning models. The results reveal that additional experimental and modeling work is needed to fully understand the behavior of SVOCs indoors and to predict exposures with greater confidence and lower uncertainty. Modeling approaches to fill some of the identified gaps are proposed. The prioritized sets of chemicals and proposed new modeling approaches will help guide future research. The inclusion of polar phases in the framework will further expand its applicability and scope. IMPACT STATEMENT: This paper compiles data on high-priority chemicals commonly found indoors and information on the availability of applicable models and model parameters to predict emission, partitioning, and subsequent exposure to these chemicals. Modeling approaches for a selection of the missing SECs (source emission categories) are proposed, to illustrate the path forward. The comprehensive data set helps inform researchers, exposure assessors, and policy makers to better understand the state of the science regarding modeling of indoor exposure to semi-volatile organic compounds (SVOCs) and per- and polyfluoroalkyl substances (PFAS).

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.

Comprehensive profile of DNA adducts as both tissue and urinary biomarkers of exposure to acrylamide and chemo-preventive effect of catechins in rats

Two representative DNA adducts from acrylamide exposure, N7-(2-carbamoyl-2-hydroxyethyl) guanine (N7-GA-Gua) and N3-(2-carbamoyl-2-hydroxyethyl) adenine (N3-GA-Ade), are important long-term exposure biomarkers for evaluating genotoxicity of acrylamide. Catechins as natural antioxidants present in tea possess multiple health benefits, and may also have the potential to protect against acrylamide-induced DNA damage. The current study developed an ultra-high performance liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) method for simultaneous analysis of N7-GA-Gua and N3-GA-Ade in tissues and urine. The validated UHPLC-MS/MS method showed high sensitivity, with limit of detection and limit of quantification ranging 0.2-0.8 and 0.5-1.5 ng/mL, respectively, and achieved qualified precision (RSD<14.0%) and spiking recovery (87.2%-110.0%) with elution within 6 min, which was suitable for the analysis of the two DNA adducts in different matrices. The levels of N7-GA-Gua and N3-GA-Ade ranged 0.9-11.9 and 0.6-3.5 μg/g creatinine in human urine samples, respectively. To investigate the interventional effects of catechins on the two DNA adducts from acrylamide exposure, rats were supplemented with three types of catechins (tea polyphenols, epigallocatechin gallate, and epicatechin) 30 min before administration with acrylamide. Our results showed that catechins effectively inhibited the formation of DNA adducts from acrylamide exposure in both urine and tissues of rats. Among three catechins, epicatechin performed the best inhibitory effect. The current study provided evidence for the chemo-preventive effect of catechins, indicating that dietary supplement of catechins may contribute to health protection against exposure to acrylamide.

Integrating toxicokinetics into toxicology studies and the human health risk assessment process for chemicals: Reduced uncertainty, better health protection

The database of practical examples where toxicokinetic (TK) data has benefitted all stages of the human health risk assessment process are increasingly being published and accepted. This review aimed to highlight and summarise notable examples and to describe the "state of the art" in this field. The overall recommendation is that for any in vivo animal study conducted, measurements of TK should be very carefully considered for inclusion as the numerous benefits this brings continues to grow, particularly during the current march towards animal free toxicology testing and ambitions to eventually conduct human health risk assessments entirely based upon non-animal methods.

Phthalates and Alternative Plasticizers Differentially affect Phenotypic Parameters in Gonadal Somatic and Germ Cell Lines

The developmental and reproductive toxicity associated with exposure to phthalates has motivated a search for alternatives. However, there is limited knowledge regarding the adverse effects of some of these chemicals. We used high-content imaging to compare the effects of mono (2-ethylhexyl) phthalate (MEHP) with six alternative plasticizers: di-2-ethylhexyl terephthalate (DEHTP); diisononyl-phthalate (DINP); di-isononylcyclohexane-1,2-dicarboxylate (DINCH); 2-ethylhexyl adipate (DEHA); 2,2,4-trimethyl 1,3-pentanediol diisobutyrate (TXIB) and di-iso-decyl-adipate (DIDA). A male germ spermatogonial cell line (C18-4), a Sertoli cell line (TM4) and two steroidogenic cell lines (MA-10 Leydig and KGN granulosa) were exposed for 48h to each chemical (0.001-100 μM). Cell images were analyzed to assess cytotoxicity and effects on phenotypic endpoints. Only MEHP (100 μM) was cytotoxic and only in C18-4 cells. However, several plasticizers had distinct phenotypic effects in all four cell lines. DINP increased Calcein intensity in C18-4 cells, whereas DIDA induced oxidative stress. In TM4 cells, MEHP, and DINCH affected lipid droplet numbers, while DEHTP and DINCH increased oxidative stress. In MA-10 cells, MEHP increased lipid droplet areas and oxidative stress; DINP decreased the number of lysosomes, while DINP, DEHA and DIDA altered mitochondrial activity. In KGN cells, MEHP, DINP and DINCH increased the number of lipid droplets, whereas DINP decreased the number of lysosomes, increased oxidative stress and affected mitochondria. The Toxicological Priority Index (ToxPi) provided a visual illustration of the cell line specificity of the effects on phenotypic parameters. The lowest administered equivalent doses were observed for MEHP. We propose that this approach may assist in screening alternative plasticizers.

An integrative approach to define chemical exposure threshold limits for endangered sea turtles

Environmental contaminants pose serious health threats to marine megafauna species, yet methods defining exposure threshold limits are lacking. Here, a three-pillar chemical risk assessment framework is presented based on (1) species- and chemical-specific lifetime bioaccumulation modelling, (2) non-destructive in vitro and in vivo toxicity threshold assessment, and (3) chemical risk quantification. We used the effects of cadmium (Cd) in green sea turtles (Chelonia mydas) as a proof of concept to evaluate the quantitative mechanistic modelling approach. A physiologically-based kinetic (PBK) model simulated Cd tissue concentrations (liver, kidney, muscle, fat, brain, scute, and 'rest of the body') in C.mydas. The validated PBK model then translated species-specific in vitro results to in vivo effects. The results showed that the resilience of C.mydas towards Cd kidney toxicity is age-dependent and differs with changing physiology and feeding ecology. Using the model in reverse mode, a steady-state exposure threshold of 0.1 µg/g dry weight Cd in forage was derived and compared to real-world exposure scenarios. Three out of the four globally distinct C.mydas populations assessed are exposed to Cd levels above this threshold limit. This approach can be adapted to other marine species and chemicals to prioritize measures for managing potentially harmful chemical exposures.

Calculating qualified non-mutagenic impurity levels: Harmonization of approaches

The presence of impurities in drugs is unavoidable. As impurities offer no direct benefit to the patient, it is critical that impurities do not compromise patient safety. Current guidelines on the derivation of acceptable impurity levels leave aspects of calculations open for interpretation, resulting in inconsistencies across industry and regulators. To understand current impurity qualification practices from a safety standpoint, regulatory expectations and the safety risk that impurities pose, the IQ DruSafe Impurities Working Group (WG) conducted a pharmaceutical industry-wide survey. Survey results highlighted areas that could benefit from harmonization, including nonclinical species/sex selection and the application of adjustment factors (i.e., body surface area). Recommendations for alignment on these topics is included in this publication. Additionally, the WG collated repeat-dose toxicity information for 181 starting materials and intermediates, reflective of pharmaceutical impurities, to understand the toxicological risks they generally pose in relation to the drug substance (DS) and the assumptions surrounding the calculation of qualified impurity levels. An evaluation of this dataset and the survey were used to harmonize how to calculate a safe limit for an impurity based on toxicology testing of the impurity when present within the DS.

Integrating endocrine-related health effects into comparative human toxicity characterization

Endocrine-disrupting chemicals have the ability to interfere with and alter functions of the hormone system, leading to adverse effects on reproduction, growth and development. Despite growing concerns over their now ubiquitous presence in the environment, endocrine-related human health effects remain largely outside of comparative human toxicity characterization frameworks as applied for example in life cycle impact assessments. In this paper, we propose a new methodological framework to consistently integrate endocrine-related health effects into comparative human toxicity characterization. We present two quantitative and operational approaches for extrapolating towards a common point of departure from both in vivo and dosimetry-adjusted in vitro endocrine-related effect data and deriving effect factors as well as corresponding characterization factors for endocrine-active/endocrine-disrupting chemicals. Following the proposed approaches, we calculated effect factors for 323 chemicals, reflecting their endocrine potency, and related characterization factors for 157 chemicals, expressing their relative endocrine-related human toxicity potential. Developed effect and characterization factors are ready for use in the context of chemical prioritization and substitution as well as life cycle impact assessment and other comparative assessment frameworks. Endocrine-related effect factors were found comparable to existing effect factors for cancer and non-cancer effects, indicating that (1) the chemicals' endocrine potency is not necessarily higher or lower than other effect potencies and (2) using dosimetry-adjusted effect data to derive effect factors does not consistently overestimate the effect of potential endocrine disruptors. Calculated characterization factors span over 8-11 orders of magnitude for different substances and emission compartments and are dominated by the range in endocrine potencies.

Assessing Human Exposure to SVOCs in Materials, Products, and Articles: A Modular Mechanistic Framework

A critical review of the current state of knowledge of chemical emissions from indoor sources, partitioning among indoor compartments, and the ensuing indoor exposure leads to a proposal for a modular mechanistic framework for predicting human exposure to semivolatile organic compounds (SVOCs). Mechanistically consistent source emission categories include solid, soft, frequent contact, applied, sprayed, and high temperature sources. Environmental compartments are the gas phase, airborne particles, settled dust, indoor surfaces, and clothing. Identified research needs are the development of dynamic emission models for several of the source emission categories and of estimation strategies for critical model parameters. The modular structure of the framework facilitates subsequent inclusion of new knowledge, other chemical classes of indoor pollutants, and additional mechanistic processes relevant to human exposure indoors. The framework may serve as the foundation for developing an open-source community model to better support collaborative research and improve access for application by stakeholders. Combining exposure estimates derived using this framework with toxicity data for different end points and toxicokinetic mechanisms will accelerate chemical risk prioritization, advance effective chemical management decisions, and protect public health.

Development and evaluation of a high throughput inhalation model for organic chemicals

Currently it is difficult to prospectively estimate human toxicokinetics (particularly for novel chemicals) in a high-throughput manner. The R software package httk has been developed, in part, to address this deficiency, and the aim of this investigation was to develop a generalized inhalation model for httk. The structure of the inhalation model was developed from two previously published physiologically based models from Jongeneelen and Berge (Ann Occup Hyg 55:841-864, 2011) and Clewell et al. (Toxicol Sci 63:160-172, 2001), while calculated physicochemical data was obtained from EPA's CompTox Chemicals Dashboard. In total, 142 exposure scenarios across 41 volatile organic chemicals were modeled and compared to published data. The slope of the regression line of best fit between log-transformed simulated and observed blood and exhaled breath concentrations was 0.46 with an r2 = 0.45 and a root mean square error (RMSE) of direct comparison between the log-transformed simulated and observed values of 1.11. Approximately 5.1% (n = 108) of the data points analyzed were >2 orders of magnitude different than expected. The volatile organic chemicals examined in this investigation represent small, generally lipophilic molecules. Ultimately this paper details a generalized inhalation component that integrates with the httk physiologically based toxicokinetic model to provide high-throughput estimates of inhalation chemical exposures.

Impact of stressors in the aviation environment on xenobiotic dosimetry in humans: physiologically based prediction of the effect of barometric pressure or altitude

Standard health risks from volatile organic compounds (VOCs) are generally interpreted at ambient environmental conditions. The aim of this study was to develop a strategy for using physiologically based pharmacokinetic (PBPK) modeling to compare known risks in the general population to calculated risks in pilots experiencing pressure-based stressors. PBPK models facilitate these comparisons by prediction of how target-tissue specific doses are altered when a stressor, such as high altitude, produces changes in physiological parameters. Cardiac output, regional blood flow, and alveolar ventilation rate following acute exposure to altitude ranging from moderate to extremely high were estimated from published data from 52 groups of human subjects. Scenarios where pilots might inhale toluene, 1,2,4-trimethylbenzene (1,2,4-TMB), or cyclohexane during routine military flight training were simulated. At the recommended Threshold Limit Values (TLV), arterial blood concentrations were predicted to be higher for exposure at 15000 ft (4572 m) than at sea level. The differences were greater for toluene and TMB, which have higher blood: air and fat: blood partition coefficients than less lipophilic cyclohexane. In summary, quantitative approaches to internal dosimetry prediction that take advantage of existing knowledge of physiological changes induced by occupational stressors possess potential as tools in performing a human health risk assessment.

Design and Rationale of the Biomarker Center of the Household Air Pollution Intervention Network (HAPIN) Trial

BACKGROUND

Biomarkers of exposure, susceptibility, and effect are fundamental for understanding environmental exposures, mechanistic pathways of effect, and monitoring early adverse outcomes. To date, no study has comprehensively evaluated a large suite and variety of biomarkers in household air pollution (HAP) studies in concert with exposure and outcome data. The Household Air Pollution Intervention Network (HAPIN) trial is a liquified petroleum gas (LPG) fuel/stove randomized intervention trial enrolling 800 pregnant women in each of four countries (i.e., Peru, Guatemala, Rwanda, and India). Their offspring will be followed from birth through 12 months of age to evaluate the role of pre- and postnatal exposure to HAP from biomass burning cookstoves in the control arm and LPG stoves in the intervention arm on growth and respiratory outcomes. In addition, up to 200 older adult women per site are being recruited in the same households to evaluate indicators of cardiopulmonary, metabolic, and cancer outcomes.

OBJECTIVES

Here we describe the rationale and ultimate design of a comprehensive biomarker plan to enable us to explore more fully how exposure is related to disease outcome.

METHODS

HAPIN enrollment and data collection began in May 2018 and will continue through August 2021. As a part of data collection, dried blood spot (DBS) and urine samples are being collected three times during pregnancy in pregnant women and older adult women. DBS are collected at birth for the child. DBS and urine samples are being collected from the older adult women and children three times throughout the child's first year of life. Exposure biomarkers that will be longitudinally measured in all participants include urinary hydroxy-polycyclic aromatic hydrocarbons, volatile organic chemical metabolites, metals/metalloids, levoglucosan, and cotinine. Biomarkers of effect, including inflammation, endothelial and oxidative stress biomarkers, lung cancer markers, and other clinically relevant measures will be analyzed in urine, DBS, or blood products from the older adult women. Similarly, genomic/epigenetic markers, microbiome, and metabolomics will be measured in older adult women samples.

DISCUSSION

Our study design will yield a wealth of biomarker data to evaluate, in great detail, the link between exposures and health outcomes. In addition, our design is comprehensive and innovative by including cutting-edge measures such as metabolomics and epigenetics. https://doi.org/10.1289/EHP5751.

Risk Assessment of Cosmetics Using Triclosan on Future Generation's Germ Cell Maturation via Lactating Mother Rats

Triclosan (TCS) is a widely used chemical in personal care and household products as an antimicrobial agent but some studies have reported it as being estrogenic. We investigated the influence of TCS on the male reproductive system of postnatal pups. Lactating mother rats (Rattus norvegicus) were given daily doses of 0 mg, 3 mg, and 5 mg/kg/day from the day of delivery until 28 days, equivalent to their natural breastfeeding duration. At 28 days, the male pups of all three groups were sacrificed and their biochemical parameters evaluated. TCS-treated pups had decreased mRNA levels for 3β hydro-hydroxysteroid dehydrogenases (3βHSD), OCT3/4, and androgen receptor (AR) (p < 0.05). The higher dose (5 mg/kg/day) male pups exhibited more significantly affected germ cell maturation and decreased body weight. In summary, TCS-treated lactating mothers passed the deleterious effects to their untreated male pups as exhibited by reduced androgens synthesis and subsequently decreased sperm count.

Integration of Food Animal Residue Avoidance Databank (FARAD) empirical methods for drug withdrawal interval determination with a mechanistic population-based interactive physiologically based pharmacokinetic (iPBPK) modeling platform: example for flunixin meglumine administration

Violative chemical residues in animal-derived food products affect food safety globally and have impact on the trade of international agricultural products. The Food Animal Residue Avoidance Databank program has been developing scientific tools to provide appropriate withdrawal interval (WDI) estimations after extralabel drug use in food animals for the past three decades. One of the tools is physiologically based pharmacokinetic (PBPK) modeling, which is a mechanistic-based approach that can be used to predict tissue residues and WDIs. However, PBPK models are complicated and difficult to use by non-modelers. Therefore, a user-friendly PBPK modeling framework is needed to move this field forward. Flunixin was one of the top five violative drug residues identified in the United States from 2010 to 2016. The objective of this study was to establish a web-based user-friendly framework for the development of new PBPK models for drugs administered to food animals. Specifically, a new PBPK model for both cattle and swine after administration of flunixin meglumine was developed. Population analysis using Monte Carlo simulations was incorporated into the model to predict WDIs following extralabel administration of flunixin meglumine. The population PBPK model was converted to a web-based interactive PBPK (iPBPK) framework to facilitate its application. This iPBPK framework serves as a proof-of-concept for further improvements in the future and it can be applied to develop new models for other drugs in other food animal species, thereby facilitating the application of PBPK modeling in WDI estimation and food safety assessment.

Evaluation of toxicokinetics of nonylphenol in the adult female Sprague-Dawley rats using a physiologically based toxicokinetic model

The physiologically based toxicokinetic (PBTK) model was firstly developed in female rats to quantitatively evaluate toxicokinetics of nonylphenol (NP). Changes in NP serum concentrations over time of single oral NP administration experiments in Sprague-Dawley rats and literature data were collected to establish and calibrate the PBTK model in the SimBiology framework. The calibrated model predicted the serum and tissue NP concentrations of repeat oral NP administration for model evaluation. NP concentrations in serum and tissues (liver, kidneys, adipose, brain, uterus and ovaries) were quantified using ultra-high performance liquid chromatography-tandem mass spectrometry (UPLC-MS/MS). The model output of the time course data (values are within the standard deviation defined for each data point) indicated proximity of predictions to reality. The coefficients of determination r² were all greater than 0.9, and the root mean squared error (RMSE) were within 0.177-2.027, which means the model predicted and observed serum NP concentrations were in excellent agreement. The results indicated that the model could contribute to a simplification of the future exposure risk assessments of NP in a more realistic scenario and provided a better understanding of the disposition process of NP in human.