Role of Physiologically Based Kinetic modelling in addressing environmental chemical mixtures - A review

Evaluation of a Proportional Response Addition Approach to Mixture Risk Assessment and Predictive Toxicology Using Data on Four Trihalomethanes from the U.S. EPA's Multiple-Purpose Design Study

In this study, proportional response addition (Prop-RA), a model for predicting response from chemical mixture exposure, is demonstrated and evaluated by statistically analyzing data on all possible binary combinations of the four regulated trihalomethanes (THMs). These THMs were the subject of a multipurpose toxicology study specifically designed to evaluate Prop-RA. The experimental design used a set of doses common to all components and mixtures, providing hepatotoxicity data on the four single THMs and the binary combinations. In Prop-RA, the contribution of each component to mixture toxicity is proportional to its fraction in the mixture based on its response at the total mixture dose. The primary analysis consisted of 160 evaluations. Statistically significant departures from the Prop-RA prediction were found for seven evaluations, with three predications that were greater than and four that were less than the predicted response; interaction magnitudes (n-fold difference in response vs. prediction) ranged from 1.3 to 1.4 for the former and 2.6 to 3.8 for the latter. These predictions support the idea that Prop-RA works best with chemicals where the effective dose ranges overlap. Prop-RA does not assume the similarity of toxic action or independence, but it can be applied to a mixture of components that affect the same organ/system, with perhaps unknown toxic modes of action.

Studying mixture effects on uptake and tissue distribution of PFAS in zebrafish (Danio rerio) using physiologically based kinetic (PBK) modelling

Per- and polyfluoroalkyl substances (PFAS) are ubiquitously distributed in the aquatic environment. They include persistent, mobile, bioaccumulative, and toxic chemicals and it is therefore critical to increase our understanding on their adsorption, distribution, metabolism, excretion (ADME). The current study focused on uptake of seven emerging PFAS in zebrafish (Danio rerio) and their potential maternal transfer. In addition, we aimed at increasing our understanding on mixture effects on ADME by developing a physiologically based kinetic (PBK) model capable of handling co-exposure scenarios of any number of chemicals. All studied chemicals were taken up in the fish to varying degrees, whereas only perfluorononanoate (PFNA) and perfluorooctanoate (PFOA) were quantified in all analysed tissues. Perfluorooctane sulfonamide (FOSA) was measured at concerningly high concentrations in the brain (Cmax over 15 μg/g) but also in the liver and ovaries. All studied PFAS were maternally transferred to the eggs, with FOSA and 6:2 perfluorooctane sulfonate (6,2 FTSA) showing significant (p < 0.02) signs of elimination from the embryos during the first 6 days of development, while perfluorobutane sulfonate (PFBS), PFNA, and perfluorohexane sulfonate (PFHxS) were not eliminated in embryos during this time-frame. The mixture PBK model resulted in >85 % of predictions within a 10-fold error and 60 % of predictions within a 3-fold error. At studied levels of PFAS exposure, competitive binding was not a critical factor for PFAS kinetics. Gill surface pH influenced uptake for some carboxylates but not the sulfonates. The developed PBK model provides an important tool in understanding kinetics under complex mixture scenarios and this use of New Approach Methodologies (NAMs) is critical in future risk assessment of chemicals and early warning systems.

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

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

The Exposome and Nutritional Pharmacology and Toxicology: A New Application for Metabolomics

The exposome refers to all of the internal and external life-long exposures that an individual experiences. These exposures, either acute or chronic, are associated with changes in metabolism that will positively or negatively influence the health and well-being of individuals. Nutrients and other dietary compounds modulate similar biochemical processes and have the potential in some cases to counteract the negative effects of exposures or enhance their beneficial effects. We present herein the concept of Nutritional Pharmacology/Toxicology which uses high-information metabolomics workflows to identify metabolic targets associated with exposures. Using this information, nutritional interventions can be designed toward those targets to mitigate adverse effects or enhance positive effects. We also discuss the potential for this approach in precision nutrition where nutrients/diet can be used to target gene-environment interactions and other subpopulation characteristics. Deriving these "nutrient cocktails" presents an opportunity to modify the effects of exposures for more beneficial outcomes in public health.

Characterisation of the risk associated with chronic lifetime exposure to mixture of chemical hazards: case study of trace elements

Risk assessment methodology, mostly commonly used, faces the complexity of the environment. Populations are exposed to multiple sources of chemicals throughout life and the chemical mixtures they are exposed change during time (lifestyle factors, regulatory decisions, etc). The risk assessment needs to consider these dynamics and the evolution of the body with age, in order to refine the exposure assessment to chemicals and to predict the health impact of these exposures. This review highlights the latest methodologies developed to improve risk assessment, especially cor heavy metals. The methodologies aim to better describe the chemical toxicokinetic and toxicodynamic as well as the exposure assessment. Human Biomonitoring (HBM) data give great opportunities to link biomarkers of exposure with an adverse effect. Physiologically-Based Toxicokinetic (PBTK) models are more and more used to simulate the evolution of biomarkers in organisms, considering the external exposures and the physiological evolutions. PBTK models may also be used to determine the routes of exposure or to predict the impacts of schemes of exposure. The major limit is the integration of several chemicals in mixture with common adverse effects and the interactions between them.

Partitioning of Persistent Organic Pollutants between Adipose Tissue and Serum in Human Studies

Blood is the most widely used matrix for biomonitoring of persistent organic pollutants (POPs). It is assumed that POPs are homogenously distributed within body lipids at steady state; however, the variability underlying the partitioning of POPs between fat compartments is poorly understood. Hence, the objective of this study was to review the state of the science about the relationships of POPs between adipose tissue and serum in humans. We conducted a narrative literature review of human observational studies reporting concentrations of POPs in paired samples of adipose tissue with other lipid-based compartments (e.g., serum lipids). The searches were conducted in SCOPUS and PUBMED. A meta-regression was performed to identify factors responsible for variability. All included studies reported high variability in the partition coefficients of POPs, mainly between adipose tissue and serum. The number of halogen atoms was the physicochemical variable most strongly and positively associated with the partition ratios, whereas body mass index was the main biological factor positively and significantly associated. To conclude, although this study provides a better understanding of partitioning of POPs to refine physiologically based pharmacokinetic and epidemiological models, further research is still needed to determine other key factors involved in the partitioning of POPs.

Evaluation of Inhalation Exposures and Potential Health Impacts of Ingredient Mixtures Using in vitro to in vivo Extrapolation

In vitro methods offer opportunities to provide mechanistic insight into bioactivity as well as human-relevant toxicological assessments compared to animal testing. One of the challenges for this task is putting in vitro bioactivity data in an in vivo exposure context, for which in vitro to in vivo extrapolation (IVIVE) translates in vitro bioactivity to clinically relevant exposure metrics using reverse dosimetry. This study applies an IVIVE approach to the toxicity assessment of ingredients and their mixtures in e-cigarette (EC) aerosols as a case study. Reported in vitro cytotoxicity data of EC aerosols, as well as in vitro high-throughput screening (HTS) data for individual ingredients in EC liquids (e-liquids) are used. Open-source physiologically based pharmacokinetic (PBPK) models are used to calculate the plasma concentrations of individual ingredients, followed by reverse dosimetry to estimate the human equivalent administered doses (EADs) needed to obtain these plasma concentrations for the total e-liquids. Three approaches (single actor approach, additive effect approach, and outcome-oriented ingredient integration approach) are used to predict EADs of e-liquids considering differential contributions to the bioactivity from the ingredients (humectant carriers [propylene glycol and glycerol], flavors, benzoic acid, and nicotine). The results identified critical factors for the EAD estimation, including the ingredients of the mixture considered to be bioactive, in vitro assay selection, and the data integration approach for mixtures. Further, we introduced the outcome-oriented ingredient integration approach to consider e-liquid ingredients that may lead to a common toxicity outcome (e.g., cytotoxicity), facilitating a quantitative evaluation of in vitro toxicity data in support of human risk assessment.

Towards a systematic use of effect biomarkers in population and occupational biomonitoring

Effect biomarkers can be used to elucidate relationships between exposure to environmental chemicals and their mixtures with associated health outcomes, but they are often underused, as underlying biological mechanisms are not understood. We aim to provide an overview of available effect biomarkers for monitoring chemical exposures in the general and occupational populations, and highlight their potential in monitoring humans exposed to chemical mixtures. We also discuss the role of the adverse outcome pathway (AOP) framework and physiologically based kinetic and dynamic (PBK/D) modelling to strengthen the understanding of the biological mechanism of effect biomarkers, and in particular for use in regulatory risk assessments. An interdisciplinary network of experts from the European chapter of the International Society for Exposure Science (ISES Europe) and the Organization for Economic Co-operation and Development (OECD) Occupational Biomonitoring activity of Working Parties of Hazard and Exposure Assessment group worked together to map the conventional framework of biomarkers and provided recommendations for their systematic use. We summarized the key aspects of this work here, and discussed these in three parts. Part I, we inventory available effect biomarkers and promising new biomarkers for the general population based on the H2020 Human Biomonitoring for Europe (HBM4EU) initiative. Part II, we provide an overview AOP and PBK/D modelling use that improved the selection and interpretation of effect biomarkers. Part III, we describe the collected expertise from the OECD Occupational Biomonitoring subtask effect biomarkers in prioritizing relevant mode of actions (MoAs) and suitable effect biomarkers. Furthermore, we propose a tiered risk assessment approach for occupational biomonitoring. Several effect biomarkers, especially for use in occupational settings, are validated. They offer a direct assessment of the overall health risks associated with exposure to chemicals, chemical mixtures and their transformation products. Promising novel effect biomarkers are emerging for biomonitoring of the general population. Efforts are being dedicated to prioritizing molecular and biochemical effect biomarkers that can provide a causal link in exposure-health outcome associations. This mechanistic approach has great potential in improving human health risk assessment. New techniques such as in silico methods (e.g. QSAR, PBK/D modelling) as well as 'omics data will aid this process. Our multidisciplinary review represents a starting point for enhancing the identification of effect biomarkers and their mechanistic pathways following the AOP framework. This may help in prioritizing the effect biomarker implementation as well as defining threshold limits for chemical mixtures in a more structured way. Several ex vivo biomarkers have been proposed to evaluate combined effects including genotoxicity and xeno-estrogenicity. There is a regulatory need to derive effect-based trigger values using the increasing mechanistic knowledge coming from the AOP framework to address adverse health effects due to exposure to chemical mixtures. Such a mechanistic strategy would reduce the fragmentation observed in different regulations. It could also stimulate a harmonized use of effect biomarkers in a more comparable way, in particular for risk assessments to chemical mixtures.

Physiologically based kinetic (PBK) modelling and human biomonitoring data for mixture risk assessment

Human biomonitoring (HBM) data can provide insight into co-exposure patterns resulting from exposure to multiple chemicals from various sources and over time. Therefore, such data are particularly valuable for assessing potential risks from combined exposure to multiple chemicals. One way to interpret HBM data is establishing safe levels in blood or urine, called Biomonitoring Equivalents (BE) or HBM health based guidance values (HBM-HBGV). These can be derived by converting established external reference values, such as tolerable daily intake (TDI) values. HBM-HBGV or BE values are so far agreed only for a very limited number of chemicals. These values can be established using physiologically based kinetic (PBK) modelling, usually requiring substance specific models and the collection of many input parameters which are often not available or difficult to find in the literature. The aim of this study was to investigate the suitability and limitations of generic PBK models in deriving BE values for several compounds with a view to facilitating the use of HBM data in the assessment of chemical mixtures at a screening level. The focus was on testing the methodology with two generic models, the IndusChemFate tool and High-Throughput Toxicokinetics package, for two different classes of compounds, phenols and phthalates. HBM data on Danish children and on Norwegian mothers and children were used to evaluate the quality of the predictions and to illustrate, by means of a case study, the overall approach of applying PBK models to chemical classes with HBM data in the context of chemical mixture risk assessment. Application of PBK models provides a better understanding and interpretation of HBM data. However, the study shows that establishing safety threshold levels in urine is a difficult and complex task. The approach might be more straightforward for more persistent chemicals that are analysed as parent compounds in blood but high uncertainties have to be considered around simulated metabolite concentrations in urine. Refining the models may reduce these uncertainties and improve predictions. Based on the experience gained with this study, the performance of the models for other chemicals could be investigated, to improve the accuracy of the simulations.