PBPK modeling to support risk assessment of pyrethroid exposure in French pregnant women.

BACKGROUND: Pyrethroids are widely used pesticides and are suspected to affect children's neurodevelopment. The characterization of pyrethroid exposure during critical windows of development, such as fetal development and prenatal life, is essential to ensure a better understanding of pyrethroids potential effects within the concept of Developmental Origins of Health and Disease. OBJECTIVE: The aim of this study was to estimate maternal exposure of French pregnant women from biomonitoring data and simulate maternal and fetal internal concentrations of 3 pyrethroids (permethrin, cypermethrin and deltamethrin) using a multi-substance pregnancy-PBPK (physiologically based pharmacokinetics) model. The estimated maternal exposures were compared to newly proposed toxicological reference values (TRV) children specific also called draft child-specific reference value to assess pyrethroid exposure risk during pregnancy i.e. during the in utero exposure period. METHODS: A pregnancy-PBPK model was developed based on an existing adult pyrethroids model. The maternal exposure to each parent compound of pregnant women of the Elfe (French Longitudinal Study since Childhood) cohort was estimated by reverse dosimetry based on urinary biomonitoring data. To identify permethrin and cypermethrin contribution to their common urinary biomarkers of exposure, an exposure ratio based on biomarkers in hair was tested. Finally, exposure estimates were compared to current and draft child-specific reference values derived from rodent prenatal and postnatal exposure studies. RESULTS: The main contributor to maternal pyrethroid diet intake is cis-permethrin. In blood, total internal concentrations main contributor is deltamethrin. In brain, the major contributors to internal pyrethroid exposure are deltamethrin for fetuses and cis-permethrin for mothers. Risk is identified only for permethrin when referring to the draft child-specific reference value. 2.5% of the population exceeded permethrin draft child-specific reference value. CONCLUSIONS: A new reverse dosimetry approach using PBPK model combined with human biomonitoring data in urine and hair was proposed to estimate Elfe pregnant population exposure to a pyrethroids mixture with common metabolites.

Development of a physiologically based toxicokinetic model for lead in pregnant women: The role of bone tissue in the maternal and fetal internal exposure.

Epidemiological studies have shown associations between prenatal exposure to lead (Pb) and neurodevelopmental effects in young children. Prenatal exposure is generally characterized by measuring the concentration in the umbilical cord at delivery or in the maternal blood during pregnancy. To assess internal Pb exposure during prenatal life, we developed a pregnancy physiologically based pharmacokinetic (p-PBPK) model that to simulates Pb levels in blood and target tissues in the fetus, especially during critical periods for brain development. An existing Pb PBPK model was adapted to pregnant women and fetuses. Using data from literature, both the additional maternal bone remodeling, that causes Pb release into the blood, and the Pb placental transfers were estimated by Bayesian inference. Additional maternal bone remodeling was estimated to start at 21.6 weeks. Placental transfers were estimated between 4.6 and 283 L.day-1 at delivery with high interindividual variability. Once calibrated, the p-PBPK model was used to simulate fetal exposure to Pb. Internal fetal exposure greatly varies over the pregnancy with two peaks of Pb levels in blood and brain at the end of the 1st and 3rd trimesters. Sensitivity analysis shows that the fetal blood lead levels are affected by the maternal burden of bone Pb via maternal bone remodeling and by fetal bone formation at different pregnancy stages. Coupling the p-PBPK model with an effect model such as an adverse outcome pathway could help to predict the effects on children's neurodevelopment.

Estimating human exposure to pyrethroids' mixtures from biomonitoring data using physiologically based pharmacokinetic modeling.

Human biomonitoring data provide evidence to exposure of environmental chemicals. Physiologically based pharmacokinetic (PBPK) modelling together with an adequate exposure scenario allows to transpose measured concentrations of chemicals or their metabolites into exposure levels, as daily intakes. In France, high levels of urinary pyrethroids metabolites have been measured in populations. Our work aims at estimating the exposure of the French ENNS cohort to mixtures of four pyrethroids (deltamethrin, permethrin, cypermethrin, and cyfluthrin) from the urinary concentrations of five pyrethroids' metabolites commonly measured in biomonitoring studies. We developed a modelling approach based on a global toxicokinetic model that accounts for the cumulative exposure to pyrethroids as some of the metabolites can be shared by several parent compounds and for human inter-individual variability in metabolism. The median of the individual daily intakes was estimated to 8.1 ng/kg bw/day for permethrin, 17.7 ng/kg bw/day for cypermethrin, 20.4 ng/kg bw/day for cyfluthrin and 34.3 ng/kg bw/day for deltamethrin leading to similar weights for the pair permethrin and cypermethrin (36%), cyfluthrin (31%) and deltamethrin (33%) to the cumulative exposure. Accounting for human variability enabled to explain some of the variations in the metabolites' levels within the cohort. The cumulative exposure was then weighted by their toxicities towards three neurotoxic effects to calculate margins of exposure (MOE). Low MOE values were always associated with high measured concentrations of metabolites in urine and the lowest MOEs were observed for the autonomic division. No risks associated with reconstructed mixtures of pyrethroids were expected for the ENNS cohort. Our approach is an asset to analyse the biomarkers of exposure to pyrethroids simultaneously and could be easily adapted to any local or national specificities in pyrethroids' exposure or populations.

Characterizing environmental geographic inequalities using an integrated exposure assessment.

BACKGROUND: At a regional or continental scale, the characterization of environmental health inequities (EHI) expresses the idea that populations are not equal in the face of pollution. It implies an analysis be conducted in order to identify and manage the areas at risk of overexposure where an increasing risk to human health is suspected. The development of methods is a prerequisite for implementing public health activities aimed at protecting populations. METHODS: This paper presents the methodological framework developed by INERIS (French National Institute for Industrial Environment and Risks) to identify a common framework for a structured and operationalized assessment of human exposure. An integrated exposure assessment approach has been developed to integrate the multiplicity of exposure pathways from various sources, through a series of models enabling the final exposure of a population to be defined. RESULTS: Measured data from environmental networks reflecting the actual contamination of the environment are used to gauge the population's exposure. Sophisticated methods of spatial analysis are applied to include additional information and take benefit of spatial and inter-variable correlation to improve data representativeness and characterize the associated uncertainty. Integrated approaches bring together all the information available for assessing the source-to-human-dose continuum using a Geographic Information System, multimedia exposure and toxicokinetic model. DISCUSSION: One of the objectives of the integrated approach was to demonstrate the feasibility of building complex realistic exposure scenarios satisfying the needs of stakeholders and the accuracy of the modelling predictions at a fine spatial-temporal resolution. A case study is presented to provide a specific application of the proposed framework and how the results could be used to identify an overexposed population. CONCLUSION: This framework could be used for many purposes, such as mapping EHI, identifying vulnerable populations and providing determinants of exposure to manage and plan remedial actions and to assess the spatial relationships between health and the environment to identify factors that influence the variability of disease patterns.

Assessing the impacts on fetal dosimetry of the modelling of the placental transfers of xenobiotics in a pregnancy physiologically based pharmacokinetic model.

The developmental origin of health and diseases theory supports the critical role of the fetal exposure to children's health. We developed a physiologically based pharmacokinetic model for human pregnancy (pPBPK) to simulate the maternal and fetal dosimetry throughout pregnancy. Four models of the placental exchanges of chemicals were assessed on ten chemicals for which maternal and fetal data were available. These models were calibrated using non-animal methods: in vitro (InV) or ex vivo (ExV) data, a semi-empirical relationship (SE), or the limitation by the placental perfusion (PL). They did not impact the maternal pharmacokinetics but provided different profiles in the fetus. The PL and InV models performed well even if the PL model overpredicted the fetal exposure for some substances. The SE and ExV models showed the lowest global performance and the SE model a tendency to underprediction. The comparison of the profiles showed that the PL model predicted an increase in the fetal exposure with the pregnancy age, whereas the ExV model predicted a decrease. For the SE and InV models, a small decrease was predicted during the second trimester. All models but the ExV one, presented the highest fetal exposure at the end of the third trimester. Global sensitivity analyses highlighted the predominant influence of the placental transfers on the fetal exposure, as well as the metabolic clearance and the fraction unbound. Finally, the four transfer models could be considered depending on the framework of the use of the pPBPK model and the availability of data or resources to inform their parametrization.

Integrative Strategy of Testing Systems for Identification of Endocrine Disruptors Inducing Metabolic Disorders-An Introduction to the OBERON Project.

Exposure to chemical substances that can produce endocrine disrupting effects represents one of the most critical public health threats nowadays. In line with the regulatory framework implemented within the European Union (EU) to reduce the levels of endocrine disruptors (EDs) for consumers, new and effective methods for ED testing are needed. The OBERON project will build an integrated testing strategy (ITS) to detect ED-related metabolic disorders by developing, improving and validating a battery of test systems. It will be based on the concept of an integrated approach for testing and assessment (IATA). OBERON will combine (1) experimental methods (in vitro, e.g., using 2D and 3D human-derived cells and tissues, and in vivo, i.e., using zebrafish at different stages), (2) high throughput omics technologies, (3) epidemiology and human biomonitoring studies and (4) advanced computational models (in silico and systems biology) on functional endpoints related to metabolism. Such interdisciplinary framework will help in deciphering EDs based on a mechanistic understanding of toxicity by providing and making available more effective alternative test methods relevant for human health that are in line with regulatory needs. Data generated in OBERON will also allow the development of novel adverse outcome pathways (AOPs). The assays will be pre-validated in order to select the test systems that will show acceptable performance in terms of relevance for the second step of the validation process, i.e., the inter-laboratory validation as ring tests. Therefore, the aim of the OBERON project is to support the organization for economic co-operation and development (OECD) conceptual framework for testing and assessment of single and/or mixture of EDs by developing specific assays not covered by the current tests, and to propose an IATA for ED-related metabolic disorder detection, which will be submitted to the Joint Research Center (JRC) and OECD community.

Investigating the interaction between melamine and cyanuric acid using a Physiologically-Based Toxicokinetic model in rainbow trout.

Following outbreaks of feed and food adulterations with a melamine and cyanuric acid mixture in 2007 and melamine in 2008 respectively, the kinetics and toxicodynamics of the mixture have been investigated particularly in sensitive species such as the rainbow trout. Tissue concentrations and intensity of the adverse effect, melamine-cyanurate crystal formation in kidney, were reported in similar experimental conditions. Here, a recent PBTK model for rainbow trout has been applied to model the kinetics of both single compounds based on residue levels in tissues. Both PBTK models for the single compounds were combined and a model of crystal formation for the mixture melamine-cyanuric acid was also added to predict the intensity of crystal formation under the assumptions that crystals formed either in urine or in kidney tissue. Modelling the kinetics of melamine and cyanuric acid provided a better understanding and prediction of intensity of crystal formation in case of sequential exposures with varying intensity or co-exposure. This study demonstrates, for the first time, how fish PBTK models can play a key role in the understanding and prediction of toxicokinetics and toxicodynamics of mixtures. This study also illustrates how adverse effects may potentially occur even when the compounds are not administered together as a mixture.

Prediction of maternal and foetal exposures to perfluoroalkyl compounds in a Spanish birth cohort using toxicokinetic modelling.

Prenatal exposures to perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) have been associated with child health outcomes, but many of these associations remain poorly characterized. The aim of this work was to provide new indicators of foetal exposure for the Spanish INMA birth cohort. First, a pregnancy and lactation physiologically based pharmacokinetic (PBPK) model was calibrated in a population framework to provide quantitative estimates for the PFOA and PFOS placental transfers in humans. The estimated distributions indicated that PFOA crosses the placental barrier at a rate three times higher than PFOS and shows a higher variability between mothers. The PBPK model was then used to back-calculate the time-varying daily intakes of the INMA mothers corrected for their individual history from a spot maternal concentration. We showed the importance of accounting for the mothers' history as different dietary intakes can result in similar measured concentrations at one time point. Finally, the foetal exposure was simulated in target organs over pregnancy using the PBPK model and the estimated maternal intakes. We showed that the pattern of PFOA and PFOS exposures varies greatly among the foetuses. About a third has levels of either one compound always higher than the levels of the other compound. The other two thirds showed different ranking of PFOA and PFOS in terms of concentrations in the target organs. Our simulated foetal exposures bring additional information to the measured maternal spot concentrations and can help to better characterize the prenatal exposure in target organs during windows of susceptibility.

Determination of maternal and foetal distribution of cis- and trans-permethrin isomers and their metabolites in pregnant rats by liquid chromatography tandem mass spectrometry (LC-MS/MS).

We developed a method to quantify cis-permethrin and trans-permethrin and their metabolites in several biological matrices in pregnant rats and foetuses using liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS). The objective was to quantify cis-permethrin and trans-permethrin in faeces, kidney, mammary gland, fat and placenta in mothers and in both maternal and foetal blood, brain and liver. The metabolites cis-3-(2,2-dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (cis-DCCA), trans-3-(2,2-dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (trans-DCCA) and 3-phenoxybenzoic acid (3-PBA) were measured in blood, liver and urine. Sample preparation was performed by liquid-liquid extraction. A purification step was not carried out except for the more complex biological samples (fat, mammary glands and faeces). Validation parameters including specificity, linearity, matrix effect, limits of quantification (LOQs), accuracy and precision were evaluated. The recoveries of target compounds ranged from 47 to 136%. LOQs were in the range 4 to 80 ng/mL for permethrin isomers and 4 to 800 ng/mL for their respective metabolites. Intra- and inter-batch precision and accuracy in matrix were better than 15%. The validated method was applied in a preliminary toxicokinetic study in pregnant rats with oral dosing of 50 mg/kg permethrin. In pregnant rats, permethrin isomers and their metabolites were quantified in all requested matrices except maternal liver and blood for trans-permethrin and cis-DCCA respectively. In foetuses, cis- and trans-permethrin were also quantified, demonstrating that the method is suitable for the analysis of foetal distribution of permethrin in toxicokinetic studies.

Physiology-based toxicokinetic modelling in the frame of the European Human Biomonitoring Initiative.

Given the opportunities provided by internal dosimetry modelling in the interpretation of human biomonitoring (HBM) data, the assessment of the links between exposure to chemicals and observed HBM data can be effectively supported by PBTK modelling. This paper gives a comprehensive review of available human PBTK models for compounds selected as a priority by the European Human Biomonitoring Initiative (HBM4EU). We highlight their advantages and deficiencies and suggest steps for advanced internal dose modelling. The review of the available PBTK models highlighted the conceptual differences between older models compared to the ones developed recently, reflecting commensurate differences in research questions. Due to the lack of coordinated strategies for deriving useful biomonitoring data for toxicokinetic properties, significant problems in model parameterisation still remain; these are further increased by the lack of human toxicokinetic data due to ethics issues. Finally, questions arise as well as to the extent they are really representative of interindividual variability. QSARs for toxicokinetic properties is a complementary approach for PBTK model parameterisation, especially for data poor chemicals. This approach could be expanded to model chemico-biological interactions such as intestinal absorption and renal clearance; this could serve the development of more complex generic PBTK models that could be applied to newly derived chemicals. Another gap identified is the framework for mixture interaction terms among compounds that could eventually interact in metabolism. From the review it was concluded that efforts should be shifted toward the development of generic multi-compartmental and multi-route models, supported by targeted biomonitoring coupled with parameterisation by both QSAR approach and experimental (in-vivo and in-vitro) data for newly developed and data poor compounds.

PBPK modeling of the cis- and trans-permethrin isomers and their major urinary metabolites in rats.

Permethrin, a pyrethroid insecticide, is suspected to induce neuronal and hormonal disturbances in humans. The widespread exposure of the populations has been confirmed by the detection of the urinary metabolites of permethrin in biomonitoring studies. Permethrin is a chiral molecule presenting two forms, the cis and the trans isomers. Because in vitro studies indicated a metabolic interaction between the trans and cis isomers of permethrin, we adapted and calibrated a PBPK model for trans- and cis-permethrin separately in rats. The model also describes the toxicokinetics of three urinary metabolites, cis- and trans-3-(2,2 dichlorovinyl)-2,2-dimethyl-(1-cyclopropane) carboxylic acid (cis- and trans-DCCA), 3-phenoxybenzoic acid (3-PBA) and 4'OH-phenoxybenzoic acid (4'-OH-PBA). In vivo experiments performed in Sprague-Dawley rats were used to calibrate the PBPK model in a Bayesian framework. The model captured well the toxicokinetics of permethrin isomers and their metabolites including the rapid absorption, the accumulation in fat, the extensive metabolism of the parent compounds, and the rapid elimination of metabolites in urine. Average hepatic clearances in rats were estimated to be 2.4 and 5.7 L/h/kg for cis- and trans-permethrin, respectively. High concentrations of the metabolite 4'-OH-PBA were measured in urine compared to cis- and trans-DCCA and 3-PBA. The confidence in the extended PBPK model was then confirmed by good predictions of published experimental data obtained using the isomers mixture. The extended PBPK model could be extrapolated to humans to predict the internal dose of exposure to permethrin from biomonitoring data in urine.

A physiologically based toxicokinetic model for the zebrafish Danio rerio.

Zebrafish (Danio rerio) is a widely used model for toxicological studies, in particular those related to investigations on endocrine disruption. The development and regulatory use of in vivo and in vitro tests based on this species can be enhanced by toxicokinetic modeling. For this reason, we propose a physiologically based toxicokinetic (PBTK) model for zebrafish describing the uptake and disposition of organic chemicals. The model is based on literature data on zebrafish, other cyprinidae and other fish families, new experimental physiological information (volumes, lipids and water contents) obtained from zebrafish, and chemical-specific parameters predicted by generic models. The relevance of available models predicting the latter parameters was evaluated with respect to gill uptake and partition coefficients in zebrafish. This evaluation benefited from the fact that the influence of confounding factors such as body weight and temperature on ventilation rate was included in our model. The predictions for six chemicals (65 data points) yielded by our PBTK model were compared to available toxicokinetics data for zebrafish and 88% of them were within a factor of 5 of the corresponding experimental values. Sensitivity analysis highlighted that the 1-octanol/water partition coefficient, the metabolism rate, and all the parameters that enable the prediction of assimilation efficiency and partitioning of chemicals need to be precisely determined in order to allow an effective toxicokinetic modeling.

Estimating the dynamic early life exposure to PFOA and PFOS of the HELIX children: Emerging profiles via prenatal exposure, breastfeeding, and diet.

In utero and children's exposure to per- and polyfluoroalkyl substances (PFAS) is a major concern in health risk assessment as early life exposures are suspected to induce adverse health effects. Our work aims to estimate children's exposure (from birth to 12 years old) to PFOA and PFOS, using a Physiologically-Based Pharmacokinetic (PBPK) modelling approach. A model for PFAS was updated to simulate the internal PFAS exposures during the in utero life and childhood, and including individual characteristics and exposure scenarios (e.g., duration of breastfeeding, weight at birth, etc.). Our approach was applied to the HELIX cohort, involving 1,239 mother-child pairs with measured PFOA and PFOS plasma concentrations at two sampling times: maternal and child plasma concentrations (6 to 12 y.o). Our model predicted an increase in plasma concentrations during fetal development and childhood until 2 y.o when the maximum concentrations were reached. Higher plasma concentrations of PFOA than PFOS were predicted until 2 y.o, and then PFOS concentrations gradually became higher than PFOA concentrations. From 2 to 8 y.o, mean concentrations decreased from 3.1 to 1.88 µg/L or ng/mL (PFOA) and from 4.77 to 3.56 µg/L (PFOS). The concentration-time profiles vary with the age and were mostly influenced by in utero exposure (on the first 4 months after birth), breastfeeding (from 5 months to 2 (PFOA) or 5 (PFOS) y.o of the children), and food intake (after 3 (PFOA) or 6 (PFOS) y.o of the children). Similar measured biomarker levels can correspond to large differences in the simulated internal exposures, highlighting the importance to investigate the children's exposure over the early life to improve exposure classification. Our approach demonstrates the possibility to simulate individual internal exposures using PBPK models when measured biomarkers are scarce, helping risk assessors in gaining insight into internal exposure during critical windows, such as early life.

Investigation of ifosfamide nephrotoxicity induced in a liver-kidney co-culture biochip.

In this article, we present a liver-kidney co-culture model in a micro fluidic biochip. The liver was modeled using HepG2/C3a and HepaRG cell lines and the kidney using MDCK cell lines. To demonstrate the synergic interaction between both organs, we investigated the effect of ifosfamide, an anticancerous drug. Ifosfamide is a prodrug which is metabolized by the liver to isophosforamide mustard, an active metabolite. This metabolism process also leads to the formation of chloroacetaldehyde, a nephrotoxic metabolite and acrolein a urotoxic one. In the biochips of MDCK cultures, we did not detect any nephrotoxic effects after 72 h of 50 µM ifosfamide exposure. However, in the liver-kidney biochips, the same 72 h exposure leads to a nephrotoxicity illustrated by a reduction of the number of MDCK cells (up to 30% in the HepaRG-MDCK) when compared to untreated co-cultures or treated MDCK monocultures. The reduction of the MDCK cell number was not related to a modification of the cell cycle repartition in ifosfamide treated cases when compared to controls. The ifosfamide biotransformation into 3-dechloroethylifosfamide, an equimolar byproduct of the chloroacetaldehyde production, was detected by mass spectrometry at a rate of apparition of 0.3 ± 0.1 and 1.1 ± 0.3 pg/h/biochips in HepaRG monocultures and HepaRG-MDCK co-cultures respectively. Any metabolite was detected in HepG2/C3a cultures. Furthermore, the ifosfamide treatment in HepaRG-MDCK co-culture system triggered an increase in the intracellular calcium release in MDCK cells on contrary to the treatment on MDCK monocultures. As 3-dechloroethylifosfamide is not toxic, we have tested the effect of equimolar choloroacetaldehyde concentration onto the MDCK cells. At this concentration, we found a quite similar calcium perturbation and MDCK nephrotoxicity via a reduction of 30% of final cell numbers such as in the ifosfamide HepaRG-MDCK co-culture experiments. Our results suggest that ifosfamide nephrotoxicity in a liver-kidney micro fluidic co-culture model using HepaRG-MDCK cells is induced by the metabolism of ifosfamide into chloroacetaldehyde whereas this pathway is not functional in HepG2/C3a-MDCK model. This study demonstrates the interest in the development of systemic organ-organ interactions using micro fluidic biochips. It also illustrated their potential in future predictive toxicity model using in vitro models as alternative methods.

Evaluation of seven drug metabolisms and clearances by cryopreserved human primary hepatocytes cultivated in microfluidic biochips.

We present characterization of the metabolic performance of human cryopreserved hepatocytes cultivated in a platform of parallelized microfluidic biochips. The RTqPCR analysis revealed that the mRNA levels of the cytochromes P450 (CYP 1A2, 2B6, 2C8, 2C9, 2C19, 2D6, 2E1, 3A4) were reduced after the adhesion period (when compared to the post-thawing step). The microfluidic perfusion played a part in stabilizing and partially recovering the levels of the HNF4α, PXR, OAPT2, CYP 1A2, 2B6, 2C19 and 3A4 mRNA on contrary to non-perfused cultures. Fluorescein diacetate staining and P-gp mRNA level illustrated the hepatocytes' polarity in the biochips. Drug metabolism was assessed using midazolam, tolbutamide, caffeine, omeprazole, dextromethorphan, acetaminophen and repaglinide as probes. Metabolite detection and quantification revealed that CYP1A2 (via the detection of paraxanthine), CYP3A4 (via 1-OH-midazolam, and omeprazole sulfone detection), CYP2C8 (via hydroxyl-repaglinide detection), CYP2C19 (via hydroxy-omeprazole detection) and CYP2D6 (via dextrorphan detection) were functional in our microfluidic configurations. Furthermore, the RTqPCR analysis showed that the drugs acted as inductors leading to overexpression of mRNA levels when compared to post-thawing values (such as for HNF4α, PXR and CYP3A4 by dextromethorpahn and omeprazole). Finally, intrinsic in vitro biochip clearances were extracted using a PBPK model for predictions. The biochip predictions were compared to literature in vitro data and in vivo situations.

A PBPK model to evaluate zebrafish eleutheroembryos' actual exposure: bisphenol A and analogs' (AF, F, and S) case studies.

The zebrafish eleutheroembryo model is increasingly used to assess the toxicity and developmental adverse effects of xenobiotics. However, the actual exposure is seldom measured (poorly accessible), while a predictive model could estimate these concentrations. The predictions with a new eleutheroembryo physiologically based pharmacokinetic (PBPK) model have been evaluated using datasets obtained from literature data for several bisphenols. The model simulated the toxicokinetics of bisphenols A (BPA), AF, F, and S through the eleutheroembryo tissues while considering the body and organ growth. We further improved the predictions by adding dynamic flows through the embryo and/or its chorion, impact of experimental temperature, metabolic clearance, and saturation of the absorption by Bayesian calibration. The model structure was determined using the BPA dataset and generalized to the other bisphenols. This model revealed the central role of the chorion in the compound uptake in the first 48 h post-fertilization. The predictions for the BPA substitutes estimated by our PBPK model were compared to available toxicokinetics data for zebrafish embryos, and 63% and 88% of them were within a twofold and fivefold error intervals of the corresponding experimental values, respectively. This model provides a tool to design new eleutheroembryo assays and evaluate the actual exposure.

PBTK modeled perfluoroalkyl acid kinetics in zebrafish eleutheroembryos suggests impacts on bioconcentrations by chorion porosity dynamics.

The zebrafish eleutheroembryo (zfe) is widely used as a model to characterize the toxicity of chemicals. However, analytical methods are still missing to measure organ concentrations. Therefore, physiologically-based toxicokinetic (PBTK) modeling may overcome current limitations to help understand the relationship between toxic effects and internal exposure in various organs. A previous PBTK model has been updated to include the chorionic transport barrier and its permeabilization, hatching dynamics within a zfe population over development, and active mediated transport mechanisms. The zfe PBTK model has been calibrated using measured time-dependent internal concentrations of PFBA, PFHxS, PFOA, and PFOS in a zfe population and evaluated using external datasets from the literature. Calibration was successful with 96% of the predictions falling within a 2-fold range of the observed concentrations. The external dataset was correctly estimated with about 50% of the predictions falling within a factor of 3 of the observed data and 10% of the predictions are out of the 10-fold error. The calibrated model suggested that active mediated transport differs between PFAS with a sulfonic and carboxylic acid functional end groups. This PBTK model predicts well the fate of PFAS with various physicochemical properties in zfe. Therefore, this model may improve the use of zfe as an alternative model in toxicokinetic-toxicodynamic studies and help to refine and reduce zfe-based experiments, while giving insights into the internal kinetics of chemicals.

Metabolomics-on-a-chip and predictive systems toxicology in microfluidic bioartificial organs.

The world faces complex challenges for chemical hazard assessment. Microfluidic bioartificial organs enable the spatial and temporal control of cell growth and biochemistry, critical for organ-specific metabolic functions and particularly relevant to testing the metabolic dose-response signatures associated with both pharmaceutical and environmental toxicity. Here we present an approach combining a microfluidic system with (1)H NMR-based metabolomic footprinting, as a high-throughput small-molecule screening approach. We characterized the toxicity of several molecules: ammonia (NH(3)), an environmental pollutant leading to metabolic acidosis and liver and kidney toxicity; dimethylsulfoxide (DMSO), a free radical-scavenging solvent; and N-acetyl-para-aminophenol (APAP, or paracetamol), a hepatotoxic analgesic drug. We report organ-specific NH(3) dose-dependent metabolic responses in several microfluidic bioartificial organs (liver, kidney, and cocultures), as well as predictive (99% accuracy for NH(3) and 94% for APAP) compound-specific signatures. Our integration of microtechnology, cell culture in microfluidic biochips, and metabolic profiling opens the development of so-called "metabolomics-on-a-chip" assays in pharmaceutical and environmental toxicology.

Predictive toxicology using systemic biology and liver microfluidic "on chip" approaches: application to acetaminophen injury.

We have analyzed transcriptomic, proteomic and metabolomic profiles of hepatoma cells cultivated inside a microfluidic biochip with or without acetaminophen (APAP). Without APAP, the results show an adaptive cellular response to the microfluidic environment, leading to the induction of anti-oxidative stress and cytoprotective pathways. In presence of APAP, calcium homeostasis perturbation, lipid peroxidation and cell death are observed. These effects can be attributed to APAP metabolism into its highly reactive metabolite, N-acetyl-p-benzoquinone imine (NAPQI). That toxicity pathway was confirmed by the detection of GSH-APAP, the large production of 2-hydroxybutyrate and 3-hydroxybutyrate, and methionine, cystine, and histidine consumption in the treated biochips. Those metabolites have been reported as specific biomarkers of hepatotoxicity and glutathione depletion in the literature. In addition, the integration of the metabolomic, transcriptomic and proteomic collected profiles allowed a more complete reconstruction of the APAP injury pathways. To our knowledge, this work is the first example of a global integration of microfluidic biochip data in toxicity assessment. Our results demonstrate the potential of that new approach to predictive toxicology.

PBPK Modeling to Simulate the Fate of Compounds in Living Organisms.

Pharmacokinetics study the fate of xenobiotics in a living organism. Physiologically based pharmacokinetic (PBPK) models provide realistic descriptions of xenobiotics' absorption, distribution, metabolism, and excretion processes. They model the body as a set of homogeneous compartments representing organs, and their parameters refer to anatomical, physiological, biochemical, and physicochemical entities. They offer a quantitative mechanistic framework to understand and simulate the time-course of the concentration of a substance in various organs and body fluids. These models are well suited for performing extrapolations inherent to toxicology and pharmacology (e.g., between species or doses) and for integrating data obtained from various sources (e.g., in vitro or in vivo experiments, structure-activity models). In this chapter, we describe the practical development and basic use of a PBPK model from model building to model simulations, through implementation with an easily accessible free software.