Characterizing the risk related to the exposure to methylmercury over a lifetime: A global approach using population internal exposure

Seafood products accumulate methylmercury throughout the food chain and are the main source of methylmercury exposure. Methylmercury may trigger a number of adverse health effects, such as neurodevelopmental or nephrotoxic effects, the risk of which cannot be ruled out for the French high consumers of seafood. The characterisation of methylmercury-related risks is generally based on short-term dietary exposure without considering changes in consumption and exposure over the lifetime. Additionally, focusing on short-term dietary exposure, the fate of methylmercury (especially its accumulation) in the organism is not considered. The present study proposes a methodology basing risk characterization on estimates of body burden over a lifetime. First, trajectories of dietary exposures throughout lifetime were constructed based on the actual concentrations of total diet studies for a fictive representative French population, taking into account the social, economic and demographic parameters of individuals. Next, the fate of methylmercury in the body was estimated, based on these trajectories, using a specific physiologically-based kinetic (PBK) model that generated a representative pool of body burden trajectories. Simulated hair mercury concentrations were closed to previously reported French representative human biomonitoring data. Results showed that at certain stages of life, concentrations of methylmercury in hair were higher than the human biomonitoring guidance value set at 2.5 μg/g of hair by JECFA. This study showed the added value, in the case of substances accumulating in the body, of estimating dietary exposure over a lifetime and using exposure biomarkers estimated by a PBK model characterize the risk.

Assessing chronic gestational exposure to environmental chemicals in pregnant women: Advancing the co-PBK model

Vulnerable populations, such as pregnant women and their fetuses, confront potential health risks due to exposure to environmental toxic compounds. Computational methods have been popular in assessing chemical exposure to populations, contrasting with traditional cohort studies for human biomonitoring. This study proposes a screening-level approach based on physiologically based kinetic (PBK) modeling to evaluate the steady-state exposure of pregnant women to environmental chemicals throughout pregnancy. To exemplify the modeling application, naphthalene was chosen. Simulation results indicated that maternal fat exhibited significant bioaccumulation potential, with the log-transformed BTF of naphthalene at 0.51 mg kg-1 per mg d-1 in the steady state. The placenta was primarily exposed to 0.83 mg/d naphthalene for a 75.2 kg pregnant woman, considering all exposure routes. In the fetal structure, single-organ fetal PBK modeling estimated a naphthalene exposure of 123.64 mg/d to the entire fetus, while multiple-organ fetal PBK modeling further revealed the bioaccumulation highest in fat tissue. The liver identified as the vital organ for metabolism, kBioT,LiverM was demonstrated with the highest sensitivity among rate constants in the maternal body. Furthermore, the first-order kinetic rate constants related to the placenta and blood were found to impact the distribution process of naphthalene in the fetus, influencing gestational exposure. In conclusion, urgent attention is needed to develop a computational biomonitoring tool for assessing toxic chemical exposure in vulnerable populations.

A systems toxicology approach for identification of disruptions in cholesterol homeostasis after aggregated exposure to mixtures of perfluorinated compounds in humans

Per- and polyfluoroalkyl substances (PFAS) are used in various household and industrial products. In humans, positive associations were reported between PFAS, including perfluorsulfonic acid and perfluorooctanoic acid, and cholesterol, a cardiometabolic risk factor. Animal studies show the opposite. Human-centered approaches are needed to better understand the effects of PFAS mixtures on cholesterol. Here, a systems toxicology approach is described, using a gene-centered cholesterol biokinetic model. PFAS exposure-gene expression relations from published data were introduced into the model. An existing PFAS physiologically based kinetic model was augmented with lung and dermal compartments and integrated with the cholesterol model to enable exposure-effect modeling. The final model was populated with data reflecting lifetime mixture exposure from: tolerable weekly intake values; the environment; high occupational exposures (ski waxing, PFAS industry). Results indicate that low level exposures (tolerable weekly intake, environmental) did not change cholesterol. In contrast, occupational exposures clearly resulted in internal PFAS exposure and disruption of cholesterol homeostasis, largely in line with epidemiological observations. Despite model limitations (eg, dynamic range, directionality), changes in cholesterol homeostasis were predicted for ski waxers, hitherto unknown from epidemiological studies. Here, future studies involving lipid metabolism could improve risk assessment.

A minimal human physiologically based kinetic model of thyroid hormones and chemical disruption of plasma thyroid hormone binding proteins

The thyroid hormones (THs), thyroxine (T4) and triiodothyronine (T3), are under homeostatic control by the hypothalamic-pituitary-thyroid axis and plasma TH binding proteins (THBPs), including thyroxine-binding globulin (TBG), transthyretin (TTR), and albumin (ALB). THBPs buffer free THs against transient perturbations and distribute THs to tissues. TH binding to THBPs can be perturbed by structurally similar endocrine-disrupting chemicals (EDCs), yet their impact on circulating THs and health risks are unclear. In the present study, we constructed a human physiologically based kinetic (PBK) model of THs and explored the potential effects of THBP-binding EDCs. The model describes the production, distribution, and metabolism of T4 and T3 in the Body Blood, Thyroid, Liver, and Rest-of-Body (RB) compartments, with explicit consideration of the reversible binding between plasma THs and THBPs. Rigorously parameterized based on literature data, the model recapitulates key quantitative TH kinetic characteristics, including free, THBP-bound, and total T4 and T3 concentrations, TH productions, distributions, metabolisms, clearance, and half-lives. Moreover, the model produces several novel findings. (1) The blood-tissue TH exchanges are fast and nearly at equilibrium especially for T4, providing intrinsic robustness against local metabolic perturbations. (2) Tissue influx is limiting for transient tissue uptake of THs when THBPs are present. (3) Continuous exposure to THBP-binding EDCs does not alter the steady-state levels of THs, while intermittent daily exposure to rapidly metabolized TBG-binding EDCs can cause much greater disruptions to plasma and tissue THs. In summary, the PBK model provides novel insights into TH kinetics and the homeostatic roles of THBPs against thyroid disrupting chemicals.

Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict In Vivo Nrf2 Activation by EGCG and Its Colonic Metabolites in Humans

(-)-Epigallocatechin gallate (EGCG) is prone to microbial metabolism when reaching the colon. This study aimed to develop a human physiologically based kinetic (PBK) model for EGCG, with sub-models for its colonic metabolites gallic acid and pyrogallol. Results show that the developed PBK model could adequately predict in vivo time-dependent blood concentrations of EGCG after either the single or repeated oral administration of EGCG under both fasting and non-fasting conditions. The predicted in vivo blood C_max of EGCG indicates that the Nrf2 activation is limited, while concentrations of its metabolites in the intestinal tract may reach levels that are higher than that of EGCG and also high enough to activate Nrf2 gene transcription. Taken together, combining in vitro data with a human PBK model allowed the prediction of a dose-response curve for EGCG-induced Nrf2-mediated gene expression in humans and provided insights into the contribution of gut microbial metabolites to this effect.

Assessment of the predictive capacity of a physiologically based kinetic model using a read-across approach

•

Potential regulatory application of PBK modelling information to assist read-across.

•

Presents workflow to read across PBK model information from data-rich to data-poor chemicals.

•

Describes appropriate analogue selection based on a set of specific criteria.

•

Uses estragole and safrole as source chemicals for a target chemical - methyleugenol.

•

Example of PBK model validation where in vivo kinetic data are lacking.

With current progress in science, there is growing interest in developing and applying Physiologically Based Kinetic (PBK) models in chemical risk assessment, as knowledge of internal exposure to chemicals is critical to understanding potential effects in vivo. In particular, a new generation of PBK models is being developed in which the model parameters are derived from in silico and in vitro methods. To increase the acceptance and use of these “Next Generation PBK models”, there is a need to demonstrate their validity. However, this is challenging in the case of data-poor chemicals that are lacking in kinetic data and for which predictive capacity cannot, therefore, be assessed. The aim of this work is to lay down the fundamental steps in using a read across framework to inform modellers and risk assessors on how to develop, or evaluate, PBK models for chemicals without in vivo kinetic data. The application of a PBK model that takes into account the absorption, distribution, metabolism and excretion characteristics of the chemical reduces the uncertainties in the biokinetics and biotransformation of the chemical of interest. A strategic flow-charting application, proposed herein, allows users to identify the minimum information to perform a read-across from a data-rich chemical to its data-poor analogue(s). The workflow analysis is illustrated by means of a real case study using the alkenylbenzene class of chemicals, showing the reliability and potential of this approach. It was demonstrated that a consistent quantitative relationship between model simulations could be achieved using models for estragole and safrole (source chemicals) when applied to methyleugenol (target chemical). When the PBK model code for the source chemicals was adapted to utilise input values relevant to the target chemical, simulation was consistent between the models. The resulting PBK model for methyleugenol was further evaluated by comparing the results to an existing, published model for methyleugenol, providing further evidence that the approach was successful. This can be considered as a “read-across” approach, enabling a valid PBK model to be derived to aid the assessment of a data poor chemical.

Gaining acceptance in next generation PBK modelling approaches for regulatory assessments - An OECD international effort

Physiologically Based Kinetic (PBK) models are valuable tools to help define safe external levels of chemicals based on internal doses at target organs in experimental animals, humans and organisms used in environmental risk assessment. As the toxicity testing paradigm shifts to alternative testing approaches, PBK model development has started to rely (mostly or entirely) on model parameters quantified using in vitro or in silico methods. Recently, the Organisation for Economic Cooperation and Development (OECD) published a guidance document (GD) describing a scientific workflow for characterising and validating PBK models developed using in vitro and in silico data. The GD provides an assessment framework for evaluating these models, with emphasis on the major uncertainties underlying model inputs and outputs. To help end-users submit or evaluate a PBK model for regulatory purposes, the GD also includes a template for documenting the model characteristics, and a checklist for evaluating the quality of a model. This commentary highlights the principles, criteria and tools laid out in the OECD PBK model GD, with the aim of facilitating the dialogue between model developers and risk assessors.