Validation of a prediction model for estimating serum concentrations of chemicals which are equivalent to toxic concentrations in vitro

Paired Liver:Plasma PFAS Concentration Ratios from Adolescents in the Teen-LABS Study and Derivation of Empirical and Mass Balance Models to Predict and Explain Liver PFAS Accumulation

Animal studies have pointed at the liver as a hotspot for per- and polyfluoroalkyl substances (PFAS) accumulation and toxicity; however, these findings have not been replicated in human populations. We measured concentrations of seven PFAS in matched liver and plasma samples collected at the time of bariatric surgery from 64 adolescents in the Teen-Longitudinal Assessment of Bariatric Surgery (Teen-LABS) study. Liver:plasma concentration ratios were perfectly explained (r2 > 0.99) in a multilinear regression (MLR) model based on toxicokinetic (TK) descriptors consisting of binding to tissue constituents and membrane permeabilities. Of the seven matched plasma and liver PFAS concentrations compared in this study, the liver:plasma concentration ratio of perfluoroheptanoic acid (PFHpA) was considerably higher than the liver:plasma concentration ratio of other PFAS congeners. Comparing the MLR model with an equilibrium mass balance model (MBM) suggested that complex kinetic transport processes are driving the unexpectedly high liver:plasma concentration ratio of PFHpA. Intratissue MBM modeling pointed to membrane lipids as the tissue constituents that drive the liver accumulation of long-chain, hydrophobic PFAS, whereas albumin binding of hydrophobic PFAS dominated PFAS distribution in plasma. The liver:plasma concentration data set, empirical MLR model, and mechanistic MBM modeling allow the prediction of liver from plasma concentrations measured in human cohort studies. Our study demonstrates that combining biomonitoring data with mechanistic modeling can identify underlying mechanisms of internal distribution and specific target organ toxicity of PFAS in humans.

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.

Update and Evaluation of a High-Throughput In Vitro Mass Balance Distribution Model: IV-MBM EQP v2.0

This study demonstrates the utility of an updated mass balance model for predicting the distribution of organic chemicals in in vitro test systems (IV-MBM EQP v2.0) and evaluates its performance with empirical data. The IV-MBM EQP v2.0 tool was parameterized and applied to four independent data sets with measured ratios of bulk medium or freely-dissolved to initial nominal concentrations (e.g., C24/C0 where C24 is the measured concentration after 24 h of exposure and C0 is the initial nominal concentration). Model performance varied depending on the data set, chemical properties (e.g., "volatiles" vs. "non-volatiles", neutral vs. ionizable organics), and model assumptions but overall is deemed acceptable. For example, the r² was greater than 0.8 and the mean absolute error (MAE) in the predictions was less than a factor of two for most neutral organics included. Model performance was not as good for the ionizable organic chemicals included but the r² was still greater than 0.7 and the MAE less than a factor of three. The IV-MBM EQP v2.0 model was subsequently applied to several hundred chemicals on Canada's Domestic Substances List (DSL) with nominal effects data (AC50s) reported for two in vitro assays. We report the frequency of chemicals with AC50s corresponding to predicted cell membrane concentrations in the baseline toxicity range (i.e., >20-60 mM) and tabulate the number of chemicals with "volatility issues" (majority of chemical in headspace) and "solubility issues" (freely-dissolved concentration greater than water solubility after distribution). In addition, the predicted "equivalent EQP blood concentrations" (i.e., blood concentration at equilibrium with predicted cellular concentration) were compared to the AC50s as a function of hydrophobicity (log octanol-water partition or distribution ratio). The predicted equivalent EQP blood concentrations exceed the AC50 by up to a factor of 100 depending on hydrophobicity and assay conditions. The implications of using AC50s as direct surrogates for human blood concentrations when estimating the oral equivalent doses using a toxicokinetic model (i.e., reverse dosimetry) are then briefly discussed.

Defining in vivo dose-response curves for kidney DNA adduct formation of aristolochic acid I in rat, mouse and human by an in vitro and physiologically based kinetic modeling approach

Aristolochic acid I (AAI) is a well-known genotoxic kidney carcinogen. Metabolic conversion of AAI into the DNA-reactive aristolactam-nitrenium ion is involved in the mode of action of tumor formation. This study aims to predict in vivo AAI-DNA adduct formation in the kidney of rat, mouse and human by translating the in vitro concentration-response curves for AAI-DNA adduct formation to the in vivo situation using physiologically based kinetic (PBK) modeling-based reverse dosimetry. DNA adduct formation in kidney proximal tubular LLC-PK1 cells exposed to AAI was quantified by liquid chromatography-electrospray ionization-tandem mass spectrometry. Subsequently, the in vitro concentration-response curves were converted to predicted in vivo dose-response curves in rat, mouse and human kidney using PBK models. Results obtained revealed a dose-dependent increase in AAI-DNA adduct formation in the rat, mouse and human kidney and the predicted DNA adduct levels were generally within an order of magnitude compared with values reported in the literature. It is concluded that the combined in vitro PBK modeling approach provides a novel way to define in vivo dose-response curves for kidney DNA adduct formation in rat, mouse and human and contributes to the reduction, refinement and replacement of animal testing.

Partition coefficients of four perfluoroalkyl acid alternatives between bovine serum albumin (BSA) and water in comparison to ten classical perfluoroalkyl acids

Perfluoroalkyl acids (PFAAs) are persistent, ubiquitous environmental contaminants and their long-chain representatives are bioaccumulative. The phase-out of these compounds (e.g. PFOA and PFOS) shifted the production to alternatives. However, little is known about the bioaccumulative behaviour of the alternatives, which are still highly fluorinated. PFAAs are predominantly detected in blood, where they bind to the transport protein serum albumin. This sorption can be described by the albumin/water partition coefficient. It is unclear whether the partition coefficients of the alternatives are lower than or in the same range as those of classical PFAAs. We determined albumin/water partition coefficients for seven perfluoroalkyl carboxylates, three perfluoroalkane sulfonates and four alternatives by dialysis experiments in a physiologically representative system. Quantification was done by LC-MS/MS and a mass balance approach. Logarithmic albumin/water partition coefficients for PFAAs range from 2.8 to 4.8 [L_water kg_albumin^-1] and increase with increasing chain length. Perfluorinated sulfonates sorb more strongly than their carboxylate counterparts. The albumin/water partition coefficients for the alternatives (HFPO-DA, DONA, 9Cl-PF3ONS and PFECHS) are in the same range as for classical PFAAs. Structural modifications such as the introduction of ether groups into the chain do not reduce sorption to albumin, whereas the chlorine atom in 9Cl-PF3ONS seems to even increase the sorption to albumin. We further investigated whether the sorption strength could be affected in the presence of medium- or long-chain fatty acids. Binding competition between medium-chain fatty acids and PFAAs appeared to be possible. However, the presence of physiologically more relevant long-chain fatty acids should not alter the albumin/water partition coefficients of PFAAs.

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.

The CompTox Chemistry Dashboard: a community data resource for environmental chemistry

Despite an abundance of online databases providing access to chemical data, there is increasing demand for high-quality, structure-curated, open data to meet the various needs of the environmental sciences and computational toxicology communities. The U.S. Environmental Protection Agency's (EPA) web-based CompTox Chemistry Dashboard is addressing these needs by integrating diverse types of relevant domain data through a cheminformatics layer, built upon a database of curated substances linked to chemical structures. These data include physicochemical, environmental fate and transport, exposure, usage, in vivo toxicity, and in vitro bioassay data, surfaced through an integration hub with link-outs to additional EPA data and public domain online resources. Batch searching allows for direct chemical identifier (ID) mapping and downloading of multiple data streams in several different formats. This facilitates fast access to available structure, property, toxicity, and bioassay data for collections of chemicals (hundreds to thousands at a time). Advanced search capabilities are available to support, for example, non-targeted analysis and identification of chemicals using mass spectrometry. The contents of the chemistry database, presently containing ~ 760,000 substances, are available as public domain data for download. The chemistry content underpinning the Dashboard has been aggregated over the past 15 years by both manual and auto-curation techniques within EPA's DSSTox project. DSSTox chemical content is subject to strict quality controls to enforce consistency among chemical substance-structure identifiers, as well as list curation review to ensure accurate linkages of DSSTox substances to chemical lists and associated data. The Dashboard, publicly launched in April 2016, has expanded considerably in content and user traffic over the past year. It is continuously evolving with the growth of DSSTox into high-interest or data-rich domains of interest to EPA, such as chemicals on the Toxic Substances Control Act listing, while providing the user community with a flexible and dynamic web-based platform for integration, processing, visualization and delivery of data and resources. The Dashboard provides support for a broad array of research and regulatory programs across the worldwide community of toxicologists and environmental scientists.

Antioxidant activity of rosmarinic acid and its principal metabolites in chemical and cellular systems: Importance of physico-chemical characteristics

Persistent accumulation of reactive oxygen species causes cellular oxidative stress which contributes strongly towards the induction and progression of various diseases. Therapeutic focus has therefore shifted towards the use of antioxidants, with recent interest in those of plant origin. In the current study, rosmarinic acid (RA) and its key metabolites were evaluated in non-cellular and cellular antioxidant assays, using quercetin (Q) as a positive control. The non-cellular assay was performed as scavenging of DPPH radical, whilst the cellular assay was performed as protection from an oxidant stress. Radical-scavenging activity of RA and two of its primary metabolites, CA and DHPLA, were comparable to that of Q, whilst FA was of lower potency and m-CoA was inactive. In the cellular assay, RA and CA were markedly less potent than Q, with DHPLA, FA and m-CoA being inactive, this being true in short-term (5-h) or long-term (20-h) exposure conditions. However, antioxidant potency of Q and methyl rosmarinate, a non-ionisable ester of RA, was similar in the non-cellular and short-term cellular assays. It is proposed that marked ionisation of organic acids such as RA and its metabolites at physiological pH greatly limits their intracellular accumulation, and so attenuates intrinsic antioxidant ability demonstrated in the non-cellular assay. This study demonstrates some of the factors that prevent well-known phytochemicals from progressing further along the drug discovery chain.

Equilibrium Sorption of Structurally Diverse Organic Ions to Bovine Serum Albumin

Reliable partitioning data are essential for assessing the bioaccumulation potential and the toxicity of chemicals. In contrast to neutral organic chemicals, the partitioning behavior of ionogenic organic chemicals (IOCs) is still a black box for environmental scientists. Partitioning to serum albumin, the major protein in blood plasma, strongly influences the freely dissolved concentration of many chemicals (including IOCs), which affects their transport and distribution in the body. Because consistent data sets for partitioning of IOCs are rarely available, bovine serum albumin-water partition coefficients (KBSA/w) were measured in this study for 45 anionic and 4 cationic organic chemicals, including various substituted benzoic and naphthoic acids, sulfonates and several pesticides and pharmaceuticals. The results of this study suggest that binding to BSA is substantially influenced by the three-dimensional structure of the chemicals and the position of substitutions on the sorbing molecules. For example, we found a difference of >1.5 log units between isomeric chemicals such as 3,4-dichlorobenzoic acid and 2,6-dichlorobenzoic acid, and 1-naphthoic acid and 2-naphthoic acid. Conventional modeling approaches (e.g., based on octanol-water partition coefficients) poorly predict log KBSA/w of organic ions (R(2) ≤ 0.5), partially because they do not capture the observed steric effects. Hence, alternative modeling strategies will be required for accurate prediction of serum albumin-water partition coefficients of organic ions.

The power of energy-resolved tandem mass spectrometry experiments for resolution of isomers: the case of drug plasma stability investigation of multidrug resistance inhibitors

RATIONALE

A series of drug plasma stability experiments were carried out to evaluate the bioavailability of three multidrug resistance inhibitors. The studied compounds are positional isomers; therefore, a chromatographic separation or taking advantage of specific collisionally activated decomposition pathways, obtained by tandem mass spectrometry (MS/MS) experiments, is necessary in order to resolve them.

METHODS

A method was developed for quantitative determination of the analytes in plasma using liquid chromatography (LC) coupled with a triple quadrupole mass spectrometer operating in MS/MS mode. Different collisional approaches were employed based on the potentiality of a triple quadrupole system. Aside from the classical product ion spectroscopy, energy-resolved MS/MS experiments and a post-processing mathematical algorithm tool (LEDA) were used to distinguish among different kinds of inhibitors present in the sample batch.

RESULTS

The developed LC/MS/MS method showed precision between 1.8-7.9%, accuracy ranging from 92.8 to 99.9% and limit of detection (LOD) values in the range 1.0-1.4 ng mL(-1) for all the analytes. The evaluation of matrix effects demonstrated that the sample preparation procedure did not affect the ionization efficiency or recovery (matrix effects and recovery larger than 88%). Finally, the LEDA tool was able to differentiate among the isomers, ensuring their proper monitoring.

CONCLUSIONS

The proposed LC/MS/MS method was suitable for evaluating the stability of the analytes in plasma samples, although small concentration variations occurred. Furthermore, the investigation on the energetics of fragmentation pathways allowed the better product ions and optimal abundance ratios to be selected for LEDA application into a multi-component analysis. Copyright © 2015 John Wiley & Sons, Ltd.

Application of mass balance models and the chemical activity concept to facilitate the use of in vitro toxicity data for risk assessment

Practical, financial, and ethical considerations related to conducting extensive animal testing have resulted in various initiatives to promote and expand the use of in vitro testing data for chemical evaluations. Nominal concentrations in the aqueous phase corresponding to an effect (or biological activity) are commonly reported and used to characterize toxicity (or biological response). However, the true concentration in the aqueous phase can be substantially different from the nominal. To support in vitro test design and aid the interpretation of in vitro toxicity data, we developed a mass balance model that can be parametrized and applied to represent typical in vitro test systems. The model calculates the mass distribution, freely dissolved concentrations, and cell/tissue concentrations corresponding to the initial nominal concentration and experimental conditions specified by the user. Chemical activity, a metric which can be used to assess the potential for baseline toxicity to occur, is also calculated. The model is first applied to a set of hypothetical chemicals to illustrate the degree to which test conditions (e.g., presence or absence of serum) influence the distribution of the chemical in the test system. The model is then applied to set of 1194 real substances (predominantly from the ToxCast chemical database) to calculate the potential range of concentrations and chemical activities under assumed test conditions. The model demonstrates how both concentrations and chemical activities can vary by orders of magnitude for the same nominal concentration.

Toxic response of HIPCO single-walled carbon nanotubes in mice and RAW264.7 macrophage cells

In this study, we identified the toxic response of pristine single-walled carbon nanotubes (P-SWCNTs) synthesized by HIPCO method in mice and RAW264.7 cells, a murine peritoneal macrophage cell line. P-SWCNT contained a large amount of Fe ion (36 wt%). In the lungs of mice 24 h after intratracheal administration, P-SWCNTs increased the secretion of IL-6 and MCP-1, and the number of total cells, the portion of neutrophils, lymphocytes, and eosinophils, also significantly increased at a 100 μg/mL of concentration. In RAW264.7 cells, cell viability and ATP production decreased in a dose-dependent manner at 24 h after exposure, whereas the generations of ROS and NO were enhanced at all concentrations together with the activation of the MAP kinase pathway. Moreover, the levels of both apoptosis- and autophagy-related proteins and ER stress-related proteins clearly increased, and the concentrations of Fe, Cu, and Zn ions, but not of Mn ions, increased in a dose-dependent manner. TEM images also revealed that P-SWCNTs induced the formation of autophagosome-like vacuoles, the dilatation of the ER, the generation of mitochondrial flocculent densities, and the separation of organelle by disappearance of the cell membrane. Taken together, we suggest that P-SWCNTs cause acute inflammatory response in the lungs of mice, and induce autophagy accompanied with apoptosis through mitochondrial dysfunction and ER stress in RAW264.7 cells. Furthermore, further study is required to elucidate how the physicochemical properties of SWCNTs determine the cell death pathway and an immune response.

Evaluation of differential cytotoxic effects of the oil spill dispersant Corexit 9500

AIMS

The British Petroleum (BP) oil spill has raised several ecological and health concerns. As the first response, BP used a chemical dispersant, Corexit-9500, to disperse the crude oil in the Gulf of Mexico to limit shoreline contamination problems. Nevertheless, portions of this oil/Corexit mixture reached the shoreline and still remain in various Gulf shore environments. The use of Corexit itself has become a significant concern since its impacts on human health and environment is unclear.

MAIN METHODS

In this study, in vitro cytotoxic effects of Corexit were evaluated using different mammalian cells.

KEY FINDINGS

Under serum free conditions, the LC50 value for Corexit in BL16/BL6 cell was 16 ppm, in 1321N1 cell was 33 ppm, in H19-7 cell was 70 ppm, in HEK293 was 93 ppm, and in HK-2 cell was 95 ppm. With regard to the mechanisms of cytotoxicity, we hypothesize that Corexit can possibly induce cytotoxicity in mammalian cells by altering the intracellular oxidative balance and inhibiting mitochondrial functions. Corexit induced increased reactive oxygen species and lipid peroxide levels; also, it depleted glutathione content and altered catalase activity in H19-7 cells. In addition, there was mitochondrial complex-I inhibition and increase in the pro-apoptotic factors including caspase-3 and BAX expression.

SIGNIFICANCE

The experimental results show changes in intracellular oxidative radicals leading to mitochondrial dysfunctions and apoptosis in Corexit treatments, possibly contributing to cell death. Our findings raise concerns about using large volumes of Corexit, a potential environmental toxin, in sensitive ocean environments.

Protein and lipid binding parameters in rainbow trout (Oncorhynchus mykiss) blood and liver fractions to extrapolate from an in vitro metabolic degradation assay to in vivo bioaccumulation potential of hydrophobic organic chemicals

Binding of hydrophobic chemicals to colloids such as proteins or lipids is difficult to measure using classical microdialysis methods due to low aqueous concentrations, adsorption to dialysis membranes and test vessels, and slow kinetics of equilibration. Here, we employed a three-phase partitioning system where silicone (polydimethylsiloxane, PDMS) serves as a third phase to determine partitioning between water and colloids and acts at the same time as a dosing device for hydrophobic chemicals. The applicability of this method was demonstrated with bovine serum albumin (BSA). Measured binding constants (K(BSAw)) for chlorpyrifos, methoxychlor, nonylphenol, and pyrene were in good agreement with an established quantitative structure-activity relationship (QSAR). A fifth compound, fluoxypyr-methyl-heptyl ester, was excluded from the analysis because of apparent abiotic degradation. The PDMS depletion method was then used to determine partition coefficients for test chemicals in rainbow trout (Oncorhynchus mykiss) liver S9 fractions (K(S9w)) and blood plasma (K(bloodw)). Measured K(S9w) and K(bloodw) values were consistent with predictions obtained using a mass-balance model that employs the octanol-water partition coefficient (K(ow)) as a surrogate for lipid partitioning and K(BSAw) to represent protein binding. For each compound, K(bloodw) was substantially greater than K(S9w), primarily because blood contains more lipid than liver S9 fractions (1.84% of wet weight vs 0.051%). Measured liver S9 and blood plasma binding parameters were subsequently implemented in an in vitro to in vivo extrapolation model to link the in vitro liver S9 metabolic degradation assay to in vivo metabolism in fish. Apparent volumes of distribution (V(d)) calculated from the experimental data were similar to literature estimates. However, the calculated binding ratios (f(u)) used to relate in vitro metabolic clearance to clearance by the intact liver were 10 to 100 times lower than values used in previous modeling efforts. Bioconcentration factors (BCF) predicted using the experimental binding data were substantially higher than the predicted values obtained in earlier studies and correlated poorly with measured BCF values in fish. One possible explanation for this finding is that chemicals bound to proteins can desorb rapidly and thus contribute to metabolic turnover of the chemicals. This hypothesis remains to be investigated in future studies, ideally with chemicals of higher hydrophobicity.

Long-term exposure to dieldrin reduces γ-aminobutyric acid type A and N-methyl-D-aspartate receptor function in primary cultures of mouse cerebellar granule cells

The organochlorine pesticide dieldrin is a persistent organic pollutant that accumulates in the fatty tissue of living organisms. In mammals, it antagonizes the GABA(A) receptor, producing convulsions after acute exposure. Although accumulation in human brain has been reported, little is known about the effects of long-term exposure to dieldrin in the nervous system. Homeostatic control of the balance between excitation and inhibition has been reported when neuronal activity is chronically altered. We hypothesized that noncytotoxic concentrations of dieldrin could decrease glutamatergic neurotransmission as a consequence of a prolonged reduction in GABA(A) receptor function. Long-term exposure of primary cerebellar granule cell cultures to 3 microM dieldrin reduced the GABA(A) receptor function to 55% of control, as measured by the GABA-induced (36)Cl(-) uptake. This exposure produced a significant reduction (approximately 35%) of the NMDA-induced increase in [Ca(2+)](i) and of the [(3)H]MK-801 binding, which was not accompanied by a reduction in the NMDA receptor subunit NR1, as determined by Western blot. Consistent with the decreased NMDA receptor function, dieldrin-treated cultures were insensitive to an excitotoxic stimulus induced by exposure to high potassium. In summary, we report that the chronic reduction of GABA(A) receptor function induced by dieldrin decreases the number of functional NMDA receptors, which may be attributable to a mechanism of synaptic scaling. These effects could underlie neural mechanisms involved in cognitive impairment produced by low-level exposure to dieldrin.