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

A novel method to derive a human safety limit for PFOA by gene expression profiling and modelling

Perfluorooctanoic acid (PFOA) is a persistent environmental contaminant that can accumulate in the human body due to its long half-life. This substance has been associated with liver, pancreatic, testicular and breast cancers, liver steatosis and endocrine disruption. PFOA is a member of a large group of substances also known as "forever chemicals" and the vast majority of substances of this group lack toxicological data that would enable their effective risk assessment in terms of human health hazards. This study aimed to derive a health-based guidance value for PFOA intake (ng/kg BW/day) from in vitro transcriptomics data. To this end, we developed an in silico workflow comprising five components: (i) sourcing in vitro hepatic transcriptomics concentration-response data; (ii) deriving molecular points of departure using BMDExpress3 and performing pathway analysis using gene set enrichment analysis (GSEA) to identify the most sensitive molecular pathways to PFOA exposure; (iii) estimating freely-dissolved PFOA concentrations in vitro using a mass balance model; (iv) estimating in vivo doses by reverse dosimetry using a PBK model for PFOA as part of a quantitative in vitro to in vivo extrapolation (QIVIVE) algorithm; and (v) calculating a tolerable daily intake (TDI) for PFOA. Fourteen percent of interrogated genes exhibited in vitro concentration-response relationships. GSEA pathway enrichment analysis revealed that "fatty acid metabolism" was the most sensitive pathway to PFOA exposure. In vitro free PFOA concentrations were calculated to be 2.9% of the nominal applied concentrations, and these free concentrations were input into the QIVIVE workflow. Exposure doses for a virtual population of 3,000 individuals were estimated, from which a TDI of 0.15 ng/kg BW/day for PFOA was calculated using the benchmark dose modelling software, PROAST. This TDI is comparable to previously published values of 1.16, 0.69, and 0.86 ng/kg BW/day by the European Food Safety Authority. In conclusion, this study demonstrates the combined utility of an "omics"-derived molecular point of departure and in silico QIVIVE workflow for setting health-based guidance values in anticipation of the acceptance of in vitro concentration-response molecular measurements in chemical risk assessment.

Rapid equilibrium dialysis, ultrafiltration or ultracentrifugation? Evaluation of methods to quantify the unbound fraction of substances in plasma

In pharmacokinetics plasma protein binding (PPB) is a well-established parameter impacting drug disposition. The unbound fraction (f_u) is arguably regarded the effective concentration at the target site. Pharmacology and toxicology, increasingly use in vitro models. The translation of in vitro concentrations to in vivo doses can be supported by toxicokinetic modelling, e.g. physiologically based toxicokinetic models (PBTK). PPB of a test substance is an input parameter for PBTK. We compared three methods to quantify f_u: rapid equilibrium dialysis (RED), ultrafiltration (UF) and ultracentrifugation (UC) using twelve substances covering a wide range of Log P_ow (-0.1 to 6.8) and molecular weights (151 and 531 g/mol): Acetaminophen, Bisphenol A, Caffeine, Colchicine, Fenarimol, Flutamide, Genistein, Ketoconazole, α-Methyltestosterone, Tamoxifen, Trenbolone and Warfarin. After RED and UF separation, three polar substances (Log P_ow < 2) were largely unbound (f_u > 70%), while more lipophilic substances were largely bound (f_u < 33%). Compared to RED or UF, UC resulted in a generally higher f_u of lipophilic substances. f_u obtained after RED and UF were more consistent with published data. For half of the substances, UC resulted in f_u higher than the reference data. UF, RED and both UF and UC, resulted in lower f_u of Flutamide, Ketoconazole and Colchicine, respectively. For f_u quantifications, the separation method should be selected according to the test substance's properties. Based on our data, RED is suitable for a broader range of substances while UC and UF are suitable for polar substances.

Derivation of a Human In Vivo Benchmark Dose for Bisphenol A from ToxCast In Vitro Concentration Response Data Using a Computational Workflow for Probabilistic Quantitative In Vitro to In Vivo Extrapolation

A computational workflow which integrates physiologically based kinetic (PBK) modelling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC), Markov Chain Monte Carlo (MCMC) simulation and the Virtual Cell Based Assay (VCBA) for the estimation of the active, free in vitro concentration of chemical in the reaction medium was developed to facilitate quantitative in vitro to in vivo extrapolation (QIVIVE). The workflow was designed to estimate parameter and model uncertainty within a computationally efficient framework. The workflow was tested using a human PBK model for bisphenol A (BPA) and high throughput screening (HTS) in vitro concentration-response data, for estrogen and pregnane X receptor activation determined in human liver and kidney cell lines, from the ToxCast/Tox21 database. In vivo benchmark dose 10% lower confidence limits (BMDL₁₀) for oral uptake of BPA (ng/kg BW/day) were calculated from the in vivo dose-responses and compared to the human equivalent dose (HED) BMDL₁₀ for relative kidney weight change in the mouse derived by European Food Safety Authority (EFSA). Three from four in vivo BMDL₁₀ values calculated in this study were similar to the EFSA values whereas the fourth was much smaller. The derivation of an uncertainty factor (UF) to accommodate the uncertainties associated with measurements using human cell lines in vitro, extrapolated to in vivo, could be useful for the derivation of Health Based Guidance Values (HBGV).

Genistein: Dual Role in Women's Health

Advanced research in recent years has revealed the important role of nutrients in the protection of women's health and in the prevention of women's diseases. Genistein is a phytoestrogen that belongs to a class of compounds known as isoflavones, which structurally resemble endogenous estrogen. Genistein is most often consumed by humans via soybeans or soya products and is, as an auxiliary medicinal, used to treat women's diseases. In this review, we focused on analyzing the geographic distribution of soybean and soya product consumption, global serum concentrations of genistein, and its metabolism and bioactivity. We also explored genistein's dual effects in women's health through gathering, evaluating, and summarizing evidence from current in vivo and in vitro studies, clinical observations, and epidemiological surveys. The dose-dependent effects of genistein, especially when considering its metabolites and factors that vary by individuals, indicate that consumption of genistein may contribute to beneficial effects in women's health and disease prevention and treatment. However, consumption and exposure levels are nuanced because adverse effects have been observed at lower concentrations in in vitro models. Therefore, this points to the duplicity of genistein as a possible therapeutic agent in some instances and as an endocrine disruptor in others.

Experimental and molecular docking studies of estrogen-like and anti-osteoporosis activity of compounds in Fructus Psoraleae

ETHNOPHARMACOLOGICAL RELEVANCE

Fructus Psoraleae (FP), dry mature fruits of Cullen corylifolium (L.) Medik., has been used clinically to treat kidney yang deficiency-induced impotence, asthma and cold pain in waist and knee caused by kidney deficiency. A study of the source of the significant kidney-enhancing effect of FP revealed that it may be due to its strong estrogen-like activity.

AIM OF THE STUDY

This study aimed to investigate the estrogen-like activity of the FP extract and 13 bioactive compounds in it, as well as the mechanisms underlying their estrogen-like and anti-osteoporosis activities.

MATERIALS AND METHODS

The estrogen-like activities of the 75% ethanol-only FP extract, and 75% ethanol plus petroleum ether, ethyl acetate, n-butanol or water FP extracts were each measured using Cell Counting Kit-8 (CCK-8) and luciferase reporter gene assays. The compounds were identified by high-performance liquid chromatography analysis. The activation of estrogen receptor signaling by the compounds was compared with that by estradiol (E₂) using the molecular docking software MOE-Dock 2008.10. The activation of the ER-Wnt-β-catenin signaling pathway was investigated using an alkaline phosphatase (ALP) assay, qPCR analysis and Western blot analysis.

RESULTS

The results revealed that the 75% ethanol plus ethyl acetate extract showed the highest estrogen-like activity among the four 75% ethanol extract fractions (further extracted with petroleum ether, ethyl acetate, n-butanol or water). Some compounds in FP showed strong estrogenic effect and anti-osteoporosis activity, and activated the Wnt-β-catenin pathway. The isoflavone compound was the most active.

CONCLUSIONS

This study demonstrated that FP has a strong estrogen-like activity and some of its component compounds have anti-osteoporosis activity by activating the ER-Wnt-β-catenin signaling pathway. Our detections provide a new insight into the mechanisms underlying the estrogen-like and anti-osteoporosis activities of FP, as well as a better understanding of structure effects.

Genistein: An Integrative Overview of Its Mode of Action, Pharmacological Properties, and Health Benefits

Genistein is an isoflavone first isolated from the brooming plant Dyer's Genista tinctoria L. and is widely distributed in the Fabaceae family. As an isoflavone, mammalian genistein exerts estrogen-like functions. Several biological effects of genistein have been reported in preclinical studies, such as the antioxidant, anti-inflammatory, antibacterial, and antiviral activities, the effects of angiogenesis and estrogen, and the pharmacological activities on diabetes and lipid metabolism. The purpose of this review is to provide up-to-date evidence of preclinical pharmacological activities with mechanisms of action, bioavailability, and clinical evidence of genistein. The literature was researched using the most important keyword "genistein" from the PubMed, Science, and Google Scholar databases, and the taxonomy was validated using The Plant List. Data were also collected from specialized books and other online resources. The main positive effects of genistein refer to the protection against cardiovascular diseases and to the decrease of the incidence of some types of cancer, especially breast cancer. Although the mechanism of protection against cancer involves several aspects of genistein metabolism, the researchers attribute this effect to the similarity between the structure of soy genistein and that of estrogen. This structural similarity allows genistein to displace estrogen from cellular receptors, thus blocking their hormonal activity. The pharmacological activities resulting from the experimental studies of this review support the traditional uses of genistein, but in the future, further investigations are needed on the efficacy, safety, and use of nanotechnologies to increase bioavailability and therapeutic efficacy.

Metabolic tracking of isoflavones in soybean products and biosamples from healthy adults after fermented soybean consumption

Fermentation may enhance the nutritional properties of foods by increasing metabolite bioactivity or bioavailability. This study explored the effect of fermentation on isoflavone bioavailability and metabolism. Isoflavone metabolites were tracked in foods and biospecimens of healthy adults after fermented soybean (FS) or non-fermented soybean (NFS) consumption in a randomized, controlled, crossover intervention study. The change in soybean isoflavones caused by fermentation resulted in faster absorption and higher bioavailability after consumption of FS. Although the urinary level of total isoflavone metabolites was similar after the consumption of the two diets, urinary genistein 7-O-sulfate was derived as a discriminant metabolite for the FS diet by partial least squares discriminant analysis. This study suggests that an isoflavone conjugate profile might be a more appropriate marker than total isoflavone levels for discriminating between the consumption of FS and NFS diets.

Potential of ToxCast Data in the Safety Assessment of Food Chemicals

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

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.

A Computational Workflow for Probabilistic Quantitative in Vitro to in Vivo Extrapolation

A computational workflow was developed to facilitate the process of quantitative in vitro to in vivo extrapolation (QIVIVE), specifically the translation of in vitro concentration-response to in vivo dose-response relationships and subsequent derivation of a benchmark dose value (BMD). The workflow integrates physiologically based pharmacokinetic (PBPK) modeling; global sensitivity analysis (GSA), Approximate Bayesian Computation (ABC) and Markov Chain Monte Carlo (MCMC) simulation. For a given set of in vitro concentration and response data the algorithm returns the posterior distribution of the corresponding in vivo, population-based dose-response values, for a given route of exposure. The novel aspect of the workflow is a rigorous statistical framework for accommodating uncertainty in both the parameters of the PBPK model (both parameter uncertainty and population variability) and in the structure of the PBPK model itself recognizing that the model is an approximation to reality. Both these sources of uncertainty propagate through the workflow and are quantified within the posterior distribution of in vivo dose for a fixed representative in vitro concentration. To demonstrate this process and for comparative purposes a similar exercise to previously published work describing the kinetics of ethylene glycol monoethyl ether (EGME) and its embryotoxic metabolite methoxyacetic acid (MAA) in rats was undertaken. The computational algorithm can be used to extrapolate from in vitro data to any organism, including human. Ultimately, this process will be incorporated into a user-friendly, freely available modeling platform, currently under development, that will simplify the process of QIVIVE.

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.

In vitro bioassays to evaluate beneficial and adverse health effects of botanicals: promises and pitfalls

This review provides an update on the promises and pitfalls when using in vitro bioassays to evaluate beneficial and adverse health effects of botanicals and botanical preparations. Important issues addressed in the paper are: (i) the type of assays and biological effects available; (ii) false-positives, false-negatives and confounding factors; (iii) matrix and combination effects; (iv) extrapolation of in vitro data to the in vivo situation; (v) when (not) to use bioassays; and (vi) identification of active constituents. It is concluded that in vitro bioassays provide models to detect beneficial as well as adverse activities, but that linking these observations to individual ingredients and extrapolations to the in vivo situation is more complicated than generally anticipated.