Utility of life stage-specific chemical risk assessments based on New Approach Methodologies (NAMs)

The safety assessments for chemicals targeted for use or expected to be exposed to specific life stages, including infancy, childhood, pregnancy and lactation, and geriatrics, need to account for extrapolation of data from healthy adults to these populations to assess their human health risk. However, often adequate and relevant toxicity or pharmacokinetic (PK) data of chemicals in specific life stages are not available. For such chemicals, New Approach Methodologies (NAMs), such as physiologically based pharmacokinetic (PBPK) modeling, biologically based dose response (BBDR) modeling, in vitro to in vivo extrapolation (IVIVE), etc. can be used to understand the variability of exposure and effects of chemicals in specific life stages and assess their associated risk. A life stage specific PBPK model incorporates the physiological and biochemical changes associated with each life stage and simulates their impact on the absorption, distribution, metabolism, and elimination (ADME) of these chemicals. In our review, we summarize the parameterization of life stage models based on New Approach Methodologies (NAMs) and discuss case studies that highlight the utility of a life stage based PBPK modeling for risk assessment. In addition, we discuss the utility of artificial intelligence (AI)/machine learning (ML) and other computational models, such as those based on in vitro data, as tools for estimation of relevant physiological or physicochemical parameters and selection of model. We also discuss existing gaps in the available toxicological datasets and current challenges that need to be overcome to expand the utility of NAMs for life stage-specific chemical risk assessment.

Ammonium perchlorate: serum dosimetry, neurotoxicity, and resilience of the neonatal rat thyroid system

The environmental contaminant perchlorate impairs the synthesis of thyroid hormones by reducing iodine uptake into the thyroid gland. Despite this known action, moderate doses of perchlorate do not significantly alter serum thyroid hormone in rat pups born to exposed dams. We examined perchlorate dosimetry and responsivity of the thyroid gland and brain in offspring following maternal exposure to perchlorate. Pregnant rat dams were delivered perchlorate in drinking water (0, 30, 100, 300, 1000 ppm) from gestational day 6 to postnatal day (PN) 21. Perchlorate was present in the placenta, milk, and serum, the latter declining in pups over the course of lactation. Serum and brain thyroid hormone were reduced in pups at birth but recovered to control levels by PN2. Dramatic upregulation of Nis was observed in the thyroid gland of the exposed pup. Despite the return of serum thyroid hormone to control levels by PN2, expression of several TH-responsive genes was altered in the PN14 pup brain. Contextual fear learning was unimpaired in the adults, supporting previous reports. Declining levels of serum perchlorate and a profound upregulation of Nis gene expression in the thyroid gland are consistent with the rapid return to the euthyroid state in the neonate. However, despite this recovery, thyroid hormone insufficiencies in serum and brain beginning in utero and present at birth appear sufficient to alter TH action in the fetus and subsequent trajectory of brain development. Biomarkers of that altered trajectory remain in the brain of the neonate, demonstrating that perchlorate is not devoid of effects on the developing brain.

Structural Malformations in the Neonatal Rat Brain Accompany Developmental Exposure to Ammonium Perchlorate

Environmental contaminants are often flagged as thyroid system disruptors due to their actions to reduce serum thyroxine (T4) in rodent models. The presence of a periventricular heterotopia (PVH), a brain malformation resulting from T4 insufficiency, has been described in response to T4 decrements induced by pharmaceuticals that reduce the hormone synthesis enzyme thyroperoxidase. In this report, we extend these observations to the environmental contaminant perchlorate, an agent that interferes with thyroid status by inhibiting iodine uptake into the thyroid gland. Pregnant rat dams were administered perchlorate in their drinking water (0, 30, 100, 300, 1000 ppm) from gestational day (GD) 6 until the weaning of pups on postnatal day (PN) 21. Serum T4 was reduced in dams and fetuses in late gestation and remained lower in lactating dams. Pup serum and brain T4, however, were not reduced beyond PN0, and small PVHs were evident in the brains of offspring when assessed on PN14. To emulate the developmental time window of the brain in humans, a second study was conducted in which pups from perchlorate-exposed dams were administered perchlorate orally from PN0 to PN6. This treatment reduced serum and brain T4 in the pup and resulted in large PVH. A third study extended the period of serum and brain TH suppression in pups by coupling maternal perchlorate exposure with maternal dietary iodine deficiency (ID). No PVHs were evident in the pups from ID dams, small PVHs were observed in the offspring of dams exposed to 300 ppm of perchlorate, and very large PVHs were present in the brains of pups born to dams receiving ID and perchlorate. These findings underscore the importance of the inclusion of serum hormone profiles in pregnant dams and fetuses in in vivo screens for thyroid-system-disrupting chemicals and indicate that chemical-induced decreases in fetal rat serum that resolve in the immediate postnatal period may still harbor considerable concern for neurodevelopment in humans.

Cross species extrapolation of the disruption of thyroid hormone synthesis by oxyfluorfen using in vitro data, physiologically based pharmacokinetic (PBPK), and thyroid hormone kinetics models

The thyroid hormones play key roles in physiological processes such as regulation of the metabolic and cardiac systems as well as the development of the brain and surrounding sympathetic nervous system. Recent efforts to screen environmental chemicals for their ability to alter thyroid hormone synthesis, transport, metabolism and/or function have identified novel chemicals that target key processes in the thyroid pathway. One newly identified chemical, oxyfluorfen, is a diphenyl-ether herbicide used for control of annual broadleaf and grassy weeds in a variety of tree fruit, nut, vine, and field crops. Using in vitro high-throughput screening (HTS) assays, oxyfluorofen was identified to be a potent inhibitor of the thyroidal sodium-iodide symporter (NIS). To quantitatively assess this inhibition mechanism in vivo, we extrapolated in vitro NIS inhibition data to in vivo disruption of thyroid hormones synthesis in rats using physiologically based pharmacokinetic (PBPK) and thyroid hormone kinetics models. The overall computational model (chemical PBPK and THs kinetic sub-models) was calibrated against in vivo data for the levels of oxyfluorfen in thyroid tissue and serum and against serum levels of thyroid hormones triiodothyronine (T3) and thyroxine (T4) in rats. The rat thyroid model was then extrapolated to humans using human in vitro HTS data for NIS inhibition and the chemical specific hepatic clearance rate in humans. The overall species extrapolated PBPK-thyroid kinetics model can be used to predict dose-response (% drop in thyroid serum levels compared to homeostasis) relationships in humans. These relationships can be used to estimate points of departure for health risks related to a drop in serum levels of TH hormones based on HTS assays in vitro to in vivo extrapolation (IVIVE), toxicokinetics, and physiological principles.

Towards translating in vitro measures of thyroid hormone system disruption to in vivo responses in the pregnant rat via a biologically based dose response (BBDR) model

Despite the number of in vitro assays that have been recently developed to identify chemicals that interfere with the hypothalamic-pituitary-thyroid axis (HPT), the translation of those in vitro results into in vivo responses (in vitro to in vivo extrapolation, IVIVE) has received limited attention from the modeling community. To help advance this field a steady state biologically based dose response (BBDR) model for the HPT axis was constructed for the pregnant rat on gestation day (GD) 20. The BBDR HPT axis model predicts plasma levels of thyroid stimulating hormone (TSH) and the thyroid hormones, thyroxine (T4) and triiodothyronine (T3). Thyroid hormones are important for normal growth and development of the fetus. Perchlorate, a potent inhibitor of thyroidal uptake of iodide by the sodium iodide symporter (NIS) protein, was used as a case study for the BBDR HPT axis model. The inhibitory blocking of the NIS by perchlorate was associated with dose-dependent steady state decreases in thyroid hormone production in the thyroid gland. The BBDR HPT axis model predictions for TSH, T3, and T4 plasma concentrations in pregnant Sprague Dawley (SD) rats were within 2-fold of observations for drinking water perchlorate exposures ranging from 10 to 30,000 μg/kg/d. In Long Evans (LE) pregnant rats, for both control and perchlorate drinking water exposures, ranging from 85 to 82,000 μg/kg/d, plasma thyroid hormone and TSH concentrations were predicted within 2 to 3.4- fold of observations. This BBDR HPT axis model provides a successful IVIVE template for thyroid hormone disruption in pregnant rats.

AOP Report: Thyroperoxidase Inhibition Leading to Altered Visual Function in Fish Via Altered Retinal Layer Structure

Thyroid hormones (THs) are involved in the regulation of many important physiological and developmental processes, including vertebrate eye development. Thyroid hormone system-disrupting chemicals (THSDCs) may have severe consequences, because proper functioning of the visual system is a key factor for survival in wildlife. However, the sequence of events leading from TH system disruption (THSD) to altered eye development in fish has not yet been fully described. The development of this adverse outcome pathway (AOP) was based on an intensive literature review of studies that focused on THSD and impacts on eye development, mainly in fish. In total, approximately 120 studies (up to the end of 2021) were used in the development of this AOP linking inhibition of the key enzyme for TH synthesis, thyroperoxidase (TPO), to effects on retinal layer structure and visual function in fish (AOP-Wiki, AOP 363). In a weight-of-evidence evaluation, the confidence levels were overall moderate, with ample studies showing the link between reduced TH levels and altered retinal layer structure. However, some uncertainties about the underlying mechanism(s) remain. Although the current weight-of-evidence evaluation is based on fish, the AOP is plausibly applicable to other vertebrate classes. Through the re-use of several building blocks, this AOP is connected to the AOPs leading from TPO and deiodinase inhibition to impaired swim bladder inflation in fish (AOPs 155-159), together forming an AOP network describing THSD in fish. This AOP network addresses the lack of thyroid-related endpoints in existing fish test guidelines for the evaluation of THSDCs. Environ Toxicol Chem 2022;41:2632-2648. © 2022 SETAC.

Gestational Exposure to Perchlorate in the Rat: Thyroid Hormones in Fetal Thyroid Gland, Serum, and Brain

Iodine is essential for the production of thyroid hormones. Perchlorate is an environmental contaminant that interferes with iodine uptake into the thyroid gland to reduce thyroid hormone synthesis. As thyroid hormones are critical for brain development, exposure to perchlorate during pregnancy is of concern for the developing fetal brain. In this study, we (1) define profiles of thyroid hormone in the maternal and fetal compartments of pregnant rats in response to inhibition of the sodium-iodide symporter (NIS) by perchlorate and (2) expand inquiry previously limited to serum to include fetal thyroid gland and brain. Perchlorate was added to the drinking water (0, 1, 30, 300, and 1000 ppm) of pregnant rat dams from gestational days (GD) 6-20. On GD20, blood, thyroid gland, and brain were collected from the fetus and dam for thyroid hormone and molecular analyses. Thyroid gland and serum thyroid hormones were dose-dependently reduced, with steeper declines evident in the fetus than in the dam. The thyroid gland revealed perturbations of thyroid hormone-action with greater sensitivity in the fetus than the dam. Thyroid hormones and thyroid hormone-responsive gene expression were reduced in the fetal cortex portending effects on brain development. These findings are the first quantitative assessments of perchlorate-induced deficits in the fetal thyroid gland and fetal brain. We provide a conceptual framework to develop a quantitative NIS adverse outcome pathway for serum thyroid hormone deficits and the potential to impact the fetal brain. Such a framework may also serve to facilitate the translation of in vitro bioactivity to the downstream in vivo consequences of NIS inhibition in the developing fetus.

Mechanistic Computational Model for Extrapolating In vitro Thyroid Peroxidase (TPO) Inhibition Data to Predict Serum Thyroid Hormone Levels in Rats

High throughput (HTP) in vitro assays are developed to screen chemicals for their potential to inhibit thyroid hormones (THs) synthesis. Some of these experiments, such as the thyroid peroxidase (TPO) inhibition assay, are based on thyroid microsomal extracts. However, the regulation of thyroid disruption chemicals (TDCs) is based on THs in vivo serum levels. This necessitates the estimation of TDCs in vivo tissue levels in the thyroid where THs synthesis inhibition by TPO takes place. The in vivo tissue levels of chemicals are controlled by pharmacokinetic determinants such as absorption, distribution, metabolism and excretion (ADME), and can be described quantitatively in physiologically based pharmacokinetic (PBPK) models. An integrative computational model including chemical specific PBPK and TH kinetics models provides a mechanistic quantitative approach to translate thyroidal HTP in vitro assays to in vivo measures of circulating THs serum levels. This computational framework is developed to quantitatively establish the linkage between applied dose, chemical thyroid tissue levels, thyroid TPO inhibition potential, and in vivo TH serum levels. Once this link is established quantitively, the overall model is used to calibrate the TH kinetics parameters using experimental data for THs levels in thyroid tissue and serum for the two drugs Propylthiouracil (PTU) and Methimazole (MMI). The calibrated quantitative framework is then evaluated against literature data for the environmental chemical ethylenethiourea (ETU). The linkage of PBPK and TH kinetics models illustrates a computational framework that can be extrapolated to humans to screen chemicals based on their exposure levels and potential to disrupt serum THs levels in vivo.

Extrapolating In Vitro Screening Assay Data for Thyroperoxidase Inhibition to Predict Serum Thyroid Hormones in the Rat

Thyroperoxidase (TPO) is an enzyme essential for thyroid hormone (TH) synthesis and a target site for a number of xenobiotics that disrupt TH homeostasis. An in vitro high-throughput screening assay for TPO inhibition, the Amplex UltraRed-TPO (AUR-TPO), has been used to screen the ToxCast chemical libraries for this action. Output from this assay would be most useful if it could be readily translated into an in vivo response, namely a reduction of TH in serum. To this end, the relationship between TPO inhibition in vitro and serum TH decreases was examined in rats exposed to 2 classic TPO inhibitors, propylthiouracil (PTU) and methimazole (MMI). Serum and gland PTU, MMI, and TH levels were quantified using tandem liquid chromatography mass spectrometry. Thyroperoxidase activity was determined in thyroid gland microsomes treated with PTU or MMI in vitro and ex vivo from thyroid gland microsomes prepared from exposed animals. A quantitative model was constructed by contrasting in vitro and ex vivo AUR-TPO results and the in vivo time-course and dose-response analysis. In vitro:ex vivo correlations of AUR-TPO outputs indicated that less than 30% inhibition of TPO in vitro was sufficient to reduce serum T4 by 20%, a degree of regulatory significance. Although further testing of model estimates using other TPO inhibitors is essential for verification of these initial findings, the results of this study provide a means to translate in vitro screening assay results into predictions of in vivo serum T4 changes to inform risk assessment.

Approaches to Dose Finding in Neonates, Illustrating the Variability between Neonatal Drug Development Programs

Drug dosing in neonates should be based on integrated knowledge concerning the disease to be treated, the physiological characteristics of the neonate, and the pharmacokinetics (PK) and pharmacodynamics (PD) of a given drug. It is critically important that all sources of information be leveraged to optimize dose selection for neonates. Sources may include data from adult studies, pediatric studies, non-clinical (juvenile) animal models, in vitro studies, and in silico models. Depending on the drug development program, each of these modalities could be used to varying degrees and with varying levels of confidence to guide dosing. This paper aims to illustrate the variability between neonatal drug development programs for neonatal diseases that are similar to those seen in other populations (meropenem), neonatal diseases related but not similar to pediatric or adult populations (clopidogrel, thyroid hormone), and diseases unique to neonates (caffeine, surfactant). Extrapolation of efficacy from older children or adults to neonates is infrequently used. Even if a disease process is similar between neonates and children or adults, such as with anti-infectives, additional dosing and safety information will be necessary for labeling, recognizing that dosing in neonates is confounded by maturational PK in addition to body size.

Neurodevelopment and Thyroid Hormone Synthesis Inhibition in the Rat: Quantitative Understanding Within the Adverse Outcome Pathway Framework

Adequate levels of thyroid hormone (TH) are needed for proper brain development, deficiencies may lead to adverse neurologic outcomes in humans and animal models. Environmental chemicals have been linked to TH disruption, yet the relationship between developmental exposures and decline in serum TH resulting in neurodevelopmental impairment is poorly understood. The present study developed a quantitative adverse outcome pathway where serum thyroxin (T4) reduction following inhibition of thyroperoxidase in the thyroid gland are described and related to deficits in fetal brain TH and the development of a brain malformation, cortical heterotopia. Pregnant rats were exposed to 6-propylthiouracil (PTU 0, 0.1, 0.5, 1, 2, or 3 parts per million [ppm]) from gestational days 6-20, sequentially increasing PTU concentrations in maternal thyroid gland and serum as well as in fetal serum. Dams exposed to 0.5 ppm PTU and higher exhibited dose-dependent decreases in thyroidal T4. Serum T4 levels in the dam were significantly decreased with exposure to 2 and 3 ppm PTU. In the fetus, T4 decrements were first observed at a lower dose of 0.5 ppm PTU. Based on these data, fetal brain T4 levels were estimated from published literature sources, and quantitatively linked to increases in the size of the heterotopia present in the brains of offspring. These data show the potential of in vivo assessments and computational descriptions of biologic responses to predict the development of this structural brain malformation and use of quantitative adverse outcome pathway approach to evaluate brain deficits that may result from exposure to other TH disruptors.

How Adverse Outcome Pathways Can Aid the Development and Use of Computational Prediction Models for Regulatory Toxicology

Efforts are underway to transform regulatory toxicology and chemical safety assessment from a largely empirical science based on direct observation of apical toxicity outcomes in whole organism toxicity tests to a predictive one in which outcomes and risk are inferred from accumulated mechanistic understanding. The adverse outcome pathway (AOP) framework provides a systematic approach for organizing knowledge that may support such inference. Likewise, computational models of biological systems at various scales provide another means and platform to integrate current biological understanding to facilitate inference and extrapolation. We argue that the systematic organization of knowledge into AOP frameworks can inform and help direct the design and development of computational prediction models that can further enhance the utility of mechanistic and in silico data for chemical safety assessment. This concept was explored as part of a workshop on AOP-Informed Predictive Modeling Approaches for Regulatory Toxicology held September 24-25, 2015. Examples of AOP-informed model development and its application to the assessment of chemicals for skin sensitization and multiple modes of endocrine disruption are provided. The role of problem formulation, not only as a critical phase of risk assessment, but also as guide for both AOP and complementary model development is described. Finally, a proposal for actively engaging the modeling community in AOP-informed computational model development is made. The contents serve as a vision for how AOPs can be leveraged to facilitate development of computational prediction models needed to support the next generation of chemical safety assessment.

Dietary Iodine Sufficiency and Moderate Insufficiency in the Lactating Mother and Nursing Infant: A Computational Perspective

The Institute of Medicine recommends that lactating women ingest 290 μg iodide/d and a nursing infant, less than two years of age, 110 μg/d. The World Health Organization, United Nations Children’s Fund, and International Council for the Control of Iodine Deficiency Disorders recommend population maternal and infant urinary iodide concentrations ≥ 100 μg/L to ensure iodide sufficiency. For breast milk, researchers have proposed an iodide concentration range of 150–180 μg/L indicates iodide sufficiency for the mother and infant, however no national or international guidelines exist for breast milk iodine concentration. For the first time, a lactating woman and nursing infant biologically based model, from delivery to 90 days postpartum, was constructed to predict maternal and infant urinary iodide concentration, breast milk iodide concentration, the amount of iodide transferred in breast milk to the nursing infant each day and maternal and infant serum thyroid hormone kinetics. The maternal and infant models each consisted of three sub-models, iodide, thyroxine (T4), and triiodothyronine (T3). Using our model to simulate a maternal intake of 290 μg iodide/d, the average daily amount of iodide ingested by the nursing infant, after 4 days of life, gradually increased from 50 to 101 μg/day over 90 days postpartum. The predicted average lactating mother and infant urinary iodide concentrations were both in excess of 100 μg/L and the predicted average breast milk iodide concentration, 157 μg/L. The predicted serum thyroid hormones (T4, free T4 (fT4), and T3) in both the nursing infant and lactating mother were indicative of euthyroidism. The model was calibrated using serum thyroid hormone concentrations for lactating women from the United States and was successful in predicting serum T4 and fT4 levels (within a factor of two) for lactating women in other countries. T3 levels were adequately predicted. Infant serum thyroid hormone levels were adequately predicted for most data. For moderate iodide deficient conditions, where dietary iodide intake may range from 50 to 150 μg/d for the lactating mother, the model satisfactorily described the iodide measurements, although with some variation, in urine and breast milk. Predictions of serum thyroid hormones in moderately iodide deficient lactating women (50 μg/d) and nursing infants did not closely agree with mean reported serum thyroid hormone levels, however, predictions were usually within a factor of two. Excellent agreement between prediction and observation was obtained for a recent moderate iodide deficiency study in lactating women. Measurements included iodide levels in urine of infant and mother, iodide in breast milk, and serum thyroid hormone levels in infant and mother. A maternal iodide intake of 50 μg/d resulted in a predicted 29–32% reduction in serum T4 and fT4 in nursing infants, however the reduced serum levels of T4 and fT4 were within most of the published reference intervals for infant. This biologically based model is an important first step at integrating the rapid changes that occur in the thyroid system of the nursing newborn in order to predict adverse outcomes from exposure to thyroid acting chemicals, drugs, radioactive materials or iodine deficiency.

Quantitative global sensitivity analysis of a biologically based dose-response pregnancy model for the thyroid endocrine system

A deterministic biologically based dose-response model for the thyroidal system in a near-term pregnant woman and the fetus was recently developed to evaluate quantitatively thyroid hormone perturbations. The current work focuses on conducting a quantitative global sensitivity analysis on this complex model to identify and characterize the sources and contributions of uncertainties in the predicted model output. The workflow and methodologies suitable for computationally expensive models, such as the Morris screening method and Gaussian Emulation processes, were used for the implementation of the global sensitivity analysis. Sensitivity indices, such as main, total and interaction effects, were computed for a screened set of the total thyroidal system descriptive model input parameters. Furthermore, a narrower sub-set of the most influential parameters affecting the model output of maternal thyroid hormone levels were identified in addition to the characterization of their overall and pair-wise parameter interaction quotients. The characteristic trends of influence in model output for each of these individual model input parameters over their plausible ranges were elucidated using Gaussian Emulation processes. Through global sensitivity analysis we have gained a better understanding of the model behavior and performance beyond the domains of observation by the simultaneous variation in model inputs over their range of plausible uncertainties. The sensitivity analysis helped identify parameters that determine the driving mechanisms of the maternal and fetal iodide kinetics, thyroid function and their interactions, and contributed to an improved understanding of the system modeled. We have thus demonstrated the use and application of global sensitivity analysis for a biologically based dose-response model for sensitive life-stages such as pregnancy that provides richer information on the model and the thyroidal system modeled compared to local sensitivity analysis.