PBTK model-based analysis of CYP3A4 induction and the toxicokinetics of the pyrrolizidine alkaloid retrorsine in man

Cytochrome P450 (CYP)3A4 induction by drugs and pesticides plays a critical role in the enhancement of pyrrolizidine alkaloid (PA) toxicity as it leads to increased formation of hepatotoxic dehydro-PA metabolites. Addressing the need for a quantitative analysis of this interaction, we developed a physiologically-based toxicokinetic (PBTK) model. Specifically, the model describes the impact of the well-characterized CYP3A4 inducer rifampicin on the kinetics of retrorsine, which is a prototypic PA and contaminant in herbal teas. Based on consumption data, the kinetics after daily intake of retrorsine were simulated with concomitant rifampicin treatment. Strongest impact on retrorsine kinetics (plasma AUC24 and C max reduced to 67% and 74% compared to the rifampicin-free reference) was predicted directly after withdrawal of rifampicin. At this time point, the competitive inhibitory effect of rifampicin stopped, while CYP3A4 induction was still near its maximum. Due to the impacted metabolism kinetics, the cumulative formation of intestinal retrorsine CYP3A4 metabolites increased to 254% (from 10 to 25 nmol), while the cumulative formation of hepatic CYP3A4 metabolites was not affected (57 nmol). Return to baseline PA toxicokinetics was predicted 14 days after stop of a 14-day rifampicin treatment. In conclusion, the PBTK model showed to be a promising tool to assess the dynamic interplay of enzyme induction and toxification pathways.

PBTK modeling of the pyrrolizidine alkaloid retrorsine to predict liver toxicity in mouse and rat

Retrorsine is a hepatotoxic pyrrolizidine alkaloid (PA) found in herbal supplements and medicines, food and livestock feed. Dose-response studies enabling the derivation of a point of departure including a benchmark dose for risk assessment of retrorsine in humans and animals are not available. Addressing this need, a physiologically based toxicokinetic (PBTK) model of retrorsine was developed for mouse and rat. Comprehensive characterization of retrorsine toxicokinetics revealed: both the fraction absorbed from the intestine (78%) and the fraction unbound in plasma (60%) are high, hepatic membrane permeation is dominated by active uptake and not by passive diffusion, liver metabolic clearance is 4-fold higher in rat compared to mouse and renal excretion contributes to 20% of the total clearance. The PBTK model was calibrated with kinetic data from available mouse and rat studies using maximum likelihood estimation. PBTK model evaluation showed convincing goodness-of-fit for hepatic retrorsine and retrorsine-derived DNA adducts. Furthermore, the developed model allowed to translate in vitro liver toxicity data of retrorsine to in vivo dose-response data. Resulting benchmark dose confidence intervals (mg/kg bodyweight) are 24.1-88.5 in mice and 79.9-104 in rats for acute liver toxicity after oral retrorsine intake. As the PBTK model was built to enable extrapolation to different species and other PA congeners, this integrative framework constitutes a flexible tool to address gaps in the risk assessment of PA.

Metabolic Pattern of Hepatotoxic Pyrrolizidine Alkaloids in Liver Cells

Contamination with 1,2-unsaturated pyrrolizidine alkaloids (PAs) is a serious problem for certain phytomedicines, foods, and animal feeds. Several of these PAs are genotoxic and carcinogenic, primarily in the liver, upon cytochrome P450 (CYP)-catalyzed activation into reactive (pyrrolic and pyrrole-like) metabolites. Here we investigated the metabolism of selected PAs (echimidine, europine, lasiocarpine, lycopsamine, retrorsine, and senecionine) in rat hepatocytes in primary culture and in human CYP3A4-transfected HepG2 cells. The open-chained diesters echimidine and lasiocarpine and the cyclic diester senecionine were extensively metabolized in rat hepatocytes into a broad spectrum of products released into the medium. A large portion of unidentified, possibly irreversibly bound, products remained in the cells while detectable amounts of reactive and other metabolites were found in the incubation media. In HepG2-CYP3A4 cells, lasiocarpine was more extensively metabolized than echimidine and senecionine which also gave rise to the release of pyrrolic metabolites. In human cells, no pyrrolic metabolites were detected in retrorsine or lycopsamine incubations, while no such metabolites were detected from europine in both cell types. Other types of metabolic changes comprised modifications such as side chain demethylation or oxygenation reactions like the formation of N-oxides. The latter, considered as a detoxification step, was a major pathway with cyclic diesters, was less distinctive for echimidine and lycopsamine and almost negligible for lasiocarpine and europine. Our data are in agreement with previously published cyto- and genotoxicity findings and suggests that the metabolic pattern may contribute substantially to the specific toxic potency of a certain congener. In addition, marked differences were found for certain congeners between rat hepatocytes and transfected human HepG2 cells, whereby a high level of bioactivation was found for lasiocarpine, whereas a very low level of bioactivation was observed for monoesters, in particular in human cells.

In vitro biotransformation of pyrrolizidine alkaloids in different species: part II-identification and quantitative assessment of the metabolite profile of six structurally different pyrrolizidine alkaloids

Pyrrolizidine alkaloids (PA) exert their toxic effects only after bioactivation. Although their toxicity has already been studied and metabolic pathways including important metabolites were described, the quantification of the latter revealed a large unknown portion of the metabolized PA. In this study, the qualitative and quantitative metabolite profiles of structurally different PAs in rat and human liver microsomes were investigated. Between five metabolites for europine and up to 48 metabolites for lasiocarpine were detected. Proposals for the chemical structure of each metabolite were derived based on fragmentation patterns using high-resolution mass spectrometry. The metabolite profiles of the diester PAs showed a relatively good agreement between both species. The metabolic reactions were summarized into three groups: dehydrogenation, oxygenation, and shortening of necic acid(s). While dehydrogenation of the necine base is considered as bioactivation, both other routes are considered as detoxification steps. The most abundant changes found for open chained diesters were dealkylations, while the major metabolic pathway for cyclic diesters was oxygenation especially at the nitrogen atom. In addition, all diester PAs formed several dehydrogenation products, via the insertion of a second double bond in the necine base, including the formation of glutathione conjugates. In rat liver microsomes, all investigated PAs formed dehydropyrrolizidine metabolites with the highest amount formed by lasiocarpine, whereas in human liver microsomes, these metabolites could only be detected for diesters. Our findings demonstrate that an extensive analysis of PA metabolism can provide the basis for a better understanding of PA toxicity and support future risk assessment.