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 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.