TKPlate 1.0: An Open-access platform for toxicokinetic and toxicodynamic modelling of chemicals to implement new approach methodologies in chemical risk assessment.

This publication is linked to the following EFSA Supporting Publications articles: http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8441/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8440/full, http://onlinelibrary.wiley.com/doi/10.2903/sp.efsa.2023.EN-8437/full.

OpenCYP: An open source database exploring human variability in activities and frequencies of polymophisms for major cytochrome P-450 isoforms across world populations.

The open source database "OpenCYP database" has been developed based on the results of extensive literature searches from the peer-reviewed literature. OpenCYP provides data on human variability on baseline of activities and polymophism frequencies for selected cytochrome P-450 isoforms (CYP1A2, CYP2A6, CYP2D6, CYP3A4/3A5 and CYP3A7) in healthy adult populations from world populations. CYP enzymatic activities were generally expressed as the metabolic ratio (MR) between an unchanged probe drug and its metabolite(s) in urine or plasma measured in healthy adults. Data on other age groups were very limited and fragmented, constituting an important data gap. Quantitative comparisons were often hampered by the different experimental conditions used. However, variability was quite limited for CYP1A2, using caffeine as a probe substrate, with a symmetrical distribution of metabolic activity values. For CYP3A4, human variability was dependent on the probe substrate itself with low variability when data considering the dextromethorphan/demethilathed metabolite MR were used and large variability when the urinary 6ß-hydroxycortisol/cortisol ratio was used. The largest variability in CYP activity was shown for CYP2D6 activity, after oral dosing of dextromethorphan, for which genetic polymorphisms are well characterised and constitute a significant source of variability. It is foreseen that the OpenCYP database can contribute to promising tools to support the further development of QIVIVE and PBK models for human risk assessment of chemicals particularly when combined with information on isoform-specific content in cells using proteomic approaches.

Application of in vitro data in physiologically-based kinetic models for quantitative in vitro-in vivo extrapolation: A case-study for baclofen.

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.

In vitro metabolism of lidocaine in subcellular post-mitochondrial fractions and precision cut slices from cattle liver.

In vitro models are widely used to study the biotransformation of xenobiotics and to provide input parameters to physiologically based kinetic models required to predict the kinetic behavior in vivo. For farm animals this is not common practice yet. The use of slaughterhouse-derived tissue material may provide opportunities to study biotransformation reactions in farm animals. The goal of the present study was to explore the potential of slaughterhouse-derived bovine liver S9 (S9) and precision cut liver slices (PCLSs) to capture observed biotransformation reactions of lidocaine in cows. The in vitro data obtained with both S9 and PCLSs confirm in vivo findings that 2,6-dimethylaniline (DMA) is an important metabolite of lidocaine in cows, being for both PCLSs and S9 the end-product. In case of S9, also conversion of lidocaine to lidocaine-N-oxide and monoethylglycinexylidine (MEXG) was observed. MEGX is considered as intermediate for DMA formation, given that this metabolite was metabolized to DMA by both PLCSs and S9. In contrast to in vivo, no in vitro conversion of DMA to 4-OH-DMA was observed. Further work is needed to explain this lack of conversion and to further evaluate the use of slaughterhouse-derived tissue materials to predict the biotransformation of xenobiotics in farm animals.

Towards a systematic method for assessing the impact of chemical pollution on ecosystem services of water systems.

Chemical pollution impinges on the quality of water systems and the ecosystem services (ESs) they provide. Expression of ESs in monetary units has become an essential tool for sustainable ecosystem management. However, the impact of chemical pollution on ESs is rarely quantified, and ES valuation often focuses on individual services without considering the total services provided by the ecosystem. The purpose of the study was to develop a stepwise approach to quantify the impact of sediment pollution on the total ES value provided by water systems. Thereby, we calculated the total ES value loss as a function of the multi-substance potentially affected fraction of species at the HC50 level (msPAF(HC50)). The function is a combination of relationships between, subsequently: the msPAF(HC50), diversity, productivity and total ES value. Regardless of the inherent differences between terrestrial and aquatic ecosystems, an increase of diversity generally corresponded to an increase in productivity with curvilinear or linear effects. A positive correlation between productivity and total values of ESs of biomes was observed. The combined relationships showed that 1% msPAF(HC50) corresponded to on average 0.5% (0.05-1.40%) of total ES value loss. The ES loss due to polluted sediments in the Waal-Meuse river estuary (the Netherlands) and Flemish waterways (Belgium) was estimated to be 0.3-5 and 0.6-10 thousand 2007$/ha/yr, respectively. Our study presents a novel methodology to assess the impact of chemical exposure on diversity, productivity, and total value that ecosystems provide. With sufficient monitoring data, our generic methodology can be applied for any chemical and region of interest and help water managers make informed decisions on cost-effective measures to remedy pollution. Acknowledging that the ES loss estimates as a function of PAF(HC50) are crude, we explicitly discuss the uncertainties in each step for further development and application of the methodology.

Inter-phenotypic differences in CYP2C9 and CYP2C19 metabolism: Bayesian meta-regression of human population variability in kinetics and application in chemical risk assessment.

Quantifying variability in pharmacokinetics (PK) and toxicokinetics (TK) provides a science-based approach to refine uncertainty factors (UFs) for chemical risk assessment. In this context, genetic polymorphisms in cytochromes P450 (CYPs) drive inter-phenotypic differences and may result in reduction or increase in metabolism of drugs or other xenobiotics. Here, an extensive literature search was performed to identify PK data for probe substrates of the human polymorphic isoforms CYP2C9 and CYP2C19. Relevant data from 158 publications were extracted for markers of chronic exposure (clearance and area under the plasma concentration-time curve) and analysed using a Bayesian meta-regression model. Enzyme function (EF), driven by inter-phenotypic differences across a range of allozymes present in extensive and poor metabolisers (EMs and PMs), and fraction metabolised (Fm), were identified as exhibiting the highest impact on the metabolism. The Bayesian meta-regression model provided good predictions for such inter-phenotypic differences. Integration of population distributions for inter-phenotypic differences and estimates for EF and Fm allowed the derivation of CYP2C9- and CYP2C19-related UFs which ranged from 2.7 to 12.7, and were above the default factor for human variability in TK (3.16) for PMs and major substrates (Fm >60%). These results provide population distributions and pathway-related UFs as conservative in silico options to integrate variability in CYP2C9 and CYP2C19 metabolism using in vitro kinetic evidence and in the absence of human data. The future development of quantitative extrapolation models is discussed with particular attention to integrating human in vitro and in vivo PK or TK data with pathway-related variability for chemical risk assessment.

Human variability in glutathione-S-transferase activities, tissue distribution and major polymorphic variants: Meta-analysis and implication for chemical risk assessment.

The input into the QIVIVE and Physiologically-Based kinetic and dynamic models of drug metabolising enzymes performance and their inter-individual differences significantly improve the modelling performance, supporting the development and integration of alternative approaches to animal testing. Bayesian meta-analyses allow generating and integrating statistical distributions with human in vitro metabolism data for quantitative in vitro-in vivo extrapolation. Such data are lacking on glutathione-S-transferases (GSTs). This paper reports for the first time results on the human variability of GST activities in healthy individuals, their tissue localisation and the frequencies of their major polymorphic variants by means of extensive literature search, data collection, data base creation and meta-analysis. A limited number of papers focussed on in vivo GST inter-individual differences in humans. Ex-vivo total GST activity without discriminating amongst isozymes is generally reported, resulting in a high inter-individual variability. The highest levels of cytosolic GSTs in humans are measured in the kidney, liver, adrenal glands and blood. The frequencies of GST polymorphisms for cytosolic isozymes in populations of different geographical ancestry were also presented. Bayesian meta-analyses to derive GST-related uncertainty factors provided uncertain estimates, due to the limited database. Considering the relevance of GST activities and their pivotal role in cellular adaptive response mechanisms to chemical stressors, further studies are needed to identify GST probe substrates for specific isozymes and quantify inter-individual differences.

Acetylcholinesterase inhibition in electric eel and human donor blood: an in vitro approach to investigate interspecies differences and human variability in toxicodynamics.

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 IC50s in human blood ranging from 0.2-29 µM and IC20s 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.