Evaluating the impact of anatomical and physiological variability on human equivalent doses using PBPK models.

Addressing human anatomical and physiological variability is a crucial component of human health risk assessment of chemicals. Experts have recommended probabilistic chemical risk assessment paradigms in which distributional adjustment factors are used to account for various sources of uncertainty and variability, including variability in the pharmacokinetic behavior of a given substance in different humans. In practice, convenient assumptions about the distribution forms of adjustment factors and human equivalent doses (HEDs) are often used. Parameters such as tissue volumes and blood flows are likewise often assumed to be lognormally or normally distributed without evaluating empirical data for consistency with these forms. In this work, we performed dosimetric extrapolations using physiologically based pharmacokinetic (PBPK) models for dichloromethane (DCM) and chloroform that incorporate uncertainty and variability to determine if the HEDs associated with such extrapolations are approximately lognormal and how they depend on the underlying distribution shapes chosen to represent model parameters. We accounted for uncertainty and variability in PBPK model parameters by randomly drawing their values from a variety of distribution types. We then performed reverse dosimetry to calculate HEDs based on animal points of departure for each set of sampled parameters. Corresponding samples of HEDs were tested to determine the impact of input parameter distributions on their central tendencies, extreme percentiles, and degree of conformance to lognormality. This work demonstrates that the measurable attributes of human variability should be considered more carefully and that generalized assumptions about parameter distribution shapes may lead to inaccurate estimates of extreme percentiles of HEDs.

A Model Template Approach for Rapid Evaluation and Application of Physiologically Based Pharmacokinetic Models: Extension to Volatile Organic Compounds.

Chemical risk assessors use physiologically based pharmacokinetic (PBPK) models to perform dosimetric calculations, including extrapolations between exposure scenarios, species, and populations of interest. Assessors should complete a thorough quality assurance (QA) review to ensure biological accuracy and correct implementation prior to using these models. This process can be time-consuming, and we developed a PBPK model template that allows for faster, more efficient QA review. The model template consists of a single model "superstructure" with equations and logic commonly found in PBPK models, allowing users to implement a wide variety of chemical-specific PBPK models. QA review can be completed more quickly than for conventional PBPK model implementations because the general model equations have already been reviewed and only parameters describing chemical-specific model and exposure scenarios need review for any given model implementation. We have expanded a previous version of the PBPK model template by adding features commonly included in PBPK models for volatile organic compounds (VOCs). We included multiple options for representing concentrations in blood, describing metabolism, and modeling gas exchange processes to allow for inhalation exposures. We created PBPK model template implementations of published models for seven VOCs: dichloromethane, methanol, chloroform, styrene, vinyl chloride, trichloroethylene, and carbon tetrachloride. Simulations performed using our template implementations matched published simulation results to a high degree of accuracy (maximum observed percent error: 1%). Thus, the model template approach can now be applied to a broader class of chemical-specific PBPK models while continuing to bolster efficiency of QA processes that should be conducted prior to using models for risk assessment applications.

Task-specific modulation of corticospinal neuron activity during motor learning in mice.

Motor skill learning relies on the plasticity of the primary motor cortex as task acquisition drives cortical motor network remodeling. Large-scale cortical remodeling of evoked motor outputs occurs during the learning of corticospinal-dependent prehension behavior, but not simple, non-dexterous tasks. Here we determine the response of corticospinal neurons to two distinct motor training paradigms and assess the role of corticospinal neurons in the execution of a task requiring precise modulation of forelimb movement and one that does not. In vivo calcium imaging in mice revealed temporal coding of corticospinal activity coincident with the development of precise prehension movements, but not more simplistic movement patterns. Transection of the corticospinal tract and optogenetic regulation of corticospinal activity show the necessity for patterned corticospinal network activity in the execution of precise movements but not simplistic ones. Our findings reveal a critical role for corticospinal network modulation in the learning and execution of precise motor movements.

Simulating toxicokinetic variability to identify susceptible and highly exposed populations.

BACKGROUND: Toxicokinetic (TK) data needed for chemical risk assessment are not available for most chemicals. To support a greater number of chemicals, the U.S. Environmental Protection Agency (EPA) created the open-source R package "httk" (High Throughput ToxicoKinetics). The "httk" package provides functions and data tables for simulation and statistical analysis of chemical TK, including a population variability simulator that uses biometrics data from the National Health and Nutrition Examination Survey (NHANES). OBJECTIVE: Here we modernize the "HTTK-Pop" population variability simulator based on the currently available data and literature. We provide explanations of the algorithms used by "httk" for variability simulation and uncertainty propagation. METHODS: We updated and revised the population variability simulator in the "httk" package with the most recent NHANES biometrics (up to the 2017-18 NHANES cohort). Model equations describing glomerular filtration rate (GFR) were revised to more accurately represent physiology and population variability. The model output from the updated "httk" package was compared with the current version. RESULTS: The revised population variability simulator in the "httk" package now provides refined, more relevant, and better justified estimations. SIGNIFICANCE: Fulfilling the U.S. EPA's mission to provide open-source data and models for evaluations and applications by the broader scientific community, and continuously improving the accuracy of the "httk" package based on the currently available data and literature.

A Model Template Approach for Rapid Evaluation and Application of Physiologically Based Pharmacokinetic Models for Use in Human Health Risk Assessments: A Case Study on Per- and Polyfluoroalkyl Substances.

Physiologically based pharmacokinetic (PBPK) models are commonly used in risk assessments to perform inter- and intraspecies extrapolations as well as to extrapolate between different dosing scenarios; however, they must first undergo quality assurance review, which can be a time-consuming process, especially when model code is not readily available. We developed and implemented (using R and MCSim) a PBPK model template capable of replicating published model results for several chemical-specific PBPK models. This model template allows for faster quality assurance review because the general model equations only need to be reviewed once, and application to a specific chemical then only requires reviewing input parameters. The model template can implement PBPK models with oral and intravenous exposure routes, varying numbers of tissue compartments, renal reabsorption, and multiple elimination pathways, including fecal, urinary, and biliary. Using the model template, we reproduced published model simulation results for perfluorohexanesulfonic acid, perfluorononanoic acid, perfluorodecanoic acid, perfluorooctanoate, and perflouorooctane sulfonate. We also show that the template can be a useful tool for identifying potential model errors. Thus, the model template allows for faster evaluation and review of published PBPK models and provides a proof of concept for using this approach with broader classes of chemical-specific PBPK models.