UDP-glucuronosyltransferases (UGTs) are involved in phase II conjugation reactions of xenobiotics and differences in their isoform activities result in interindividual kinetic differences of UGT probe substrates. Here, extensive literature searches were performed to identify probe substrates (14) for various UGT isoforms (UGT1A1, UGT1A3, UGT1A4, UGT1A6, UGT1A9, UGT2B7 and UGT2B15) and frequencies of human polymorphisms. Chemical-specific pharmacokinetic data were collected in a database to quantify interindividual differences in markers of acute (Cmax) and chronic (area under the curve, clearance) exposure. Using this database, UGT-related uncertainty factors were derived and compared to the default factor (i.e. 3.16) allowing for interindividual differences in kinetics. Overall, results show that pharmacokinetic data are predominantly available for Caucasian populations and scarce for other populations of different geographical ancestry. Furthermore, the relationships between UGT polymorphisms and pharmacokinetic parameters are rarely addressed in the included studies. The data show that UGT-related uncertainty factors were mostly below the default toxicokinetic uncertainty factor of 3.16, with the exception of five probe substrates (1-OH-midazolam, ezetimibe, raltegravir, SN38 and trifluoperazine), with three of these substrates being metabolised by the polymorphic isoform 1A1. Data gaps and future work to integrate UGT-related variability distributions with in vitro data to develop quantitative in vitro–in vivo extrapolations in chemical risk assessment are discussed.

With an increasing need to incorporate new approach methodologies (NAMs) in chemical risk assessment and the concomitant need to phase out animal testing, the interpretation of in vitro assay readouts for quantitative hazard characterisation becomes more important. Physiologically based kinetic (PBK) models, which simulate the fate of chemicals in tissues of the body, play an essential role in extrapolating in vitro effect concentrations to in vivo bioequivalent exposures. As PBK-based testing approaches evolve, it will become essential to standardise PBK modelling approaches towards a consensus approach that can be used in quantitative in vitro-to-in vivo extrapolation (QIVIVE) studies for regulatory chemical risk assessment based on in vitro assays. Based on results of an ECETOC expert workshop, steps are recommended that can improve regulatory adoption: (1) define context and implementation, taking into consideration model complexity for building fit-for-purpose PBK models, (2) harmonise physiological input parameters and their distribution and define criteria for quality chemical-specific parameters, especially in the absence of in vivo data, (3) apply Good Modelling Practices (GMP) to achieve transparency and design a stepwise approach for PBK model development for risk assessors, (4) evaluate model predictions using alternatives to in vivo PK data including read-across approaches, (5) use case studies to facilitate discussions between modellers and regulators of chemical risk assessment. Proof-of-concepts of generic PBK modelling approaches are published in the scientific literature at an increasing rate. Working on the previously proposed steps is, therefore, needed to gain confidence in PBK modelling approaches for regulatory use.

In chemical risk assessment, default uncertainty factors are used to account for interspecies and interindividual differences, and differences in toxicokinetics and toxicodynamics herein. However, these default factors come with little scientific support. Therefore, our aim was to develop an in vitro method, using acetylcholinesterase (AChE) inhibition as a proof of principle, to assess both interspecies and interindividual differences in toxicodynamics. Electric eel enzyme and human blood of 20 different donors (12 men/8 women) were exposed to eight different compounds (chlorpyrifos, chlorpyrifos-oxon, phosmet, phosmet-oxon, diazinon, diazinon-oxon, pirimicarb, rivastigmine) and inhibition of AChE was measured using the Ellman method. The organophosphate parent compounds, chlorpyrifos, phosmet and diazinon, did not show inhibition of AChE. All other compounds showed concentration-dependent inhibition of AChE, with IC₅₀s in human blood ranging from 0.2–29 µM and IC₂₀s ranging from 0.1–18 µM, indicating that AChE is inhibited at concentrations relevant to the in vivo human situation. The oxon analogues were more potent inhibitors of electric eel AChE compared to human AChE. The opposite was true for carbamates, pointing towards interspecies differences for AChE inhibition. Human interindividual variability was low and ranged from 5–25%, depending on the concentration. This study provides a reliable in vitro method for assessing human variability in AChE toxicodynamics. The data suggest that the default uncertainty factor of ~ 3.16 may overestimate human variability for this toxicity endpoint, implying that specific toxicodynamic-related adjustment factors can support quantitative in vitro to in vivo extrapolations that link kinetic and dynamic data to improve chemical risk assessment.

Emerging advances in the field of in vitro toxicity testing attempt to meet the need for reliable human-based safety assessment in drug development. Intrahepatic cholangiocyte organoids (ICOs) are described as a donor-derived in vitro model for disease modelling and regenerative medicine. Here, we explored the potential of hepatocyte-like ICOs (HL-ICOs) in in vitro toxicity testing by exploring the expression and activity of genes involved in drug metabolism, a key determinant in drug-induced toxicity, and the exposure of HL-ICOs to well-known hepatotoxicants. The current state of drug metabolism in HL-ICOs showed levels comparable to those of PHHs and HepaRGs for CYP3A4; however, other enzymes, such as CYP2B6 and CYP2D6, were expressed at lower levels. Additionally, EC50 values were determined in HL-ICOs for acetaminophen (24.0–26.8 mM), diclofenac (475.5–>500 µM), perhexiline (9.7–>31.5 µM), troglitazone (23.1–90.8 µM), and valproic acid (>10 mM). Exposure to the hepatotoxicants showed EC50s in HL-ICOs comparable to those in PHHs and HepaRGs; however, for acetaminophen exposure, HL-ICOs were less sensitive. Further elucidation of enzyme and transporter activity in drug metabolism in HL-ICOs and exposure to a more extensive compound set are needed to accurately define the potential of HL-ICOs in in vitro toxicity testing.

Amphibian populations are undergoing a global decline worldwide. Such decline has been attributed to their unique physiology, ecology, and exposure to multiple stressors including chemicals, temperature, and biological hazards such as fungi of the Batrachochytrium genus, viruses such as Ranavirus, and habitat reduction. There are limited toxicity data for chemicals available for amphibians and few quantitative structure-activity relationship (QSAR) models have been developed and are publicly available. Such QSARs provide important tools to assess the toxicity of chemicals particularly in a data poor context. QSARs provide important tools to assess the toxicity of chemicals particularly when no toxicological data are available. This manuscript provides a description and validation of a regression-based QSAR model to predict, in a quantitative manner, acute lethal toxicity of aromatic chemicals in tadpoles of the Japanese brown frog (Rana japonica). QSAR models for acute median lethal molar concentrations (LC50-12 h) of waterborne chemicals using the Monte Carlo method were developed. The statistical characteristics of the QSARs were described as average values obtained from five random distributions into training and validation sets. Predictions from the model gave satisfactory results for the overall training set (R² = 0.72 and RMSE = 0.33) and were even more robust for the validation set (R² = 0.96 and RMSE = 0.11). Further development of QSAR models in amphibians, particularly for other life stages and species, are discussed.

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.