Moving from external exposure concentration to internal dose: duration extrapolation based on physiologically based pharmacokinetic derived estimates of internal dose

Recent Progress on Physiologically Based Pharmacokinetic (PBPK) Model: A Review Based on Bibliometrics

Physiologically based pharmacokinetic/toxicokinetic (PBPK/PBTK) models are designed to elucidate the mechanism of chemical compound action in organisms based on the physiological, biochemical, anatomical, and thermodynamic properties of organisms. After nearly a century of research and practice, good results have been achieved in the fields of medicine, environmental science, and ecology. However, there is currently a lack of a more systematic review of progress in the main research directions of PBPK models, especially a more comprehensive understanding of the application in aquatic environmental research. In this review, a total of 3974 articles related to PBPK models from 1996 to 24 March 2024 were collected. Then, the main research areas of the PBPK model were categorized based on the keyword co-occurrence maps and cluster maps obtained by CiteSpace. The results showed that research related to medicine is the main application area of PBPK. Four major research directions included in the medical field were "drug assessment", "cross-species prediction", "drug-drug interactions", and "pediatrics and pregnancy drug development", in which "drug assessment" accounted for 55% of the total publication volume. In addition, bibliometric analyses indicated a rapid growth trend in the application in the field of environmental research, especially in predicting the residual levels in organisms and revealing the relationship between internal and external exposure. Despite facing the limitation of insufficient species-specific parameters, the PBPK model is still an effective tool for improving the understanding of chemical-biological effectiveness and will provide a theoretical basis for accurately assessing potential risks to ecosystems and human health. The combination with the quantitative structure-activity relationship model, Bayesian method, and machine learning technology are potential solutions to the previous research gaps.

Dietary exposure and accumulation of per- and polyfluoroalkyl substances alters growth and reduces body condition of post-metamorphic salamanders

Per- and polyfluoroalkyl substances (PFAS) are contaminants of concern due to their persistence, potential to bioaccumulate, and toxicity. While dietary exposure is the primary route of exposure for terrestrial species, data on dietary PFAS uptake and adverse effects are largely restricted to mammals. As such, substantial data gaps exist that hinder ecological risk assessment, including environmentally relevant exposure levels and taxa. Using a 30-d laboratory experiment, we examined the effects of dietary PFAS-exposure on post-metamorphic tiger salamanders (Ambystoma tigrinum). We fed salamanders crickets exposed to perfluorooctanoic acid (PFOA), perfluorooctane sulfonate (PFOS), perfluorohexane sulfonate (PFHxS), or 6:2 fluorotelomer sulfonate (6:2 FTS) at low (<1.0), medium (2-5), or high (16-62) ng PFAS/g/d (wet weight) dose rates. We found that only PFOS resulted in substantial biomagnification. Despite limited evidence for biomagnification, PFAS altered growth and generally reduced body condition. Salamanders with the highest burdens of PFOS grew less and had lower body conditions, while burdens of PFHxS and PFOA were only associated with reduced growth. There was no evidence that environmentally relevant doses of PFAS increase liver size in salamanders. Our results demonstrate that dietary exposure and accumulation of PFAS can impact fitness-related traits in amphibians and contribute to trophic transfer in terrestrial food webs.

Toxicity and accumulation of 6-OH-BDE-47 and newly synthesized 6,6'-diOH-BDE-47 in early life-stages of Zebrafish (Danio rerio)

Dihydroxylated polybrominated diphenyl ethers (diOH-PBDEs) appear to be natural products or metabolites of PBDEs in some marine organisms, yet its toxicity is still largely unknown. With a newly lab-synthesized diOH-PBDE, 6,6'-dihydroxy-2,2',4'4'-tetrabromodiphenyl ether (6,6'-diOH-BDE-47) in hand, the present study has provided the first data set to compare 6-hydroxy-2,2',4'4'- tetrabromodiphenyl ether (6-OH-BDE-47) and 6,6'-diOH-BDE-47 for their acute toxicity and accumulation, and thyroid hormone levels in treated zebrafish larvae. By real time-PCR technique, transcripts of hypothalamic-pituitary-thyroid axis associated genes were also investigated in developing larvae at 96 h post fertilization (96 hpf). Apparently, 6,6'-diOH-BDE-47 was less toxic than that of 6-OH-BDE-47: 1) the 96-h LC₅₀ (96-h median lethal concentration) of 6-OH-BDE-47 and 6,6'-diOH-BDE-47 were 235 nM and 516 nM, respectively; 2) although severe developmental delays and morphological deformities were observed in zebrafish larvae in high exposure doses, at the exposure concentration of 1-50 nM, the accumulated 6-OH-BDE-47 and 6,6'-diOH-BDE-47 is ranged between 226-2279 nmol/g and 123-539 nmol/g in treated larvae; and 3) for 6-OH-BDE-47, its bioconcentration factor (BCF) were 1.83- to 4.30-fold more than that of 6,6'-diOH-BDE-47, suggesting that the lower internal exposure concentration of 6,6'-diOH-BDE-47 may lead to lower toxicity. The increased thyroid hormone levels were recorded for 1 nM of 6-OH-BDE-47 and 20 nM of 6,6'-diOH-BDE-47, and the exposures both significantly increased thyroid gland-specific transcription of thyroglobulin gene, indicating an adverse effect associated with the HPT axis. Therefore, 6,6'-diOH-BDE-47, with lower toxicity compared to that of 6-OH-BDE-47, still possesses hazards and environmental risk.

Nonlinear dose-time-response functions and health-protective exposure limits for inflammation-mediated diseases

Why have occupational safety regulations in the United States not been more successful in protecting worker health from mesothelioma risks, while apparently succeeding relatively well in reducing silicosis risks? This paper briefly discusses biological bases for thresholds and nonlinearities in exposure-response functions for respirable crystalline silica (RCS) and asbestos, based on modeling a chronic inflammation mode of action (mediated by activation of the NLRP3 inflammasome, for both RCS and asbestos). It applies previously published physiologically based pharmacokinetic (PBPK) models to perform computational experiments illuminating how different time courses of exposure with the same time-weighted average (TWA) concentration affect internal doses in target tissues (lung for RCS and mesothelium for asbestos). Key conclusions are that (i) For RCS, but not asbestos, limiting average (TWA) exposure concentrations also tightly constrains internal doses and ability to trigger chronic inflammation and resulting increases in disease risks (ii) For asbestos, excursions (i.e., spikes in concentrations); and especially the times between them are crucial drivers of internal doses and time until chronic inflammation; and hence (iii) These dynamic aspects of exposure, which are not addressed by current occupational safety regulations, should be constrained to better protect worker health. Adjusting permissible average exposure concentration limits (PELs) and daily excursion limits (ELs) is predicted to have little impact on reducing mesothelioma risks, but increasing the number of days between successive excursions is predicted to be relatively effective in reducing worker risks, even if it has little or no impact on TWA average concentrations.

The role of chemical interactions between thorium, cerium, and lanthanum in lymphocyte toxicity

Thorium, cerium, and lanthanum are metals present in several types of minerals, the most common of which is monazite. Cerium and lanthanum are elements in the lanthanides series. Thorium, an actinide metal, is a hazardous element due to its radioactive characteristics. There is a lack of information describing the possible chemical interactions among these elements and the effects they may have on humans. Toxicological analyses were performed using cell viability, cell death, and DNA damage assays. Chemical interactions were evaluated based on the Loewe additivity model. The results indicate that thorium and cerium individually have no toxic effects on lymphocytes. However, thorium associated with lanthanum increases the toxicity of this element, thereby reducing the viability of lymphocytes at low concentrations of metals in the mixture.

Blood transcriptomics: applications in toxicology

The number of new chemicals that are being synthesized each year has been steadily increasing. While chemicals are of immense benefit to mankind, many of them have a significant negative impact, primarily owing to their inherent chemistry and toxicity, on the environment as well as human health. In addition to chemical exposures, human exposures to numerous non-chemical toxic agents take place in the environment and workplace. Given that human exposure to toxic agents is often unavoidable and many of these agents are found to have detrimental human health effects, it is important to develop strategies to prevent the adverse health effects associated with toxic exposures. Early detection of adverse health effects as well as a clear understanding of the mechanisms, especially at the molecular level, underlying these effects are key elements in preventing the adverse health effects associated with human exposure to toxic agents. Recent developments in genomics, especially transcriptomics, have prompted investigations into this important area of toxicology. Previous studies conducted in our laboratory and elsewhere have demonstrated the potential application of blood gene expression profiling as a sensitive, mechanistically relevant and practical surrogate approach for the early detection of adverse health effects associated with exposure to toxic agents. The advantages of blood gene expression profiling as a surrogate approach to detect early target organ toxicity and the molecular mechanisms underlying the toxicity are illustrated and discussed using recent studies on hepatotoxicity and pulmonary toxicity. Furthermore, the important challenges this emerging field in toxicology faces are presented in this review article.

Dynamic exposure of organisms and passive samplers to hydrophobic chemicals

An insight into the dynamic aspects of the accumulation process is essential for understanding bioaccumulation as well as effect studies of hydrophobic organic chemicals. This review presents an overview of kinetic studies with organisms (fish, bivalve, crustacean, insect, worm, algae, and protozoan) as well as passive samplers (solid and liquid phase microextraction, semipermeable membrane device, polymer sheet, solid-phase extraction, Chemcatcher, etc.) for the uptake of neutral nonpolar chemicals from the aqueous phase. Information about uptake rates, elimination rates, and 95% equilibration times was collected and analyzed with diffusion based models. The present literature review suggests that the surface to volume ratio appears to be a critical parameter for the uptake rate of the more hydrophobic chemicals both for samplers and organisms. In addition, as a very first approximation, the combination of the first-order kinetic model with the assumption that diffusion through the aqueous boundary layers is rate limiting, gives a reasonable description of the experimental kinetic data. In this way, the presented model might be used to estimate uptake and elimination rate constants of chemicals by organisms or passive samplers.

Role of exhaled breath biomarkers in environmental health science

As a discipline of public health, environmental health science is the study of the linkage from environmental pollution sources to eventual adverse health outcome. This progression may be divided into two components, (1) "exposure assessment," which deals with the source terms, environmental transport, human exposure routes, and internal dose, and (2) "health effects," which deals with metabolism, cell damage, DNA changes, pathology, and onset of disease. The primary goal of understanding the linkage from source to health outcome is to provide the most effective and efficient environmental intervention methods to reduce health risk to the population. Biomarker measurements address an individual response to a common external environmental stressor. Biomarkers are substances within an individual and are subdivided into chemical markers, exogenous metabolites, endogenous response chemicals, and complex adducts (e.g., proteins, DNA). Standard biomarker measurements are performed in blood, urine, or other biological media such as adipose tissue and lavage fluid. In general, sample collection is invasive, requires medical personnel and a controlled environment, and generates infectious waste. Exploiting exhaled breath as an alternative or supplement to established biomarker measurements is attractive primarily because it allows a simpler collection procedure in the field for numerous individuals. Furthermore, because breath is a gas-phase matrix, volatile biomarkers become more readily accessible to analysis. This article describes successful environmental health applications of exhaled breath and proposes future research directions from the perspective of U.S. Environmental Protection Agency (EPA) human exposure research.

Approaches for applications of physiologically based pharmacokinetic models in risk assessment

Physiologically based pharmacokinetic (PBPK) models are particularly useful for simulating exposures to environmental toxicants for which, unlike pharmaceuticals, there is often little or no human data available to estimate the internal dose of a putative toxic moiety in a target tissue or an appropriate surrogate. This article reviews the current state of knowledge and approaches for application of PBPK models in the process of deriving reference dose, reference concentration, and cancer risk estimates. Examples drawn from previous U.S. Environmental Protection Agency (EPA) risk assessments and human health risk assessments in peer-reviewed literature illustrate the ways and means of using PBPK models to quantify the pharmacokinetic component of the interspecies and intraspecies uncertainty factors as well as to conduct route to route, high dose to low dose and duration extrapolations. The choice of the appropriate dose metric is key to the use of the PBPK models for the various applications in risk assessment. Issues related to whether uncertainty factors are most appropriately applied before or after derivation of human equivalent dose (or concentration) continue to be explored. Scientific progress in the understanding of life stage and genetic differences in dosimetry and their impacts on variability in susceptibility, as well as ongoing development of analytical methods to characterize uncertainty in PBPK models, will make their use in risk assessment increasingly likely. As such, it is anticipated that when PBPK models are used to express adverse tissue responses in terms of the internal target tissue dose of the toxic moiety rather than the external concentration, the scientific basis of, and confidence in, risk assessments will be enhanced.

Integrating epidemiology and toxicology in neurotoxicity risk assessment

Neurotoxicity risk assessments depend on the best available scientific information, including data from animal toxicity studies, human experimental studies and human epidemiology studies. There are several factors to consider when evaluating the comparability of data from studies. Regarding the epidemiology literature, issues include choice of study design, use of appropriate controls, methods of exposure assessment, subjective or objective evaluation of neurological status, and assessment and statistical control of potential confounding factors, including co-exposure to other agents. Animal experiments must be evaluated regarding factors such as dose level and duration, procedures used to assess neurological or behavioural status, and appropriateness of inference from the animal model to human neurotoxicity. Major factors that may explain apparent differences between animal and human studies include: animal neurological status may be evaluated with different procedures than those used in humans; animal studies may involve shorter exposure durations and higher dose levels; and most animal studies evaluate a single substance whereas humans typically are exposed to multiple agents. The comparability of measured outcomes in animals and humans may be improved by considering functional domains rather than individual test measures. The application of predictive models, weight of evidence considerations and meta-analysis can help evaluate the consistency of outcomes across studies. An appropriate blend of scientific information from toxicology and epidemiology studies is necessary to evaluate potential human risks of exposure to neurotoxic substances.

Whole-body physiologically based pharmacokinetic models

This review summarizes the most recent developments in and applications of physiologically based pharmacokinetic (PBPK) modeling methodology originating from both the pharmaceutical and environmental toxicology areas. It focuses on works published in the last 5 years, although older seminal papers have also been referenced. After a brief introduction to the field and several essential definitions, the main body of the text is structured to follow the major steps of a typical PBPK modeling exercise. Various applications of the methodology are briefly described. The major future trends and perspectives are outlined. The main conclusion from the review of the available literature is that PBPK modeling, despite its obvious potential and recent incremental developments, has not taken the place it deserves, especially in pharmaceutical and drug development sciences.

Acute Toluene Exposure and Rat Visual Function in Proportion to Momentary Brain Concentration

Acute exposure to toluene was assessed in two experiments to determine the relationship between brain toluene concentration and changes in neurophysiological function. The concentration of toluene in brain tissue at the time of assessment was estimated using a physiologically based pharmacokinetic model. Brain neurophysiological function was measured using pattern-elicited visual evoked potentials (VEP) recorded from electrodes located over visual cortex of adult male Long-Evans rats. In the first experiment, VEPs were recorded before and during exposure to control air or toluene at 1000 ppm for 4 h, 2000 ppm for 2 h, 3000 ppm for 1.3 h, or 4000 ppm for 1 h. In the second experiment, VEPs were recorded during and after exposure to clean air or 3000 or 4000 ppm toluene. In both experiments, the response amplitude of the major spectral component of the VEP (F2 at twice the stimulus rate in steady-state responses) was reduced by toluene. A logistic function was fit to baseline-adjusted F2 amplitudes from the first experiment that described a significant relationship between brain toluene concentration and VEP amplitude deficits. In the second experiment, 3000 ppm caused equivalent VEP deficits during or after exposure as a function of estimated brain concentration, but 4000 ppm showed a rapid partial adaptation to the acute effects of toluene after exposure. In general, however, the neurophysiological deficits caused by acute toluene exposure could be described by estimates of the momentary concentration of toluene in the brain at the time of VEP evaluation.

Model studies for evaluating the neurobehavioral effects of complex hydrocarbon solvents

As part of a project designed to develop a framework for extrapolating acute central nervous system (CNS) effects of hydrocarbon solvents in animals to humans, experimental studies were conducted in rats and human volunteers in which acute CNS effects were measured and toxicokinetic data were collected. A complex hydrocarbon solvent, white spirit (WS) was used as a model solvent and two marker compounds for WS, 1,2,4-trimethyl benzene (TMB) and n-decane (NDEC), were analyzed to characterize internal exposure after WS inhalation. Toxicokinetic data on blood and brain concentrations of the two marker compounds in the rat, together with in vitro partition coefficients were used to develop physiologically based pharmacokinetic (PBPK) models for TMB and NDEC. The rat models were then allometrically scaled to obtain models for inhalatory exposure for man. The human models were validated with blood and alveolar air kinetics of TMB and NDEC, measured in human volunteers. Using these models, it was predicted that external exposures to WS in the range of 344-771mg/m(3) would produce brain concentrations similar to those in rats exposed to 600mg/m(3) WS, the no effect level (NOEL) for acute CNS effects. Assuming similar brain concentration-effect relations for humans and rats, the NOEL for acute CNS effects in humans should be in this range. The prediction was consistent with data from a human volunteer study in which the only statistically significant finding was a small change in the simple reaction time test following 4h exposure to approximately 570mg/m(3) WS. Thus, the data indicated that the results of animal studies could be used to predict a no effect level for acute CNS depression in humans, consistent with the framework described above.

Momentary Brain Concentration of Trichloroethylene Predicts the Effects on Rat Visual Function

The relationship between the concentration of trichloroethylene (TCE) in the brain and changes in brain function, indicated by the amplitude of steady-state pattern-elicited visual evoked potentials (VEP), was evaluated in Long-Evans rats. VEPs were recorded from visual cortex following stimulation of the eyes and, thus, reflect the function of the afferent visual pathway and, in broad terms, may be indicative of overall brain function. The concentration of TCE in the brain at the time of VEP testing (i.e., momentary brain concentration) was hypothesized to predict the amplitude of the VEP across a range of inhalation concentrations, both during and after exposure. Awake restrained rats were exposed to clean air or TCE in the following combinations of concentration and duration: 500 ppm (4 h), 1000 ppm (4 h), 2000 (2 h), 3000 ppm (1.3 h), 4000 ppm (1 h), and 5000 ppm (0.8 h). VEPs were recorded several times during the exposure session, and afterward for experimental sessions of less than 4 h total duration (i.e., concentrations from 2000 to 5000 ppm). The sample collection time for each VEP was about 1 min. Brain concentrations of TCE were predicted using a physiologically based pharmacokinetic (PBPK) model. VEP waveforms were submitted to spectral analysis, and the amplitude of the largest response component, occurring at twice the temporal stimulation rate (F2), was measured. Exposure to all air concentrations of TCE in the study reduced F2 amplitude. The reduction of F2 amplitude was proportional to momentary brain TCE concentration during and after exposure. A logistical function fit to the combined data from all exposure conditions described a statistically significant relationship with 95% confidence limits between brain TCE concentration and F2 amplitude. The results support the hypothesis that momentary brain concentration of TCE is an appropriate dose metric to describe the effects of acute TCE inhalation exposure on rat VEPs. The combination of the PBPK model predicting brain TCE concentration from the exposure conditions with the logistical function predicting F2 amplitude from the brain TCE concentration constitute a quantitative exposure-dose-response model describing an acute change in neurological function following exposure to an important hazardous air pollutant.