Aloe-emodin, a natural hydroxyanthraquinone, exerts both adverse and protective effects. This study aimed at investigating these potential effects of aloe-emodin in humans upon the use of food supplements and herbal medicines using a physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach. For this, PBK models in rats and humans were established for aloe-emodin including its active metabolite rhein and used to convert in vitro data on hepatotoxicity, nephrotoxicity, reactive oxidative species (ROS) generation, and Nrf2 induction to corresponding in vivo dose-response curves, from which points of departure (PODs) were derived by BMD analysis. The derived PODs were subsequently compared to the estimated daily intakes (EDIs) resulting from the use of food supplements or herbal medicines. It is concluded that the dose levels of aloe-emodin from food supplements or herbal medicines are unlikely to induce toxicity, ROS generation, or Nrf2 activation in liver and kidney.

Worldwide use of organophosphate pesticides as agricultural chemicals aims to maintain a stable food supply, while their toxicity remains a major public health concern. A common mechanism of acute neurotoxicity following organophosphate pesticide exposure is the inhibition of acetylcholinesterase (AChE). To support Next Generation Risk Assessment for public health upon acute neurotoxicity induced by organophosphate pesticides, physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach was employed in this study, with fenitrothion (FNT) as an exemplary organophosphate pesticide. Rat and human PBK models were parametrized with data derived from in silico predictions and in vitro incubations. Then, PBK model-based QIVIVE was performed to convert species-specific concentration-dependent AChE inhibition obtained from in vitro blood assays to corresponding in vivo dose-response curves, from which points of departure (PODs) were derived. The obtained values for rats and humans were comparable with reported no-observed-adverse-effect levels (NOAELs). Humans were found to be more susceptible than rats toward erythrocyte AChE inhibition induced by acute FNT exposure due to interspecies differences in toxicokinetics and toxicodynamics. The described approach adequately predicts toxicokinetics and acute toxicity of FNT, providing a proof-of-principle for applying this approach in a 3R-based chemical risk assessment paradigm.

Current climate trends are likely to expand the geographic distribution of the toxigenic microalgae and concomitant phycotoxins, making intoxications by such toxins a global phenomenon. Among various phycotoxins, saxitoxin (STX) acts as a neurotoxin that might cause severe neurological symptoms in mammals following consumptions of contaminated seafood. To derive a point of departure (POD) for human health risk assessment upon acute neurotoxicity induced by oral STX exposure, a physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach was employed. The PBK models for rats, mice, and humans were built using parameters from the literature, in vitro experiments, and in silico predictions. Available in vitro toxicity data for STX were converted to in vivo dose-response curves via the PBK models established for these three species, and POD values were derived from the predicted curves and compared to reported in vivo toxicity data. Interspecies differences in acute STX toxicity between rodents and humans were found, and they appeared to be mainly due to differences in toxicokinetics. The described approach resulted in adequate predictions for acute oral STX exposure, indicating that new approach methodologies, when appropriately integrated, can be used in a 3R-based chemical risk assessment paradigm.

Perfluoroalkyl substances (PFAS) are persistent and pose a risk to human health. High throughput screening (HTS) cell-based bioassays may inform risk assessment of PFAS provided that quantitative in vitro to in vivo extrapolation (QIVIVE) can be developed. The QIVIVE ratio is the ratio of nominal (Cnom) or freely dissolved concentration (Cfree) in human blood to Cnom or Cfree in the bioassays. Considering that the concentrations of PFAS in human plasma and in vitro bioassays may vary by orders of magnitude, we tested the hypothesis that anionic PFAS bind to proteins concentration-dependently and therefore the binding differs substantially between human plasma and bioassays, which has an impact on QIVIVE. Solid phase microextraction (SPME) with C18-coated fibers served to quantify the Cfree of four anionic PFAS (perfluorobutanoate (PFBA), perfluorooctanoate (PFOA), perfluorohexane sulfonate (PFHxS) and perfluorooctane sulfonate (PFOS)) in the presence of proteins and lipid, medium components, cells and human plasma over five orders of magnitude in concentrations. The C18-SPME method was used to quantify the non-linear binding to proteins, human plasma and medium, and the partition constants to cells. These binding parameters were used to predict Cfree of PFAS in cell bioassays and human plasma by a concentration-dependent mass balance model (MBM). The approach was illustrated with a reporter gene assay indicating activation of the peroxisome proliferator-activated receptor gamma (PPARγ-GeneBLAzer). Blood plasma levels were collected from literature for occupational exposure and the general population. The QIVIVEnom ratios were higher than the QIVIVEfree ratios due to the strong affinity to proteins and large differences in protein contents between human blood and bioassays. For human health risk assessment, the QIVIVEfree ratios of many in vitro assays need to be combined to cover all health relevant endpoints. If Cfree cannot be measured, they can be estimated with the MBM and concentration-dependent distribution ratios.

New approach methodologies predicting human cardiotoxicity are of interest to support or even replace in vivo-based drug safety testing. The present study presents an in vitro-in silico approach to predict the effect of inter-individual and inter-ethnic kinetic variations in the cardiotoxicity of R- and S-methadone in the Caucasian and the Chinese population. In vitro cardiotoxicity data, and metabolic data obtained from two approaches, using either individual human liver microsomes or recombinant cytochrome P450 enzymes (rCYPs), were integrated with physiologically based kinetic (PBK) models and Monte Carlo simulations to predict inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. Chemical specific adjustment factors were defined and used to derive dose-response curves for the sensitive individuals. Our simulations indicated that Chinese are more sensitive towards methadone-induced cardiotoxicity with Margin of Safety values being generally two-fold lower than those for Caucasians for both methadone enantiomers. Individual PBK models using microsomes and PBK models using rCYPs combined with Monte Carlo simulations predicted similar inter-individual and inter-ethnic variations in methadone-induced cardiotoxicity. The present study illustrates how inter-individual and inter-ethnic variations in cardiotoxicity can be predicted by combining in vitro toxicity and metabolic data, PBK modelling and Monte Carlo simulations. The novel methodology can be used to enhance cardiac safety evaluations and risk assessment of chemicals.

Physiologically-based kinetic (PBK) models can simulate concentrations of chemicals in tissues over time without animal experiments. Nevertheless, in vivo data are often used to parameterise PBK models. This study aims to illustrate that a combination of kinetic and dynamic readouts from in vitro assays can be used to parameterise PBK models simulating neurologically-active concentrations of xenobiotics. Baclofen, an intrathecally administered drug to treat spasticity, was used as a proof-of-principle xenobiotic. An in vitro blood-brain barrier (BBB) model was used to determine the BBB permeability of baclofen needed to simulate plasma and cerebrospinal concentrations. Simulated baclofen concentrations in individuals and populations of adults and children generally fall within 2-fold of measured clinical study concentrations. Further, in vitro micro-electrode array recordings were used to determine the effect of baclofen on neuronal activity (cell signalling). Using quantitative in vitro-in vivo extrapolations (QIVIVE) corresponding doses of baclofen were estimated. QIVIVE showed that up to 4600 times lower intrathecal doses than oral and intravenous doses induce comparable neurological effects. Most simulated doses were in the range of administered doses. This show that PBK models predict concentrations in the central nervous system for various routes of administration accurately without the need for additional in vivo data.

The development of non-animal-based new approach methodologies (NAMs) for chemical risk assessment and safety evaluation is urgently needed. The aim of the present study was to investigate the applicability of an in vitro-in silico approach to predict human cardiotoxicity of the herbal alkaloid ibogaine and its metabolite noribogaine, which are promising anti-addiction drugs. Physiologically based kinetic (PBK) models were developed using in silico-derived parameters and biokinetic data obtained from in vitro liver microsomal incubations and Caco-2 transport studies. Human induced pluripotent stem cell-derived cardiomyocytes combined with a multi-electrode array (MEA) assay were used to determine in vitro concentration-dependent cardiotoxicity reflected by prolongation of field potential duration, which was subsequently translated to in vivo dose-dependent prolongation of the QTc (heart rate corrected duration from ventricular depolarization to repolarization) using PBK modeling-based reverse dosimetry. Results showed that the predictions matched well with in vivo kinetic data and QTc data for ibogaine and noribogaine available in the literature, indicating a good performance of the NAM. Benchmark dose analysis of the predicted dose response curves adequately predicted the onset of in vivo cardiotoxicity detected by QTc prolongation upon oral exposure to ibogaine and noribogaine. The present study provides an additional proof-of-principle of using PBK modeling-based reverse dosimetry as a NAM to predict human cardiotoxicity.

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.