Integrating in vitro data and physiologically based kinetic (PBK) modelling to assess the in vivo potential developmental toxicity of a series of phenols

Predicting the in vivo developmental toxicity of fenarimol from in vitro toxicity data using PBTK modelling-facilitated reverse dosimetry approach

In vitro methods are widely used in modern toxicological testing; however, the data cannot be directly employed for risk assessment. In vivo toxicity of chemicals can be predicted from in vitro data using physiologically based toxicokinetic (PBTK) modelling-facilitated reverse dosimetry (PBTK-RD). In this study, a minimal-PBTK model was constructed to predict the in-vivo kinetic profile of fenarimol (FNL) in rats and humans. The model was verified by comparing the observed and predicted pharmacokinetics of FNL for rats (calibrator) and further applied to humans. Using the PBTK-RD approach, the reported in vitro developmental toxicity data for FNL was translated to in vivo dose-response data to predict the assay equivalent oral dose in rats and humans. The predicted assay equivalent rat oral dose (36.46 mg/kg) was comparable to the literature reported in vivo BMD10 value (22.8 mg/kg). The model was also employed to derive the chemical-specific adjustment factor (CSAF) for interspecies toxicokinetics variability of FNL. Further, Monte Carlo simulations were performed to predict the population variability in the plasma concentration of FNL and to derive CSAF for intersubject human kinetic differences. The comparison of CSAF values for interspecies and intersubject toxicokinetic variability with their respective default values revealed that the applied uncertainty factors were adequately protective.

Predicting acute paraquat toxicity using physiologically based kinetic modelling incorporating in vitro active renal excretion via the OCT2 transporter

Including active renal excretion in physiologically based kinetic (PBK) models can improve their use in quantitative in vitro- in vivo extrapolation (QIVIVE) as a new approach methodology (NAM) for predicting the acute toxicity of organic cation transporter 2 (OCT2) substrates like paraquat (PQ). To realise this NAM, kinetic parameters Vmax and Km for in vitro OCT2 transport of PQ were obtained from the literature. Appropriate scaling factors were applied to translate the in vitro Vmax to an in vivo Vmax. in vitro cytotoxicity data were defined in the rat RLE-6TN and L2 cell lines and the human A549 cell line. The developed PQ PBK model was used to apply reverse dosimetry for QIVIVE translating the in vitro cytotoxicity concentration-response curves to predicted in vivo toxicity dose-response curves after which the lower and upper bound benchmark dose (BMD) for 50% lethality (BMDL50 and BMDU50) were derived by applying BMD analysis. Comparing the predictions to the in vivo reported LD50 values resulted in a conservative prediction for rat and a comparable prediction for human showing proof of principle on the inclusion of active renal excretion and prediction of PQ acute toxicity for the developed NAM.

Advancing the use of new approach methodologies for assessing teratogenicity: Building a tiered approach

Many New Approach Methodologies (NAMs) have been developed for the safety assessment of new ingredients. Research into reproductive toxicity and teratogenicity is a particularly high priority, especially given their mechanistic complexity. Forty-six non-teratogenic and 39 teratogenic chemicals were screened for teratogenic potential using the in silico DART model from the OECD QSAR Toolbox; the devTox quickPredict™ (devTox assay) test and the Zebrafish Embryotoxicity Test (ZET). The sensitivity and specificity were 94.7% and 84.1%, respectively, for the DART tree (83 chemicals), 86.1% and 35.6% for the devTox (81 chemicals) and 77.8% and 76.7% for the ZET (57 chemicals). Fifty-three chemicals were tested in all three assays and when results were combined and based on a "2 out of 3 rule", the sensitivity and specificity were 96.0% and 71.4%, respectively. The specificity of the devTox assay for a sub-set of 43 chemicals was increased from 26.1% to 82.6% by incorporating human plasma concentrations into the assay interpretation. When all 85 chemicals were assessed in a decision tree approach, there was an excellent predictivity and assay robustness of 90%. In conclusion, all three models exhibited a good sensitivity and specificity, especially when outcomes from all three were combined or used in "2 out of 3" or a tiered decision tree approach. The latter is an interesting predictive approach for evaluating the teratogenic potential of new chemicals. Future investigations will extend the number of chemicals tested, as well as explore ways to refine the results and obtain a robust Integrated Testing Strategy to evaluate teratogenic potential.

Establishing relationships between particle-induced in vitro and in vivo inflammation endpoints to better extrapolate between in vitro markers and in vivo fibrosis

BACKGROUND

Toxicity assessment for regulatory purposes is starting to move away from traditional in vivo methods and towards new approach methodologies (NAM) such as high-throughput in vitro models and computational tools. For materials with limited hazard information, utilising quantitative Adverse Outcome Pathways (AOPs) in a testing strategy involving NAM can produce information relevant for risk assessment. The aim of this work was to determine the feasibility of linking in vitro endpoints to in vivo events, and moreover to key events associated with the onset of a chosen adverse outcome to aid in the development of NAM testing strategies. To do this, we focussed on the adverse outcome pathway (AOP) relating to the onset of pulmonary fibrosis.

RESULTS

We extracted in vivo and in vitro dose-response information for particles known to induce this pulmonary fibrosis (crystalline silica, specifically α-quartz). To test the in vivo-in vitro extrapolation (IVIVE) determined for crystalline silica, cerium dioxide nanoparticles (nano-CeO₂) were used as a case study allowing us to evaluate our findings with a less studied substance. The IVIVE methodology outlined in this paper is formed of five steps, which can be more generally summarised into two categories (i) aligning the in vivo and in vitro dosimetry, (ii) comparing the dose-response curves and derivation of conversion factors.

CONCLUSION

Our analysis shows promising results with regards to correlation of in vitro cytokine secretion to in vivo acute pulmonary inflammation assessed by polymorphonuclear leukocyte influx, most notable is the potential of using IL-6 and IL-1β cytokine secretion from simple in vitro submerged models as a screening tool to assess the likelihood of lung inflammation at an early stage in product development, hence allowing a more targeted investigation using either a smaller, more targeted in vivo study or in the future a more complex in vitro protocol. This paper also highlights the strengths and limitations as well as the current difficulties in performing IVIVE assessment and suggestions for overcoming these issues.

Advancing New Approach Methodologies (NAMs) for Tobacco Harm Reduction: Synopsis from the 2021 CORESTA SSPT-NAMs Symposium

New approach methodologies (NAMs) are emerging chemical safety assessment tools consisting of in vitro and in silico (computational) methodologies intended to reduce, refine, or replace (3R) various in vivo animal testing methods traditionally used for risk assessment. Significant progress has been made toward the adoption of NAMs for human health and environmental toxicity assessment. However, additional efforts are needed to expand their development and their use in regulatory decision making. A virtual symposium was held during the 2021 Cooperation Centre for Scientific Research Relative to Tobacco (CORESTA) Smoke Science and Product Technology (SSPT) conference (titled "Advancing New Alternative Methods for Tobacco Harm Reduction"), with the goals of introducing the concepts and potential application of NAMs in the evaluation of potentially reduced-risk (PRR) tobacco products. At the symposium, experts from regulatory agencies, research organizations, and NGOs shared insights on the status of available tools, strengths, limitations, and opportunities in the application of NAMs using case examples from safety assessments of chemicals and tobacco products. Following seven presentations providing background and application of NAMs, a discussion was held where the presenters and audience discussed the outlook for extending the NAMs toxicological applications for tobacco products. The symposium, endorsed by the CORESTA In Vitro Tox Subgroup, Biomarker Subgroup, and NextG Tox Task Force, illustrated common ground and interest in science-based engagement across the scientific community and stakeholders in support of tobacco regulatory science. Highlights of the symposium are summarized in this paper.

The Evolution of Regulatory Toxicology: Where is the Gardener?

There is a need for paradigm change in the methodology employed for toxicological testing and assessment. It could be said that this change is well on its way, through an evolutionary progress analogous to that of natural selection. Darwin's Theory of Evolution has defined the idea of evolution and descendancy since the last third of the 19th century. Increasingly, this concept of 'evolution' is being applied beyond the field of biology. This Comment article discusses the progress of toxicological testing in the context of 'evolutionary pressure' and deliberates how this process can help foster the development, implementation and acceptance of mechanistic and human-relevant methods in this field. By comparing the current regulatory landscape in toxicity testing and assessment to specific elements in Charles Darwin's evolutionary theory, we aim to better understand the needs and requirements for the future.

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.

Use of physiologically based kinetic modeling-facilitated reverse dosimetry to predict in vivo acute toxicity of tetrodotoxin in rodents

In this study, the ability of a new in vitro/in silico quantitative in vitro–in vivo extrapolation (QIVIVE) methodology was assessed to predict the in vivo neurotoxicity of tetrodotoxin (TTX) in rodents. In vitro concentration–response data of TTX obtained in a multielectrode array assay with primary rat neonatal cortical cells and in an effect study with mouse neuro-2a cells were quantitatively extrapolated into in vivo dose–response data, using newly developed physiologically based kinetic (PBK) models for TTX in rats and mice. Incorporating a kidney compartment accounting for active renal excretion in the PBK models proved to be essential for its performance. To evaluate the predictions, QIVIVE-derived dose–response data were compared with in vivo data on neurotoxicity in rats and mice upon oral and parenteral dosing. The results revealed that for both rats and mice the predicted dose–response data matched the data from available in vivo studies well. It is concluded that PBK modeling-based reserve dosimetry of in vitro TTX effect data can adequately predict the in vivo neurotoxicity of TTX in rodents, providing a novel proof-of-principle for this methodology.

Integrating in vitro chemical transplacental passage into a generic PBK model: A QIVIVE approach

With the increasing application of cell culture models as primary tools for predicting chemical safety, the quantitative extrapolation of the effective dose from in vitro to in vivo (QIVIVE) is of increasing importance. For developmental toxicity this requires scaling the in vitro observed dose-response characteristics to in vivo fetal exposure, while integrating maternal in vivo kinetics during pregnancy, in particular transplacental transfer. Here the transfer of substances across the placental barrier, has been studied using the in vitro BeWo cell assay and six embryotoxic compounds of different kinetic complexity. The BeWo assay results were incorporated in an existing generic Physiologically Based Kinetic (PBK) model which for this purpose was extended with rat pregnancy. Finally, as a "proof of principle", the BeWo PBK model was used to perform a QIVIVE based on developmental toxicity as observed in various different in vitro toxicity assays. The BeWo results illustrated different transport profiles of the chemicals across the BeWo monolayer, allocating the substances into two distinct groups: the 'quickly-transported' and the 'slowly-transported'. BeWo PBK exposure simulations during gestation were compared to experimentally measured maternal blood and fetal concentrations and a reverse dosimetry approach was applied to translate in vitro observed embryotoxicity into equivalent in vivo dose-response curves. This approach allowed for a direct comparison of the in vitro dose-response characteristics as observed in the Whole Embryo Culture (WEC), and the Embryonic Stem Cell test (cardiac:ESTc and neural:ESTn) with in vivo rat developmental toxicity data. Overall, the in vitro to in vivo comparisons suggest a promising future for the application of such QIVIVE methodologies for screening and prioritization purposes of developmental toxicants. Nevertheless, the clear need for further improvements is acknowledged for a wider application of the approach in chemical safety assessment.

A review of pregnancy-induced changes in opioid pharmacokinetics, placental transfer, and fetal exposure: Towards fetomaternal physiologically-based pharmacokinetic modeling to improve the treatment of neonatal opioid withdrawal syndrome

Physiologically-based pharmacokinetic (PBPK) modeling has emerged as a useful tool to study pharmacokinetics (PK) in special populations, such as pregnant women, fetuses, and newborns, where practical hurdles severely limit the study of drug behavior. PK in pregnant women is variable and everchanging, differing greatly from that in their nonpregnant female and male counterparts typically enrolled in clinical trials. PBPK models can accommodate pregnancy-induced physiological and metabolic changes, thereby providing mechanistic insights into maternal drug disposition and fetal exposure. Fueled by the soaring opioid epidemic in the United States, opioid use during pregnancy continues to rise, leading to an increased incidence of neonatal opioid withdrawal syndrome (NOWS). The severity of NOWS is influenced by a complex interplay of extrinsic and intrinsic factors, and varies substantially between newborns, but the extent of prenatal opioid exposure is likely the primary driver. Fetomaternal PBPK modeling is an attractive approach to predict in utero opioid exposure. To facilitate the development of fetomaternal PBPK models of opioids, this review provides a detailed overview of pregnancy-induced changes affecting the PK of commonly used opioids during gestation. Moreover, the placental transfer of these opioids is described, along with their disposition in the fetus. Lastly, the implementation of these factors into PBPK models is discussed. Fetomaternal PBPK modeling of opioids is expected to provide improved insights in fetal opioid exposure, which allows for prediction of postnatal NOWS severity, thereby opening the way for precision postnatal treatment of these vulnerable infants.

Predicting the in vivo developmental toxicity of benzo[a]pyrene (BaP) in rats by an in vitro-in silico approach

Developmental toxicity testing is an animal-intensive endpoints in toxicity testing and calls for animal-free alternatives. Previous studies showed the applicability of an in vitro–in silico approach for predicting developmental toxicity of a range of compounds, based on data from the mouse embryonic stem cell test (EST) combined with physiologically based kinetic (PBK) modelling facilitated reverse dosimetry. In the current study, the use of this approach for predicting developmental toxicity of polycyclic aromatic hydrocarbons (PAHs) was evaluated, using benzo[a]pyrene (BaP) as a model compound. A rat PBK model of BaP was developed to simulate the kinetics of its main metabolite 3-hydroxybenzo[a]pyrene (3-OHBaP), shown previously to be responsible for the developmental toxicity of BaP. Comparison to in vivo kinetic data showed that the model adequately predicted BaP and 3-OHBaP blood concentrations in the rat. Using this PBK model and reverse dosimetry, a concentration–response curve for 3-OHBaP obtained in the EST was translated into an in vivo dose–response curve for developmental toxicity of BaP in rats upon single or repeated dose exposure. The predicted half maximal effect doses (ED₅₀) amounted to 67 and 45 mg/kg bw being comparable to the ED₅₀ derived from the in vivo dose–response data reported for BaP in the literature, of 29 mg/kg bw. The present study provides a proof of principle of applying this in vitro–in silico approach for evaluating developmental toxicity of BaP and may provide a promising strategy for predicting the developmental toxicity of related PAHs, without the need for extensive animal testing.

Establishment of an in vitro method of rabbit embryo toxicity with toxicokinetics study

This report introduces a novel method, rabbit whole embryo culture (WEC) combined with toxicokinetics (TK), for toxicity testing. Rodent WEC has been extensively used for in vitro screening of developmental toxicity. To improve the reliability of in vitro data, it is important to consider TK and species specificity. To test the utility and effectiveness of this method, we investigated the toxic effect of thalidomide on rabbit embryos and its behavior in test systems both in vitro and in vivo under the same experimental condition. The data showed that thalidomide induced embryo malformations such as embryonic brain hypoplasia, short limb buds, and declined embryonic growth both in vitro and in vivo. The toxic effect increased with the increasing exposure of the embryo to thalidomide. In addition, we observed similar toxic effects and exposure-effect relationships in vivo and in vitro. Therefore, we preliminarily conclude that this new method can effectively predict and explain thalidomide toxicity. Furthermore, we investigated the behavior of test compounds in the WEC system for the first time, and this method is expected to be an important technique for in vitro toxicity study after extensive verification.

Development of a Generic Physiologically Based Kinetic Model to Predict In Vivo Uterotrophic Responses Induced by Estrogenic Chemicals in Rats Based on In Vitro Bioassays

The present study assessed the potential of a generic physiologically based kinetic (PBK) model to convert in vitro data for estrogenicity to predict the in vivo uterotrophic response in rats for diethylstibestrol (DES), ethinylestradiol (EE2), genistein (GEN), coumestrol (COU), and methoxychlor (MXC). PBK models were developed using a generic approach and in vitro concentration-response data from the MCF-7 proliferation assay and the yeast estrogen screening assay were translated into in vivo dose-response data. Benchmark dose analysis was performed on the predicted data and available in vivo uterotrophic data to evaluate the model predictions. The results reveal that the developed generic PBK model adequate defines the in vivo kinetics of the estrogens. The predicted dose-response data of DES, EE2, GEN, COU, and MXC matched the reported in vivo uterus weight response in a qualitative way, whereas the quantitative comparison was somewhat hampered by the variability in both in vitro and in vivo data. From a safety perspective, the predictions based on the MCF-7 proliferation assay would best guarantee a safe point of departure for further risk assessment although it may be conservative. The current study indicates the feasibility of using a combination of in vitro toxicity data and a generic PBK model to predict the relative in vivo uterotrophic response for estrogenic chemicals.

Derivation of a Human In Vivo Benchmark Dose for Perfluorooctanoic Acid From ToxCast In Vitro Concentration-Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation

A computational workflow which integrates physiologically based kinetic (PBK) modeling, global sensitivity analysis (GSA), approximate Bayesian computation (ABC), and Markov Chain Monte Carlo (MCMC) simulation was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow accounts for parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for perfluorooctanoic acid (PFOA) and high throughput screening (HTS) in vitro concentration–response data, determined in a human liver cell line, from the ToxCast/Tox21 database. In vivo benchmark doses (BMDs) for PFOA intake (ng/kg BW/day) and drinking water exposure concentrations (µg/L) were calculated from the in vivo dose responses and compared to intake values derived by the European Food Safety Authority (EFSA). The intake benchmark dose lower confidence limit (BMDL₅) of 0.82 was similar to 0.86 ng/kg BW/day for altered serum cholesterol levels derived by EFSA, whereas the intake BMDL₅ of 6.88 was six-fold higher than the value of 1.14 ng/kg BW/day for altered antibody titer also derived by the EFSA. Application of a chemical-specific adjustment factor (CSAF) of 1.4, allowing for inter-individual variability in kinetics, based on biological half-life, gave an intake BMDL₅ of 0.59 for serum cholesterol and 4.91 (ng/kg BW/day), for decreased antibody titer, which were 0.69 and 4.31 the EFSA-derived values, respectively. The corresponding BMDL₅ for drinking water concentrations, for estrogen receptor binding activation associated with breast cancer, pregnane X receptor binding associated with altered serum cholesterol levels, thyroid hormone receptor α binding leading to thyroid disease, and decreased antibody titer (pro-inflammation from cytokines) were 0.883, 0.139, 0.086, and 0.295 ng/ml, respectively, with application of no uncertainty factors. These concentrations are 5.7-, 36-, 58.5-, and 16.9-fold lower than the median measured drinking water level for the general US population which is approximately, 5 ng/ml.

Physiologically based kinetic modelling based prediction of in vivo rat and human acetylcholinesterase (AChE) inhibition upon exposure to diazinon

The present study predicts in vivo human and rat red blood cell (RBC) acetylcholinesterase (AChE) inhibition upon diazinon (DZN) exposure using physiological based kinetic (PBK) modelling-facilitated reverse dosimetry. Due to the fact that both DZN and its oxon metabolite diazoxon (DZO) can inhibit AChE, a toxic equivalency factor (TEF) was included in the PBK model to combine the effect of DZN and DZO when predicting in vivo AChE inhibition. The PBK models were defined based on kinetic constants derived from in vitro incubations with liver fractions or plasma of rat and human, and were used to translate in vitro concentration-response curves for AChE inhibition obtained in the current study to predicted in vivo dose-response curves. The predicted dose-response curves for rat matched available in vivo data on AChE inhibition, and the benchmark dose lower confidence limits for 10% inhibition (BMDL10 values) were in line with the reported BMDL10 values. Humans were predicted to be 6-fold more sensitive than rats in terms of AChE inhibition, mainly because of inter-species differences in toxicokinetics. It is concluded that the TEF-coded DZN PBK model combined with quantitative in vitro to in vivo extrapolation (QIVIVE) provides an adequate approach to predict RBC AChE inhibition upon acute oral DZN exposure, and can provide an alternative testing strategy for derivation of a point of departure (POD) in risk assessment.

Incorporating renal excretion via the OCT2 transporter in physiologically based kinetic modelling to predict in vivo kinetics of mepiquat in rat

The present study aimed at incorporating active renal excretion via the organic cation transporter 2 (OCT2) into a generic rat physiologically based kinetic (PBK) model using an in vitro human renal proximal tubular epithelial cell line (SA7K) and mepiquat chloride (MQ) as the model compound. The Vmax (10.5 pmol/min/mg protein) and Km (20.6 μM) of OCT2 transport of MQ were determined by concentration-dependent uptake in SA7K cells using doxepin as inhibitor. PBK model predictions incorporating these values in the PBK model were 6.7-8.4-fold different from the reported in vivo data on the blood concentration of MQ in rat. Applying an overall scaling factor that also corrects for potential differences in OCT2 activity in the SA7K cells and in vivo kidney cortex and species differences resulted in adequate predictions for in vivo kinetics of MQ in rat (2.3-3.2-fold). The results indicate that using SA7K cells to define PBK parameters for active renal OCT2 mediated excretion with adequate scaling enables incorporation of renal excretion via the OCT2 transporter in PBK modelling to predict in vivo kinetics of mepiquat in rat. This study demonstrates a proof-of-principle on how to include active renal excretion into generic PBK models.

Novel testing strategy for prediction of rat biliary excretion of intravenously administered estradiol-17β glucuronide

The aim of the present study was to develop a generic rat physiologically based kinetic (PBK) model that includes a novel testing strategy where active biliary excretion is incorporated using estradiol-17β glucuronide (E217βG) as the model substance. A major challenge was the definition of the scaling factor for the in vitro to in vivo conversion of the PBK-model parameter Vmax. In vitro values for the Vmax and Km for transport of E217βG were found in the literature in four different studies based on experiments with primary rat hepatocytes. The required scaling factor was defined based on fitting the PBK model-based predicted values to reported experimental data on E217βG blood levels and cumulative biliary E217βG excretion. This resulted in a scaling factor of 129 mg protein/g liver. With this scaling factor the PBK model predicted the in vivo data for blood and cumulative biliary E217βG levels with on average of less than 1.8-fold deviation. The study provides a proof of principle on how biliary excretion can be included in a generic PBK model using primary hepatocytes to define the kinetic parameters that describe the biliary excretion.

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.

Integrating biokinetics and in vitro studies to evaluate developmental neurotoxicity induced by chlorpyrifos in human iPSC-derived neural stem cells undergoing differentiation towards neuronal and glial cells

•

Human iPSC-derived NSCs undergoing differentiation possess some metabolic competence.

•

CPF entered the cells and was biotrasformed into its two main metabolites (CPFO and TCP).

•

After repeated exposure, very limited bioaccumulation of CPF was observed.

•

Treatment with CPF decreased neurite outgrowth, synapse number and electrical activity.

•

Treatment with CPF increased BDNF levels and the percentage of astrocytes.

For some complex toxicological endpoints, chemical safety assessment has conventionally relied on animal testing. Apart from the ethical issues, also scientific considerations have been raised concerning the traditional approach, highlighting the importance for considering real life exposure scenario. Implementation of flexible testing strategies, integrating multiple sources of information, including in vitro reliable test methods and in vitro biokinetics, would enhance the relevance of the obtained results. Such an approach could be pivotal in the evaluation of developmental neurotoxicity (DNT), especially when applied to human cell-based models, mimicking key neurodevelopmental processes, relevant to human brain development.

Here, we integrated the kinetic behaviour with the toxicodynamic alterations of chlorpyrifos (CPF), such as in vitro endpoints specific for DNT evaluation, after repeated exposure during differentiation of human neural stem cells into a mixed culture of neurons and astrocytes. The upregulation of some cytochrome P450 and glutathione S-transferase genes during neuronal differentiation and the formation of the two major CPF metabolites (due to bioactivation and detoxification) supported the metabolic competence of the used in vitro model. The alterations in the number of synapses, neurite outgrowth, brain derived neurotrophic factor, the proportion of neurons and astrocytes, as well as spontaneous electrical activity correlated well with the CPF ability to enter the cells and be bioactivated to CPF-oxon. Overall, our results confirm that combining in vitro biokinetics and assays to evaluate effects on neurodevelopmental endpoints in human cells should be regarded as a key strategy for a quantitative characterization of DNT effects.

Predicting the Acute Liver Toxicity of Aflatoxin B1 in Rats and Humans by an In Vitro-In Silico Testing Strategy

SCOPE

High-level exposure to aflatoxin B1 (AFB1) is known to cause acute liver damage and fatality in animals and humans. The intakes actually causing this acute toxicity have so far been estimated based on AFB1 levels in contaminated foods or biomarkers in serum. The aim of the present study is to predict the doses causing acute liver toxicity of AFB1 in rats and humans by an in vitro-in silico testing strategy.

METHODS AND RESULTS

Physiologically based kinetic (PBK) models for AFB1 in rats and humans are developed. The models are used to translate in vitro concentration-response curves for cytotoxicity in primary rat and human hepatocytes to in vivo dose-response curves using reverse dosimetry. From these data, the dose levels at which toxicity would be expected are obtained and compared to toxic dose levels from available rat and human case studies on AFB1 toxicity. The results show that the in vitro-in silico testing strategy can predict dose levels causing acute toxicity of AFB1 in rats and human.

CONCLUSIONS

Quantitative in vitro in vivo extrapolation (QIVIVE) using PBK modeling-based reverse dosimetry can predict AFB1 doses that cause acute liver toxicity in rats and human.

Source apportionment of phenolic compounds based on a simultaneous monitoring of surface water and emission sources: A case study in a typical region adjacent to Taihu Lake watershed

A comprehensive characterisation of four phenolic compounds in surface water and various emission sources was conducted simultaneously in a typical industrial city upstream of the Taihu Lake watershed. The overall concentrations of the target compounds ranged from 0.15 to 2.75 μg/L in the surface water and from 0.43 to 377 μg/L in the emission sources. Relatively high concentration levels were observed in August, during which the rainy season typically occurred in the study area, indicating seasonal emission sources. The spatial distribution revealed severe phenolic-compound contamination in the northeast part of the study area. According to the relationships between the surface water and emission sources, combined-sewer overflow and surface runoff from agricultural activities were identified as seasonal emission sources that were responsible for the increased concentration levels in wet seasons. Indirect discharge of industrial wastewater was proposed as a persistent emission source responsible for the severe contamination level in the northeast region over the entire sampling period, accounting for 16.5% of the phenolic-compound load in a demonstration river section. The findings of this study are useful for identifying the sources of phenolic compounds and controlling the contamination from the main sources in a typical industrial city.

Integrating in vitro data and physiologically based kinetic modeling-facilitated reverse dosimetry to predict human cardiotoxicity of methadone

Development of novel testing strategies to detect adverse human health effects is of interest to replace in vivo-based drug and chemical safety testing. The aim of the present study was to investigate whether physiologically based kinetic (PBK) modeling-facilitated conversion of in vitro toxicity data is an adequate approach to predict in vivo cardiotoxicity in humans. To enable evaluation of predictions made, methadone was selected as the model compound, being a compound for which data on both kinetics and cardiotoxicity in humans are available. A PBK model for methadone in humans was developed and evaluated against available kinetic data presenting an adequate match. Use of the developed PBK model to convert concentration–response curves for the effect of methadone on human-induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CM) in the so-called multi electrode array (MEA) assay resulted in predictions for in vivo dose–response curves for methadone-induced cardiotoxicity that matched the available in vivo data. The results also revealed differences in protein plasma binding of methadone to be a potential factor underlying variation between individuals with respect to sensitivity towards the cardiotoxic effects of methadone. The present study provides a proof-of-principle of using PBK modeling-based reverse dosimetry of in vitro data for the prediction of cardiotoxicity in humans, providing a novel testing strategy in cardiac safety studies.

Liver toxicity of anthraquinones: A combined in vitro cytotoxicity and in silico reverse dosimetry evaluation

Anthraquinones are found in a variety of consumer products such as dietary supplements, traditional Chinese medicines, and drugs. Along with their widespread use, potential safety concerns have emerged, especially liver toxicity. Therefore, there is a need to conduct rapid and inexpensive safety assessment for anthraquinones due to a lack of animal and human toxicological data. Here, a combined in vitro cytotoxicity and in silico reverse dosimetry approach was adopted to consider the potential human liver toxicity of 16 anthraquinones and derivatives. First, cytotoxicity (EC₅₀) in two human liver cell lines (HepG2/C3A and HuH-7) was measured under two conditions (single and repeated dosing, 72 h). Second, toxic doses (D_tox) required to yield plasma steady-state concentrations (C_ss) equal to in vitro EC₅₀ values were predicted by reverse dosimetry simulation using a PBPK model. Finally, D_tox was compared to literature-derived estimated daily intake (EDI) of anthraquinones to assess safety. Among the 16 anthraquinones, rhein was identified as a potential hepatotoxicant due to a combination of cytotoxicity, plasma concentration, and daily intake level. These in vitro and in silico findings provide preliminary data and guidance for further animal and clinical studies to confirm liver toxicity of anthraquinones.

Use of Physiologically Based Kinetic Modeling to Predict Rat Gut Microbial Metabolism of the Isoflavone Daidzein to S-Equol and Its Consequences for ERα Activation

SCOPE

To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds.

METHODS AND RESULTS

Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported.

CONCLUSION

The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.

Integrating in vitro testing and physiologically-based pharmacokinetic (PBPK) modelling for chemical liver toxicity assessment-A case study of troglitazone

In vitro to in vivo extrapolation (IVIVE) for next-generation risk assessment (NGRA) of chemicals requires computational modeling and faces unique challenges. Using mitochondria-related toxicity data of troglitazone (TGZ), a prototype drug known for liver toxicity, from HepaRG, HepG2, HC-04, and primary human hepatocytes, we explored inherent uncertainties in IVIVE, including cell models, cellular response endpoints, and dose metrics. A human population physiologically-based pharmacokinetic (PBPK) model for TGZ was developed to predict in vivo doses from in vitro point-of-departure (POD) concentrations. Compared to the 200-800 mg/d dose range of TGZ where liver injury was observed clinically, the predicted POD doses for the mean and top one percentile of the PBPK population were 28-372 and 15-178 mg/d respectively based on C_max dosimetry, and 185-2552 and 83-1010 mg/d respectively based on AUC. In conclusion, although with many uncertainties, integrating in vitro assays and PBPK modeling is promising in informing liver toxicity-inducing TGZ doses.

Combining In Vitro Data and Physiologically Based Kinetic Modeling Facilitates Reverse Dosimetry to Define In Vivo Dose-Response Curves for Bixin- and Crocetin-Induced Activation of PPARγ in Humans

SCOPE

It is investigated whether at realistic dietary intake bixin and crocetin could induce peroxisome proliferator-activated receptor γ (PPARγ)-mediated gene expression in humans using a combined in vitro-in silico approach.

METHODS AND RESULTS

Concentration-response curves obtained from in vitro PPARγ-reporter gene assays are converted to in vivo dose-response curves using physiologically based kinetic modeling-facilitated reverse dosimetry, from which the benchmark dose levels resulting in a 50% effect above background level (BMD50 ) are predicted and subsequently compared to dietary exposure levels. Bixin and crocetin activated PPARγ-mediated gene transcription in a concentration-dependent manner with similar potencies. Due to differences in kinetics, the predicted BMD50 values for in vivo PPARγ activation are about 30-fold different, amounting to 115 and 3505 mg kg bw-1 for crocetin and bixin, respectively. Human dietary and/or supplemental estimated daily intakes may reach these BMD50 values for crocetin but not for bixin, pointing at better possibilities for in vivo PPARγ activation by crocetin.

CONCLUSION

Based on a combined in vitro-in silico approach, it is estimated whether at realistic dietary intakes plasma concentrations of bixin and crocetin are likely to reach concentrations that activate PPARγ-mediated gene expression, without the need for a human intervention study.

Biotransformation and bioactivation reactions - 2018 literature highlights

In the past three decades, ADME sciences have become an integral component of the drug discovery and development process. At the same time, the field has continued to evolve, thus, requiring ADME scientists to be knowledgeable of and engage with diverse aspects of drug assessment: from pharmacology to toxicology, and from in silico modeling to in vitro models and finally in vivo models. Progress in this field requires deliberate exposure to different aspects of ADME; however, this task can seem daunting in the current age of mass information. We hope this review provides a focused and brief summary of a wide array of critical advances over the past year and explains the relevance of this research ( Table 1 ). We divided the articles into categories of (1) drug optimization, (2) metabolites and drug metabolizing enzymes, and (3) bioactivation. This annual review is the fourth of its kind (Baillie et al. 2016 ; Khojasteh et al. 2017 , 2018 ). We have followed the same format we used in previous years in terms of the selection of articles and the authoring of each section. This effort in itself also continues to evolve. I am pleased that Rietjens, Miller, and Mitra have again contributed to this annual review. We would like to welcome Namandjé N. Bumpus, James P. Driscoll, and Donglu Zhang as authors for this year's issue. We strive to maintain a balance of authors from academic and industry settings. We would be pleased to hear your opinions of our commentary, and we extend an invitation to anyone who would like to contribute to a future edition of this review. Cyrus Khojasteh, on behalf of the authors.

The in vivo developmental toxicity of diethylstilbestrol (DES) in rat evaluated by an alternative testing strategy

In the present study, we evaluated an alternative testing strategy to quantitatively predict the in vivo developmental toxicity of the synthetic hormone diethylstilbestrol (DES). To this end, a physiologically based kinetic (PBK) model was defined that was subsequently used to translate concentration-response data for the in vitro developmental toxicity of DES, obtained in the ES-D3 cell differentiation assay, into predicted in vivo dose-response data for developmental toxicity. The previous studies showed that the PBK model-facilitated reverse dosimetry approach is a useful approach to quantitatively predict the developmental toxicity of several developmental toxins. The results obtained in the present study show that the PBK model adequately predicted DES blood concentrations in rats. Further studies revealed that DES tested positive in the ES-D3 differentiation assay and that DES-induced inhibition of the ES-D3 cell differentiation could be counteracted by the estrogen receptor alpha (ERα) antagonist fulvestrant, indicating that the in vitro ES-D3 cell differentiation assay was able to mimic the role of ERα reported in the mode of action underlying the developmental toxicity of DES in vivo. In spite of this, combining these in vitro data with the PBK model did not adequately predict the in vivo developmental toxicity of DES in a quantitative way. It is concluded that although the EST qualifies DES as a developmental toxin and detects the role of ERα in this process, the ES-D3 cell differentiation assay of the EST apparently does not adequately capture the processes underlying DES-induced developmental toxicity in vivo.

A Computational Workflow for Probabilistic Quantitative in Vitro to in Vivo Extrapolation

A computational workflow was developed to facilitate the process of quantitative in vitro to in vivo extrapolation (QIVIVE), specifically the translation of in vitro concentration-response to in vivo dose-response relationships and subsequent derivation of a benchmark dose value (BMD). The workflow integrates physiologically based pharmacokinetic (PBPK) modeling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC) simulation. For a given set of in vitro concentration and response data the algorithm returns the posterior distribution of the corresponding in vivo, population-based dose-response values, for a given route of exposure. The novel aspect of the workflow is a rigorous statistical framework for accommodating uncertainty in both the parameters of the PBPK model (both parameter uncertainty and population variability) and in the structure of the PBPK model itself recognizing that the model is an approximation to reality. Both these sources of uncertainty propagate through the workflow and are quantified within the posterior distribution of in vivo dose for a fixed representative in vitro concentration. To demonstrate this process and for comparative purposes a similar exercise to previously published work describing the kinetics of ethylene glycol monoethyl ether (EGME) and its embryotoxic metabolite methoxyacetic acid (MAA) in rats was undertaken. The computational algorithm can be used to extrapolate from in vitro data to any organism, including human. Ultimately, this process will be incorporated into a user-friendly, freely available modeling platform, currently under development, that will simplify the process of QIVIVE.

Toxicokinetics in Risk Evaluations

Quantitative predictions of in vivo chemical levels based on in vitro data will become a cornerstone of next generation nonanimal risk evaluations. Both regulatory and scientific experience with quantitative toxicokinetics must increase now for this transition to happen.

Use of physiologically based kinetic modelling-facilitated reverse dosimetry to convert in vitro cytotoxicity data to predicted in vivo liver toxicity of lasiocarpine and riddelliine in rat

Lasiocarpine and riddelliine are pyrrolizidine alkaloids (PAs) present in food and able to cause liver toxicity. The aim of this study was to investigate whether physiologically based kinetic (PBK) modelling-facilitated reverse dosimetry can adequately translate in vitro concentration-response curves for toxicity of lasiocarpine and riddelliine to in vivo liver toxicity data for the rat. To this purpose, PBK models were developed for lasiocarpine and riddelliine, and predicted blood concentrations were compared to available literature data to evaluate the models. Concentration-response curves obtained from in vitro cytotoxicity assays in primary rat hepatocytes were converted to in vivo dose-response curves from which points of departure (PODs) were derived and that were compared to available literature data on in vivo liver toxicity. The results showed that the predicted PODs fall well within the range of PODs derived from available in vivo toxicity data. To conclude, this study shows the proof-of-principle for a method to predict in vivo liver toxicity for PAs by an alternative testing strategy integrating in vitro cytotoxicity assays with in silico PBK modelling-facilitated reverse dosimetry. The approach may facilitate prediction of acute liver toxicity for the large number of PAs for which in vivo toxicity data are lacking.

Towards a generic physiologically based kinetic model to predict in vivo uterotrophic responses in rats by reverse dosimetry of in vitro estrogenicity data

Physiologically based kinetic (PBK) modelling-based reverse dosimetry is a promising tool for the prediction of in vivo developmental toxicity using in vitro concentration–response data. In the present study, the potential of this approach to predict the dose-dependent increase of uterus weight in rats upon exposure to estrogenic chemicals was assessed. In vitro concentration–response data of 17β-estradiol (E2) and bisphenol A (BPA) obtained in the MCF-7/BOS proliferation assay, the U2OS ER-CALUX assay and the yeast estrogen screen (YES) assay, were translated into in vivo dose–response data in rat, using a PBK model with a minimum number of in vitro and in silico determined parameter values. To evaluate the predictions made, benchmark dose (BMD) analysis was performed on the predicted dose–response data and the obtained BMDL₁₀ values were compared with BMDL₁₀ values derived from data on the effects of E2 and BPA in the uterotrophic assay reported in the literature. The results show that predicted dose–response data of E2 and BPA matched with the data from in vivo studies when predictions were made based on YES assay data. The YES assay-based predictions of the BMDL₁₀ values differed 3.9-fold (E2) and 4.7- to 13.4-fold (BPA) from the BMDL₁₀ values obtained from the in vivo data. The present study provides the proof-of-principle that PBK modelling-based reverse dosimetry of YES assay data using a minimum PBK model can predict dose-dependent in vivo uterus growth caused by estrogenic chemicals. In future studies, the approach should be extended to include other estrogens.