The comprehensive understanding of the precise mode of action and/or adverse outcome pathway (MoA/AOP) of chemicals has become a key step toward the development of a new generation of predictive toxicology tools. One of the challenges of this process is to test the feasibility of the molecular modelling approaches to explore key molecular initiating events (MIE) within the integrated strategy of MoA/AOP characterisation. The description of MoAs leading to toxicity and liver damage has been the focus of much interest. Growing evidence underlines liver PPARγ ligand-dependent activation as a key MIE in the elicitation of liver steatosis. Synthetic PPARγ full agonists are of special concern, since they may trigger a number of adverse effects not observed with partial agonists. In this study, molecular modelling was performed based on the PPARγ complexes with full agonists extracted from the Protein Data Bank. The receptor binding pocket was analysed, and the specific ligand-receptor interactions were identified for the most active ligands. A pharmacophore model was derived, and the most important pharmacophore features were outlined and characterised in relation to their specific role for PPARγ activation. The results are useful for the characterisation of the chemical space of PPARγ full agonists and could facilitate the development of preliminary filtering rules for the effective virtual ligand screening of compounds with PPARγ full agonistic activity.

Silymarin is the standardized extract from the fruits of Silybum marianum (L.) Gaertn., a well-known hepatoprotectant and antioxidant. Recently, bioactive compounds of silymarin, i.e., silybins and their 2,3-dehydro derivatives, have been shown to exert anticancer activities, yet with unclear mechanisms. This study combines in silico and in vitro methods to reveal the potential interactions of optically pure silybins and dehydrosilybins with novel protein targets. The shape and chemical similarity with approved drugs were evaluated in silico, and the potential for interaction with the Hedgehog pathway receptor Smoothened (SMO) and BRAF kinase was confirmed by molecular docking. In vitro studies on SMO and BRAF V600E kinase activity and in BRAF V600E A-375 human melanoma cell lines were further performed to examine their effects on these proteins and cancer cell lines and to corroborate computational predictions. Our in silico results direct to new potential targets of silymarin constituents as dual inhibitors of BRAF and SMO, two major targets in anticancer therapy. The experimental studies confirm that BRAF kinase and SMO may be involved in mechanisms of anticancer activities, demonstrating dose-dependent profiles, with dehydrosilybins showing stronger effects than silybins. The results of this work outline the dual SMO/BRAF effect of flavonolignans from Silybum marianum with potential clinical significance. Our approach can be applied to other natural products to reveal their potential targets and mechanism of action.