Toxicokinetics of the Antidepressant Fluoxetine and Its Active Metabolite Norfluoxetine in Caenorhabditis elegans and Their Comparative Potency.

The nematode Caenorhabditis elegans is a valuable model for ecotoxicological research, yet limited attention has been given to understanding how it absorbs, distributes, metabolizes, and excretes chemicals. This is crucial for C. elegans because the organism is known to have strong uptake barriers that are known to be susceptible to potential confounding effects of the presence of Escherichia coli as a food source. One frequently studied compound in C. elegans is the antidepressant fluoxetine, which has an active metabolite norfluoxetine. In this study, we evaluated the toxicokinetics and relative potency of norfluoxetine and fluoxetine in chemotaxis and activity tests. Toxicokinetics experiments were conducted with varying times, concentrations of fluoxetine, and in the absence or presence of E. coli, simulated with a one-compartment model. Our findings demonstrate that C. elegans can take up fluoxetine and convert it into norfluoxetine. Norfluoxetine proved slightly more potent and had a longer elimination half-life. The bioconcentration factor, uptake, and elimination rate constants depended on exposure levels, duration, and the presence of E. coli in the exposure medium. These findings expand our understanding of toxicokinetic modeling in C. elegans for different exposure scenarios, underlining the importance of considering norfluoxetine formation in exposure and bioactivity assessments of fluoxetine.

Integrating In Vitro Data and Physiologically Based Kinetic Modeling to Predict and Compare Acute Neurotoxic Doses of Saxitoxin in Rats, Mice, and Humans.

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.

Acetylcholinesterase Inhibition in Rats and Humans Following Acute Fenitrothion Exposure Predicted by Physiologically Based Kinetic Modeling-Facilitated Quantitative In Vitro to In Vivo Extrapolation.

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.

Inter-individual variation in chlorpyrifos toxicokinetics characterized by physiologically based kinetic (PBK) and Monte Carlo simulation comparing human liver microsome and Supersome™ cytochromes P450 (CYP)-specific kinetic data as model input.

The present study compares two approaches to evaluate the effects of inter-individual differences in the biotransformation of chlorpyrifos (CPF) on the sensitivity towards in vivo red blood cell (RBC) acetylcholinesterase (AChE) inhibition and to calculate a chemical-specific adjustment factor (CSAF) to account for inter-individual differences in kinetics (HKAF). These approaches included use of a Supersome™ cytochromes P450 (CYP)-based and a human liver microsome (HLM)-based physiologically based kinetic (PBK) model, both combined with Monte Carlo simulations. The results revealed that bioactivation of CPF exhibits biphasic kinetics caused by distinct differences in the Km of CYPs involved, which was elucidated by Supersome™ CYP rather than by HLM. Use of Supersome™ CYP-derived kinetic data was influenced by the accuracy of the intersystem extrapolation factors (ISEFs) required to scale CYP isoform activity of Supersome™ to HLMs. The predicted dose-response curves for average, 99th percentile and 1st percentile sensitive individuals were found to be similar in the two approaches when biphasic kinetics was included in the HLM-based approach, resulting in similar benchmark dose lower confidence limits for 10% inhibition (BMDL10) and HKAF values. The variation in metabolism-related kinetic parameters resulted in HKAF values at the 99th percentile that were slightly higher than the default uncertainty factor of 3.16. While HKAF values up to 6.9 were obtained when including also the variability in other influential PBK model parameters. It is concluded that the Supersome™ CYP-based approach appeared most adequate for identifying inter-individual variation in biotransformation of CPF and its resulting RBC AChE inhibition.

Physiologically based kinetic modelling predicts the in vivo relative potency of riddelliine N-oxide compared to riddelliine in rat to be dose dependent.

Pyrrolizidine alkaloids (PAs) are toxic plant constituents occurring often in their N-oxide form. This raises the question on the relative potency (REP) values of PA-N-oxides compared to the corresponding parent PAs. The present study aims to quantify the in vivo REP value of riddelliine N-oxide compared to riddelliine using physiologically based kinetic (PBK) modelling, taking into account that the toxicity of riddelliine N-oxide depends on its conversion to riddelliine by intestinal microbiota and in the liver. The models predicted a lower Cmax and higher Tmax for the blood concentration of riddelliine upon oral administration of riddelliine N-oxide compared to the Cmax and Tmax predicted for an equimolar oral dose of riddelliine. Comparison of the area under the riddelliine concentration-time curve (AUCRID) obtained upon dosing either the N-oxide or riddelliine itself revealed a ratio of 0.67, which reflects the in vivo REP for riddelliine N-oxide compared to riddelliine, and appeared to closely match the REP value derived from available in vivo data. The models also predicted that the REP value will decrease with increasing dose level, because of saturation of riddelliine N-oxide reduction by the intestinal microbiota and of riddelliine clearance by the liver. It is concluded that PBK modeling provides a way to define in vivo REP values of PA-N-oxides as compared to their parent PAs, without a need for animal experiments.

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.

Molecular Dynamics and In Vitro Quantification of Safrole DNA Adducts Reveal DNA Adduct Persistence Due to Limited DNA Distortion Resulting in Inefficient Repair.

The formation and repair of N2-(trans-isosafrol-3'-yl)-2'-deoxyguanosine (S-3'-N2-dG) DNA adduct derived from the spice and herbal alkenylbenzene constituent safrole were investigated. DNA adduct formation and repair were studied in vitro and using molecular dynamics (MD) simulations. DNA adduct formation was quantified using liquid chromatography-mass spectrometry (LCMS) in wild type and NER (nucleotide excision repair) deficient CHO cells and also in HepaRG cells and primary rat hepatocytes after different periods of repair following exposure to safrole or 1'-hydroxysafrole (1'-OH safrole). The slower repair of the DNA adducts found in NER deficient cells compared to that in CHO wild type cells indicates a role for NER in repair of S-3'-N2-dG DNA adducts. However, DNA repair in liver cell models appeared to be limited, with over 90% of the adducts remaining even after 24 or 48 h recovery. In our further studies, MD simulations indicated that S-3'-N2-dG adduct formation causes only subtle changes in the DNA structure, potentially explaining inefficient activation of NER. Inefficiency of NER mediated repair of S-3'-N2-dG adducts points at persistence and potential bioaccumulation of safrole DNA adducts upon daily dietary exposure.

Monocrotaline-induced liver toxicity in rat predicted by a combined in vitro physiologically based kinetic modeling approach.

The aim of the present study was to use an in vitro-in silico approach to predict the in vivo acute liver toxicity of monocrotaline and to characterize the influence of its metabolism on its relative toxic potency compared to lasiocarpine and riddelliine. In the absence of data on acute liver toxicity of monocrotaline upon oral exposure, the predicted dose-response curve for acute liver toxicity in rats and the resulting benchmark dose lower and upper confidence limits for 10% effect (BMDL10 and BMDU10) were compared to data obtained in studies with intraperitoneal or subcutaneous dosing regimens. This indicated the predicted BMDL10 value to be in line with the no-observed-adverse-effect levels (NOAELs) derived from availabe in vivo studies. The predicted BMDL10-BMDU10 of 1.1-4.9 mg/kg bw/day also matched the oral dose range of 1-3 mg PA/kg bw/day at which adverse effects in human are reported. A comparison to the oral toxicity of the related pyrrolizidine alkaloids (PAs) lasiocarpine and riddelliine revealed that, although in the rat hepatocytes monocrotaline was less toxic than lasiocarpine and riddelliine, due to its relatively inefficient clearance, its in vivo acute liver toxicity was predicted to be comparable. It is concluded that the combined in vitro-PBK modeling approach can provide insight in monocrotaline-induced acute liver toxicity in rats, thereby filling existing gaps in the database on PA toxicity. Furthermore, the results reveal that the kinetic and metabolic properties of PAs can vary substantially and should be taken into account when considering differences in relative potency between different PAs.

Cellular levels and molecular dynamics simulations of estragole DNA adducts point at inefficient repair resulting from limited distortion of the double-stranded DNA helix.

Estragole, naturally occurring in a variety of herbs and spices, can form DNA adducts after bioactivation. Estragole DNA adduct formation and repair was studied in in vitro liver cell models, and a molecular dynamics simulation was used to investigate the conformation dependent (in)efficiency of N2-(trans-isoestragol-3'-yl)-2'-deoxyguanosine (E-3'-N2-dG) DNA adduct repair. HepG2, HepaRG cells, primary rat hepatocytes and CHO cells (including CHO wild-type and three NER-deficient mutants) were exposed to 50 µM estragole or 1'-hydroxyestragole and DNA adduct formation was quantified by LC-MS immediately following exposure and after a period of repair. Results obtained from CHO cell lines indicated that NER plays a role in repair of E-3'-N2-dG adducts, however, with limited efficiency since in the CHO wt cells 80% DNA adducts remained upon 24 h repair. Inefficiency of DNA repair was also found in HepaRG cells and primary rat hepatocytes. Changes in DNA structure resulting from E-3'-N2-dG adduct formation were investigated by molecular dynamics simulations. Results from molecular dynamics simulations revealed that conformational changes in double-stranded DNA by E-3'-N2-dG adduct formation are small, providing a possible explanation for the restrained repair, which may require larger distortions in the DNA structure. NER-mediated enzymatic repair of E-3'-N2-dG DNA adducts upon exposure to estragole will be limited, providing opportunities for accumulation of damage upon repeated daily exposure. The inability of this enzymatic repair is likely due to a limited distortion of the DNA double-stranded helix resulting in inefficient activation of nucleotide excision repair.

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.

Use of an in vitro-in silico testing strategy to predict inter-species and inter-ethnic human differences in liver toxicity of the pyrrolizidine alkaloids lasiocarpine and riddelliine.

Lasiocarpine and riddelliine are pyrrolizidine alkaloids (PAs) known to cause liver toxicity. The aim of this study was to predict the inter-species and inter-ethnic human differences in acute liver toxicity of lasiocarpine and riddelliine using physiologically based kinetic (PBK) modelling based reverse dosimetry of in vitro toxicity data. The concentration-response curves of in vitro cytotoxicity of lasiocarpine and riddelliine defined in pooled human hepatocytes were translated to in vivo dose-response curves by PBK models developed using kinetic data obtained from incubations with pooled tissue fractions from Chinese and Caucasian individuals, providing PBK models for the average Chinese and average Caucasian, respectively. From the predicted in vivo dose-response curves, the benchmark dose lower and upper confidence limits for 5% effect (BMDL5 and BMDU5) were derived and subsequently compared to those previously obtained in rat to evaluate inter-species differences. The inter-species differences amounted to 2.0-fold for lasiocarpine and 8.2-fold for riddelliine with humans being more sensitive than rats. The inter-ethnic human differences varied 2.0-fold for lasiocarpine and 5.0-fold for riddelliine with the average Caucasian being more sensitive than the average Chinese. In conclusion, the present study provides the proof-of-principle to predict inter-species and inter-ethnic differences in in vivo liver toxicity for PAs by an alternative testing strategy integrating in vitro cytotoxicity data with PBK modelling-based reverse dosimetry.

Study on inter-ethnic human differences in bioactivation and detoxification of estragole using physiologically based kinetic modeling.

Considering the rapid developments in food safety in the past decade in China, it is of importance to obtain insight into what extent safety and risk assessments of chemicals performed for the Caucasian population apply to the Chinese population. The aim of the present study was to determine physiologically based kinetic (PBK) modeling-based predictions for differences between Chinese and Caucasians in terms of metabolic bioactivation and detoxification of the food-borne genotoxic carcinogen estragole. The PBK models were defined based on kinetic constants for hepatic metabolism derived from in vitro incubations using liver fractions of the two ethnic groups, and used to evaluate the inter-ethnic differences in metabolic activation and detoxification of estragole. The models predicted that at realistic dietary intake levels, only 0.02% of the dose was converted to the ultimate carcinogenic metabolite 1'-sulfooxyestragole in Chinese subjects, whereas this amounted to 0.09% of the dose in Caucasian subjects. Detoxification of 1'-hydroxyestragole, mainly via conversion to 1'-oxoestragole, was similar within the two ethnic groups. The 4.5-fold variation in formation of the ultimate carcinogenic metabolite of estragole accompanied by similar rates of detoxification may indicate a lower risk of estragole for the Chinese population at similar levels of exposure. The study provides a proof of principle for how PBK modeling can identify differences in ethnic sensitivity and provide a more refined risk assessment for a specific ethnic group for a compound of concern.

Physiologically based kinetic modeling of the bioactivation of myristicin.

The present study describes physiologically based kinetic (PBK) models for the alkenylbenzene myristicin that were developed by extension of the PBK models for the structurally related alkenylbenzene safrole in rat and human. The newly developed myristicin models revealed that the formation of the proximate carcinogenic metabolite 1'-hydroxymyristicin in liver is at most 1.8 fold higher in rat than in human and limited for the ultimate carcinogenic metabolite 1'-sulfoxymyristicin to (2.8-4.0)-fold higher in human. In addition, a comparison was made between the relative importance of bioactivation for myristicin and safrole. Model predictions indicate that for these related compounds, the formation of the 1'-sulfoxy metabolites in rat and human liver is comparable with a difference of <2.2-fold over a wide dose range. The results from this PBK analysis support that risk assessment of myristicin may be based on the BMDL10 derived for safrole of 1.9-5.1 mg/kg bw per day. Using an estimated daily intake of myristicin of 0.0019 mg/kg bw per day resulting from the use of herbs and spices, this results in MOE values for myristicin that amount to 1000-2700, indicating a priority for risk management. The results obtained illustrate that PBK modeling provides insight into possible species differences in the metabolic activation of myristicin. Moreover, they provide an example of how PBK modeling can facilitate a read-across in risk assessment from a compound for which in vivo toxicity studies are available to a related compound for which tumor data are not reported, thus contributing to alternatives in animal testing.

Determination and risk assessment of naturally occurring genotoxic and carcinogenic alkenylbenzenes in nutmeg-based plant food supplements.

A risk assessment of nutmeg-based plant food supplements (PFS) containing different alkenylbenzenes was performed based on the alkenylbenzene levels quantified in a series of PFS collected via the online market. The estimated daily intake (EDI) of the alkenylbenzenes amounted to 0.3 to 312 µg kg-1 body weight (bw) for individual alkenylbenzenes, to 1.5 to 631 µg kg-1 bw when adding up the alkenylbenzene levels assuming equal potency, and to 0.4 to 295 µg kg-1 bw when expressed in safrole equivalents using toxic equivalency factors (TEFs). The margin of exposure approach (MOE) was used to evaluate the potential risks. Independent of the method used for the intake estimate, the MOE values obtained were generally lower than 10000 indicating a priority for risk management. When taking into account that PFS may be used for shorter periods of time and using Haber's rule to correct for shorter than lifetime exposure it was shown that limiting exposure to only 1 or 2 weeks would result in MOE values that would be, with the presently determined levels of alkenylbenzenes and proposed uses of the PFS, of low priority for risk management (MOE > 10000). It is concluded that the results of the present paper reveal that nutmeg-based PFS consumption following recommendations for daily intake especially for longer periods of time raise a concern. Copyright © 2017 John Wiley & Sons, Ltd.

DL-propargylglycine reduces blood pressure and renal injury but increases kidney weight in angiotensin-II infused rats.

Hydrogen sulfide (H2S), carbon monoxide (CO) and nitric oxide (NO) share signaling and vasorelaxant properties and are involved in proliferation and apoptosis. Inhibiting NO production or availability induces hypertension and proteinuria, which is prevented by concomitant blockade of the H2S producing enzyme cystathionine γ-lyase (CSE) by d,l-propargylglycine (PAG). We hypothesized that blocking H2S production ameliorates Angiotensin II (AngII)-induced hypertension and renal injury in a rodent model. Effects of concomitant administration of PAG or saline were therefore studied in healthy (CON) and AngII hypertensive rats. In CON rats, PAG did not affect systolic blood pressure (SBP), but slightly increased proteinuria. In AngII rats PAG reduced SBP, proteinuria and plasma creatinine (180 ± 12 vs. 211 ± 19 mmHg; 66 ± 35 vs. 346 ± 92 mg/24 h; 24 ± 6 vs. 47 ± 15 µmol/L, respectively; p < 0.01). Unexpectedly, kidney to body weight ratio was increased in all groups by PAG (p < 0.05). Renal injury induced by AngII was reduced by PAG (p < 0.001). HO-1 gene expression was increased by PAG alone (p < 0.05). PAG increased inner cortical tubular cell proliferation after 1 week and decreased outer cortical tubular nucleus number/field after 4 weeks. In vitro proximal tubular cell size increased after exposure to PAG. In summary, blocking H2S production with PAG reduced SBP and renal injury in AngII infused rats. Independent of the cardiovascular and renal effects, PAG increased HO-1 gene expression and kidney weight. PAG alone increased tubular cell size and proliferation in-vivo and in-vitro. Our results are indicative of a complex interplay of gasotransmitter signaling/action of mutually compensatory nature in the kidney.

In vitro-in silico study on the influence of dose, fraction bioactivated and endpoint used on the relative potency value of pyrrolizidine alkaloid N-oxides compared to parent pyrrolizidine alkaloids.

Pyrrolizidine alkaloids (PAs) and their N-oxides (PA-N-oxides) are phytotoxins found in food, feed and the environment. Yet, limited data exist from which the relative potency of a PA-N-oxide relative to its corresponding PA (REPPANO to PA) can be defined. This study aims to investigate the influence of dose, fraction bioactivated and endpoint on the REPPANO to PA of a series of pyrrolizidine N-oxides using in vitro-in silico data and physiologically based kinetic (PBK) modeling. The first endpoint used to calculate the REPPANO to PA was the ratio of the area under the concentration-time curve of PA resulting from an oral dose of PA-N-oxide divided by that from an equimolar dose of PA (Method 1). The second endpoint was the ratio of the amount of pyrrole-protein adducts formed under these conditions (Method 2). REPPANO to PA values appeared to decrease with increasing dose, with the decrease for Method 2 already starting at lower dose level than for Method 1. At dose levels as low as estimated daily human intakes, REPPANO to PA values amounted to 0.92, 0.81, 0.78, and 0.68 for retrorsine N-oxide, seneciphylline N-oxide, riddelliine N-oxide and senecivernine N-oxide, respectively, and became independent of the dose or fraction bioactivated, because no GSH depletion, saturation of PA clearance or PA-N-oxide reduction occurs. Overall, the results demonstrate the strength of using PBK modeling in defining REPPANO to PA values, thereby substantiating the use of the same approach for other PA-N-oxides for which in vivo data are lacking.

Physiologically based kinetic modeling of senecionine N-oxide in rats as a new approach methodology to define the effects of dose and endpoint used on relative potency values of pyrrolizidine alkaloid N-oxides.

Over 1,000 pyrrolizidine alkaloids (PAs) and their N-oxides (PA-N-oxides) occur in 3% of all flowering plants. PA-N-oxides are toxic when reduced to their parent PAs, which are bioactivated into pyrrole intermediates that generate protein- and DNA-adducts resulting in liver toxicity and carcinogenicity. Literature data for senecionine N-oxide in rats indicate that the relative potency (REP) value of this PA-N-oxide compared to its parent PA senecionine varies with the endpoint used. The first endpoint was the ratio between the area under the concentration-time curve (AUC) for senecionine upon dosing senecionine N-oxide or an equimolar dose of senecionine, while the second endpoint was the ratio between the amount for pyrrole-protein adducts formed under these conditions. This study aimed to investigate the mode of action underlying this endpoint dependent REP value for senecionine N-oxide with physiologically based kinetic (PBK) modeling. Results obtained reveal that limitation of 7-GS-DHP adduct formation due to GSH depletion, resulting in increased pyrrole-protein adduct formation, occurs more likely upon high dose oral PA administration than upon an equimolar dose of PA-N-oxide. At high dose levels, this results in a lower REP value when based on pyrrole-protein adduct levels than when based on PA concentrations. At low dose levels, the difference no longer exists. Altogether, the results of the study show how the REP value for senecionine N-oxide depends on dose and endpoint used, and that PBK modeling provides a way to characterize REP values for PA-N-oxides at realistic low dietary exposure levels, thus reducing the need for animal experiments.

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.

Physiologically-Based Kinetic Modeling Predicts Similar In Vivo Relative Potency of Senecionine N-Oxide for Rat and Human at Realistic Low Exposure Levels.

SCOPE: This study aims to determine if previously developed physiologically-based kinetic (PBK) model in rat can be modified for senecionine (SEN) and its N-oxide (SENO), and be used to investigate potential species differences between rat and human in relative potency (REP) of the N-oxide relative to the parent pyrrolizidine alkaloid (PA). METHODS AND RESULTS: In vitro derived kinetic parameters including the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for SENO reduction and SEN clearance are used to define the PBK models. The rat model is validated with published animal data, and the toxicokinetic profiles of SEN from either orally-administered SENO or SEN are simulated. REP values of SENO relative to SEN amount to 0.84 and 0.89 in rat and human, respectively. CONCLUSION: The REP value can be dose- and species-dependent, with the values for rat and human being comparable at low realistic exposure scenarios. In summary, PBK modeling serves as a valuable New Approach Methodology (NAM) tool for predicting REP values of PA-N-oxides and may actually result in more accurate REP values for human risk assessment than what would be defined using in vivo animal experiments.