Physiologically Based Kinetic Models Research Articles

Modeling chemical bioaccumulation in snakes, part 1: Model development

Environmental chemical emission influences ecological health to some extent. Predators (e.g., snakes) could bioaccumulate chemicals along the food chain, which also leaves potential health implications on their reproduction. For the difficulty of collecting related biomatrices for exposure assessment, part 1 of this study proposed a modeling method relying on physiologically based kinetic (PBK) theory to estimate snake chronic exposure to environmental chemicals. In the steady state, the biotransfer factors of chemicals produced by the PBK model can indicate a snake’s chronic internal exposure to environmental chemicals and their potential for bioaccumulation at this level of the food web. Specifically, 3074 organic chemicals were compelled into the dataset for PBK modeling (part 2 of the study). The modeling framework covered the physiological process of the skin to consider shed snakeskin as a potential biomarker for future study. The proposed modeling approach was integrated into a spreadsheet, enabling the modification of input values to simulate outcomes for a wide range of chemical and snake species. The proposed model can help assess the ecological risks of environmental chemicals and quantify their behavior in the food web.

Use of physiologically based kinetic modeling to predict neurotoxicity and genotoxicity of methylglyoxal in humans

This study aimed to evaluate human neurotoxicity and genotoxicity risks from dietary and endogenous methylglyoxal (MGO), utilizing physiologically based kinetic (PBK) modeling-facilitated reverse dosimetry as a new approach methodology (NAM) to extrapolate in vitro toxicity data to in vivo dose-response predictions. A human PBK model was defined based on a newly developed and evaluated mouse model enabling the translation of in vitro toxicity data for MGO from human stem cell-derived neurons and WM-266-4 melanoma cells into quantitative human in vivo toxicity data and subsequent risk assessment by the margin of exposure (MOE) approach. The results show that the MOEs resulting from daily dietary intake did not raise a concern for endpoints for neurotoxicity including mitochondrial function, cytotoxicity, and apoptosis, while those for DNA adduct formation could not exclude a concern over genotoxicity. Endogenous MGO formation, especially under diabetic conditions, resulted in MOEs that raised concern not only for genotoxicity but also for some of the neurotoxicity endpoints evaluated. Thus, the results also point to the importance of taking the endogenous levels into account in the risk assessment of MGO.

Species specific kinetics of imidacloprid and carbendazim in mouse and rat and consequences for biomonitoring

This study aimed to develop physiologically based kinetic (PBK) models to predict the blood concentrations of imidacloprid and carbendazim and their primary metabolites 5-hydroxy-imidacloprid and 2-aminobenzimidazole after single or repeated oral exposure in mouse (Mus musculus), and compare this to corresponding kinetic data in rat (Rattus norvegicus). PBK model constants for conversion of imidacloprid and carbendazim and formation and clearance of their selected primary metabolites were quantified by in vitro mouse liver microsomal and S9 incubations. The performance of the newly developed PBK models was evaluated, based on a comparison to available literature data, showing that the models performed well. Predictions made were also compared to results from PBK model simulations for rats reported previously to obtain insight in species dependent differences in kinetics of these pesticides. The results thus obtained revealed substantial species differences in kinetics for these two pesticides between mouse and rat, especially for imidacloprid and to a lesser extent for carbendazim. Repeated dose PBK model simulations revealed that the models can facilitate estimation of external exposure levels under wildlife conditions based on internal blood concentrations of the parent compound. The rate of conversion and liver volume fraction were shown to influence the accuracy of these predictions with lower values providing less variable outcomes. It is concluded that PBK modeling provides a new approach methodology of use for wildlife biomonitoring studies and that results of the present study facilitate benchmarking of the species and compounds for which kinetics enable this with sufficient accuracy.

Physiologically based Kinetic Modeling-Facilitated Quantitative In Vitro to In Vivo Extrapolation to Predict the Effects of Aloe-Emodin in Rats and Humans.

Aloe-emodin, a natural hydroxyanthraquinone, exerts both adverse and protective effects. This study aimed at investigating these potential effects of aloe-emodin in humans upon the use of food supplements and herbal medicines using a physiologically based kinetic (PBK) modeling-facilitated quantitative in vitro to in vivo extrapolation (QIVIVE) approach. For this, PBK models in rats and humans were established for aloe-emodin including its active metabolite rhein and used to convert in vitro data on hepatotoxicity, nephrotoxicity, reactive oxidative species (ROS) generation, and Nrf2 induction to corresponding in vivo dose-response curves, from which points of departure (PODs) were derived by BMD analysis. The derived PODs were subsequently compared to the estimated daily intakes (EDIs) resulting from the use of food supplements or herbal medicines. It is concluded that the dose levels of aloe-emodin from food supplements or herbal medicines are unlikely to induce toxicity, ROS generation, or Nrf2 activation in liver and kidney.

Advancing PFAS risk assessment: Integrative approaches using agent-based modelling and physiologically-based kinetic for environmental and health safety

Per- and polyfluoroalkyl substances (PFAS), ubiquitous in a myriad of consumer and industrial products, and depending on the doses of exposure represent a hazard to both environmental and public health, owing to their persistent, mobile, and bio accumulative properties. These substances exhibit long half-lives in humans and can induce potential immunotoxic effects at low exposure levels, sparking growing concerns. While the European Food Safety Authority (EFSA) has assessed the risk to human health related to the presence of PFAS in food, in which a reduced antibody response to vaccination in infants was considered as the most critical human health effect, a comprehensive grasp of the molecular mechanisms spearheading PFAS-induced immunotoxicity is yet to be attained. Leveraging modern computational tools, including the Agent-Based Model (ABM) Universal Immune System Simulator (UISS) and Physiologically Based Kinetic (PBK) models, a deeper insight into the complex mechanisms of PFAS was sought. The adapted UISS serves as a vital tool in chemical risk assessments, simulating the host immune system's reactions to diverse stimuli and monitoring biological entities within specific adverse health contexts. In tandem, PBK models unravelling PFAS' biokinetics within the body i.e. absorption, distribution, metabolism, and elimination, facilitating the development of time-concentration profiles from birth to 75 years at varied dosage levels, thereby enhancing UISS-TOX's predictive abilities. The integrated use of these computational frameworks shows promises in leveraging new scientific evidence to support risk assessments of PFAS. This innovative approach not only allowed to bridge existing data gaps but also unveiled complex mechanisms and the identification of unanticipated dynamics, potentially guiding more informed risk assessments, regulatory decisions, and associated risk mitigations measures for the future.

PBK models to predict internal and external dose levels following oral exposure of rats to imidacloprid and carbendazim

Monitoring oral exposure to pesticides in wildlife is crucial for assessing environmental risks and preventing adverse effects on non-target species. Traditionally, this requires invasive tissue sampling, raising ethical, regulatory, and economic concerns. To address this gap, our study aims to develop a method for assessing external oral dose levels in rats using physiologically-based kinetic (PBK) modeling based on blood concentration levels of two pesticides, imidacloprid and carbendazim, and one of their primary metabolites. We utilized in vitro metabolic kinetic data from hepatic microsomal and S9 incubations to inform our models. These models were then evaluated by comparing their predictions with existing in vivo experimental data from the literature. Our results demonstrate that the models provide accurate predictions, presenting a novel in vitro and in silico approach for environmental exposure and risk assessment of pesticides. This methodology has the potential for application in wildlife species, advancing the frontier of knowledge in non-invasive pesticide exposure assessment.

Physiologically based kinetic (PBK) modeling of propiconazole using a machine learning-enhanced read-across approach for interspecies extrapolation

A significant challenge in the traditional human health risk assessment of agrochemicals is the uncertainty in quantifying the interspecies differences between animal models and humans. To work toward a more accurate and animal-free risk determination, new approaches such as physiologically based kinetic (PBK) modeling have been used to perform dosimetry extrapolation from animals to humans. However, the regulatory use and acceptance of PBK modeling is limited for chemicals that lack in vivo animal pharmacokinetic (PK) data, given the inability to evaluate models. To address these challenges, this study developed PBK models in the absence of in vivo PK data for the fungicide propiconazole, an activator of constitutive androstane receptor (CAR)/pregnane X receptor (PXR). A fit-for-purpose read-across approach was integrated with hierarchical clustering − an unsupervised machine learning algorithm, to bridge the knowledge gap. The integration allowed the incorporation of a broad spectrum of attributes for analog consideration, and enabled the analog selection in a simple, reproducible, and objective manner. The applicability was evaluated and demonstrated using penconazole (source) and three pseudo-unknown target chemicals (epoxiconazole, tebuconazole and triadimefon). Applying this machine learning-enhanced read-across approach, difenoconazole was selected as the most appropriate analog for propiconazole. A mouse PBK model was developed and evaluated for difenoconazole (source), with the mode of action of CAR/PXR activation incorporated to simulate the in vivo autoinduction of metabolism. The difenoconazole mouse model then served as a template for constructing the propiconazole mouse model. A parallelogram approach was subsequently applied to develop the propiconazole rat and human models, enabling a quantitative assessment of interspecies differences in dosimetry. This integrated approach represents a substantial advancement toward refining risk assessment of propiconazole within the framework of animal alternative safety assessment strategies.

From hazard to risk prioritization: a case study to predict drug-induced cholestasis using physiologically based kinetic modeling

Cholestasis is characterized by hepatic accumulation of bile acids. Clinical manifestation of cholestasis only occurs in a small proportion of exposed individuals. The present study aims to develop a new approach methodology (NAM) to predict drug-induced cholestasis as a result of drug-induced hepatic bile acid efflux inhibition and the resulting bile acid accumulation. To this end, hepatic concentrations of a panel of drugs were predicted by a generic physiologically based kinetic (PBK) drug model. Their effects on hepatic bile acid efflux were incorporated in a PBK model for bile acids. The predicted bile acid accumulation was used as a measure for a drug’s cholestatic potency. The selected drugs were known to inhibit hepatic bile acid efflux in an assay with primary suspension-cultured hepatocytes and classified as common, rare, or no for cholestasis incidence. Common cholestasis drugs included were atorvastatin, chlorpromazine, cyclosporine, glimepiride, ketoconazole, and ritonavir. The cholestasis incidence of the drugs appeared not to be adequately predicted by their Ki for inhibition of hepatic bile acid efflux, but rather by the AUC of the PBK model predicted internal hepatic drug concentration at therapeutic dose level above this Ki. People with slower drug clearance, a larger bile acid pool, reduced bile salt export pump (BSEP) abundance, or given higher than therapeutic dose levels were predicted to be at higher risk to develop drug-induced cholestasis. The results provide a proof-of-principle of using a PBK-based NAM for cholestasis risk prioritization as a result of transporter inhibition and identification of individual risk factors.

Systematic evaluation of high-throughput PBK modelling strategies for the prediction of intravenous and oral pharmacokinetics in humans

Physiologically based kinetic (PBK) modelling offers a mechanistic basis for predicting the pharmaco-/toxicokinetics of compounds and thereby provides critical information for integrating toxicity and exposure data to replace animal testing with in vitro or in silico methods. However, traditional PBK modelling depends on animal and human data, which limits its usefulness for non-animal methods. To address this limitation, high-throughput PBK modelling aims to rely exclusively on in vitro and in silico data for model generation. Here, we evaluate a variety of in silico tools and different strategies to parameterise PBK models with input values from various sources in a high-throughput manner. We gather 2000 + publicly available human in vivo concentration–time profiles of 200 + compounds (IV and oral administration), as well as in silico, in vitro and in vivo determined compound-specific parameters required for the PBK modelling of these compounds. Then, we systematically evaluate all possible PBK model parametrisation strategies in PK-Sim and quantify their prediction accuracy against the collected in vivo concentration–time profiles. Our results show that even simple, generic high-throughput PBK modelling can provide accurate predictions of the pharmacokinetics of most compounds (87% of Cmax and 84% of AUC within tenfold). Nevertheless, we also observe major differences in prediction accuracies between the different parameterisation strategies, as well as between different compounds. Finally, we outline a strategy for high-throughput PBK modelling that relies exclusively on freely available tools. Our findings contribute to a more robust understanding of the reliability of high-throughput PBK modelling, which is essential to establish the confidence necessary for its utilisation in Next-Generation Risk Assessment.

Physiologically based kinetic (PBK) modeling as a new approach methodology (NAM) for predicting systemic levels of gut microbial metabolites

There is a clear need to develop new approach methodologies (NAMs) that combine in vitro and in silico testing to reduce and replace animal use in chemical risk assessment. Physiologically based kinetic (PBK) models are gaining popularity as NAMs in toxico/pharmacokinetics, but their coverage of complex metabolic pathways occurring in the gut are incomplete. Chemical modification of xenobiotics by the gut microbiome plays a critical role in the host response, for example, by prolonging exposure to harmful metabolites, but there is not a comprehensive approach to quantify this impact on human health. There are examples of PBK models that have implemented gut microbial biotransformation of xenobiotics with the gut as a dedicated metabolic compartment. However, the integration of microbial metabolism and parameterization of PBK models is not standardized and has only been applied to a few chemical transformations. A challenge in this area is the measurement of microbial metabolic kinetics, for which different fermentation approaches are used. Without a standardized method to measure gut microbial metabolism ex vivo/in vitro, the kinetic constants obtained will lead to conflicting conclusions drawn from model predictions. Nevertheless, there are specific cases where PBK models accurately predict systemic concentrations of gut microbial metabolites, offering potential solutions to the challenges outlined above. This review focuses on models that integrate gut microbial bioconversions and use ex vivo/in vitro methods to quantify metabolic constants that accurately represent in vivo conditions.

Parameter grouping and co-estimation in physiologically based kinetic models using genetic algorithms.

Physiologically based kinetic (PBK) models are widely used in pharmacology and toxicology for predicting the internal disposition of substances upon exposure, voluntarily or not. Due to their complexity, a large number of model parameters need to be estimated, either through in silico tools, in vitro experiments, or by fitting the model to in vivo data. In the latter case, fitting complex structural models on in vivo data can result in overparameterization and produce unrealistic parameter estimates. To address these issues, we propose a novel parameter grouping approach, which reduces the parametric space by co-estimating groups of parameters across compartments. Grouping of parameters is performed using genetic algorithms and is fully automated, based on a novel goodness-of-fit metric. To illustrate the practical application of the proposed methodology, two case studies were conducted. The first case study demonstrates the development of a new PBK model, while the second focuses on model refinement. In the first case study, a PBK model was developed to elucidate the biodistribution of titanium dioxide (TiO2) nanoparticles in rats following intravenous injection. A variety of parameter estimation schemes were employed. Comparative analysis based on goodness-of-fit metrics demonstrated that the proposed methodology yields models that outperform standard estimation approaches, while utilizing a reduced number of parameters. In the second case study, an existing PBK model for perfluorooctanoic acid (PFOA) in rats was extended to incorporate additional tissues, providing a more comprehensive portrayal of PFOA biodistribution. Both models were validated through independent in vivo studies to ensure their reliability.

Nano- and microplastic PBK modeling in the context of human exposure and risk assessment

Insufficient data on nano- and microplastics (NMP) hinder robust evaluation of their potential health risks. Methodological disparities and the absence of established toxicity thresholds impede the comparability and practical application of research findings. The diverse attributes of NMP, such as variations in sizes, shapes, and compositions, complicate human health risk assessment. Although probability density functions (PDFs) show promise in capturing this diversity, their integration into risk assessment frameworks is limited. Physiologically based kinetic (PBK) models offer a potential solution to bridge the gap between external exposure and internal dosimetry for risk evaluation. However, the heterogeneity of NMP poses challenges for accurate biodistribution modeling.A literature review, encompassing both experimental and modeling studies, was conducted to examine biodistribution studies of monodisperse micro- and nanoparticles. The literature search in PubMed and Scopus databases yielded 39 studies that met the inclusion criteria. Evaluation criteria were adapted from previous Quality Assurance and Quality Control (QA-QC) studies, best practice guidelines from WHO (2010), OECD guidance (2021), and additional criteria specific to NMP risk assessment. Subsequently, a conceptual framework for a comprehensive NMP-PBK model was developed, addressing the multidimensionality of NMP particles.Parameters for an NMP-PBK model are presented. QA-QC evaluations revealed that most experimental studies scored relatively well (>0) in particle characterizations and environmental settings but fell short in criteria application for biodistribution modeling. The evaluation of modeling studies revealed that information regarding the model type and allometric scaling requires improvement. Three potential applications of PDFs in PBK modeling of NMP are identified: capturing the multidimensionality of the NMP continuum, quantifying the probabilistic definition of external exposure, and calculating the bio-accessibility fraction of NMP in the human body. A framework for an NMP-PBK model is proposed, integrating PDFs to enhance the assessment of NMP’s impact on human health.

Assessing chronic gestational exposure to environmental chemicals in pregnant women: Advancing the co-PBK model

Vulnerable populations, such as pregnant women and their fetuses, confront potential health risks due to exposure to environmental toxic compounds. Computational methods have been popular in assessing chemical exposure to populations, contrasting with traditional cohort studies for human biomonitoring. This study proposes a screening-level approach based on physiologically based kinetic (PBK) modeling to evaluate the steady-state exposure of pregnant women to environmental chemicals throughout pregnancy. To exemplify the modeling application, naphthalene was chosen. Simulation results indicated that maternal fat exhibited significant bioaccumulation potential, with the log-transformed BTF of naphthalene at 0.51 mg kg−1 per mg d−1 in the steady state. The placenta was primarily exposed to 0.83 mg/d naphthalene for a 75.2 kg pregnant woman, considering all exposure routes. In the fetal structure, single-organ fetal PBK modeling estimated a naphthalene exposure of 123.64 mg/d to the entire fetus, while multiple-organ fetal PBK modeling further revealed the bioaccumulation highest in fat tissue. The liver identified as the vital organ for metabolism, kBioT,LiverM was demonstrated with the highest sensitivity among rate constants in the maternal body. Furthermore, the first-order kinetic rate constants related to the placenta and blood were found to impact the distribution process of naphthalene in the fetus, influencing gestational exposure. In conclusion, urgent attention is needed to develop a computational biomonitoring tool for assessing toxic chemical exposure in vulnerable populations.

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.

Using Read-Across to build Physiologically-Based Kinetic models: Part 2. Case studies for atenolol and flumioxazin

Read-across, wherein information from a data-rich chemical is used to make a prediction for a similar chemical that lacks the relevant data, is increasingly being accepted as an alternative to animal testing. Identifying chemicals that can be considered as similar (analogues) is crucial to the process. Two resources have been developed previously to address the issue of analogue selection and facilitate physiologically-based kinetic (PBK) model development, using read-across. Chemical-specific PBK models, available in the literature, were collated to form a PBK model dataset (PMD) of over 7,500 models. A KNIME workflow was created to accompany this PMD that can aid the selection of appropriate chemical analogues from chemicals within this dataset (i.e. chemicals that are similar to a target of interest and are known to have an existing PBK model). Information from the PBK model for the source chemical can then be used in a read-across approach to inform the development of a new PBK model for the target. The application of these resources is tested here using two case studies (i) for the drug atenolol and (ii) for the plant protection product, flumioxazin. New PBK models were constructed for these two target chemicals using data obtained from source chemicals, identified by the workflow as being similar (analogues). In each case, the published PBK model for the source chemical was initially reproduced, as accurately as possible, before being adapted and used as a template for the target chemical. The performance of the new PBK models was assessed by comparing simulation outputs to existing data on key kinetic properties for the targets. The results demonstrate that a read-across approach can be successfully applied to develop new PBK models for data-poor chemicals, thus enabling their deployment during early-stage risk assessment. This assists prediction of internal exposure whilst reducing reliance on animal testing.

Using read-across to build physiologically-based kinetic models: Part 1. Development of a KNIME workflow to assist analogue selection for PBK modelling

Read-across refers to the process by which information from one (source) chemical is used to infer information about another similar (target) chemical. This method can be used to fill data gaps and so inform safety assessment where data are lacking for chemicals of interest. As one chemical cannot be considered as absolutely similar to another, only similar with respect to a given property, it is essential to justify the selection of similar chemicals (analogues) for the purposes of read-across. A previously created dataset of available physiologically-based kinetic (PBK) models (referred to as the PBK modelling dataset or PMD) was used in the development of a KNIME workflow. KNIME is a freely-available, open-source analytics platform that allows users to create workflows to analyse and visualise data. The KNIME workflow described here was designed to identify chemical analogues with a corresponding model in the PMD. The PMD combined with the KWAAS enables PBK model information from source chemical(s) to be used in a read-across approach to help develop new PBK models for target chemicals. This KNIME workflow was applied to six chemicals, representing different types of chemical classes (drugs, cosmetics, botanicals, industrial chemicals, pesticides, and food additives) to assess its applicability across various industries. Information acquired from these PBK models can be used to support safety assessment of chemicals and reduce reliance on animal testing.

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.

ADME characterization and PBK model development of 3 highly protein-bound UV filters through topical application.

Estimating human exposure in the safety assessment of chemicals is crucial. Physiologically based kinetic (PBK) models which combine information on exposure, physiology, and chemical properties, describing the absorption, distribution, metabolism, and excretion (ADME) processes of a chemical, can be used to calculate internal exposure metrics such as maximum concentration and area under the concentration-time curve in plasma or tissues of a test chemical in next-generation risk assessment. This article demonstrates the development of PBK models for 3 UV filters, specifically octyl methoxycinnamate, octocrylene, and 4-methylbenzylidene camphor. The models were parameterized entirely based on data obtained from in vitro and/or in silico methods in a bottom-up modeling approach and then validated based on human dermal pharmacokinetic (PK) data. The 3 UV filters are "difficult to test" in in vitro test systems due to high lipophilicity, high binding affinity for proteins, and nonspecific binding, for example, toward plastic. This research work presents critical considerations in ADME data generation, interpretation, and parameterization to assure valid PBK model development to increase confidence in using PBK modeling to help make safety decisions in the absence of human PK data. The developed PBK models of the 3 chemicals successfully simulated the plasma concentration profiles of clinical PK data following dermal application, indicating the reliability of the ADME data generated and the parameters determined. The study also provides insights and lessons learned for characterizing ADME and developing PBK models for highly lipophilic and protein-bound chemicals in the future.

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.

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.