Toxins in Botanical Drugs and Plant-derived Food and Feed - from Science to Regulation: A Workshop Review.

In September 2022, the 3rd International Workshop on pyrrolizidine alkaloids (PAs) and related phytotoxins was held on-line, entitled 'Toxins in botanical drugs and plant-derived food and feed - from science to regulation'. The workshop focused on new findings about the occurrence, exposure, toxicity, and risk assessment of PAs. In addition, new scientific results related to the risk assessment of alkenylbenzenes, a distinct class of herbal constituents, were presented. The presence of PAs and alkenylbenzenes in plant-derived food, feed, and herbal medicines has raised health concerns with respect to their acute and chronic toxicity but mainly related to the genotoxic and carcinogenic properties of several congeners. The compounds are natural constituents of a variety of plant families and species widely used in medicinal, food, and feed products. Their individual occurrence, levels, and toxic properties, together with the broad range of congeners present in nature, represent a striking challenge to modern toxicology. This review tries to provide an overview of the current knowledge on these compounds and indicates needs and perspectives for future research.

A Computational Inter-Species Study on Safrole Phase I Metabolism-Dependent Bioactivation: A Mechanistic Insight into the Study of Possible Differences among Species.

Safrole, a 162.2 Da natural compound belonging to the alkenylbenzenes class, is classified as a possible carcinogen to humans by IARC (group IIB) and has proven to be genotoxic and carcinogenic to rodents. Despite its use as a food or feed additive, it is forbidden in many countries due to its documented toxicity; yet, it is still broadly present within food and feed and is particularly abundant in spices, herbs and essential oils. Specifically, safrole may exert its toxicity upon bioactivation to its proximate carcinogen 1'-hydroxy-safrole via specific members of the cytochrome P450 protein family with a certain inter/intra-species variability. To investigate this variability, an in-silico workflow based on molecular modelling, docking and molecular dynamics has been successfully applied. This work highlighted the mechanistic basis underpinning differences among humans, cats, chickens, goats, sheep, dogs, mice, pigs, rats and rabbits. The chosen metric to estimate the likeliness of formation of 1'-hydroxy-safrole by the species-specific cytochrome P450 under investigation allowed for the provision of a knowledge-based ground to rationally design and prioritise further experiments and deepen the current understanding of alkenylbenzenes bioactivation and CYPs mechanics. Both are crucial for a more informed framework of analysis for safrole toxicity.

Perfluoroalkyl substances (PFASs) are substrates of the renal human organic anion transporter 4 (OAT4).

Poly- and perfluoroalkyl substances (PFASs) are omnipresent in the environment and have been shown to accumulate in humans. Most PFASs are not biotransformed in animals and humans, so that elimination is largely dependent on non-metabolic clearance via bile and urine. Accumulation of certain PFASs in humans may relate to their reabsorption from the pre-urine by transporter proteins in the proximal tubules of the kidney, such as URAT1 and OAT4. The present study assessed the in vitro transport of 7 PFASs (PFHpA, PFOA, PFNA, PFDA, PFBS, PFHxS and PFOS) applying URAT1- or OAT4-transfected human embryonic kidney (HEK) cells. Virtually no transport of PFASs could be measured in URAT1-transfected HEK cells. All PFASs, except PFBS, showed clear uptake in OAT4-transfected HEK cells. In addition, these in vitro results were further supported by in silico docking and molecular dynamic simulation studies assessing transporter-ligand interactions. Information on OAT4-mediated transport may provide insight into the accumulation potential of PFASs in humans, but other kinetic aspects may play a role and should also be taken into account. Quantitative information on all relevant kinetic processes should be integrated in physiologically based kinetic (PBK) models, to predict congener-specific accumulation of PFASs in humans in a more accurate manner.

A congener-specific modelling approach for the transfer of polychlorinated dibenzo-p-dioxins and dibenzofurans and dioxin-like polychlorinated biphenyls from feed to eggs of laying hens.

Calibration of a kinetic model for the transfer of PCDD/Fs and dl-PCBs from feed to the hen's body and eggs was thus far restricted to the total TEQ concentration, i.e. the summed concentrations of PCDD/Fs and dl-PCBs expressed in terms of equivalents of 2,3,7,8-TCDD. However, this approach may lead to over- or underestimation of the transfer if the mixture contains congeners with kinetic characteristics which differ considerably from those used in such a model. This paper extends a previous transfer model of PCDD/Fs and dl-PCBs from feed to egg yolk fat and abdominal fat of high production laying hens, based on the total TEQ approach, to the level of individual congeners. Both modelling approaches are compared and the new approach is presented as a webtool application. This congener-specific approach enabled the calibration of 25 of the 29 relevant PCDD/F and dl-PCB congeners with respect to their individual transfer characteristics to body fat and egg yolk fat and their clearance from the body. Limitations of the available experimental data prevented the calibration of 1,2,3,4,6,7,8-HpCDD, OCDD, OCDF and PCB 123. The fraction transferred to egg yolk fat after long-term daily intake of contaminated feed was found to be at least 0.78 for 2,3,7,8-TCDD, 0.75 for PeCDD, 0.42-0.61 for HxCDDs, 0.70 for 2,3,7,8-TCDF, 0.71 for PeCDF, 0.54-0.60 for HxCDFs, 0.18-0.24 for HpCDFs and 0.89-1.00 for dl-PCBs. Various experimental and feed incident mixtures were used to compare the total TEQ- model with the congener-specific approach. An overestimation of the transfer by the total TEQ method was shown in particular for mixtures with a substantial contribution of hexa-, hepta- and octa-PCDD/Fs to the total TEQ level.

Impact of in vitro experimental variation in kinetic parameters on physiologically based kinetic (PBK) model simulations.

In vitro toxicokinetic data are critical in meeting an increased regulatory need to improve chemical safety evaluations towards a better understanding of internal human chemical exposure and toxicity. In vitro intrinsic hepatic clearance (CLint), the fraction unbound in plasma (fup), and the intestinal apparent permeability (Papp) are important parameters as input in a physiologically based kinetic (PBK) model to make first estimates of internal exposure after oral dosing. In the present study we explored the experimental variation in the values for these parameters as reported in the literature. Furthermore, the impact that this experimental variation has on PBK model predictions of maximum plasma concentration (Cmax) and the area under the concentration time curve (AUC0-24h) was determined. As a result of the experimental variation in CLint, Papp, and fup, the predicted variation in Cmax for individual compounds ranged between 1.4- to 28-fold, and the predicted variation in AUC0-24h ranged between 1.4- and 23-fold. These results indicate that there are still some important steps to take to achieve robust data that can be used in regulatory applications. To gain regulatory acceptance of in vitro kinetic data and PBK models based on in vitro input data, the boundaries in experimental conditions as well as the applicability domain and the use of different in vitro kinetic models need to be described in guidance documents.

Towards best use and regulatory acceptance of generic physiologically based kinetic (PBK) models for in vitro-to-in vivo extrapolation (IVIVE) in chemical risk assessment.

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.

Predictive performance of next generation human physiologically based kinetic (PBK) models based on in vitro and in silico input data.

The goal of the present study was to assess the predictive performance of a generic human physiologically based kinetic (PBK) model based on in vitro and in silico input data and the effect of using different input approaches for chemical parameterization on those predictions. For this purpose, a dataset was created of 38,772 Cmax predictions for 44 compounds by applying different combinations of in vitro and in silico approaches for chemical parameterization, and these predicted Cmax values were compared to reported in vivo data. Best results were achieved when the hepatic clearance was parameterized based on in vitro (i.e., hepatocytes or liver S9) measured intrinsic clearance values, the method of Rodgers and Rowland for calculating tissue:plasma partition coefficients, and the method of Lobell and Sivarajah for calculating the fraction unbound in plasma. With these parameters, the median Cmax values of 34 out of the 44 compounds were predicted within 5-fold of the observed Cmax, and the Cmax values of 19 compounds were predicted within 2-fold. The median Cmax values of 10 compounds were more than 5-fold overestimated. Underestimations (> 5-fold) did not occur. A comparison of the current generic PBK model structure with chemical-specific PBK models available in literature was made to identify possible kinetic processes not included in the generic PBK model that might explain the overestimations. Overall, the results provide crucial insights into the predictive performance of PBK models based on in vitro and in silico input and the influence of different input approaches on the model predictions.

Predictive Performance of Next Generation Physiologically Based Kinetic (PBK) Model Predictions in Rats Based on In Vitro and In Silico Input Data.

The goal of the present study was to assess the predictive performance of a minimal generic rat physiologically based kinetic (PBK) model based on in vitro and in silico input data to predict peak plasma concentrations (Cmax) upon single oral dosing. To this purpose, a dataset was generated of 3960 Cmax predictions for 44 compounds, applying different combinations of in vitro and in silico approaches for chemical parameterization, and comparison of the predictions to reported in vivo data. Best performance was obtained when (1) the hepatic clearance was parameterized based on in vitro measured intrinsic clearance values, (2) the method of Rodgers and Rowland for calculating partition coefficients, and (3) in silico calculated fraction unbound plasma and Papp values (the latter especially for very lipophilic compounds). Based on these input data, the median Cmax of 32 compounds could be predicted within 10-fold of the observed Cmax, with 22 out of these 32 compounds being predicted within 5-fold, and 8 compounds within 2-fold. Overestimations of more than 10-fold were observed for 12 compounds, whereas no underestimations of more than 10-fold occurred. Median Cmax predictions were frequently found to be within 10-fold of the observed Cmax when the scaled unbound hepatic intrinsic clearance (Clint,u) was either higher than 20 l/h or lower than 1 l/h. Similar findings were obtained with a test set of 5 in-house BASF compounds. Overall, this study provides relevant insights in the predictive performance of a minimal PBK model based on in vitro and in silico input 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.

Development of a Web-Based Toolbox to Support Quantitative In-Vitro-to-In-Vivo Extrapolations (QIVIVE) within Nonanimal Testing Strategies.

The goal of the present study was to develop an online web-based toolbox that contains generic physiologically based kinetic (PBK) models for rats and humans, including underlying calculation tools to predict plasma protein binding and tissue:plasma distribution, to be used for quantitative in-vitro-to-in-vivo extrapolations (QIVIVE). The PBK models within the toolbox allow first estimations of internal plasma and tissue concentrations of chemicals to be made, based on the logP and pKa of the chemicals and values for intestinal uptake and intrinsic hepatic clearance. As a case study, the toolbox was used to predict oral equivalent doses of in vitro ToxCast bioactivity data for the food additives methylparaben, propyl gallate, octyl gallate, and dodecyl gallate. These oral equivalent doses were subsequently compared with human exposure estimates, as a low tier assessment allowing prioritization for further assessment. The results revealed that daily intake levels of especially propyl gallate can lead to internal plasma concentrations that are close to in vitro biological effect concentrations, particularly with respect to the inhibition of human thyroid peroxidase (TPO). Estrogenic effects were not considered likely to be induced by the food additives, as daily exposure levels of the different compounds remained 2 orders of magnitude below the oral equivalent doses for in vitro estrogen receptor activation. Overall, the results of the study show how the toolbox, which is freely accessible through www.qivivetools.wur.nl, can be used to obtain initial internal dose estimates of chemicals and to prioritize chemicals for further assessment, based on the comparison of oral equivalent doses of in vitro biological activity data with human exposure levels.

Exploring the Potential of ToxCast Data in Supporting Read-Across for Evaluation of Food Chemical Safety.

The intention of this study was to determine the utility of high-throughput screening (HTS) data, as exemplified by ToxCast and Tox21, for application in toxicological read-across in food-relevant chemicals. Key questions were addressed on the extent to which the HTS data could provide information enabling (1) the elucidation of underlying bioactivities associated with apical toxicological outcomes, (2) the closing of existing toxicological data gaps, and (3) the definition of the boundaries of chemical space across which bioactivity could reliably be extrapolated. Results revealed that many biological targets apparently activated within the chemical groupings lack, at this time, validated toxicity pathway associations. Therefore, as means of providing proof-of-principle, a comparatively well-characterized end point-estrogenicity-was selected for evaluation. This was facilitated through the preparation of two exploratory case studies, focusing upon groupings of paraben-gallates and pyranone-type compounds (notably flavonoids). Within both, the HTS data were seen to reflect estrogenic potencies in a manner which broadly corresponded to established structure-activity group relationships, with parabens and flavonoids displaying greater estrogen receptor affinity than benzoate esters and alternative pyranone-containing molecules, respectively. As such, utility in the identification of out-of-domain compounds was demonstrated, indicating potential for application in addressing point (3) as detailed above.

Defining in vivo dose-response curves for kidney DNA adduct formation of aristolochic acid I in rat, mouse and human by an in vitro and physiologically based kinetic modeling approach.

Aristolochic acid I (AAI) is a well-known genotoxic kidney carcinogen. Metabolic conversion of AAI into the DNA-reactive aristolactam-nitrenium ion is involved in the mode of action of tumor formation. This study aims to predict in vivo AAI-DNA adduct formation in the kidney of rat, mouse and human by translating the in vitro concentration-response curves for AAI-DNA adduct formation to the in vivo situation using physiologically based kinetic (PBK) modeling-based reverse dosimetry. DNA adduct formation in kidney proximal tubular LLC-PK1 cells exposed to AAI was quantified by liquid chromatography-electrospray ionization-tandem mass spectrometry. Subsequently, the in vitro concentration-response curves were converted to predicted in vivo dose-response curves in rat, mouse and human kidney using PBK models. Results obtained revealed a dose-dependent increase in AAI-DNA adduct formation in the rat, mouse and human kidney and the predicted DNA adduct levels were generally within an order of magnitude compared with values reported in the literature. It is concluded that the combined in vitro PBK modeling approach provides a novel way to define in vivo dose-response curves for kidney DNA adduct formation in rat, mouse and human and contributes to the reduction, refinement and replacement of animal testing.

New approach methodologies (NAMs) for human-relevant biokinetics predictions. Meeting the paradigm shift in toxicology towards an animal-free chemical risk assessment

For almost fifteen years, the availability and regulatory acceptance of new approach methodologies (NAMs) to assess the absorption, distribution, metabolism and excretion (ADME/biokinetics) in chemical risk evaluations are a bottleneck. To enhance the field, a team of 24 experts from science, industry, and regulatory bodies, including new generation toxicologists, met at the Lorentz Centre in Leiden, The Netherlands. A range of possibilities for the use of NAMs for biokinetics in risk evaluations were formulated (for example to define species differences and human variation or to perform quantitative in vitro-in vivo extrapolations). To increase the regulatory use and acceptance of NAMs for biokinetics for these ADME considerations within risk evaluations, the development of test guidelines (protocols) and of overarching guidance documents is considered a critical step. To this end, a need for an expert group on biokinetics within the Organisation of Economic Cooperation and Development (OECD) to supervise this process was formulated. The workshop discussions revealed that method development is still required, particularly to adequately capture transporter mediated processes as well as to obtain cell models that reflect the physiology and kinetic characteristics of relevant organs. Developments in the fields of stem cells, organoids and organ-on-a-chip models provide promising tools to meet these research needs in the future.

Potential of ToxCast Data in the Safety Assessment of Food Chemicals.

Tox21 and ToxCast are high-throughput in vitro screening programs coordinated by the U.S. National Toxicology Program and the U.S. Environmental Protection Agency, respectively, with the goal of forecasting biological effects in vivo based on bioactivity profiling. The present study investigated whether mechanistic insights in the biological targets of food-relevant chemicals can be obtained from ToxCast results when the chemicals are grouped according to structural similarity. Starting from the 556 direct additives that have been identified in the ToxCast database by Karmaus et al. [Karmaus, A. L., Trautman, T. D., Krishan, M., Filer, D. L., and Fix, L. A. (2017). Curation of food-relevant chemicals in ToxCast. Food Chem. Toxicol. 103, 174-182.], the results showed that, despite the limited number of assays in which the chemical groups have been tested, sufficient results are available within so-called "DNA binding" and "nuclear receptor" target families to profile the biological activities of the defined chemical groups for these targets. The most obvious activity identified was the estrogen receptor-mediated actions of the chemical group containing parabens and structurally related gallates, as well the chemical group containing genistein and daidzein (the latter 2 being particularly active toward estrogen receptor ß as a potential health benefit). These group effects, as well as the biological activities of other chemical groups, were evaluated in a series of case studies. Overall, the results of the present study suggest that high-throughput screening data could add to the evidence considered for regulatory risk assessment of food chemicals and to the evaluation of desirable effects of nutrients and phytonutrients. The data will be particularly useful for providing mechanistic information and to fill data gaps with read-across.

Towards harmonization of test methods for in vitro hepatic clearance studies.

Non-animal methods for toxicokinetics, such as in vitro hepatic metabolic clearance studies, play an important role in chemical risk evaluations. To gain regulatory acceptance of such clearance data, the development of a test guideline for performing in vitro hepatic clearance studies is crucial. The aim of the present study was to obtain insight in the experimental conditions of clearance studies that influence obtained intrinsic clearance (CLint) values. To that end, in vitro hepatic CLint data obtained with rat or human hepatocytes and methodological aspects of the experiments, were collected from 42 different suitable studies published between 1995 and 2018. The CLint values for the majority of chemicals differed by more than one order of magnitude. We estimated the systematic effect of different experimental setups on the CLint values using a random forest regression analysis, revealing that 'hepatocyte concentration', 'species' (rat or human hepatocytes) and 'culture medium' have the largest impact. Calculating unbound CLint (CLint,u) values slightly reduced the variation for most chemicals. Given that in vivo clearance is in general underpredicted based on in vitro clearance data, a harmonized protocol is preferably based on a protocol that provides relatively high in vitro CLint values.

Quantitative in vitro-to-in vivo extrapolation (QIVIVE) of estrogenic and anti-androgenic potencies of BPA and BADGE analogues.

The goal of the present study was to obtain an in vivo relevant prioritization method for the endocrine potencies of different polycarbonate monomers, by combining in vitro bioassay data with physiologically based kinetic (PBK) modelling. PBK models were developed for a selection of monomers, including bisphenol A (BPA), two bisphenol F (BPF) isomers and four different bisphenol A diglycidyl ethers (BADGEs), using in vitro input data. With these models, the plasma concentrations of the compounds were simulated, providing means to estimate the dose levels at which the in vitro endocrine effect concentrations are reached. The results revealed that, whereas the in vitro relative potencies of different BADGEs (predominantly anti-androgenic effects) can be up to fourfold higher than BPA, the estimated in vivo potencies based on the oral equivalent doses are one to two orders of magnitude lower than BPA because of fast detoxification of the BADGEs. In contrast, the relative potencies of 2,2-BPF and 4,4-BPF increase when accounting for the in vivo availability. 4,4-BPF is estimated to be fivefold more potent than BPA in humans in vivo in inducing estrogenic effects and both 2,2-BPF and 4,4-BPF are estimated to be, respectively, 7 and 11-fold more potent in inducing anti-androgenic effects. These relative potencies were considered to be first-tier estimates, particularly given that the potential influence of intestinal metabolism on the in vivo availability was not accounted for. Overall, it can be concluded that both 2,2-BPF and 4,4-BPF are priority compounds.

Impact of nanoparticle surface functionalization on the protein corona and cellular adhesion, uptake and transport.

BACKGROUND: Upon ingestion, nanoparticles can interact with the intestinal epithelial barrier potentially resulting in systemic uptake of nanoparticles. Nanoparticle properties have been described to influence the protein corona formation and subsequent cellular adhesion, uptake and transport. Here, we aimed to study the effects of nanoparticle size and surface chemistry on the protein corona formation and subsequent cellular adhesion, uptake and transport. Caco-2 intestinal cells, were exposed to negatively charged polystyrene nanoparticles (PSNPs) (50 and 200 nm), functionalized with sulfone or carboxyl groups, at nine nominal concentrations (15-250 µg/ml) for 10 up to 120 min. The protein coronas were analysed by LC-MS/MS. RESULTS: Subtle differences in the protein composition of the two PSNPs with different surface chemistry were noted. High-content imaging analysis demonstrated that sulfone PSNPs were associated with the cells to a significantly higher extent than the other PSNPs. The apparent cellular adhesion and uptake of 200 nm PSNPs was not significantly increased compared to 50 nm PSNPs with the same surface charge and chemistry. Surface chemistry outweighs the impact of size on the observed PSNP cellular associations. Also transport of the sulfone PSNPs through the monolayer of cells was significantly higher than that of carboxyl PSNPs. CONCLUSIONS: The results suggest that the composition of the protein corona and the PSNP surface chemistry influences cellular adhesion, uptake and monolayer transport, which might be predictive of the intestinal transport potency of NPs.

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.

Data Analyses and Modelling for Risk Based Monitoring of Mycotoxins in Animal Feed.

Following legislation, European Member States should have multi-annual control programs for contaminants, such as for mycotoxins, in feed and food. These programs need to be risk based implying the checks are regular and proportional to the estimated risk for animal and human health. This study aimed to prioritize feed products in the Netherlands for deoxynivalenol and aflatoxin B1 monitoring. Historical mycotoxin monitoring results from the period 2007-2016 were combined with data from other sources. Based on occurrence, groundnuts had high priority for aflatoxin B1 monitoring; some feed materials (maize and maize products and several oil seed products) and complete/complementary feed excluding dairy cattle and young animals had medium priority; and all other animal feeds and feed materials had low priority. For deoxynivalenol, maize by-products had a high priority, complete and complementary feed for pigs had a medium priority and all other feed and feed materials a low priority. Also including health consequence estimations showed that feed materials that ranked highest for aflatoxin B1 included sunflower seed and palmkernel expeller/extracts and maize. For deoxynivalenol, maize products were ranked highest, followed by various small grain cereals (products); all other feed materials were of lower concern. Results of this study have proven to be useful in setting up the annual risk based control program for mycotoxins in animal feed and feed materials.

Expert opinions on the acceptance of alternative methods in food safety evaluations: Formulating recommendations to increase acceptance of non-animal methods for kinetics.

Inclusion of alternative methods that replace, reduce, or refine (3R) animal testing within regulatory safety evaluations of chemicals generally faces many hurdles. The goal of the current work is to i) collect responses from key stakeholders involved in food safety evaluations on what they consider the most relevant factors that influence the acceptance and use of 3R methods and to ii) use these responses to formulate activities needed to increase the acceptance and use of 3R methods, particularly for kinetics. The stakeholders were contacted by e-mail for their opinions, asking the respondents to write down three barriers and/or drivers and scoring these by distributing 5 points over the three factors. The main barriers that obtained the highest aggregated scores were i) uncertain predictability 3R methods/lack of validation, ii) insufficient guidance regulators/industry and iii) insufficient harmonization of legislation. The major driver identified was the possibility of 3R methods to provide more mechanistic information. Based on the results, recommendations are given to enhance the acceptance and application of 3R toxicokinetic methods in food safety evaluations. These include steering of regulatory data requirements as well as creating (funding) opportunities for development and validation of alternative methods for kinetics and development of guidances.