Toxicokinetic Modeling of Per- and Polyfluoroalkyl Substance Concentrations within Developing Zebrafish (Danio rerio) Populations.

Per- and polyfluoroalkyl substances (PFAS) are pervasive environmental contaminants, and their relative stability and high bioaccumulation potential create a challenging risk assessment problem. Zebrafish (Danio rerio) data, in principle, can be synthesized within a quantitative adverse outcome pathway (qAOP) framework to link molecular activity with individual or population level hazards. However, even as qAOP models are still in their infancy, there is a need to link internal dose and toxicity endpoints in a more rigorous way to further not only qAOP models but adverse outcome pathway frameworks in general. We address this problem by suggesting refinements to the current state of toxicokinetic modeling for the early development zebrafish exposed to PFAS up to 120 h post-fertilization. Our approach describes two key physiological transformation phenomena of the developing zebrafish: dynamic volume of an individual and dynamic hatching of a population. We then explore two different modeling strategies to describe the mass transfer, with one strategy relying on classical kinetic rates and the other incorporating mechanisms of membrane transport and adsorption/binding potential. Moving forward, we discuss the challenges of extending this model in both timeframe and chemical class, in conjunction with providing a conceptual framework for its integration with ongoing qAOP modeling efforts.

Case study on the impact of the source of metabolism parameters in next generation physiologically based pharmacokinetic models: Implications for occupational exposures to trimethylbenzenes.

Physiologically based pharmacokinetic (PBPK) models are a means of making important linkages between exposure assessment and in vitro toxicity. A key constraint on rapid application of PBPK models in risk assessment is traditional reliance on substance-specific in vivo toxicokinetic data to evaluate model quality. Bounding conditions, in silico, in vitro, and chemical read-across approaches have been proposed as alternative sources for metabolic clearance estimates. A case study to test consistency of predictive ability across these approaches was conducted using trimethylbenzenes (TMB) as prototype chemicals. Substantial concordance was found among TMB isomers with respect to accuracy (or inaccuracy) of approaches to estimating metabolism; for example, the bounding conditions never reproduced the human in vivo toxicokinetic data within two-fold. Using only approaches that gave acceptable prediction of in vivo toxicokinetics for the source compound (1,2,4-TMB) substantially narrowed the range of plausible internal doses for a given external dose for occupational, emergency response, and environmental/community health risk assessment scenarios for TMB isomers. Thus, risk assessments developed using the target compound models with a constrained subset of metabolism estimates (determined for source chemical models) can be used with greater confidence that internal dosimetry will be estimated with accuracy sufficient for the purpose at hand.

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.

Physiologically based pharmacokinetic (PBPK) modeling of perfluorohexane sulfonate (PFHxS) in humans.

Per- and polyfluoroalkyl substances (PFAS) are persistent, man-made compounds prevalent in the environment and consistently identified in human biomonitoring samples. In particular, perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and perfluorohexane sulfonic acid (PFHxS) have been identified at U.S. Air Force installations. The study of human toxicokinetics and physiologically based pharmacokinetic (PBPK) modeling of PFHxS has been less robust and has been limited in scope and application as compared to PFOS and PFOA. The primary goal of the current effort was to develop a PBPK model describing PFHxS disposition in humans that can be applied to retrospective, current, and future human health risk assessment of PFHxS. An existing model developed for PFOS and PFOA was modified and key parameter values for exposure and toxicokinetics were calibrated for PFHxS prediction based on human biomonitoring data, particularly general population serum levels from the U.S. Centers for Disease Prevention and Control (CDC) National Health and Nutrition Examination Survey (NHANES). Agreement between the model and the calibration and evaluation data was excellent and recapitulated observed trends across sex, age, and calendar years. Confidence in the model is greatest for application to adults in the 2000-2018 time frame and for shorter-term future projections.

Probabilistic pharmacokinetic modeling of airborne lead corresponding to toxicologically relevant blood lead levels in workers.

The Department of Defense (DOD) commissioned the National Research Council (NRC) to assess the potential health effects associated with exposure of DOD personnel to lead (Pb) at firing ranges. In that report, NRC concluded that the current Occupational Safety and Health Administration permissible exposure limit and the blood lead levels (BLLs) on which it was based were not sufficiently protective of worker populations covered under the general industry standard. In support of future selection of an occupational exposure limit, the relationship of airborne Pb levels to BLLs is of interest to the DOD. A subset of the BLLs identified as relevant to the management of health risks of exposed workers was selected as targets for extrapolation to equivalent airborne Pb values. The existing O'Flaherty physiologically based pharmacokinetic model for Pb in humans was modified to facilitate probabilistic predictions of DOD worker population BLLs, including 95th percentile values, based on current worker characteristics. Workplace airborne Pb 8-h time-weighted average concentrations of 1.1, 4.0, 6.8, or 9.8 µg/m3 are anticipated to maintain BLLs below 5, 10, 15, or 20 µg/dl, respectively, in the vast majority of DOD workers exposed to Pb under full-time working lifetime occupational exposure.

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.

Impact of stressors in the aviation environment on xenobiotic dosimetry in humans: physiologically based prediction of the effect of +Gz-forces.

The application of physiologically based modeling approaches in evaluating health risks in diverse environments is limited by scarcity of comprehensive reviews detailing how physiological parameters are altered due to stressors. A modern high-performance aviation environment in particular has the potential for simultaneous exposure to chemical and non-chemical stressors which may interact via non-chemical stressor-mediated pharmacokinetic alterations. To support physiologically based pharmacokinetic (PBPK) modeling of in-flight disposition inhaled chemicals, literature review, and synthesis was conducted to determine the impact of gravitational (+Gz) forces on PBPK modeling inputs. Specifically, changes in cardiac output and related parameters heart rate and stroke volume, breathing frequency, tidal volume, and pulmonary and alveolar ventilation rate in vivo were extracted from 36 publications and related mathematically to +Gz intensity. A scenario was simulated where a pilot performing test flights might inhale organic chemicals at the occupational exposure guideline level while experiencing sustained, elevated +Gz. Peak arterial blood concentrations of 1,2,4-trimethylbenzene during a 1 h-flight at +4 Gz were predicted to increase 2-fold relative to would occur on the ground under baseline conditions. This case study demonstrates the potential value of scenario-specific physiological information in assessing changes in risk-relevant internal dosimetry, providing better information for potential risk management actions.

Considerations for Development of Exposure Limits for Chemicals Encountered During Aircraft Operation.

BACKGROUND: Military aircrews' health status is critical to their mission readiness, as they perform physically and cognitively demanding tasks in nontraditional work environments. Research Objectives: Our objective is to develop a broad operational risk assessment framework and demonstrate its applicability to health risks to aircrews because of airborne chemical exposure, considering stressors such as heat and exertion. METHODS: Extrapolation of generic exposure standards to military aviation-specific conditions can include computation of risk-relevant internal dosimetry estimates by incorporating changes in breathing patterns and blood flow distribution because of aspects of the in-flight environment. We provide an example of the effects of exertion on peak blood concentrations of 1,2,4-trimethylbenzene computed using a physiologically based pharmacokinetic model. RESULTS: Existing published collections on the effects of flight-related stressors on breathing patterns and blood flow address only a limited number of stressors. Although data exist that can be used to develop operational exposure limits specific to military aircrew activities, efforts to integrate this information in specific chemical assessments have been limited. CONCLUSIONS: Efforts to develop operational exposure limits would benefit from guidance on how to make use of existing assessments and expanded databases of the impact of environmental stressors on adult human physiology.

Route-to-route extrapolation of 1,2-dichloroethane studies from the oral route to inhalation using physiologically based pharmacokinetic models.

To help develop a comprehensive, quantitative understanding of the hazards of 1,2-dichloroethane (ethylene dichloride, EDC, CAS No. 107-06-2) exposure by the inhalation route, the results of existing subchronic studies and an extended one-generation reproductive toxicity (EOGRT) study recently conducted by the oral route in rats were extrapolated using a physiologically based pharmacokinetic (PBPK) model. The no observed adverse effects levels (NOAELs) for the endpoints of neurotoxicity and reproductive/developmental toxicity were the highest tested doses of 169 and 155 mg/kg-day, respectively. These NOAELs were equivalent to continuous exposure of rats to minimums of 76 ppm and 62 ppm EDC, respectively, using total metabolism of EDC as the dose metric that is equivalent in the oral and inhalation scenarios. In contrast, the subchronic study NOAEL of 37.5 mg/kg-day corresponded to continuous inhalation of 4.4 ppm EDC, based on equivalent extrahepatic metabolism. The selection of the internal metric which serves to establish route-to-route equivalency was found to profoundly influence the NOAEL-equivalent inhalation exposure concentration and thus will be a key determinant of inhalation toxicity reference criteria developed on the basis of EDC studies conducted by the oral route.

Acute toxicity when concentration varies with time: A case study with carbon monoxide inhalation by rats.

Exposure to time-varying concentrations of toxic compounds is the norm in both occupational settings and daily human life, but little has been done to investigate the impact of variations in concentration on toxic outcomes; this case study with carbon monoxide helps fill that gap. Median acute lethality of 10-, 20-, 40-, and 60-min continuous exposures of rats to carbon monoxide was well described by the toxic load model (k = C(n) × t; k is constant, C = test concentration, n = toxic load exponent, and t = exposure duration) with n = 1.74. Dose response-relationships for 1-h exposures including a recovery period between 10- or 20-min pulses showed greater similarity (in both median lethality and steepness of dose-response curve) to continuous exposures with equivalent pulse duration and concentration, rather than a 60-min exposure with equivalent time-weighted average concentrations or toxic load. When pulses were of unequal concentration (3:1 ratio), only the high concentration pulse contributed to lethality. These findings show that fluctuations or interruptions in exposure over a short time scale (60 min or less) can have a substantial impact on outcomes (when n > 1), and thus high-resolution monitoring data are needed to aid interpretation of resulting outcomes.

Advances in Inhalation Dosimetry Models and Methods for Occupational Risk Assessment and Exposure Limit Derivation.

The purpose of this article is to provide an overview and practical guide to occupational health professionals concerning the derivation and use of dose estimates in risk assessment for development of occupational exposure limits (OELs) for inhaled substances. Dosimetry is the study and practice of measuring or estimating the internal dose of a substance in individuals or a population. Dosimetry thus provides an essential link to understanding the relationship between an external exposure and a biological response. Use of dosimetry principles and tools can improve the accuracy of risk assessment, and reduce the uncertainty, by providing reliable estimates of the internal dose at the target tissue. This is accomplished through specific measurement data or predictive models, when available, or the use of basic dosimetry principles for broad classes of materials. Accurate dose estimation is essential not only for dose-response assessment, but also for interspecies extrapolation and for risk characterization at given exposures. Inhalation dosimetry is the focus of this paper since it is a major route of exposure in the workplace. Practical examples of dose estimation and OEL derivation are provided for inhaled gases and particulates.

Bayesian evaluation of a physiologically-based pharmacokinetic (PBPK) model of long-term kinetics of metal nanoparticles in rats.

Biomathematical modeling quantitatively describes the disposition of metal nanoparticles in lungs and other organs of rats. In a preliminary model, adjustable parameters were calibrated to each of three data sets using a deterministic approach, with optimal values varying among the different data sets. In the current effort, Bayesian population analysis using Markov chain Monte Carlo (MCMC) simulation was used to recalibrate the model while improving assessments of parameter variability and uncertainty. The previously-developed model structure and some physiological parameter values were modified to improve physiological realism. The data from one of the three previously-identified studies and from two other studies were used for model calibration. The data from the one study that adequately characterized mass balance were used to generate parameter distributions. When data from a second study of the same nanomaterial (iridium) were added, the level of agreement was still acceptable. Addition of another data set (for silver nanoparticles) led to substantially lower precision in parameter estimates and large discrepancies between the model predictions and experimental data for silver nanoparticles. Additional toxicokinetic data are needed to further evaluate the model structure and performance and to reduce uncertainty in the kinetic processes governing in vivo disposition of metal nanoparticles.

Toxicokinetic Model Development for the Insensitive Munitions Component 2,4-Dinitroanisole.

The Armed Forces are developing new explosives that are less susceptible to unintentional detonation (insensitive munitions [IMX]). 2,4-Dinitroanisole (DNAN) is a component of IMX. Toxicokinetic data for DNAN are required to support interpretation of toxicology studies and refinement of dose estimates for human risk assessment. Male Sprague-Dawley rats were dosed by gavage (5, 20, or 80 mg DNAN/kg), and blood and tissue samples were analyzed to determine the levels of DNAN and its metabolite 2,4-dinitrophenol (DNP). These data and data from the literature were used to develop preliminary physiologically based pharmacokinetic (PBPK) models. The model simulations indicated saturable metabolism of DNAN in rats at higher tested doses. The PBPK model was extrapolated to estimate the toxicokinetics of DNAN and DNP in humans, allowing the estimation of human-equivalent no-effect levels of DNAN exposure from no-observed adverse effect levels determined in laboratory animals, which may guide the selection of exposure limits for DNAN.

Toxicokinetic Model Development for the Insensitive Munitions Component 3-Nitro-1,2,4-Triazol-5-One.

3-Nitro-1,2,4-triazol-5-one (NTO) is a component of insensitive munitions that are potential replacements for conventional explosives. Toxicokinetic data can aid in the interpretation of toxicity studies and interspecies extrapolation, but only limited data on the toxicokinetics and metabolism of NTO are available. To supplement these limited data, further in vivo studies of NTO in rats were conducted and blood concentrations were measured, tissue distribution of NTO was estimated using an in silico method, and physiologically based pharmacokinetic models of the disposition of NTO in rats and macaques were developed and extrapolated to humans. The model predictions can be used to extrapolate from designated points of departure identified from rat toxicology studies to provide a scientific basis for estimates of acceptable human exposure levels for NTO.

Risk assessments for chronic exposure of children and prospective parents to ethylbenzene (CAS No. 100-41-4).

Potential chronic health risks for children and prospective parents exposed to ethylbenzene were evaluated in response to the Voluntary Children's Chemical Evaluation Program. Ethylbenzene exposure was found to be predominately via inhalation with recent data demonstrating continuing decreases in releases and both outdoor and indoor concentrations over the past several decades. The proportion of ethylbenzene in ambient air that is attributable to the ethylbenzene/styrene chain of commerce appears to be relatively very small, less than 0.1% based on recent relative emission estimates. Toxicity reference values were derived from the available data, with physiologically based pharmacokinetic models and benchmark dose methods used to assess dose-response relationships. An inhalation non-cancer reference concentration or RfC of 0.3 parts per million (ppm) was derived based on ototoxicity. Similarly, an oral non-cancer reference dose or RfD of 0.5 mg/kg body weight/day was derived based on liver effects. For the cancer assessment, emphasis was placed upon mode of action information. Three of four rodent tumor types were determined not to be relevant to human health. A cancer reference value of 0.48 ppm was derived based on mouse lung tumors. The risk characterization for ethylbenzene indicated that even the most highly exposed children and prospective parents are not at risk for non-cancer or cancer effects of ethylbenzene.

Evaluation of submarine atmospheres: effects of carbon monoxide, carbon dioxide and oxygen on general toxicology, neurobehavioral performance, reproduction and development in rats. II. Ninety-day study.

Carbon monoxide (CO), carbon dioxide (CO2) and low-level oxygen (O2) (hypoxia) are submarine atmosphere components of highest concern because of a lack of toxicological data available to address the potential effects from long-duration, combined exposures on female reproductive and developmental health. In this study, subchronic toxicity of mixed atmospheres of these three submarine air components was evaluated in rats. Male and female rats were exposed via inhalation to clean air (0.4 ppm CO; 0.13% CO2; 20.6% O2) (control), a low-dose (5.0 ppm CO; 0.41% CO2; 17.1% O2), a mid-dose (13.9 ppm CO; 1.19 or 1.20% CO2; 16.1% O2) and a high-dose (89.9 ppm CO; 2.5% CO2; 15.0% O2) gas mixture for 23 h per day for 70 d premating and a 14-d mating period. Impregnated dams continued exposure to gestation day 19. Adverse reproductive effects were not identified in exposed parents (P0) or first (F1) and second generation (F2) offspring during mating, gestation or parturition. No adverse changes to the estrous cycle or in reproductive hormone concentrations were identified. The exposure-related effects were reduced weight gains and adaptive up-regulation of erythropoiesis in male rats from the high-dose group. No adverse, dose-related health effects on clinical data or physiological data were observed. Neurobehavioral tests identified no apparent developmental deficits at the tested levels of exposure. In summary, subchronic exposures to the submarine atmosphere gases did not affect the ability of the exposed rats or their offspring to reproduce and did not appear to have any significant adverse health effects.

Evaluating the validity and applicable domain of the toxic load model: impact of concentration vs. time profile on inhalation lethality of hydrogen cyanide.

The ten Berge model (or "toxic load" model) is often used to estimate the acute toxicity for varying combinations of inhaled concentration and duration. Expressed as C(n) × t = toxic load (TL), TLs are assumed constant for various combinations of concentration (C) and time (t). Experimental data in a recent acute inhalation study of rats exposed to time-varying concentrations of hydrogen cyanide (HCN) supported the validity of the toxic load model except under very brief, discontinuous, high concentration exposures. In the present investigation, experiments were conducted to extend the evaluation of the applicable domain of the model for acute lethality of HCN in the rat (cumulative exposure range of 2900-11,000 ppm min). The lethality of HCN over very short (< 5 min) durations of high concentrations did not conform to the toxic load model. A value of n=1.57 was determined for uninterrupted exposures ⩾ 5 min. For 30-min exposures, the presence or absence of a gap between two exposure pulses of different concentrations, the relative duration, relative height, and the ordering of the pulses (low then high, vs. high then low) did not appear to have a meaningful impact on the toxic load required for median lethality.

Evaluation of submarine atmospheres: effects of carbon monoxide, carbon dioxide and oxygen on general toxicology, neurobehavioral performance, reproduction and development in rats. I. Subacute exposures.

The inhalation toxicity of submarine contaminants is of concern to ensure the health of men and women aboard submarines during operational deployments. Due to a lack of adequate prior studies, potential general, neurobehavioral, reproductive and developmental toxicity was evaluated in male and female rats exposed to mixtures of three critical submarine atmospheric components: carbon monoxide (CO) and carbon dioxide (CO2; levels elevated above ambient), and oxygen (O2; levels decreased below ambient). In a 14-day, 23 h/day, whole-body inhalation study of exposure to clean air (0.4 ppm CO, 0.1% CO2 and 20.6% O2), low-dose, mid-dose and high-dose gas mixtures (high dose of 88.4 ppm CO, 2.5% CO2 and 15.0% O2), no adverse effects on survival, body weight or histopathology were observed. Reproductive, developmental and neurobehavioral performance were evaluated after a 28-day exposure in similar atmospheres. No adverse effects on estrus phase, mating, gestation or parturition were observed. No developmental or functional deficits were observed in either exposed parents or offspring related to motor activity, exploratory behavior or higher-level cognitive functions (learning and memory). Only minimal effects were discovered in parent-offspring emotionality tests. While statistically significant increases in hematological parameters were observed in the offspring of exposed parents compared to controls, these parameters remained within normal clinical ranges for blood cells and components and were not considered adverse. In summary, subacute exposures to elevated concentrations of the submarine atmosphere gases did not affect the ability of rats to reproduce and did not appear to have any significant adverse health effects.

Impact of non-constant concentration exposure on lethality of inhaled hydrogen cyanide.

The ten Berge model, also known as the toxic load model, is an empirical approach in hazard assessment modeling for estimating the relationship between the inhalation toxicity of a chemical and the exposure duration. The toxic load (TL) is normally expressed as a function of vapor concentration (C) and duration (t), with TL equaling C(n) × t being a typical form. Hypothetically, any combination of concentration and time that yields the same "toxic load" will give a constant biological response. These formulas have been developed and tested using controlled, constant concentration animal studies, but the validity of applying these assumptions to time-varying concentration profiles has not been tested. Experiments were designed to test the validity of the model under conditions of non-constant acute exposure. Male Sprague-Dawley rats inhaled constant or pulsed concentrations of hydrogen cyanide (HCN) generated in a nose-only exposure system for 5, 15, or 30 min. The observed lethality of HCN for the 11 different C versus t profiles was used to evaluate the ability of the model to adequately describe the lethality of HCN under the conditions of non-constant inhalation exposure. The model was found to be applicable under the tested conditions, with the exception of the median lethality of very brief, high concentration, discontinuous exposures.

Subacute effects of inhaled Jet Fuel-A (Jet A) on airway and immune function in female rats.

Two studies were conducted to assess the potential airway and immune effects following subacute (14 d) exposure of female rats to 500, 1000 or 2000 mg/m³ of Jet-A for 4 h/d. The first study used Sprague-Dawley rats; the second study included both Fischer 344 (F344) and Sprague-Dawley rats. In the first study, exposure to 2000 mg/m³ jet fuel may have caused significant upper airway inflammation on day 7 post-exposure, as indicated by elevated protein and lactate dehydrogenase in nasal lavage fluid, but any inflammation resolved by day 14 post-exposure. No significant impact on immune cell populations in the spleens was observed. The histological examination showed no evidence of infectious or toxic effect. In the second study, body weights of the F344 rats in the 2000 mg/m³ group were depressed, as compared to the controls, at the end of the exposure. Some lung lavage fluid markers were increased at 24 h after the final exposure, however, no test article-induced histological changes were observed in the lungs, nasal cavities, or any other tissue of any of the jet fuel exposed animals. Overall, these studies demonstrated limited evidence of effects of 14 d of exposure to Jet A on the airways, immune system, or any other organ or system of female Sprague-Dawley and F344 rats, with no remarkable differences between strains. The lack of identified significant airway or immune effects was in contrast to previous examinations of jet fuel for pulmonary toxicity in mice and rats and for immunotoxicity in mice.