Prediction of in vivo developmental toxicity of all-trans-retinoic acid based on in vitro toxicity data and in silico physiologically based kinetic modeling

New approach methodologies to confirm developmental toxicity of pharmaceuticals based on weight of evidence

The aim of embryo-fetal developmental toxicity assessments for pharmaceuticals is to inform potential risk of adverse pregnancy outcome, which has traditionally relied on studies in pregnant animals. Recent updates to international safety guidelines (ICH S5R3) have incorporated information on how to use weight of evidence and alternative assays to reduce animal use while still informing risk of fetal harm. Uptake of these alternative approaches has been slow due to limitations in understanding how alternative assays translate to in vivo effects and then relevance to human exposure. To understand the predictivity of new approach methodologies for developmental toxicity (DevTox NAMs), we used two pharmaceutical examples (glasdegib and lorlatinib) to illustrate the value of DevTox NAMs to complement weight of evidence (WoE) assessments while considering the relationship of concentration-effect levels in NAMs to in vivo studies. The in vitro results generated in a battery of assays (mEST, rWEC, zebrafish, and human based stem cells) confirmed the WoE based on literature and further confirmed by preliminary embryo-fetal development data. The data generated for these two compounds supports integrating DevTox NAMs into the developmental toxicity assessment for advanced cancer indications.

Bayesian analysis of physiologically based toxicokinetic (PBTK) modeling for pentachlorophenol exposure in pregnant women

Pentachlorophenol (PCP) is a persistent organic compound that is widely present in the environment. The estimation of internal exposure levels for a given external exposure using toxicokinetic models is key to the human health risk assessment of PCP. The present study developed a physiologically based multicompartmental pharmacokinetic (PBTK) model to describe and predict the behavior of pentachlorophenol (PCP) in an organism. The model consists of stomach, intestines, adipose tissue, kidneys and fast- and poorly perfused tissues that are interconnected via blood circulation. We constructed a PBTK model of PCP in rats and extrapolated it to human dietary PCP exposure. The toxicokinetic data of PCP in human tissues and excreta were obtained from the published literature. Based on the collected PCP dietary survey and internal exposure data of pregnant women in Shanghai, Bayesian statistical analysis was performed for the model using Markov chain Monte Carlo (MCMC) simulation. The posterior distributions of the sensitive parameters were estimated, and the model was parameter optimized and validated using the pregnant women's test dataset. The results showed that the root mean square error (RMSE) improved 37.3% compared to the original model, and a systematic literature search revealed that the optimized model achieved acceptable prediction results for other datasets in China. A PCP metabolism model based on the exposure characteristics of pregnant women in China was constructed in the present study. The model provides a theoretical basis for the study of PCP toxicity and risk assessment.

Perspectives on chick embryo models in developmental and reproductive toxicity screening

Teratology, the study of congenital anomalies and their causative factors intersects with developmental and reproductive toxicology, employing innovative methodologies. Evaluating the potential impacts of teratogens on fetal development and assessing human risk is an essential prerequisite in preclinical research. The chicken embryo model has emerged as a powerful tool for understanding human embryonic development due to its remarkable resemblance to humans. This model offers a unique platform for investigating the effects of substances on developing embryos, employing techniques such as ex ovo and in ovo assays, chorioallantoic membrane assays, and embryonic culture techniques. The advantages of chicken embryonic models include their accessibility, cost-effectiveness, and biological relevance to vertebrate development, enabling efficient screening of developmental toxicity. However, these models have limitations, such as the absence of a placenta and maternal metabolism, impacting the study of nutrient exchange and hormone regulation. Despite these limitations, understanding and mitigating the challenges posed by the absence of a placenta and maternal metabolism are critical for maximizing the utility of the chick embryo model in developmental toxicity testing. Indeed, the insights gained from utilizing these assays and their constraints can significantly contribute to our understanding of the developmental impacts of various agents. This review underscores the utilization of chicken embryonic models in developmental toxicity testing, highlighting their advantages and disadvantages by addressing the challenges posed by their physiological differences from mammalian systems.

Alternatives to animal testing in toxicity testing: Current status and future perspectives in food safety assessments

The development of alternative methods to animal testing has gained great momentum since Russel and Burch introduced the "3Rs" concept of Reduction, Refinement, and Replacement of animals in safety testing in 1959. Several alternatives to animal testing have since been introduced, including but not limited to in vitro and in chemico test systems, in silico models, and computational models (e.g., [quantitative] structural activity relationship models, high-throughput screens, organ-on-chip models, and genomics or bioinformatics) to predict chemical toxicity. Furthermore, several agencies have developed robust integrated testing strategies to determine chemical toxicity. The cosmetics sector is pioneering the adoption of alternative methodologies for safety evaluations, and other sectors are aiming to completely abandon animal testing by 2035. However, beyond the use of in vitro genetic testing, agencies regulating the food industry have been slow to implement alternative methodologies into safety evaluations compared with other sectors; setting health-based guidance values for food ingredients requires data from systemic toxicity, and to date, no standalone validated alternative models to assess systemic toxicity exist. The abovementioned models show promise for assessing systemic toxicity with further research. In this paper, we review the current alternatives and their applicability and limitations in food safety evaluations.

Anticancer Activities of Tetrasubstituted Imidazole-Pyrimidine-Sulfonamide Hybrids as Inhibitors of EGFR Mutants

A new series of tetrasubstituted imidazole derivatives carrying pyrimidine sulfonamide pharmacophores has been synthesized and evaluated for their anticancer activities. In-vitro screening of these hybrids against a full 60-cell-line panel at a single dose of 10 μM showed significant growth inhibition of up to 95 %. The most active compound showed in-vitro anticancer activities against (i) abnormal HER2 and (ii) two mutants for EGFR. Apoptotic gene expression revealed that lead compounds induced MCF-7 cell line apoptosis together with considerable change in the Bax/Bcl-2 expression ratio. One lead compound led to a significant cell-cycle S-phase arrest, while another blocked the cell cycle at G1/S-phase causing the accumulation of cells. Docking analysis of these two hybrids adopted the orientation and binding interactions with a higher liability to enter the active side pocket of HER2, L858R, and T790 M, preferable to that of co-crystallized ligands. Modelling simulation was consistent with the acquired biological evaluation.

New approach methodologies: A quantitative in vitro to in vivo extrapolation case study with PFASs

PER: and polyfluoroalkyl substances (PFASs) have been associated with increased blood lipids in humans. Perfluorooctanoic acid (PFOA) has been also linked with elevated alanine transferase (ALT) serum levels in humans, and in rodents the liver is a main target organ for many PFASs. With the focus on New Approach Methodologies, the chronic oral equivalent effect doses were calculated for PFOA, PFNA (perfluorononanoic acid), PFHxS (perfluorohexanesulfonic acid) and PFOS (perfluorooctane sulfonic acid) based on in vitro effects measured in the HepaRG cell line. Selected in vitro readouts were considered biomarkers for lipid disturbances and hepatotoxicity. Concentration-response data obtained from HepaRG cells on triglyceride (TG) accumulation and expression changes of 12 selected genes (some involved in cholesterol homeostasis) were converted into corresponding human dose-response data, using physiologically based kinetic (PBK) model-facilitated reverse dosimetry. Next to this, the biokinetics of the chemicals were studied in the cell system. The current European dietary PFASs exposure overlaps with the calculated oral equivalent effect doses, indicating that the latter may lead to interference with hepatic gene expression and lipid metabolism. These findings illustrate an in vitro-in silico methodology, which can be applied for more PFASs, to select those that should be prioritized for further hazard characterization.

New [1,2,4]triazolo[4,3-c]quinazolines as intercalative Topo II inhibitors: Design, synthesis, biological evaluation, and in silico studies

Fifteen quinazoline derivatives were designed and synthesized as DNA intercalators. The cytotoxicity of the designed members was assessed against HCT-116 and HepG2 cancer cell lines. In addition, the topoisomerase II (Topo II) inhibitory effect was assessed. Compound 16 was the most cytotoxic and Topo II inhibitor with low cytotoxicity against Vero cells. Compounds 16, 17, and 18 showed significant DNA binding affinities. Compound 16 showed Topo II catalytic inhibitory effect at a concentration of 10 μM. Further mechanistic investigations revealed the capability of compound 16 to induce apoptosis in HCT-116 cells and arrest the growth at the S and G2/M phases. Also, compound 16 showed a significant increase in the level of BAX (2.18-fold) and a marked decrease in the level of Bcl-2 (1.9-fold) compared to the control cells. In silico studies revealed the ability of the synthesized members to bind to the DNA-Topo II complex.

Advancing New Approach Methodologies (NAMs) for Tobacco Harm Reduction: Synopsis from the 2021 CORESTA SSPT-NAMs Symposium

New approach methodologies (NAMs) are emerging chemical safety assessment tools consisting of in vitro and in silico (computational) methodologies intended to reduce, refine, or replace (3R) various in vivo animal testing methods traditionally used for risk assessment. Significant progress has been made toward the adoption of NAMs for human health and environmental toxicity assessment. However, additional efforts are needed to expand their development and their use in regulatory decision making. A virtual symposium was held during the 2021 Cooperation Centre for Scientific Research Relative to Tobacco (CORESTA) Smoke Science and Product Technology (SSPT) conference (titled "Advancing New Alternative Methods for Tobacco Harm Reduction"), with the goals of introducing the concepts and potential application of NAMs in the evaluation of potentially reduced-risk (PRR) tobacco products. At the symposium, experts from regulatory agencies, research organizations, and NGOs shared insights on the status of available tools, strengths, limitations, and opportunities in the application of NAMs using case examples from safety assessments of chemicals and tobacco products. Following seven presentations providing background and application of NAMs, a discussion was held where the presenters and audience discussed the outlook for extending the NAMs toxicological applications for tobacco products. The symposium, endorsed by the CORESTA In Vitro Tox Subgroup, Biomarker Subgroup, and NextG Tox Task Force, illustrated common ground and interest in science-based engagement across the scientific community and stakeholders in support of tobacco regulatory science. Highlights of the symposium are summarized in this paper.

Safety assessment of natural products in Malaysia: current practices, challenges, and new strategies

The belief that natural products are inherently safe is a primary reason for consumers to choose traditional medicines and herbal supplements for health maintenance and disease prevention. Unfortunately, some natural products on the market have been found to contain toxic compounds, such as heavy metals and microbes, as well as banned ingredients such as aristolochic acids. It shows that the existing regulatory system is inadequate and highlights the importance of thorough safety evaluations. In Malaysia, the National Pharmaceutical Regulatory Agency is responsible for the regulatory control of medicinal products and cosmetics, including natural products. For registration purpose, the safety of natural products is primarily determined through the review of documents, including monographs, research articles and scientific reports. One of the main factors hampering safety evaluations of natural products is the lack of toxicological data from animal studies. However, international regulatory agencies such as the European Food Safety Authority and the United States Food and Drug Administration are beginning to accept data obtained using alternative strategies such as non-animal predictive toxicological tools. Our paper discusses the use of state-of-the-art techniques, including chemometrics, in silico modelling and omics technologies and their applications to the safety assessments of natural products.

IVIVE: Facilitating the Use of In Vitro Toxicity Data in Risk Assessment and Decision Making

During the past few decades, the science of toxicology has been undergoing a transformation from observational to predictive science. New approach methodologies (NAMs), including in vitro assays, in silico models, read-across, and in vitro to in vivo extrapolation (IVIVE), are being developed to reduce, refine, or replace whole animal testing, encouraging the judicious use of time and resources. Some of these methods have advanced past the exploratory research stage and are beginning to gain acceptance for the risk assessment of chemicals. A review of the recent literature reveals a burst of IVIVE publications over the past decade. In this review, we propose operational definitions for IVIVE, present literature examples for several common toxicity endpoints, and highlight their implications in decision-making processes across various federal agencies, as well as international organizations, including those in the European Union (EU). The current challenges and future needs are also summarized for IVIVE. In addition to refining and reducing the number of animals in traditional toxicity testing protocols and being used for prioritizing chemical testing, the goal to use IVIVE to facilitate the replacement of animal models can be achieved through their continued evolution and development, including a strategic plan to qualify IVIVE methods for regulatory acceptance.

Use of physiologically based kinetic modeling-facilitated reverse dosimetry to predict in vivo acute toxicity of tetrodotoxin in rodents

In this study, the ability of a new in vitro/in silico quantitative in vitro–in vivo extrapolation (QIVIVE) methodology was assessed to predict the in vivo neurotoxicity of tetrodotoxin (TTX) in rodents. In vitro concentration–response data of TTX obtained in a multielectrode array assay with primary rat neonatal cortical cells and in an effect study with mouse neuro-2a cells were quantitatively extrapolated into in vivo dose–response data, using newly developed physiologically based kinetic (PBK) models for TTX in rats and mice. Incorporating a kidney compartment accounting for active renal excretion in the PBK models proved to be essential for its performance. To evaluate the predictions, QIVIVE-derived dose–response data were compared with in vivo data on neurotoxicity in rats and mice upon oral and parenteral dosing. The results revealed that for both rats and mice the predicted dose–response data matched the data from available in vivo studies well. It is concluded that PBK modeling-based reserve dosimetry of in vitro TTX effect data can adequately predict the in vivo neurotoxicity of TTX in rodents, providing a novel proof-of-principle for this methodology.

Integrating in vitro chemical transplacental passage into a generic PBK model: A QIVIVE approach

With the increasing application of cell culture models as primary tools for predicting chemical safety, the quantitative extrapolation of the effective dose from in vitro to in vivo (QIVIVE) is of increasing importance. For developmental toxicity this requires scaling the in vitro observed dose-response characteristics to in vivo fetal exposure, while integrating maternal in vivo kinetics during pregnancy, in particular transplacental transfer. Here the transfer of substances across the placental barrier, has been studied using the in vitro BeWo cell assay and six embryotoxic compounds of different kinetic complexity. The BeWo assay results were incorporated in an existing generic Physiologically Based Kinetic (PBK) model which for this purpose was extended with rat pregnancy. Finally, as a "proof of principle", the BeWo PBK model was used to perform a QIVIVE based on developmental toxicity as observed in various different in vitro toxicity assays. The BeWo results illustrated different transport profiles of the chemicals across the BeWo monolayer, allocating the substances into two distinct groups: the 'quickly-transported' and the 'slowly-transported'. BeWo PBK exposure simulations during gestation were compared to experimentally measured maternal blood and fetal concentrations and a reverse dosimetry approach was applied to translate in vitro observed embryotoxicity into equivalent in vivo dose-response curves. This approach allowed for a direct comparison of the in vitro dose-response characteristics as observed in the Whole Embryo Culture (WEC), and the Embryonic Stem Cell test (cardiac:ESTc and neural:ESTn) with in vivo rat developmental toxicity data. Overall, the in vitro to in vivo comparisons suggest a promising future for the application of such QIVIVE methodologies for screening and prioritization purposes of developmental toxicants. Nevertheless, the clear need for further improvements is acknowledged for a wider application of the approach in chemical safety assessment.

Discovery of new VEGFR-2 inhibitors based on bis([1, 2, 4]triazolo)[4,3-a:3',4'-c]quinoxaline derivatives as anticancer agents and apoptosis inducers

Herein, a new wave of bis([1, 2, 4]triazolo)[4,3-a:3',4'-c]quinoxaline derivatives have been successfully designed and synthesised. The synthesised derivatives were biologically investigated for their cytotoxic activities against HepG2 and MCF-7. Also, the tested compounds were further examined in vitro for their VEGFR-2 inhibitory activity. The most promising derivative 23j was further investigated for its apoptotic behaviour in HepG2 cell lines using flow cytometric and western-plot analyses. Additional in-silico studies were performed to predict how the synthesised compounds can bind to VEGFR-2 and to determine the drug-likeness profiling of these derivatives. The results revealed that compounds 23a, 23i, 23j, 23l, and 23n displayed the highest antiproliferative activities against the two cell lines with IC50 values ranging from 6.4 to 19.4 µM. Furthermore, compounds 23a, 23d, 23h, 23i, 23j, 23l, 23 m, and 23n showed the highest VEGFR-2 inhibitory activities with IC50 values ranging from 3.7 to 11.8 nM, comparing to sorafenib (IC50 = 3.12 nM). Moreover, compound 23j arrested the HepG2 cell growth at the G2/M phase and induced apoptosis by 40.12% compared to the control cells (7.07%). As well, such compound showed a significant increase in the level of caspase-3 (1.36-fold), caspase-9 (2.80-fold), and BAX (1.65-fold), and exhibited a significant decrease in Bcl-2 level (2.63-fold).

PBK Model-Based Prediction of Intestinal Microbial and Host Metabolism of Zearalenone and Consequences for its Estrogenicity

Scope

The aim of the present study is to develop physiologically‐based kinetic (PBK) models for rat and human that include intestinal microbial and hepatic metabolism of zearalenone (ZEN) in order to predict systemic concentrations of ZEN and to obtain insight in the contribution of metabolism by the intestinal microbiota to the overall metabolism of ZEN.

Methods and Results

In vitro derived kinetic parameters, apparent maximum velocities (V_max) and Michaelis–Menten constants (K_m) for liver and intestinal microbial metabolism of ZEN are included in the PBK models. The models include a sub‐model for the metabolite, α‐zearalenol (α‐ZEL), a metabolite known to be 60‐times more potent as an estrogen than ZEN. Integrating intestinal microbial ZEN metabolism into the PBK models revealed that hepatic metabolism drives the formation of α‐ZEL. Furthermore, the models predicted that at the tolerable daily intake (TDI) of 0.25 µg kg⁻¹ bw the internal concentration of ZEN and α‐ZEL are three‐orders of magnitude below concentrations reported to induce estrogenicity in vitro.

Conclusion

It is concluded that combining kinetic data on liver and intestinal microbial metabolism in a PBK model facilitates a holistic view on the role of the intestinal microbiota in the overall metabolism of the foodborne xenobiotic ZEN and its bioactivation to α‐ZEL.

Design and Synthesis of New Quinoxaline Derivatives as Potential Histone Deacetylase Inhibitors Targeting Hepatocellular Carcinoma: In Silico, In Vitro, and SAR Studies

Guided by the structural optimization principle and the promising anticancer effect of the quinoxaline nucleus, a new series of novel HDAC inhibitors were designed and synthesized. The synthesized compounds were designed to bear the reported pharmacophoric features of the HDAC inhibitors in addition to an extra moiety to occupy the non-used vacant deep pocket of the HDAC receptor. The newly prepared compounds were evaluated for their in vitro anti-proliferative activities against HepG-2 and HuH-7 liver cancer cell lines. The tested compounds showed promising anti-proliferative activities against both cell lines. The most active ten candidates (6 c , 6 d , 6 f , 6 g , 6 k , 6 l , 7 b , 8, 10 h , and 12) were further evaluated for their effect on the gene expression levels of Bax as an apoptotic marker and Bcl-2 as an anti-apoptotic one. Moreover, they were evaluated for their ability to inhibit histone deacetylase (HDAC1, HDAC4, and HDAC6) activities. Compound 6 c achieved the best cytotoxic activities on both HepG-2 and HuH-7 cell lines with IC50 values of 1.53 and 3.06 µM, respectively, and also it showed the most inhibitory activities on HDAC1, HDAC4, and HDAC6 with IC50 values of 1.76, 1.39, and 3.46 µM, respectively, compared to suberoylanilide hydroxamic acid (SAHA) as a reference drug (IC50 = 0.86, 0.97, and 0.93 µM, respectively). Furthermore, it achieved a more characteristic arrest in the growth of cell population of HepG-2 at both G0/G1 and S phases with 1.23-, and 1.18-fold, respectively, compared to that of the control, as determined by cell cycle analysis. Also, compound 6 c showed a marked elevation in the AnxV-FITC apoptotic HepG-2 cells percentage in both early and late phases increasing the total apoptosis percentage by 9.98-, and 10.81-fold, respectively, compared to the control. Furthermore, docking studies were carried out to identify the proposed binding mode of the synthesized compounds towards the prospective target (HDAC4). In silico ADMET and toxicity studies revealed that most of the synthesized compounds have accepted profiles of drug-likeness with low toxicity. Finally, an interesting SAR analysis was concluded to help the future design of more potent HDACIs in the future by medicinal chemists.

Predicting the in vivo developmental toxicity of benzo[a]pyrene (BaP) in rats by an in vitro-in silico approach

Developmental toxicity testing is an animal-intensive endpoints in toxicity testing and calls for animal-free alternatives. Previous studies showed the applicability of an in vitro–in silico approach for predicting developmental toxicity of a range of compounds, based on data from the mouse embryonic stem cell test (EST) combined with physiologically based kinetic (PBK) modelling facilitated reverse dosimetry. In the current study, the use of this approach for predicting developmental toxicity of polycyclic aromatic hydrocarbons (PAHs) was evaluated, using benzo[a]pyrene (BaP) as a model compound. A rat PBK model of BaP was developed to simulate the kinetics of its main metabolite 3-hydroxybenzo[a]pyrene (3-OHBaP), shown previously to be responsible for the developmental toxicity of BaP. Comparison to in vivo kinetic data showed that the model adequately predicted BaP and 3-OHBaP blood concentrations in the rat. Using this PBK model and reverse dosimetry, a concentration–response curve for 3-OHBaP obtained in the EST was translated into an in vivo dose–response curve for developmental toxicity of BaP in rats upon single or repeated dose exposure. The predicted half maximal effect doses (ED₅₀) amounted to 67 and 45 mg/kg bw being comparable to the ED₅₀ derived from the in vivo dose–response data reported for BaP in the literature, of 29 mg/kg bw. The present study provides a proof of principle of applying this in vitro–in silico approach for evaluating developmental toxicity of BaP and may provide a promising strategy for predicting the developmental toxicity of related PAHs, without the need for extensive animal testing.

New bis([1,2,4]triazolo)[4,3-a:3',4'-c]quinoxaline derivatives as VEGFR-2 inhibitors and apoptosis inducers: Design, synthesis, in silico studies, and anticancer evaluation

A new series of bis([1,2,4]triazolo)[4,3-a:3',4'-c]quinoxaline derivatives were designed and synthesized to have the main essential pharmacophoric features of VEGFR-2 inhibitors. VEGFR-2 inhibitory activities were assessed for the designed compounds. In addition, cytotoxic activity was evaluated for all derivatives against two human cancer cell lines namely, HepG-2 and MCF-7. The most cytotoxic compound 20 h was subjected to further biological investigations including cell cycle, apoptosis, caspase-3, caspase-9, BAX, and Bcl-2 analyses. Different in silico studies as docking, ADMET and toxicity were carried out. The results exhibited that compounds 20b, 20e, 20h and20mshowed promising VEGFR-2 inhibitory activities with IC₅₀values of 5.7, 6.7, 3.2, and 3.1 µM, respectively. Moreover, these promising members exhibited the highest antiproliferative activities against the two cell lines with IC₅₀values ranging from 3.3 to 14.2 µM, comparing to sorafenib (IC₅₀ = 2.17 and 3.43 µM against HepG2 and MCF-7, respectively). Additionally, compound 20h induced cell cycle arrest of HepG2 cells at G2/M phase. Also, such compound increased the progress of apoptosis by 3.5-fold compared to the control. As well, compound 20h showed a significant increase in the level of caspase-3 (2.07-fold), caspase-9 (1.72-fold), and BAX (1.83-fold), and a significant decrease in Bcl-2 level (1.92-fold). The in silico studies revealed that the synthesized compounds have binding pattern like that of sorafenib.

Discovery of thieno[2,3-d]pyrimidine-based derivatives as potent VEGFR-2 kinase inhibitors and anti-cancer agents

Vascular endothelial growth factor-2 (VEGFR-2) is considered one of the most important factors in tumor angiogenesis, and consequently a number of anticancer therapeutics have been developed to inhibit VEGFR-2 signaling. Accordingly, eighteen derivatives of thieno[2,3-d]pyrimidines having structural characteristics similar to VEGFR-2 inhibitors were designed and synthesized. Anticancer activities of the new derivatives were assessed against three human cancer cell lines (HCT-116, HepG2, and MCF-7) using MTT. Sorafenib was used as positive control. Compounds 17c-i, and 20b showed excellent anticancer activities against HCT-116 and HepG2 cell lines, while compounds 17i and 17g was found to be active against MCF-7 cell line. Compound 17f exhibited the highest cytotoxic activities against the examined cell lines, HCT-116 and HepG2, with IC₅₀ values of 2.80 ± 0.16 and 4.10 ± 0.45 µM, respectively. Aiming at exploring the mechanism of action of these compounds, the most active cytotoxic derivatives were in vitro tested for their VEGFR-2 inhibitory activity. Compound 17f showed high activity against VEGFR-2 with an IC₅₀ value of 0.23 ± 0.03 µM, that is equal to that of reference, sorafenib (IC₅₀ = 0.23 ± 0.04 µM). Molecular docking studies also were performed to investigate the possible binding interactions of the target compounds with the active sites of VEGFR-2. The synthesized compounds were analyzed for their ADMET and toxicity properties. Results showed that most of the compounds have low to very low BBB penetration levels and they have non-inhibitory effect against CYP2D6. All compounds were predicted to be non-toxic against developmental toxicity potential model except compounds 17b and 20b.

New quinoxaline derivatives as VEGFR-2 inhibitors with anticancer and apoptotic activity: Design, molecular modeling, and synthesis

New series of [1,2,4]triazolo[4,3-a]quinoxalin-4(5H)-one and [1,2,4]triazolo[4,3-a]quinoxaline derivatives have been designed, synthesized, and biologically assessed for their anti-proliferative activities against two selected tumor cell lines MCF-7 and HepG2. Comparing to sorafenib (IC₅₀ = 2.17 ± 0.13 and 3.51 ± 0.21 µM against MCF-7 and HepG2, respectively), compound 25d, 25e, 25i, and 27e exhibited the highest activities against the examined cell lines with IC₅₀ values extending from 4.1 ± 0.4 to 11.7 ± 1.1 µM. Furthermore, VEGFR-2 inhibitory activities were assessed for all the synthesized compounds as potential mechanisms for their anti-proliferative activities. Compounds 25d, 25e, 25i, and 27e displayed prominent inhibitory efficiency versus VEGFR-2 kinase with IC₅₀ value ranging from 3.4 ± 0.3 to 6.8 ± 0.5 nM. Fascinatingly, the results of VEGFR-2 inhibitory assays were matched with that of the cytotoxicity data, where the most potent anti-proliferative derivatives exhibited promising VEGFR-2 inhibitory activities. Further studies displayed the ability of compound 25d to induce apoptosis in HepG2 cells and can arrest the growth of such cells at the G2/M phase. Also, compound 25d produced a significant increase in the level of BAX/Bcl-2 ratio (3.8-fold), caspase- 3 (1.8-fold), and caspase-9 (1.9-fold) compared to the control cells. Molecular docking studies were carried out to investigate the possible binding interaction inside the active site of the VEGFR-2.

Incorporating renal excretion via the OCT2 transporter in physiologically based kinetic modelling to predict in vivo kinetics of mepiquat in rat

The present study aimed at incorporating active renal excretion via the organic cation transporter 2 (OCT2) into a generic rat physiologically based kinetic (PBK) model using an in vitro human renal proximal tubular epithelial cell line (SA7K) and mepiquat chloride (MQ) as the model compound. The Vmax (10.5 pmol/min/mg protein) and Km (20.6 μM) of OCT2 transport of MQ were determined by concentration-dependent uptake in SA7K cells using doxepin as inhibitor. PBK model predictions incorporating these values in the PBK model were 6.7-8.4-fold different from the reported in vivo data on the blood concentration of MQ in rat. Applying an overall scaling factor that also corrects for potential differences in OCT2 activity in the SA7K cells and in vivo kidney cortex and species differences resulted in adequate predictions for in vivo kinetics of MQ in rat (2.3-3.2-fold). The results indicate that using SA7K cells to define PBK parameters for active renal OCT2 mediated excretion with adequate scaling enables incorporation of renal excretion via the OCT2 transporter in PBK modelling to predict in vivo kinetics of mepiquat in rat. This study demonstrates a proof-of-principle on how to include active renal excretion into generic PBK models.

New quinoxaline-2(1H)-ones as potential VEGFR-2 inhibitors: design, synthesis, molecular docking, ADMET profile and anti-proliferative evaluations

Eleven new quinoxaline derivatives were designed and synthesized as modified VEGFR-2 inhibitors of our previous work.

Novel testing strategy for prediction of rat biliary excretion of intravenously administered estradiol-17β glucuronide

The aim of the present study was to develop a generic rat physiologically based kinetic (PBK) model that includes a novel testing strategy where active biliary excretion is incorporated using estradiol-17β glucuronide (E217βG) as the model substance. A major challenge was the definition of the scaling factor for the in vitro to in vivo conversion of the PBK-model parameter Vmax. In vitro values for the Vmax and Km for transport of E217βG were found in the literature in four different studies based on experiments with primary rat hepatocytes. The required scaling factor was defined based on fitting the PBK model-based predicted values to reported experimental data on E217βG blood levels and cumulative biliary E217βG excretion. This resulted in a scaling factor of 129 mg protein/g liver. With this scaling factor the PBK model predicted the in vivo data for blood and cumulative biliary E217βG levels with on average of less than 1.8-fold deviation. The study provides a proof of principle on how biliary excretion can be included in a generic PBK model using primary hepatocytes to define the kinetic parameters that describe the biliary excretion.

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.

Identification, isolation, structural characterization, in silico toxicity prediction and in vitro cytotoxicity assay of simeprevir acidic and oxidative degradation products

Simeprevir is a new direct-acting antiviral drug used for the treatment of chronic hepatitis C. In this work, a simple, fast and economical chromatographic method was developed for the determination of simeprevir in the presence of its acidic and oxidative degradation products. The stress studies performed herein showed that simeprevir degraded under acidic and oxidative conditions but was stable under thermal and alkaline conditions. Chromatographic separation was achieved on a reversed-phase Eclipse XDB C₁₈ column (4.6 × 150 mm, 5 μm). The mobile phase consisted of methanol-0.05 M ammonium acetate (pH 4) (90 : 10, v/v) and was used at a flow rate of 1 mL min⁻¹. The column effluent was monitored at 237 nm. The calibration curve was linear over the concentration range of 0.1–20 μg mL⁻¹. The relative standard deviations for the intra-day and inter-day precision were less than 2%, and good percentage recoveries that met the acceptance criteria of the International Conference on Harmonization (ICH) guidelines were obtained. The robustness was assessed using the Plackett–Burman design. The simeprevir degradation products were isolated by flash chromatography and confirmed by ¹H NMR and LC-MS/MS techniques. The fully validated chromatographic method can be applied as a stability-indicating method for simeprevir and for routine analysis during quality control. Additionally, in silico toxicity prediction of the degradation products demonstrated a hepatotoxicity alert for DP 1, DP 2, DP 4 and DP 5 and a carcinogenicity alert for DP 3. In view of safety aspects, an in vitro cytotoxicity assay was carried out for simeprevir degradation products. They were found to be non-toxic in vitro at the tested concentrations.

Simeprevir is a new direct-acting antiviral drug used for the treatment of chronic hepatitis C.

Predicting the Acute Liver Toxicity of Aflatoxin B1 in Rats and Humans by an In Vitro-In Silico Testing Strategy

SCOPE

High-level exposure to aflatoxin B1 (AFB1) is known to cause acute liver damage and fatality in animals and humans. The intakes actually causing this acute toxicity have so far been estimated based on AFB1 levels in contaminated foods or biomarkers in serum. The aim of the present study is to predict the doses causing acute liver toxicity of AFB1 in rats and humans by an in vitro-in silico testing strategy.

METHODS AND RESULTS

Physiologically based kinetic (PBK) models for AFB1 in rats and humans are developed. The models are used to translate in vitro concentration-response curves for cytotoxicity in primary rat and human hepatocytes to in vivo dose-response curves using reverse dosimetry. From these data, the dose levels at which toxicity would be expected are obtained and compared to toxic dose levels from available rat and human case studies on AFB1 toxicity. The results show that the in vitro-in silico testing strategy can predict dose levels causing acute toxicity of AFB1 in rats and human.

CONCLUSIONS

Quantitative in vitro in vivo extrapolation (QIVIVE) using PBK modeling-based reverse dosimetry can predict AFB1 doses that cause acute liver toxicity in rats and human.

Use of Physiologically Based Kinetic Modeling to Predict Rat Gut Microbial Metabolism of the Isoflavone Daidzein to S-Equol and Its Consequences for ERα Activation

SCOPE

To predict gut microbial metabolism of xenobiotics and the resulting plasma concentrations of metabolites formed, an in vitro-in silico-based testing strategy is developed using the isoflavone daidzein and its gut microbial metabolite S-equol as model compounds.

METHODS AND RESULTS

Anaerobic rat fecal incubations are optimized and performed to derive the apparent maximum velocities (Vmax ) and Michaelis-Menten constants (Km ) for gut microbial conversion of daidzein to dihydrodaidzein, S-equol, and O-desmethylangolensin, which are input as parameters for a physiologically based kinetic (PBK) model. The inclusion of gut microbiota in the PBK model allows prediction of S-equol concentrations and slightly reduced predicted maximal daidzein concentrations from 2.19 to 2.16 µm. The resulting predicted concentrations of daidzein and S-equol are comparable to in vivo concentrations reported.

CONCLUSION

The optimized in vitro approach to quantify kinetics for gut microbial conversions, and the newly developed PBK model for rats that includes gut microbial metabolism, provide a unique tool to predict the in vivo consequences of daidzein microbial metabolism for systemic exposure of the host to daidzein and its metabolite S-equol. The predictions reveal a dominant role for daidzein in ERα-mediated estrogenicity despite the higher estrogenic potency of its microbial metabolite S-equol.

Integrating in vitro testing and physiologically-based pharmacokinetic (PBPK) modelling for chemical liver toxicity assessment-A case study of troglitazone

In vitro to in vivo extrapolation (IVIVE) for next-generation risk assessment (NGRA) of chemicals requires computational modeling and faces unique challenges. Using mitochondria-related toxicity data of troglitazone (TGZ), a prototype drug known for liver toxicity, from HepaRG, HepG2, HC-04, and primary human hepatocytes, we explored inherent uncertainties in IVIVE, including cell models, cellular response endpoints, and dose metrics. A human population physiologically-based pharmacokinetic (PBPK) model for TGZ was developed to predict in vivo doses from in vitro point-of-departure (POD) concentrations. Compared to the 200-800 mg/d dose range of TGZ where liver injury was observed clinically, the predicted POD doses for the mean and top one percentile of the PBPK population were 28-372 and 15-178 mg/d respectively based on C_max dosimetry, and 185-2552 and 83-1010 mg/d respectively based on AUC. In conclusion, although with many uncertainties, integrating in vitro assays and PBPK modeling is promising in informing liver toxicity-inducing TGZ doses.

Combining In Vitro Data and Physiologically Based Kinetic Modeling Facilitates Reverse Dosimetry to Define In Vivo Dose-Response Curves for Bixin- and Crocetin-Induced Activation of PPARγ in Humans

SCOPE

It is investigated whether at realistic dietary intake bixin and crocetin could induce peroxisome proliferator-activated receptor γ (PPARγ)-mediated gene expression in humans using a combined in vitro-in silico approach.

METHODS AND RESULTS

Concentration-response curves obtained from in vitro PPARγ-reporter gene assays are converted to in vivo dose-response curves using physiologically based kinetic modeling-facilitated reverse dosimetry, from which the benchmark dose levels resulting in a 50% effect above background level (BMD50 ) are predicted and subsequently compared to dietary exposure levels. Bixin and crocetin activated PPARγ-mediated gene transcription in a concentration-dependent manner with similar potencies. Due to differences in kinetics, the predicted BMD50 values for in vivo PPARγ activation are about 30-fold different, amounting to 115 and 3505 mg kg bw-1 for crocetin and bixin, respectively. Human dietary and/or supplemental estimated daily intakes may reach these BMD50 values for crocetin but not for bixin, pointing at better possibilities for in vivo PPARγ activation by crocetin.

CONCLUSION

Based on a combined in vitro-in silico approach, it is estimated whether at realistic dietary intakes plasma concentrations of bixin and crocetin are likely to reach concentrations that activate PPARγ-mediated gene expression, without the need for a human intervention study.

Biotransformation and bioactivation reactions - 2018 literature highlights

In the past three decades, ADME sciences have become an integral component of the drug discovery and development process. At the same time, the field has continued to evolve, thus, requiring ADME scientists to be knowledgeable of and engage with diverse aspects of drug assessment: from pharmacology to toxicology, and from in silico modeling to in vitro models and finally in vivo models. Progress in this field requires deliberate exposure to different aspects of ADME; however, this task can seem daunting in the current age of mass information. We hope this review provides a focused and brief summary of a wide array of critical advances over the past year and explains the relevance of this research ( Table 1 ). We divided the articles into categories of (1) drug optimization, (2) metabolites and drug metabolizing enzymes, and (3) bioactivation. This annual review is the fourth of its kind (Baillie et al. 2016 ; Khojasteh et al. 2017 , 2018 ). We have followed the same format we used in previous years in terms of the selection of articles and the authoring of each section. This effort in itself also continues to evolve. I am pleased that Rietjens, Miller, and Mitra have again contributed to this annual review. We would like to welcome Namandjé N. Bumpus, James P. Driscoll, and Donglu Zhang as authors for this year's issue. We strive to maintain a balance of authors from academic and industry settings. We would be pleased to hear your opinions of our commentary, and we extend an invitation to anyone who would like to contribute to a future edition of this review. Cyrus Khojasteh, on behalf of the authors.

The in vivo developmental toxicity of diethylstilbestrol (DES) in rat evaluated by an alternative testing strategy

In the present study, we evaluated an alternative testing strategy to quantitatively predict the in vivo developmental toxicity of the synthetic hormone diethylstilbestrol (DES). To this end, a physiologically based kinetic (PBK) model was defined that was subsequently used to translate concentration-response data for the in vitro developmental toxicity of DES, obtained in the ES-D3 cell differentiation assay, into predicted in vivo dose-response data for developmental toxicity. The previous studies showed that the PBK model-facilitated reverse dosimetry approach is a useful approach to quantitatively predict the developmental toxicity of several developmental toxins. The results obtained in the present study show that the PBK model adequately predicted DES blood concentrations in rats. Further studies revealed that DES tested positive in the ES-D3 differentiation assay and that DES-induced inhibition of the ES-D3 cell differentiation could be counteracted by the estrogen receptor alpha (ERα) antagonist fulvestrant, indicating that the in vitro ES-D3 cell differentiation assay was able to mimic the role of ERα reported in the mode of action underlying the developmental toxicity of DES in vivo. In spite of this, combining these in vitro data with the PBK model did not adequately predict the in vivo developmental toxicity of DES in a quantitative way. It is concluded that although the EST qualifies DES as a developmental toxin and detects the role of ERα in this process, the ES-D3 cell differentiation assay of the EST apparently does not adequately capture the processes underlying DES-induced developmental toxicity in vivo.

Next generation physiologically based kinetic (NG-PBK) models in support of regulatory decision making

The fields of toxicology and chemical risk assessment seek to reduce, and eventually replace, the use of animals for the prediction of toxicity in humans. In this context, physiologically based kinetic (PBK) modelling based on in vitro and in silico kinetic data has the potential to a play significant role in reducing animal testing, by providing a methodology capable of incorporating in vitro human data to facilitate the development of in vitro to in vivo extrapolation of hazard information. In the present article, we discuss the challenges in: 1) applying PBK modelling to support regulatory decision making under the toxicology and risk-assessment paradigm shift towards animal replacement; 2) constructing PBK models without in vivo animal kinetic data, while relying solely on in vitro or in silico methods for model parameterization; and 3) assessing the validity and credibility of PBK models built largely using non-animal data. The strengths, uncertainties, and limitations of PBK models developed using in vitro or in silico data are discussed in an effort to establish a higher degree of confidence in the application of such models in a regulatory context. The article summarises the outcome of an expert workshop hosted by the European Commission Joint Research Centre (EC-JRC) - European Union Reference Laboratory for Alternatives to Animal Testing (EURL ECVAM), on "Physiologically-Based Kinetic modelling in risk assessment - reaching a whole new level in regulatory decision-making" held in Ispra, Italy, in November 2016, along with results from an international survey conducted in 2017 and recently reported activities occurring within the PBK modelling field. The discussions presented herein highlight the potential applications of next generation (NG)-PBK modelling, based on new data streams.

Use of physiologically based kinetic modelling-facilitated reverse dosimetry to convert in vitro cytotoxicity data to predicted in vivo liver toxicity of lasiocarpine and riddelliine in rat

Lasiocarpine and riddelliine are pyrrolizidine alkaloids (PAs) present in food and able to cause liver toxicity. The aim of this study was to investigate whether physiologically based kinetic (PBK) modelling-facilitated reverse dosimetry can adequately translate in vitro concentration-response curves for toxicity of lasiocarpine and riddelliine to in vivo liver toxicity data for the rat. To this purpose, PBK models were developed for lasiocarpine and riddelliine, and predicted blood concentrations were compared to available literature data to evaluate the models. Concentration-response curves obtained from in vitro cytotoxicity assays in primary rat hepatocytes were converted to in vivo dose-response curves from which points of departure (PODs) were derived and that were compared to available literature data on in vivo liver toxicity. The results showed that the predicted PODs fall well within the range of PODs derived from available in vivo toxicity data. To conclude, this study shows the proof-of-principle for a method to predict in vivo liver toxicity for PAs by an alternative testing strategy integrating in vitro cytotoxicity assays with in silico PBK modelling-facilitated reverse dosimetry. The approach may facilitate prediction of acute liver toxicity for the large number of PAs for which in vivo toxicity data are lacking.

Development and validation of a physiology-based model for the prediction of pharmacokinetics/toxicokinetics in rabbits

The environmental fates of pharmaceuticals and the effects of crop protection products on non-target species are subjects that are undergoing intense review. Since measuring the concentrations and effects of xenobiotics on all affected species under all conceivable scenarios is not feasible, standard laboratory animals such as rabbits are tested, and the observed adverse effects are translated to focal species for environmental risk assessments. In that respect, mathematical modelling is becoming increasingly important for evaluating the consequences of pesticides in untested scenarios. In particular, physiologically based pharmacokinetic/toxicokinetic (PBPK/TK) modelling is a well-established methodology used to predict tissue concentrations based on the absorption, distribution, metabolism and excretion of drugs and toxicants. In the present work, a rabbit PBPK/TK model is developed and evaluated with data available from the literature. The model predictions include scenarios of both intravenous (i.v.) and oral (p.o.) administration of small and large compounds. The presented rabbit PBPK/TK model predicts the pharmacokinetics (Cmax, AUC) of the tested compounds with an average 1.7-fold error. This result indicates a good predictive capacity of the model, which enables its use for risk assessment modelling and simulations.

Towards a generic physiologically based kinetic model to predict in vivo uterotrophic responses in rats by reverse dosimetry of in vitro estrogenicity data

Physiologically based kinetic (PBK) modelling-based reverse dosimetry is a promising tool for the prediction of in vivo developmental toxicity using in vitro concentration–response data. In the present study, the potential of this approach to predict the dose-dependent increase of uterus weight in rats upon exposure to estrogenic chemicals was assessed. In vitro concentration–response data of 17β-estradiol (E2) and bisphenol A (BPA) obtained in the MCF-7/BOS proliferation assay, the U2OS ER-CALUX assay and the yeast estrogen screen (YES) assay, were translated into in vivo dose–response data in rat, using a PBK model with a minimum number of in vitro and in silico determined parameter values. To evaluate the predictions made, benchmark dose (BMD) analysis was performed on the predicted dose–response data and the obtained BMDL₁₀ values were compared with BMDL₁₀ values derived from data on the effects of E2 and BPA in the uterotrophic assay reported in the literature. The results show that predicted dose–response data of E2 and BPA matched with the data from in vivo studies when predictions were made based on YES assay data. The YES assay-based predictions of the BMDL₁₀ values differed 3.9-fold (E2) and 4.7- to 13.4-fold (BPA) from the BMDL₁₀ values obtained from the in vivo data. The present study provides the proof-of-principle that PBK modelling-based reverse dosimetry of YES assay data using a minimum PBK model can predict dose-dependent in vivo uterus growth caused by estrogenic chemicals. In future studies, the approach should be extended to include other estrogens.

In vitro to in vivo extrapolation of effective dosimetry in developmental toxicity testing: Application of a generic PBK modelling approach

Incorporation of kinetics to quantitative in vitro to in vivo extrapolations (QIVIVE) is a key step for the realization of a non-animal testing paradigm, in the sphere of regulatory toxicology. The use of Physiologically-Based Kinetic (PBK) modelling for determining systemic doses of chemicals at the target site is accepted to be an indispensable element for such purposes. Nonetheless, PBK models are usually designed for a single or a group of compounds and are considered demanding, with respect to experimental data needed for model parameterization. Alternatively, we evaluate here the use of a more generic approach, i.e. the so-called IndusChemFate model, which is based on incorporated QSAR model parametrization. The model was used to simulate the in vivo kinetics of three diverse classes of developmental toxicants: triazoles, glycol ethers' alkoxyacetic acid metabolites and phthalate primary metabolites. The model required specific input per each class of compounds. These compounds were previously tested in three alternative assays: the whole-embryo culture (WEC), the zebrafish embryo test (ZET), and the mouse embryonic stem cell test (EST). Thereafter, the PBK-simulated blood levels at toxic in vivo doses were compared to the respective in vitro effective concentrations. Comparisons pertaining to relative potency and potency ranking with integration of kinetics were similar to previously obtained comparisons. Additionally, all three in vitro systems produced quite comparable results, and hence, a combination of alternative tests is still preferable for predicting the endpoint of developmental toxicity in vivo. This approach is put forward as biologically more plausible since plasma concentrations, rather than external administered doses, constitute the most direct in vivo dose metric.

In vitro-in silico-based analysis of the dose-dependent in vivo oestrogenicity of the soy phytoestrogen genistein in humans

BACKGROUND AND PURPOSE

The in vivo oestrogenicity of genistein and its glycoside genistin is still under debate. The present study aimed to develop a physiologically based kinetic (PBK) model that provides insight in dose-dependent plasma concentrations of genistein aglycone and its metabolites and enables prediction of in vivo oestrogenic effective dose levels of genistein and genistin in humans.

EXPERIMENTAL APPROACH

A PBK model for genistein and genistin in humans was developed based on in vitro metabolic parameters. The model obtained was used to translate in vitro oestrogenic concentration-response curves of genistein to in vivo oestrogenic dose-response curves for intake of genistein and genistin.

KEY RESULTS

The model predicted that genistein-7-O-glucuronide was the major circulating metabolite and that levels of the free aglycone were generally low [0.5-17% of total plasma genistein at oral doses from 0.01 to 50 mg (kg·bw)^-1 ]. The predicted in vivo benchmark dose for 5% response values for oestrogenicity varied between 0.06 and 4.39 mg kg^-1 genistein. For genistin, these values were 1.3-fold higher. These values are in line with reported human data and show that oestrogenic responses can be expected at an Asian dietary and a supplementary intake, while intake resulting from a Western diet may not be effective.

CONCLUSIONS AND IMPLICATIONS

The present study shows how plasma concentrations of genistein and its metabolites and oestrogenic dose levels of genistein in humans can be predicted by combining in vitro oestrogenicity with PBK model-based reverse dosimetry, eliminating the need for human intervention studies.

Physiologically based kinetic modeling of hesperidin metabolism and its use to predict in vivo effective doses in humans

SCOPE

To develop a physiologically based kinetic (PBK) model that describes the absorption, distribution, metabolism, and excretion of hesperidin in humans, enabling the translation of in vitro concentration-response curves to in vivo dose-response curves.

METHODS AND RESULTS

The PBK model for hesperidin in humans was developed based on in vitro metabolic parameters. Hesperidin was predicted to mainly occur in the systemic circulation as different monoglucuronides. The plasma concentrations of hesperidin aglycone (hesperetin) was predicted to be <0.02 mg/L at an oral dose of 50 mg/kg bw. The developed PBK model allowed conversion of in vitro concentration-response curves for different effects to in vivo dose-response curves. The BMD₀₅ (benchmark dose for 5% response) values for protein kinase A inhibition ranged between 135 and 529 mg/kg bw hesperidin, and for inhibition of endothelial cell migration and prostaglandin E₂ and nitric oxide production ranged between 2.19 and 44 mg/kg bw hesperidin. These values are in line with reported human data showing in vivo effects by hesperidin and show that these effects may occur at Western dietary and supplementary intake of hesperidin.

CONCLUSIONS

The developed PBK model adequately predicts absorption, distribution, metabolism, and excretion of hesperidin in humans and allows to evaluate the human in vivo situation without the need for human intervention studies.

Physiologically Based Pharmacokinetic Model of All-trans-Retinoic Acid with Application to Cancer Populations and Drug Interactions

All-trans retinoic acid (atRA) is a front-line treatment of acute promyelocytic leukemia (APL). Due to its activity in regulating the cell cycle, it has also been evaluated for the treatment of other cancers. However, the efficacy of atRA has been limited by atRA inducing its own metabolism during therapy, resulting in a decrease of atRA exposure during continuous dosing. Frequent relapse occurs in patients receiving atRA monotherapy. In an attempt to combat therapy resistance, inhibitors of atRA metabolism have been developed. Of these, ketoconazole and liarozole have shown some benefits, but their usage is limited by side effects and low potency toward the cytochrome P450 26A1 isoform (CYP26A1), the main atRA hydroxylase. We determined the pharmacokinetic basis of therapy resistance to atRA and tested whether the complex disposition kinetics of atRA could be predicted in healthy subjects and in cancer patients in the presence and absence of inhibitors of atRA metabolism using physiologically based pharmacokinetic (PBPK) modeling. A PBPK model of atRA disposition was developed and verified in healthy individuals and in cancer patients. The population-based PBPK model of atRA disposition incorporated saturable metabolic clearance of atRA, induction of CYP26A1 by atRA, and the absorption and distribution kinetics of atRA. It accurately predicted the changes in atRA exposure after continuous dosing and when coadministered with ketoconazole and liarozole. The developed model will be useful in interpretation of atRA disposition and efficacy, design of novel dosing strategies, and development of next-generation atRA metabolism inhibitors.

Use of Physiologically Based Kinetic Modeling-Based Reverse Dosimetry to Predict in Vivo Toxicity from in Vitro Data

The development of reliable nonanimal based testing strategies, such as in vitro bioassays, is the holy grail in current human safety testing of chemicals. However, the use of in vitro toxicity data in risk assessment is not straightforward. One of the main issues is that concentration-response curves from in vitro models need to be converted to in vivo dose-response curves. These dose-response curves are needed in toxicological risk assessment to obtain a point of departure to determine safe exposure levels for humans. Recent scientific developments enable this translation of in vitro concentration-response curves to in vivo dose-response curves using physiologically based kinetic (PBK) modeling-based reverse dosimetry. The present review provides an overview of the examples available in the literature on the prediction of in vivo toxicity using PBK modeling-based reverse dosimetry of in vitro toxicity data, showing that proofs-of-principle are available for toxicity end points ranging from developmental toxicity, nephrotoxicity, hepatotoxicity, and neurotoxicity to DNA adduct formation. This review also discusses the promises and pitfalls, and the future perspectives of the approach. Since proofs-of-principle available so far have been provided for the prediction of toxicity in experimental animals, future research should focus on the use of in vitro toxicity data obtained in human models to predict the human situation using human PBK models. This would facilitate human- instead of experimental animal-based approaches in risk assessment.

Integrating in vitro data and physiologically based kinetic (PBK) modelling to assess the in vivo potential developmental toxicity of a series of phenols

Toxicity outcomes derived in vitro do not always reflect in vivo toxicity values, which was previously observed for a series of phenols tested in the embryonic stem cell test (EST). Translation of in vitro data to the in vivo situation is therefore an important, but still limiting step for the use of in vitro toxicity outcomes in the safety assessment of chemicals. The aim of the present study was to translate in vitro embryotoxicity data for a series of phenols to in vivo developmental toxic potency values for the rat by physiologically based kinetic (PBK) modelling-based reverse dosimetry. To this purpose, PBK models were developed for each of the phenols. The models were parameterised with in vitro-derived values defining metabolism and transport of the compounds across the intestinal and placental barrier and with in silico predictions and data from the literature. Using PBK-based reverse dosimetry, in vitro concentration–response curves from the EST were translated into in vivo dose–response curves from which points of departure (PoDs) were derived. The predicted PoDs differed less than 3.6-fold from PoDs derived from in vivo toxicity data for the phenols available in the literature. Moreover, the in vitro PBK-based reverse dosimetry approach could overcome the large disparity that was observed previously between the in vitro and the in vivo relative potency of the series of phenols. In conclusion, this study shows another proof-of-principle that the in vitro PBK approach is a promising strategy for non-animal-based safety assessment of chemicals.

Predicting points of departure for risk assessment based on in vitro cytotoxicity data and physiologically based kinetic (PBK) modeling: The case of kidney toxicity induced by aristolochic acid I

Aristolochic acids are naturally occurring nephrotoxins. This study aims to investigate whether physiologically based kinetic (PBK) model-based reverse dosimetry could convert in vitro concentration-response curves of aristolochic acid I (AAI) to in vivo dose response-curves for nephrotoxicity in rat, mouse and human. To achieve this extrapolation, PBK models were developed for AAI in these different species. Subsequently, concentration-response curves obtained from in vitro cytotoxicity models were translated to in vivo dose-response curves using PBK model-based reverse dosimetry. From the predicted in vivo dose-response curves, points of departure (PODs) for risk assessment could be derived. The PBK models elucidated species differences in the kinetics of AAI with the overall catalytic efficiency for metabolic conversion of AAI to aristolochic acid Ia (AAIa) being 2-fold higher for rat and 64-fold higher for mouse than human. Results show that the predicted PODs generally fall within the range of PODs derived from the available in vivo studies. This study provides proof of principle for a new method to predict a POD for in vivo nephrotoxicity by integrating in vitro toxicity testing with in silico PBK model-based reverse dosimetry.