Abstract

Thyroid stimulating hormone receptor (TSHR) is an integral part of the hypothalamic-pituitary-thyroid axis. Notably, dysregulation in TSHR activation in humans can lead to adverse effects such as Graves disease, hypothyroidism and Hashimotos disease. Moreover, animal studies have shown that binding of endocrine disrupting chemicals (EDCs) with TSHR can lead to developmental toxicity. Several such chemicals have also been screened for their adverse physiological effects in human cell lines through various high-throughput assays under the ToxCast project. The vast resource of data generated through ToxCast has enabled the development of different toxicity predictors, but they can be limited in their predictive ability due to the heterogeneity in structure-activity relationships among chemicals. In an attempt to explore this heterogeneity, we systematically investigated structure-activity and structure-mechanism relationships among the TSHR binding chemicals from ToxCast. By employing structure-activity similarity (SAS) map, we identified 79 activity cliffs among 509 chemicals in the TSHR agonist dataset and 69 activity cliffs among 650 chemicals in the TSHR antagonist dataset. Further, by using the matched molecular pair (MMP) approach, we find that the resultant activity cliffs (MMP-cliffs) are a subset of activity cliffs identified via the SAS map approach. Moreover, by leveraging ToxCast mechanism of action (MOA) annotations for chemicals common to both TSHR agonist and antagonist datasets, we identified 3 chemical pairs as Strong MOA-cliffs and 19 chemical pairs as Weak MOA-cliffs. In conclusion, the insights from this systematic analysis of the structure-activity as well as the structure-mechanism relationships of TSHR binding chemicals are likely to inform ongoing efforts towards development of better predictive toxicity models for characterizing the chemical exposome.