Computational modeling of serum-binding proteins and clearance in extrapolations across life stages and species for endocrine active compounds

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

Comparative estrogenicity of endogenous, environmental and dietary estrogens in pregnant women II: Total estrogenicity calculations accounting for competitive protein and receptor binding and potency

Evaluating the biological significance of human-relevant exposures to environmental estrogens involves assessing the individual and total estrogenicity of endogenous and exogenous estrogens found in serum, for example from biomonitoring studies. We developed a method for this assessment by integrating approaches for (i) measuring total hormone concentrations by mass spectrometry (Fleck et al., 2018), (ii) calculating hormone bioavailable concentrations in serum and, (iii) solving multiple equilibria between estrogenic ligands and receptors, and (iv) quantitatively describing key elements of estrogen potency. The approach was applied to endogenous (E1, E2, E3, E4), environmental (BPA), and dietary Genistein (GEN), Daidzein (DDZ) estrogens measured in the serum of thirty pregnant women. Fractional receptor occupancy (FRO) based estrogenicity was dominated by E1, E2 and E3 (ER-α, 94.4-99.2% (median: 97.3%), ER-β, 82.7-97.7% (median: 92.8%), as was the total response (TR), which included ligand specific differences in recruitment of co-activator proteins (RCA). The median FRO for BPA was at least five orders of magnitude lower than E1, E2 and E3, and three orders of magnitude lower than the fetal derived E4 and GEN and DDZ. BPA contributed less than 1/1000th of the normal daily variability in total serum estrogenicity in this cohort of pregnant women.

Evaluation and Optimization of Pharmacokinetic Models for in Vitro to in Vivo Extrapolation of Estrogenic Activity for Environmental Chemicals

BACKGROUND

To effectively incorporate in vitro data into regulatory use, confidence must be established in the quantitative extrapolation of in vitro activity to relevant end points in animals or humans.

OBJECTIVE

Our goal was to evaluate and optimize in vitro to in vivo extrapolation (IVIVE) approaches using in vitro estrogen receptor (ER) activity to predict estrogenic effects measured in rodent uterotrophic studies.

METHODS

We evaluated three pharmacokinetic (PK) models with varying complexities to extrapolate in vitro to in vivo dosimetry for a group of 29 ER agonists, using data from validated in vitro [U.S. Environmental Protection Agency (U.S. EPA) ToxCast™ ER model] and in vivo (uterotrophic) methods. In vitro activity values were adjusted using mass-balance equations to estimate intracellular exposure via an enrichment factor (EF), and steady-state model calculations were adjusted using fraction of unbound chemical in the plasma ([Formula: see text]) to approximate bioavailability. Accuracy of each model-adjustment combination was assessed by comparing model predictions with lowest effect levels (LELs) from guideline uterotrophic studies.

RESULTS

We found little difference in model predictive performance based on complexity or route-specific modifications. Simple adjustments, applied to account for in vitro intracellular exposure (EF) or chemical bioavailability ([Formula: see text]), resulted in significant improvements in the predictive performance of all models.

CONCLUSION

Computational IVIVE approaches accurately estimate chemical exposure levels that elicit positive responses in the rodent uterotrophic bioassay. The simplest model had the best overall performance for predicting both oral (PPK_EF) and injection (PPK_[Formula: see text]) LELs from guideline uterotrophic studies, is freely available, and can be parameterized entirely using freely available in silico tools. https://doi.org/10.1289/EHP1655.

Sex differences, endogenous sex-hormone hormones, sex-hormone binding globulin, and exogenous disruptors in diabetes and related metabolic outcomes

In assessing clinical and pathophysiological development of type 2 diabetes (T2D), the critical role of the sex steroids axis is underappreciated, particularly concerning the sex-specific relationships with many relevant cardiometabolic outcomes. In this issue of the Journal of Diabetes, we provide a comprehensive overview of these significant associations of germline variants in the genes governing the sex steroid pathways, plasma levels of steroid hormones, and sex hormone-binding globulin (SHBG) with T2D risk that have been observed in many clinical and high-quality large prospective cohorts of men and women across ethnic populations. Together, this body of evidence indicates that sex steroids and SHBG should be routinely incorporated into clinical characterization of T2D patients, particularly in screening prediabetic patients, such as those with metabolic syndrome, using plasma levels of SHBG. Given that several germline mutations in the SHBG gene have also been directly related to both plasma concentrations of SHBG and clinical manifestation of T2D, targeting signals in the sex steroid axis, particularly SHBG, may have significant utility in the prediction and treatment of T2D. Further, many of the environmental endocrine disrupting chemicals may exert their potential adverse effects on cardiometabolic outcomes via either estrogenic or androgenic signaling pathways, highlighting the importance of using the sex steroids and SHBG as important biochemical markers in both clinical and population studies in studying sex-specific mechanisms in the pathogenesis of T2D and its complications, as well as the need to equitably allocate resources in studying both men and women.

Integration of a plasma protein binding factor to the Chemical-Specific Adjustment Factor (CSAF) for facilitating the estimation of uncertainties in interspecies extrapolations when deriving health-based exposure limits for active pharmaceutical ingredients: Investigation of recent drug datasets

The objective was to challenge cross-species extrapolation factors with which to scale animal doses to human by any route for non-carcinogenic endpoints. The conventional hypothesis of the toxicokinetics (TK)-toxicodynamics (TD) relationship was equal toxicity at equal plasma level of the total drug moiety in each species, but this should also follow the free drug assumption, which states that only the unbound drug moiety in plasma may elicit a TD effect in tissue. Therefore, a protein binding factor (PBF) was combined with the Chemical-Specific Adjustment Factor (CSAF) (i.e., CSAF x PBF). The value of PBF of each drug was set equal to the ratio between human and animals of the unbound fraction in plasma (fu_p). Recent drug datasets were investigated. Our results indicate that any CSAF value would be increased or decreased while PBF deviates to the unity, and this required more attention. Accordingly, further testing indicated that the CSAF values set equal to basic allometric uncertainty factors according to the conventional hypothesis (dog∼2, monkey∼3.1, rat∼7, mouse∼12) would increase by including PBF for 30% of the drugs tested that showed a superior fu_p value in human compared to animals. However, default uncertainty factors in the range of 10-100 were less frequently exceeded. Overall, PBF could be combined with any other uncertainty factor to get reliable estimate of CSAF for each bound drug in deriving health-based exposure limits.

The virtual cell based assay: Current status and future perspectives

In order to replace the use of animals in toxicity testing, there is a need to predict in vivo toxic doses from concentrations that cause toxicological effects in relevant in vitro systems. The Virtual Cell Based Assay (VCBA) estimates time-dependent concentration of a test chemical in the cell and cell culture for a given in vitro system. The concentrations in the different compartments of the cell and test system are derived from ordinary differential equations, physicochemical parameters of the test chemical and properties of the cell line. The VCBA has been developed for a range of cell lines including BALB/c 3T3 cells, HepG2, HepaRG, lung A459 cells, and cardiomyocytes. The model can be used to design and refine in vitro experiments and extrapolate in vitro effective concentrations to in vivo doses that can be applied in risk assessment. In this paper, we first discuss potential applications of the VCBA: i) design of in vitro High Throughput Screening (HTS) experiments; ii) hazard identification (based on acute systemic toxicity); and iii) risk assessment. Further extension of the VCBA is discussed in the second part, exploring potential application to i) manufactured nanomaterials, ii) additional cell lines and endpoints, and considering iii) other opportunities.

•

VCBA as an alternative approach can be applied in the domain of nanotoxicology.

•

VCBA can support better testing strategies in acute toxicity.

•

Refinement of the VCBA taking into account biological oscillators could improve toxicity prediction.

•

Extensions of the VCBA can capture effects related to additional subcellular compartments.

Systems Biology to Support Nanomaterial Grouping

The assessment of potential health risks of engineered nanomaterials (ENMs) is a challenging task due to the high number and great variety of already existing and newly emerging ENMs. Reliable grouping or categorization of ENMs with respect to hazards could help to facilitate prioritization and decision making for regulatory purposes. The development of grouping criteria, however, requires a broad and comprehensive data basis. A promising platform addressing this challenge is the systems biology approach. The different areas of systems biology, most prominently transcriptomics, proteomics and metabolomics, each of which provide a wealth of data that can be used to reveal novel biomarkers and biological pathways involved in the mode-of-action of ENMs. Combining such data with classical toxicological data would enable a more comprehensive understanding and hence might lead to more powerful and reliable prediction models. Physico-chemical data provide crucial information on the ENMs and need to be integrated, too. Overall statistical analysis should reveal robust grouping and categorization criteria and may ultimately help to identify meaningful biomarkers and biological pathways that sufficiently characterize the corresponding ENM subgroups. This chapter aims to give an overview on the different systems biology technologies and their current applications in the field of nanotoxicology, as well as to identify the existing challenges.

Adverse effects of bisphenol A (BPA) on the dopamine system in two distinct cell models and corpus striatum of the Sprague-Dawley rat

The aim of this study was to examine the effects of bisphenol A (BPA) on the brain dopamine (DA) system utilizing both in vitro models (GH3 cells, a rat pituitary cell line, and SH-SY5Y cells, a human neuroblastoma cell line) and an animal model such as Sprague-Dawley (SD) rats. First, cellular DA uptake was measured 2 or 8 h following BPA exposure (0.1-400 μM) in SH-SY5Y cells, where a significant increase in DA uptake was noted. BPA exerted no marked effect on dopamine active transporter levels in GH3 cells exposed for 8 or 24 h. However, SH-SY5Y cells displayed an increase in dopamine transporter (DAT) levels following 24 h of exposure to BPA. In contrast to DAT levels, BPA exposure produced no marked effect on DA D1 receptor levels in SH-SY5Y cells, yet a significant decrease in GH3 cells following both 8- and 24-h exposure periods was noted, suggesting that BPA exerts differential effects dependent upon cell type. BPA produced no significant effects on prolactin levels at 2 h, but a marked fall occurred at 24 h of exposure in GH3 cells. Finally, to examine the influence of dietary developmental exposure to BPA on brain DA levels in F1 offspring, SD rats were exposed to BPA (0.5-20 mg/kg) through maternal transfer and/or diet and striatal DA levels were measured on postnatal day (PND) 60 using high-performance liquid chromatography (HPLC). Data demonstrated that chronic exposure to BPA did not significantly alter striatal DA levels in the SD rat.

Comparative iron oxide nanoparticle cellular dosimetry and response in mice by the inhalation and liquid cell culture exposure routes

Background

Toxicity testing the rapidly growing number of nanomaterials requires large scale use of in vitro systems under the presumption that these systems are sufficiently predictive or descriptive of responses in in vivo systems for effective use in hazard ranking. We hypothesized that improved relationships between in vitro and in vivo models of experimental toxicology for nanomaterials would result from placing response data in vitro and in vivo on the same dose scale, the amount of material associated with cells.

Methods

Balb/c mice were exposed nose-only to an aerosol (68.6 nm CMD, 19.9 mg/m³, 4 hours) generated from of 12.8 nm superparamagnetic iron oxide particles (SPIO). Target cell doses were calculated, histological evaluations conducted, and biomarkers of response were identified by global transcriptomics. Representative murine epithelial and macrophage cell types were exposed in vitro to the same material in liquid suspension for four hours and levels of nanoparticle regulated cytokine transcripts identified in vivo were quantified as a function of measured nanoparticle cellular dose.

Results

Target tissue doses of 0.009-0.4 μg SPIO/cm² in lung led to an inflammatory response in the alveolar region characterized by interstitial inflammation and macrophage infiltration. In vitro, higher target tissue doses of ~1.2-4 μg SPIO/ cm² of cells were required to induce transcriptional regulation of markers of inflammation, CXCL2 & CCL3, in C10 lung epithelial cells. Estimated in vivo macrophage SPIO nanoparticle doses ranged from 1-100 pg/cell, and induction of inflammatory markers was observed in vitro in macrophages at doses of 8-35 pg/cell.

Conclusions

Application of target tissue dosimetry revealed good correspondence between target cell doses triggering inflammatory processes in vitro and in vivo in the alveolar macrophage population, but not in the epithelial cells of the alveolar region. These findings demonstrate the potential for target tissue dosimetry to enable the more quantitative comparison of in vitro and in vivo systems and advance their use for hazard assessment and extrapolation to humans. The mildly inflammogentic cellular doses experienced by mice were similar to those calculated for humans exposed to the same material at the existing permissible exposure limit of 10 mg/m³ iron oxide (as Fe).

Electronic supplementary material

The online version of this article (doi:10.1186/s12989-014-0046-4) contains supplementary material, which is available to authorized users.

Comparison of life-stage-dependent internal dosimetry for bisphenol A, ethinyl estradiol, a reference estrogen, and endogenous estradiol to test an estrogenic mode of action in Sprague Dawley rats

Bisphenol A (BPA) was administered by gavage (2.5-300,000 μg/kg body weight (bw)/day) to pregnant Sprague Dawley dams, newborn pups, and continuing into adulthood. Aglycone (i.e., unconjugated and active) and conjugated (i.e., inactive) BPA were evaluated by liquid chromatography electrospray tandem mass spectrometry (LC-ES/MS/MS) in serum to better interpret toxicological endpoints measured in the study. Ethinyl estradiol (EE2, 0.5 and 5 μg/kg bw/day) and the endogenous hormones, 17β-estradiol (E2) and testosterone, were similarly evaluated. Mean BPA aglycone levels in vehicle and naïve control rat serum (0.02-0.5 ng/ml) indicated sample processing artifact, consistent with literature reports of a propensity for postexposure blood contamination by BPA. Direct measurements of BPA-glucuronide in vehicle and naïve control serum (2-10nM) indicated unintentional exposure and metabolism at levels similar to those produced by 2.5 μg/kg bw/day BPA (7-10nM), despite careful attention to potential BPA inputs (diet, drinking water, vehicle, cages, bedding, and dust) and rigorous dosing solution certification and delivery. The source of this exposure could not be identified, but interpretation of the toxicological effects, observed only at the highest BPA doses, was not compromised. Internal exposures to BPA and EE2 aglycones were highest in young rats. When maximal serum concentrations from the two highest BPA doses and both EE2 doses were compared with concurrent levels of endogenous E2, the ERα binding equivalents were similar to or above those of endogenous E2 in male and female rats of all ages tested. Such evaluations of estrogenic internal dosimetry and comprehensive evaluation of contamination impact should aid in extrapolating risks from human BPA exposures.

A computational model of the hypothalamic: pituitary: gonadal axis in female fathead minnows (Pimephales promelas) exposed to 17α-ethynylestradiol and 17β-trenbolone

Background

Endocrine disrupting chemicals (e.g., estrogens, androgens and their mimics) are known to affect reproduction in fish. 17α-ethynylestradiol is a synthetic estrogen used in birth control pills. 17β-trenbolone is a relatively stable metabolite of trenbolone acetate, a synthetic androgen used as a growth promoter in livestock. Both 17α-ethynylestradiol and 17β-trenbolone have been found in the aquatic environment and affect fish reproduction. In this study, we developed a physiologically-based computational model for female fathead minnows (FHM, Pimephales promelas), a small fish species used in ecotoxicology, to simulate how estrogens (i.e., 17α-ethynylestradiol) or androgens (i.e., 17β-trenbolone) affect reproductive endpoints such as plasma concentrations of steroid hormones (e.g., 17β-estradiol and testosterone) and vitellogenin (a precursor to egg yolk proteins).

Results

Using Markov Chain Monte Carlo simulations, the model was calibrated with data from unexposed, 17α-ethynylestradiol-exposed, and 17β-trenbolone-exposed FHMs. Four Markov chains were simulated, and the chains for each calibrated model parameter (26 in total) converged within 20,000 iterations. With the converged parameter values, we evaluated the model's predictive ability by simulating a variety of independent experimental data. The model predictions agreed with the experimental data well.

Conclusions

The physiologically-based computational model represents the hypothalamic-pituitary-gonadal axis in adult female FHM robustly. The model is useful to estimate how estrogens (e.g., 17α-ethynylestradiol) or androgens (e.g., 17β-trenbolone) affect plasma concentrations of 17β-estradiol, testosterone and vitellogenin, which are important determinants of fecundity in fish.

Differential accumulation levels in the brain of rats exposed to the endocrine disruptor 4-tert-octylphenol (OP)

Octylphenol (OP) is an endocrine-disrupting chemical that accumulates in various organs. It has also been shown to exert noxious effects on the central nervous system. In the present study, we measured in Sprague-Dawley rats the degree of OP accumulation in different areas of the brain and investigated the effect of OP in pain modulation. Two groups of male Sprague-Dawley rats were treated for 20 days with 50mg/kg BW/day of OP (group 1) or vehicle (group 2). At the end of the treatment, the formalin test was performed to evaluate the effect of OP exposure on pain. Soon after, rats were sacrificed, and the accumulation of OP in the cerebral cortex, hippocampus, hypothalamus, cerebellum, thalamus, striatum, mesencephalus and ventral hindbrain was measured by HPLC analysis. The results showed a greater accumulation of OP in the cerebral cortex compared to all the other areas; there was also more accumulation in the cerebellum compared to the mesencephalus and thalamus. No accumulation was found in the striatum. These results suggest that there is a preferential accumulation of OP in different areas of the brain with consequences to neural behaviour. On the contrary, experiments on facial grooming did not show significant effects of OP on pain.

ISDD: A computational model of particle sedimentation, diffusion and target cell dosimetry for in vitro toxicity studies

BACKGROUND

The difficulty of directly measuring cellular dose is a significant obstacle to application of target tissue dosimetry for nanoparticle and microparticle toxicity assessment, particularly for in vitro systems. As a consequence, the target tissue paradigm for dosimetry and hazard assessment of nanoparticles has largely been ignored in favor of using metrics of exposure (e.g. μg particle/mL culture medium, particle surface area/mL, particle number/mL). We have developed a computational model of solution particokinetics (sedimentation, diffusion) and dosimetry for non-interacting spherical particles and their agglomerates in monolayer cell culture systems. Particle transport to cells is calculated by simultaneous solution of Stokes Law (sedimentation) and the Stokes-Einstein equation (diffusion).

RESULTS

The In vitro Sedimentation, Diffusion and Dosimetry model (ISDD) was tested against measured transport rates or cellular doses for multiple sizes of polystyrene spheres (20-1100 nm), 35 nm amorphous silica, and large agglomerates of 30 nm iron oxide particles. Overall, without adjusting any parameters, model predicted cellular doses were in close agreement with the experimental data, differing from as little as 5% to as much as three-fold, but in most cases approximately two-fold, within the limits of the accuracy of the measurement systems. Applying the model, we generalize the effects of particle size, particle density, agglomeration state and agglomerate characteristics on target cell dosimetry in vitro.

CONCLUSIONS

Our results confirm our hypothesis that for liquid-based in vitro systems, the dose-rates and target cell doses for all particles are not equal; they can vary significantly, in direct contrast to the assumption of dose-equivalency implicit in the use of mass-based media concentrations as metrics of exposure for dose-response assessment. The difference between equivalent nominal media concentration exposures on a μg/mL basis and target cell doses on a particle surface area or number basis can be as high as three to six orders of magnitude. As a consequence, in vitro hazard assessments utilizing mass-based exposure metrics have inherently high errors where particle number or surface areas target cells doses are believed to drive response. The gold standard for particle dosimetry for in vitro nanotoxicology studies should be direct experimental measurement of the cellular content of the studied particle. However, where such measurements are impractical, unfeasible, and before such measurements become common, particle dosimetry models such as ISDD provide a valuable, immediately useful alternative, and eventually, an adjunct to such measurements.

Physiologically based modeling of p-tert-octylphenol kinetics following intravenous, oral or subcutaneous exposure in male and female Sprague-Dawley rats

The objective of this study was to develop a physiologically based pharmacokinetic (PBPK) model for p-tert-octylphenol (OP) for understanding the qualitative and quantitative determinants of its kinetics in Sprague-Dawley rats. Compartments of the PBPK model included the liver, richly perfused tissues, poorly perfused tissues, reproductive tissues, adipose tissue and subcutaneous space, in which OP uptake was described as a blood flow- or a membrane diffusion-limited process. The PBPK model successfully simulated previously published data on blood and tissue OP concentrations in Sprague-Dawley rats following oral, intravenous (i.v.) or subcutaneous (s.c.) routes. The model predicted that OP concentrations would reach 6.8, 13.8 and 27.9 ng ml(-1) (male) and 7.2, 14.7 and 31.4 ng ml(-1) (female), 4 h after a single i.v. dose of 2, 4 and 8 mg kg(-1), respectively. The model also predicted that OP concentrations would reach 53.3, 134.8 and 271.2 ng ml(-1) (male) and 87.4, 221.4 and 449.7 ng ml(-1) (female) 4 h after a single oral dose (50, 125 and 250 mg kg(-1)) and that, 4 h after a single s.c. dose (125 mg kg(-1)), OP concentrations would reach 111.3 ng ml(-1) (male) and 121.6 ng ml(-1). A marked sex difference was seen in blood and tissue OP concentrations. This was reflected in the model by a gender-specific maximal velocity of metabolism (V(max)) that was higher (1.77 x) in male than in female rats. Further studies are required to elucidate the mechanism underlying the gender differences and to evaluate whether that is also observed in humans.

Comparing single and repeated dosimetry data for perfluorooctane sulfonate in rats

Perfluorooctane sulfonate (PFOS) is a member of a class of perfluorinated chemicals used in a variety of consumer and industrial applications because of their oleophobic and hydrophobic properties. It has been shown to cause toxicity in adult and developing laboratory animals. Because PFOS has also been shown to be widely distributed throughout the environment, there have been concerns about its potential health risk to humans. Limited pharmacokinetic data for PFOS are available in rodents and humans, while epidemiological studies of workers and extensive toxicity studies in rodents have been performed. The existing pharmacokinetic and toxicity database in rodents can be useful in the cross-species extrapolations needed to evaluate and interpret internal dosimetry in humans. A mathematical model that describes the disposition of PFOS in adult rats following intravenous, oral, and chronic dietary exposures was developed to gain a better understanding of the pharmacokinetics of PFOS and to determine whether single-dose kinetics are predictive of repeated-dose kinetics. In order to characterize existing time-course data, time-dependent and concentration-dependent changes in the pharmacokinetic parameters for urinary and biliary clearance and liver distribution were needed. Whether these time-dependent changes represent inconsistencies across experiments, effects of aging in the rats, or chemically induced changes in pharmacokinetics remains to be determined.

In vitro molecular mechanisms of bisphenol A action

Bisphenol A (BPA, 2,2-bis (4-hydroxyphenyl) propane; CAS# 80-05-7) is a chemical used primarily in the manufacture of polycarbonate plastic, epoxy resins and as a non-polymer additive to other plastics. Recent evidence has demonstrated that human and wildlife populations are exposed to levels of BPA which cause adverse reproductive and developmental effects in a number of different wildlife species and laboratory animal models. However, there are major uncertainties surrounding the spectrum of BPA's mechanisms of action, the tissue-specific impacts of exposures, and the critical windows of susceptibility during which target tissues are sensitive to BPA exposures. As a foundation to address some of those uncertainties, this review was prepared by the "In vitro" expert sub-panel assembled during the "Bisphenol A: An Examination of the Relevance of Ecological, In vitro and Laboratory Animal Studies for Assessing Risks to Human Health" workshop held in Chapel Hill, NC, Nov 28-29, 2006. The specific charge of this expert panel was to review and assess the strength of the published literature pertaining to the mechanisms of BPA action. The resulting document is a detailed review of published studies that have focused on the mechanistic basis of BPA action in diverse experimental models and an assessment of the strength of the evidence regarding the published BPA research.

Particokinetics In Vitro: Dosimetry Considerations for In Vitro Nanoparticle Toxicity Assessments

The rapid growth in the use of in vitro methods for nanoparticle toxicity assessment has proceeded with limited consideration of the unique kinetics of these materials in solution. Particles in general and nanoparticles specifically, diffuse, settle, and agglomerate in cell culture media as a function of systemic and particle properties: media density and viscosity and particle size, shape, charge and density, for example. Cellular dose then is also a function of these factors as they determine the rate of transport of nanoparticles to cells in culture. Here we develop and apply the principles of dosimetry in vitro and outline an approach for simulation of nanoparticle particokinetics in cell culture systems. We illustrate that where equal mass concentrations (mug/ml) imply equal doses for dissimilar materials, the corresponding particle number or surface area concentration doses differ by orders of magnitude. More importantly, when rates of diffusional and gravitational particle delivery are accounted for, trends and magnitude of the cellular dose as a function of particle size and density differ significantly from those implied by "concentration" doses. For example, 15-nm silver nanoparticles appear approximately 4000 times more potent than micron-sized cadmium oxide particles on a cm(2)/ml media basis, but are only approximately 50 times more potent when differences in delivery to adherent cells are considered. We conclude that simple surrogates of dose can cause significant misinterpretation of response and uptake data for nanoparticles in vitro. Incorporating particokinetics and principles of dosimetry would significantly improve the basis for nanoparticle toxicity assessment, increasing the predictive power and scalability of such assays.

The acquisition and application of absorption, distribution, metabolism, and excretion (ADME) data in agricultural chemical safety assessments

A proposal has been developed by the Agricultural Chemical Safety Assessment (ACSA) Technical Committee of the ILSI Health and Environmental Sciences Institute (HESI) for an improved approach to assessing the safety of crop protection chemicals. The goal is to ensure that studies are scientifically appropriate and necessary without being redundant, and that tests emphasize toxicological endpoints and exposure durations that are relevant for risk assessment. Incorporation of pharmacokinetic studies describing absorption, distribution, metabolism, and excretion is an essential tool for improving the design and interpretation of toxicity studies and their application for safety assessment. A tiered approach is described in which basic pharmacokinetic studies, similar to those for pharmaceuticals, are conducted for regulatory submission. Subsequent tiers provide additional information in an iterative manner, depending on pharmacokinetic properties, toxicity study results, and the intended uses of the compound.

Measurement of the free concentration of octylphenol in biological samples with negligible depletion-solid phase microextraction (nd-SPME): Analysis of matrix effects

A negligible depletion-solid phase microextraction (nd-SPME) method is presented to measure free concentrations of octylphenol in biological samples. Potential confounding factors, such as matrix effects, are studied as well. Fouling of the fibre appears to occur, but it does not seem to reduce or enhance the measured uptake of octylphenol. In the setup applied here, without any agitation, it has also been found that there is a large effect of protein presence on the kinetics of octylphenol uptake. In addition, an apparent affinity constant of octylphenol for bovine serum albumin was determined.

Physiologically-based pharmacokinetic modeling of genistein in rats, Part I: Model development

Genistein is a phytoestrogen-a plant-derived compound that binds to and activates the estrogen receptor-occurring at high levels in soy beans and food products, leading to widespread human exposure. The numerous scientific publications available describing genistein's dosimetry, mechanisms of action, and identified or putative health effects in both experimental animals and humans make it ideal for examination as an example of endocrine-active compound (EAC). We developed a physiologically-based pharmacokinetic (PBPK) model to quantify the internal, target-tissue dosimetry of genistein in adult rats. Complexities of the model include enterohepatic circulation, binding of both genistein and its conjugates to plasma proteins, and the multiple compartments used to describe transport through the bile duct and gastrointestinal tract. Other aspects of the model are simple perfusion-limited transport to the tissue groups and first-order rates of metabolism, uptake, and excretion. We describe here the model structure and initial calibration of the model by fitting to a large data set for Wistar rats. The model structure can be readily extrapolated to describe genistein dosimetry in humans or modified to describe the dosimetry of other phytoestrogens and phenolic EACs. The model does a fair job of capturing the pharmacokinetics. Although it does not describe the interindividual variability and we have not identified a single set of parameters that provide a good fit to the data for both oral and intravenous exposures, we believe it provides a good initial attempt at PBPK modeling for genistein, which can serve as a template for other phytoestrogens and in the design of future experiments and research that can be used to fill data gaps and better estimate model parameters.

Computational pharmacokinetics during developmental windows of susceptibility

Computational modeling has an increasing role in analyses of biological effects, including how the body handles chemicals (i.e., pharmacokinetics or toxicokinetics) and how the body responds to chemicals (i.e., pharmacodynamics or toxicodynamics). Pharmacokinetic models increasingly describe not just adult humans and animals, but also changes with age and life stage (e.g., pregnancy and fetal exposures, lactational exposures, and childhood growth). Physiologically based pharmacokinetic models provide an important route to estimate the potential changes in internal dose that may occur throughout the life cycle. These models require inputs describing changes in physiology, metabolism, and exposure with age and life stage. A particular challenge exists when the "equivalent" developmental period in the rodents and humans differs (e.g., early postnatal in rats and in utero in humans) such that the "equivalent" window of susceptibility to toxic effects of the chemical may involve substantially different exposures (e.g., lactational versus placental transfer). Pharmacodynamic modeling could similarly address changes with age, but few such models currently exist. The growth of systems biology is anticipated to change this over the coming decade.

Evaluation of Oral and Intravenous Route Pharmacokinetics, Plasma Protein Binding, and Uterine Tissue Dose Metrics of Bisphenol A: A Physiologically Based Pharmacokinetic Approach

Bisphenol A (BPA) is a weakly estrogenic monomer used in the production of polycarbonate plastic and epoxy resins, both of which are used in food contact and other applications. A physiologically based pharmacokinetic (PBPK) model of BPA pharmacokinetics in rats and humans was developed to provide a physiological context in which the processes controlling BPA pharmacokinetics (e.g., plasma protein binding, enterohepatic recirculation of the glucuronide [BPAG]) could be incorporated. A uterine tissue compartment was included to allow the correlation of simulated estrogen receptor (ER) binding of BPA with increases in uterine wet weight (UWW) in rats. Intravenous- and oral-route blood kinetics of BPA in rats and oral-route plasma and urinary elimination kinetics in humans were well described by the model. Simulations of rat oral-route BPAG pharmacokinetics were less exact, most likely the result of oversimplification of the GI tract compartment. Comparison of metabolic clearance rates derived from fitting rat i.v. and oral-route data implied that intestinal glucuronidation of BPA is significant. In rats, but not humans, terminal elimination rates were strongly influenced by enterohepatic recirculation. In the absence of BPA binding to plasma proteins, simulations showed high ER occupancy at doses without uterine effects. Restricting free BPA to the measured unbound amount demonstrated the importance of including plasma binding in BPA kinetic models: the modeled relationship between ER occupancy and UWW increases was consistent with expectations for a receptor-mediated response with low ER occupancy at doses with no response and increasing occupancy with larger increases in UWW.

Framework for evaluation of physiologically-based pharmacokinetic models for use in safety or risk assessment

Proposed applications of increasingly sophisticated biologically-based computational models, such as physiologically-based pharmacokinetic models, raise the issue of how to evaluate whether the models are adequate for proposed uses, including safety or risk assessment. A six-step process for model evaluation is described. It relies on multidisciplinary expertise to address the biological, toxicological, mathematical, statistical, and risk assessment aspects of the modeling and its application. The first step is to have a clear definition of the purpose(s) of the model in the particular assessment; this provides critical perspectives on all subsequent steps. The second step is to evaluate the biological characterization described by the model structure based on the intended uses of the model and available information on the compound being modeled or related compounds. The next two steps review the mathematical equations used to describe the biology and their implementation in an appropriate computer program. At this point, the values selected for the model parameters (i.e., model calibration) must be evaluated. Thus, the fifth step is a combination of evaluating the model parameterization and calibration against data and evaluating the uncertainty in the model outputs. The final step is to evaluate specialized analyses that were done using the model, such as modeling of population distributions of parameters leading to population estimates for model outcomes or inclusion of early pharmacodynamic events. The process also helps to define the kinds of documentation that would be needed for a model to facilitate its evaluation and implementation.