The structural and dynamical properties of the peroxisome proliferator-activated receptor γ (PPARγ) nuclear receptor have been broadly studied in its agonist state but little is known about the key features required for the receptor antagonistic activity. Here we report a series of molecular dynamics (MD) simulations in combination with free energy estimation of the recently discovered class of non-covalent PPARγ antagonists. Their binding modes and dynamical behavior are described in details. Two key interactions have been detected within the cavity between helices H3, H11 and the activation helix H12, as well as with H12. The strength of the ligand-amino acid residues interactions has been analyzed in relation to the specificity of the ligand dynamical and antagonistic features. According to our results, the PPARγ activation helix does not undergo dramatic conformational changes, as seen in other nuclear receptors, but rather perturbations that occur through a significant ligand-induced reshaping of the ligand-receptor and the receptor-coactivator binding pockets. The H12 residue Tyr473 and the charge clamp residue Glu471 play a central role for the receptor transformations. Our results also demonstrate that MD can be a helpful tool for the compound phenotype characterization (full agonists, partial agonists or antagonists) when insufficient experimental data are available.
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.
Comprehensive understanding of the precise mode of action/adverse outcome pathway (MoA/AOP) of chemicals becomes a key step towards superseding the current repeated dose toxicity testing methodology with new generation predictive toxicology tools. The description and characterization of the toxicological MoA leading to non-alcoholic fatty liver disease (NAFLD) are of specific interest, due to its increasing incidence in the modern society. Growing evidence stresses on the PPARγ ligand-dependent dysregulation as a key molecular initiating event (MIE) for this adverse effect. The aim of this work was to analyze and systematize the numerous scientific data about the steatogenic role of PPARγ. Over 300 papers were ranked according to preliminary defined criteria and used as reliable and significant sources of data about the PPARγ-dependent prosteatotic MoA. A detailed analysis was performed regarding proteins which PPARγ-mediated expression changes had been confirmed to be prosteatotic by most experimental evidence. Two probable toxicological MoAs from PPARγ ligand binding to NAFLD were described according to the Organisation for Economic Cooperation and Development (OECD) concepts: (i) PPARγ activation in hepatocytes and (ii) PPARγ inhibition in adipocytes. The possible events at different levels of biological organization starting from the MIE to the organ response and the connections between them were described in details.
The aim of this paper was to provide a proof of concept demonstrating that molecular modelling methodologies can be employed as a part of an integrated strategy to support toxicity prediction consistent with the mode of action/adverse outcome pathway (MoA/AOP) framework.To illustrate the role of molecular modelling in predictive toxicology, a case study was undertaken in which molecular modelling methodologies were employed to predict the activation of the peroxisome proliferator-activated nuclear receptor γ (PPARγ) as a potential molecular initiating event (MIE) for liver steatosis. A stepwise procedure combining different in silico approaches (virtual screening based on docking and pharmacophore filtering, and molecular field analysis) was developed to screen for PPARγ full agonists and to predict their transactivation activity (EC50). The performance metrics of the classification model to predict PPARγ full agonists were balanced accuracy = 81%, sensitivity = 85% and specificity = 76%. The 3D QSAR model developed to predict EC50 of PPARγ full agonists had the following statistical parameters: q
Graphical abstract Highlights
Silymarin, the active constituent of Silybum marianum (milk thistle), and its main component, silybin, are products with well-known hepatoprotective, cytoprotective, antioxidant, and chemopreventative properties. Despite substantial in vitro and in vivo investigations of these flavonolignans, their mechanisms of action and potential toxic effects are not fully defined. In this study we explored important ADME/Tox properties and biochemical interactions of selected flavonolignans using in silico methods. A quantitative structure-activity relationship (QSAR) model based on data from a parallel artificial membrane permeability assay (PAMPA) was used to estimate bioavailability after oral administration. Toxic effects and metabolic transformations were predicted using the knowledge-based expert systems Derek Nexus and Meteor Nexus (Lhasa Ltd). Potential estrogenic activity of the studied silybin congeners was outlined. To address further the stereospecificity of this effect the stereoisomeric forms of silybin were docked into the ligand-binding domain of the human estrogen receptor alpha (ERɑ) (MOE software, CCG). According to our results both stereoisomers can be accommodated into the ERɑ active site, but different poses and interactions were observed for silybin A and silybin B.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.