The application of molecular modelling in the safety assessment of chemicals: A case study on ligand-dependent PPARγ dysregulation

Computational Strategies for Assessing Adverse Outcome Pathways: Hepatic Steatosis as a Case Study

The evolving landscape of chemical risk assessment is increasingly focused on developing tiered, mechanistically driven approaches that avoid the use of animal experiments. In this context, adverse outcome pathways have gained importance for evaluating various types of chemical-induced toxicity. Using hepatic steatosis as a case study, this review explores the use of diverse computational techniques, such as structure-activity relationship models, quantitative structure-activity relationship models, read-across methods, omics data analysis, and structure-based approaches to fill data gaps within adverse outcome pathway networks. Emphasizing the regulatory acceptance of each technique, we examine how these methodologies can be integrated to provide a comprehensive understanding of chemical toxicity. This review highlights the transformative impact of in silico techniques in toxicology, proposing guidelines for their application in evidence gathering for developing and filling data gaps in adverse outcome pathway networks. These guidelines can be applied to other cases, advancing the field of toxicological risk assessment.

Identification of molecular initiating events (MIE) using chemical database analysis and nuclear receptor activity assays for screening potential inhalation toxicants

An adverse outcome pathway (AOP) framework can facilitate the use of alternative assays in chemical regulations by providing scientific evidence. Previously, an AOP, peroxisome proliferative-activating receptor gamma (PPARγ) antagonism that leads to pulmonary fibrosis, was developed. Based on a literature search, PPARγ inactivation has been proposed as a molecular initiating event (MIE). In addition, a list of candidate chemicals that could be used in the experimental validation was proposed using toxicity database and deep learning models. In this study, the screening of environmental chemicals for MIE was conducted using in silico and in vitro tests to maximize the applicability of this AOP for screening inhalation toxicants. Initially, potential inhalation exposure chemicals that are active in three or more key events were selected, and in silico molecular docking was performed. Among the chemicals with low binding energy to PPARγ, nine chemicals were selected for validation of the AOP using in vitro PPARγ activity assay. As a result, rotenone, triorthocresyl phosphate, and castor oil were proposed as PPARγ antagonists and stressor chemicals of the AOP. Overall, the proposed tiered approach of the database-in silico-in vitro can help identify the regulatory applicability and assist in the development and experimental validation of AOP.

Multi-step structure-activity relationship screening efficiently predicts diverse PPARγ antagonists

In discovering the potential antagonist of peroxisome proliferator-activated receptor gamma (PPARγ), the structure-activity relationship (SAR) is a useful in silico method. However, it is difficult for conventional SAR approaches to predict the activities of antagonists owing to the large structural diversity of antagonistic compounds. This study provides evidence that multi-step SAR screening is applicable for predicting PPARγ antagonists by combining different complementary methodologies. We constructed three models: read-across-like SAR, docking-simulation-interpreting SAR, and deep-learning-based SAR. To provide user-customized prediction results, our multi-step SAR screening model combined the three SAR models in a stepwise manner, which subdivided them according to potential levels of the PPARγ antagonist. The read-across-like SAR, which considered specific antagonist scaffolds, revealed the highest positive predictive value (PPV). The docking-simulation-interpreting SAR, which considered the molecular surface features, revealed high statistics for the PPV and the true-positive rate (TPR). The deep-learning-based SAR showed the highest TPR at the last classification step. This multi-step SAR screening covered the antagonists of high reliability provided by a read-across-like SAR, as well as the antagonists of diverse scaffolds provided by docking-simulation-interpreting SAR and deep-learning-based SAR. Therefore, to predict PPARγ antagonists, multi-step SAR screening could be as a useful tool.

The GOLIATH Project: Towards an Internationally Harmonised Approach for Testing Metabolism Disrupting Compounds

The purpose of this project report is to introduce the European "GOLIATH" project, a new research project which addresses one of the most urgent regulatory needs in the testing of endocrine-disrupting chemicals (EDCs), namely the lack of methods for testing EDCs that disrupt metabolism and metabolic functions. These chemicals collectively referred to as "metabolism disrupting compounds" (MDCs) are natural and anthropogenic chemicals that can promote metabolic changes that can ultimately result in obesity, diabetes, and/or fatty liver in humans. This project report introduces the main approaches of the project and provides a focused review of the evidence of metabolic disruption for selected EDCs. GOLIATH will generate the world's first integrated approach to testing and assessment (IATA) specifically tailored to MDCs. GOLIATH will focus on the main cellular targets of metabolic disruption-hepatocytes, pancreatic endocrine cells, myocytes and adipocytes-and using an adverse outcome pathway (AOP) framework will provide key information on MDC-related mode of action by incorporating multi-omic analyses and translating results from in silico, in vitro, and in vivo models and assays to adverse metabolic health outcomes in humans at real-life exposures. Given the importance of international acceptance of the developed test methods for regulatory use, GOLIATH will link with ongoing initiatives of the Organisation for Economic Development (OECD) for test method (pre-)validation, IATA, and AOP development.

High-throughput screening and chemotype-enrichment analysis of ToxCast phase II chemicals evaluated for human sodium-iodide symporter (NIS) inhibition

In support of the Endocrine Disruptor Screening Program (EDSP), the U.S.EPA's Office of Research and Development (ORD) is developing high-throughput screening (HTS) approaches to identify chemicals that alter target sites in the thyroid hormone (TH) pathway. The sodium iodide symporter (NIS) is a transmembrane glycoprotein that mediates iodide uptake into the thyroid as the initial step of TH biosynthesis. Previously, we screened 293 ToxCast chemicals (ph1v2) using a HEK293T cell line expressing human NIS in parallel radioactive iodide uptake (RAIU) and cell viability assays to identify potential environmental NIS inhibitors. Here, we expanded NIS inhibitor screening for a set of 768 ToxCast Phase II (ph2) chemicals, and applied a novel computational toxicology approach based on the ToxPrint chemotype to identify chemical substructures associated with NIS inhibition. Following single-concentration screening (at 1 × 10-4 M with a 20% inhibition cutoff), 235 samples (228 chemicals) were further tested in multiple-concentration (1 × 10-9 - 1 × 10-4 M) format in both RAIU and cell viability assays. The 167 chemicals that exhibited significant RAIU inhibition were then prioritized using combined RAIU and cell viability responses that were normalized relative to the known NIS inhibitor sodium perchlorate. Some of the highest ranked chemicals, such as PFOS, tributyltin chloride, and triclocarban, have been previously reported to be thyroid disruptors. In addition, several novel chemicals were identified as potent NIS inhibitors. The present results were combined with the previous ph1v2 screening results to produce two sets of binary hit-calls for 1028 unique chemicals, consisting of 273 positives exhibiting significant RAIU inhibition, and 63 positives following application of a cell viability filter. A ToxPrint chemotype-enrichment analysis identified >20 distinct chemical substructural features, represented in >60% of the active chemicals, as significantly enriched in each NIS inhibition hit-call space. A shared set of 9 chemotypes enriched in both hit-call sets indicates stable chemotype signals (insensitive to cytotoxicity filters) that can help guide structure-activity relationship (SAR) investigations and inform future research.

Agonism activities of lyso-phosphatidylcholines (LPC) Ligands binding to peroxisome proliferator-activated receptor gamma (PPARγ)

PPARγ is an isoform of peroxisome proliferator-activated receptor (PPAR) belonging to a super family of nuclear receptors and is a primary target of the effective drug to treat the type II diabetes. The experiments found that Lyso-phosphatidylcholines (LPC) could bind to PPARγ, but the binding modes remain unknown. We used the Molecular Docking and Molecular Dynamic (MD) simulations to study the binding of four LPC ligands (LPC16:0, LPC18:0, LPC18:1-1 and LPC18:1-2) to PPARγ. The two-step MD simulations were employed to determine the final binding modes. The 20 ns MD simulations for four final LPC-PPARγ complexes were performed to analyze their structures, the binding key residues, and agonism activities. The results reveal that three LPC ligands (LPC16:0, LPC18:0 and LPC18:1-1) bind to Arm II and III regions of the Ligand Binding Domain (LBD) pocket, whereas they do not interact with Tyr473 of Helix 12 (H12). In contrast, LPC18:1-2 can form the hydrogen bonds with Tyr473 and bind into Arm I and II regions. Comparing with the paradigm systems of the full agonist (Rosiglitazone-PPARγ) and the partial agonist (MRL24-PPARγ), our results indicate that LPC16:0, LPC18:0 and LPC18:1-1 could be the potential partial agonists and LPC18:1-2 could be a full agonist. The in-depth analysis of the residue fluctuations and structure alignment confirm the present prediction of the LPC agonism activities.Communicated by Ramaswamy H. Sarma.

In Silico Molecular Docking and In Vivo Validation with Caenorhabditis elegans to Discover Molecular Initiating Events in Adverse Outcome Pathway Framework: Case Study on Endocrine-Disrupting Chemicals with Estrogen and Androgen Receptors

Molecular docking is used to analyze structural complexes of a target with its ligand for understanding the chemical and structural basis of target specificity. This method has the potential to be applied for discovering molecular initiating events (MIEs) in the Adverse Outcome Pathway framework. In this study, we aimed to develop in silico–in vivo combined approach as a tool for identifying potential MIEs. We used environmental chemicals from Tox21 database to identify potential endocrine-disrupting chemicals (EDCs) through molecular docking simulation, using estrogen receptor (ER), androgen receptor (AR) and their homology models in the nematode Caenorhabditis elegans (NHR-14 and NHR-69, respectively). In vivo validation was conducted on the selected EDCs with C. elegans reproductive toxicity assay using wildtype N2, nhr-14, and nhr-69 loss-of-function mutant strains. The chemicals showed high binding affinity to tested receptors and showed the high in vivo reproductive toxicity, and this was further confirmed using the mutant strains. The present study demonstrates that the binding affinity from the molecular docking potentially correlates with in vivo toxicity. These results prove that our in silico–in vivo combined approach has the potential to be applied for identifying MIEs. This study also suggests the potential of C. elegans as useful in the in vivo model for validating the in silico approach.

Targets and mechanisms of chemically induced aneuploidy. Part 1 of the report of the 2017 IWGT workgroup on assessing the risk of aneugens for carcinogenesis and hereditary diseases

An aneuploidy workgroup was established as part of the 7th International Workshops on Genotoxicity Testing. The workgroup conducted a review of the scientific literature on the biological mechanisms of aneuploidy in mammalian cells and methods used to detect chemical aneugens. In addition, the current regulatory framework was discussed, with the objective to arrive at consensus statements on the ramifications of exposure to chemical aneugens for human health risk assessment. As part of these efforts, the workgroup explored the use of adverse outcome pathways (AOPs) to document mechanisms of chemically induced aneuploidy in mammalian somatic cells. The group worked on two molecular initiating events (MIEs), tubulin binding and binding to the catalytic domain of aurora kinase B, which result in several adverse outcomes, including aneuploidy. The workgroup agreed that the AOP framework provides a useful approach to link evidence for MIEs with aneuploidy on a cellular level. The evidence linking chemically induced aneuploidy with carcinogenicity and hereditary disease was also reviewed and is presented in two companion papers. In addition, the group came to the consensus that the current regulatory test batteries, while not ideal, are sufficient for the identification of aneugens and human risk assessment. While it is obvious that there are many different MIEs that could lead to the induction of aneuploidy, the most commonly observed mechanisms involving chemical aneugens are related to tubulin binding and, to a lesser extent, inhibition of mitotic kinases. The comprehensive review presented here should help with the identification and risk management of aneugenic agents.

Molecular determinants of PPARγ partial agonism and related in silico/in vivo studies of natural saponins as potential type 2 diabetes modulators

The metabolic syndrome, which includes hypertension, type 2 diabetes (T2D) and obesity, has reached an epidemic-like scale. Saponins and sapogenins are considered as valuable natural products for ameliorating this pathology, possibly through the nuclear receptor PPARγ activation. The aims of this study were: to look for in vivo antidiabetic effects of a purified saponins' mixture (PSM) from Astragalus corniculatus Bieb; to reveal by in silico methods the molecular determinants of PPARγ partial agonism, and to investigate the potential PPARγ participation in the PSM effects. In the in vivo experiments spontaneously hypertensive rats (SHRs) with induced T2D were treated with PSM or pioglitazone as a referent PPARγ full agonist, and pathology-relevant biochemical markers were analysed. The results provided details on the PSM modulation of the glucose homeostasis and its potential mechanism. The in silico studies focused on analysis of the protein-ligand interactions in crystal structures of human PPARγ-partial agonist complexes, pharmacophore modelling and molecular docking. They outlined key pharmacophoric features, typical for the PPARγ partial agonists, which were used for pharmacophore-based docking of the main PSM sapogenin. The in silico studies, strongly suggest possible involvement of PPARγ-mediated mechanisms in the in vivo antidiabetic and antioxidant effects of PSM from A. corniculatus.

Ab initio chemical safety assessment: A workflow based on exposure considerations and non-animal methods

We describe and illustrate a workflow for chemical safety assessment that completely avoids animal testing. The workflow, which was developed within the SEURAT-1 initiative, is designed to be applicable to cosmetic ingredients as well as to other types of chemicals, e.g. active ingredients in plant protection products, biocides or pharmaceuticals. The aim of this work was to develop a workflow to assess chemical safety without relying on any animal testing, but instead constructing a hypothesis based on existing data, in silico modelling, biokinetic considerations and then by targeted non-animal testing. For illustrative purposes, we consider a hypothetical new ingredient x as a new component in a body lotion formulation. The workflow is divided into tiers in which points of departure are established through in vitro testing and in silico prediction, as the basis for estimating a safe external dose in a repeated use scenario. The workflow includes a series of possible exit (decision) points, with increasing levels of confidence, based on the sequential application of the Threshold of Toxicological (TTC) approach, read-across, followed by an "ab initio" assessment, in which chemical safety is determined entirely by new in vitro testing and in vitro to in vivo extrapolation by means of mathematical modelling. We believe that this workflow could be applied as a tool to inform targeted and toxicologically relevant in vitro testing, where necessary, and to gain confidence in safety decision making without the need for animal testing.