In Silico Prediction of Endocrine Disrupting Chemicals Using Single-Label and Multilabel Models

Recognizing high-priority disinfection byproducts based on experimental and predicted endocrine disrupting data: Virtual screening and in vitro study

So far, about 130 disinfection by-products (DBPs) and several DBPs-groups have had their potential endocrine-disrupting effects tested on some endocrine endpoints. However, it is still not clear which specific DBPs, DBPs-groups/subgroups may be the most toxic substances or groups/subgroups for any given endocrine endpoint. In this study, we attempt to address this issue. First, a list of relevant DBPs was updated, and 1187 DBPs belonging to 4 main-groups (aliphatic, aromatic, alicyclic, heterocyclic) and 84 subgroups were described. Then, the high-priority endocrine endpoints, DBPs-groups/subgroups, and specific DBPs were determined from 18 endpoints, 4 main-groups, 84 subgroups, and 1187 specific DBPs by a virtual-screening method. The results demonstrate that most of DBPs could not disturb the endocrine endpoints in question because the proportion of active compounds associated with the endocrine endpoints ranged from 0 (human thyroid receptor beta) to 32% (human transthyretin (hTTR)). All the endpoints with a proportion of active compounds greater than 10% belonged to the thyroid system, highlighting that the potential disrupting effects of DBPs on the thyroid system should be given more attention. The aromatic and alicyclic DBPs may have higher priority than that of aliphatic and heterocyclic DBPs by considering the activity rate and potential for disrupting effects. There were 2 (halophenols and estrogen DBPs), 12, and 24 subgroups that belonged to high, moderate, and low priority classes, respectively. For individual DBPs, there were 23 (2%), 193 (16%), and 971 (82%) DBPs belonging to the high, moderate, and low priority groups, respectively. Lastly, the hTTR binding affinity of 4 DBPs was determined by an in vitro assay and all the tested DBPs exhibited dose-dependent binding potency with hTTR, which was consistent with the predicted result. Thus, more efforts should be performed to reveal the potential endocrine disruption of those high research-priority main-groups, subgroups, and individual DBPs.

On the additive artificial intelligence-based discovery of nanoparticle neurodegenerative disease drug delivery systems

Neurodegenerative diseases are characterized by slowly progressing neuronal cell death. Conventional drug treatment strategies often fail because of poor solubility, low bioavailability, and the inability of the drugs to effectively cross the blood-brain barrier. Therefore, the development of new neurodegenerative disease drugs (NDDs) requires immediate attention. Nanoparticle (NP) systems are of increasing interest for transporting NDDs to the central nervous system. However, discovering effective nanoparticle neuronal disease drug delivery systems (N2D3Ss) is challenging because of the vast number of combinations of NP and NDD compounds, as well as the various assays involved. Artificial intelligence/machine learning (AI/ML) algorithms have the potential to accelerate this process by predicting the most promising NDD and NP candidates for assaying. Nevertheless, the relatively limited amount of reported data on N2D3S activity compared to assayed NDDs makes AI/ML analysis challenging. In this work, the IFPTML technique, which combines information fusion (IF), perturbation theory (PT), and machine learning (ML), was employed to address this challenge. Initially, we conducted the fusion into a unified dataset comprising 4403 NDD assays from ChEMBL and 260 NP cytotoxicity assays from journal articles. Through a resampling process, three new working datasets were generated, each containing 500,000 cases. We utilized linear discriminant analysis (LDA) along with artificial neural network (ANN) algorithms, such as multilayer perceptron (MLP) and deep learning networks (DLN), to construct linear and non-linear IFPTML models. The IFPTML-LDA models exhibited sensitivity (Sn) and specificity (Sp) values in the range of 70% to 73% (>375,000 training cases) and 70% to 80% (>125,000 validation cases), respectively. In contrast, the IFPTML-MLP and IFPTML-DLN achieved Sn and Sp values in the range of 85% to 86% for both training and validation series. Additionally, IFPTML-ANN models showed an area under the receiver operating curve (AUROC) of approximately 0.93 to 0.95. These results indicate that the IFPTML models could serve as valuable tools in the design of drug delivery systems for neurosciences.

deepFPlearn +: enhancing toxicity prediction across the chemical universe using graph neural networks

SUMMARY

Sophisticated approaches for the in silico prediction of toxicity are required to support the risk assessment of chemicals. The number of chemicals on the global chemical market and the speed of chemical innovation stand in massive contrast to the capacity for regularizing chemical use. We recently proved our ready-to-use application deepFPlearn as a suitable approach for this task. Here, we present its extension deepFPlearn+ incorporating (i) a graph neural network to feed our AI with a more sophisticated molecular structure representation and (ii) alternative train-test splitting strategies that involve scaffold structures and the molecular weights of chemicals. We show that the GNNs outperform the previous model substantially and that our models can generalize on unseen data even with a more robust and challenging test set. Therefore, we highly recommend the application of deepFPlearn+ on the chemical inventory to prioritize chemicals for experimental testing or any chemical subset of interest in monitoring studies.

AVAILABILITY AND IMPLEMENTATION

The software is compatible with python 3.6 or higher, and the source code can be found on our GitHub repository: https://github.com/yigbt/deepFPlearn. The data underlying this article are available in Zenodo, and can be accessed with the link below: https://zenodo.org/record/8146252. Detailed installation guides via Docker, Singularity, and Conda are provided within the repository for operability across all operating systems.

EDC-Predictor: A Novel Strategy for Prediction of Endocrine-Disrupting Chemicals by Integrating Pharmacological and Toxicological Profiles

Identification of endocrine-disrupting chemicals (EDCs) is crucial in the reduction of human health risks. However, it is hard to do so because of the complex mechanisms of the EDCs. In this study, we propose a novel strategy named EDC-Predictor to integrate pharmacological and toxicological profiles for the prediction of EDCs. Different from conventional methods that only focus on a few nuclear receptors (NRs), EDC-Predictor considers more targets. It uses computational target profiles from network-based and machine learning-based methods to characterize compounds, including both EDCs and non-EDCs. The best model constructed by these target profiles outperformed those models by molecular fingerprints. In a case study to predict NR-related EDCs, EDC-Predictor showed a wider applicability domain and higher accuracy than four previous tools. Another case study further demonstrated that EDC-Predictor could predict EDCs targeting other proteins rather than NRs. Finally, a free web server was developed to make EDC prediction easier (http://lmmd.ecust.edu.cn/edcpred/). In summary, EDC-Predictor would be a powerful tool in EDC prediction and drug safety assessment.

Modeling and insights into the structural characteristics of endocrine-disrupting chemicals

Endocrine-disrupting chemicals (EDCs) can cause serious harm to human health and the environment; therefore, it is important to rapidly and correctly identify EDCs. Different computational models have been proposed for the prediction of EDCs over the past few decades, but the reported models are not always easily available, and few studies have investigated the structural characteristics of EDCs. In the present study, we have developed a series of artificial intelligence models targeting EDC receptors: the androgen receptor (AR); estrogen receptor (ER); and pregnane X receptor (PXR). The consensus models achieved good predictive results for validation sets with balanced accuracy values of 87.37%, 90.13%, and 79.21% for AR, ER, and PXR binding assays, respectively. Analysis of the physical-chemical properties suggested that several chemical properties were significantly (p < 0.05) different between EDCs and non-EDCs. We also identified structural alerts that can indicate an EDC, which were integrated into the web server SApredictor. These models and structural characteristics can provide useful tools and information in the discrimination and mechanistic understanding of EDCs in drug discovery and environmental risk assessment.

Prediction of inotropic effect based on calcium transients in human iPSC-derived cardiomyocytes and machine learning

Functional changes to cardiomyocytes are undesirable during drug discovery and identifying the inotropic effects of compounds is hence necessary to decrease the risk of cardiovascular adverse effects in the clinic. Recently, approaches leveraging calcium transients in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) have been developed to detect contractility changes, induced by a variety of mechanisms early during drug discovery projects. Although these approaches have been able to provide some predictive ability, we hypothesised that using additional waveform parameters could offer improved insights, as well as predictivity. In this study, we derived 25 parameters from each calcium transient waveform and developed a modified Random Forest method to predict the inotropic effects of the compounds. In total annotated data for 48 compounds were available for modelling, out of which 31 were inotropes. The results show that the Random Forest model with a modified purity criterion performed slightly better than an unmodified algorithm in terms of the Area Under the Curve, giving values of 0.84 vs 0.81 in a cross-validation, and outperformed the ToxCast Pipeline model, for which the highest value was 0.76 when using the best-performing parameter, PW10. Our study hence demonstrates that more advanced parameters derived from waveforms, in combination with additional machine learning methods, provide improved predictivity of cardiovascular risk associated with inotropic effects.

Bioremediation of Endocrine Disrupting Chemicals- Advancements and Challenges

Endocrine Disrupting Chemicals (EDCs), major group of recalcitrant compounds, poses a serious threat to the health and future of millions of human beings, and other flora and fauna for years to come. A close analysis of various xenobiotics undermines the fact that EDC is structurally diverse chemical compounds generated as a part of anthropogenic advancements as well as part of their degradation. Regardless of such structural diversity, EDC is common in their ultimate drastic effect of impeding the proper functioning of the endocrinal system, basic physiologic systems, resulting in deregulated growth, malformations, and cancerous outcomes in animals as well as humans. The current review outlines an overview of various EDCs, their toxic effects on the ecosystem and its inhabitants. Conventional remediation methods such as physico-chemical methods and enzymatic approaches have been put into action as some form of mitigation measures. However, the last decade has seen the hunt for newer technologies and methodologies at an accelerated pace. Genetically engineered microbial degradation, gene editing strategies, metabolic and protein engineering, and in-silico predictive approaches - modern day's additions to our armamentarium in combating the EDCs are addressed. These additions have greater acceptance socially with lesser dissonance owing to reduced toxic by-products, lower health trepidations, better degradation, and ultimately the prevention of bioaccumulation. The positive impact of such new approaches on controlling the menace of EDCs has been outlaid. This review will shed light on sources of EDCs, their impact, significance, and the different remediation and bioremediation approaches, with a special emphasis on the recent trends and perspectives in using sustainable approaches for bioremediation of EDCs. Strict regulations to prevent the release of estrogenic chemicals to the ecosystem, adoption of combinatorial methods to remove EDC and prevalent use of bioremediation techniques should be followed in all future endeavors to combat EDC pollution. Moreover, the proper development, growth and functioning of future living forms relies on their non-exposure to EDCs, thus remediation of such chemicals present even in nano-concentrations should be addressed gravely.

Profiling prediction of nuclear receptor modulators with multi-task deep learning methods: toward the virtual screening

Nuclear receptors (NRs) are ligand-activated transcription factors, which constitute one of the most important targets for drug discovery. Current computational strategies mainly focus on a single target, and the transfer of learned knowledge among NRs was not considered yet. Herein we proposed a novel computational framework named NR-Profiler for prediction of potential NR modulators with high affinity and specificity. First, we built a comprehensive NR data set including 42 684 interactions to connect 42 NRs and 31 033 compounds. Then, we used multi-task deep neural network and multi-task graph convolutional neural network architectures to construct multi-task multi-classification models. To improve the predictive capability and robustness, we built a consensus model with an area under the receiver operating characteristic curve (AUC) = 0.883. Compared with conventional machine learning and structure-based approaches, the consensus model showed better performance in external validation. Using this consensus model, we demonstrated the practical value of NR-Profiler in virtual screening for NRs. In addition, we designed a selectivity score to quantitatively measure the specificity of NR modulators. Finally, we developed a freely available standalone software for users to make profiling predictions for their compounds of interest. In summary, our NR-Profiler provides a useful tool for NR-profiling prediction and is expected to facilitate NR-based drug discovery.

AI for predicting chemical-effect associations at the chemical universe level-deepFPlearn

Many chemicals are present in our environment, and all living species are exposed to them. However, numerous chemicals pose risks, such as developing severe diseases, if they occur at the wrong time in the wrong place. For the majority of the chemicals, these risks are not known. Chemical risk assessment and subsequent regulation of use require efficient and systematic strategies. Lab-based methods—even if high throughput—are too slow to keep up with the pace of chemical innovation. Existing computational approaches are designed for specific chemical classes or sub-problems but not usable on a large scale. Further, the application range of these approaches is limited by the low amount of available labeled training data. We present the ready-to-use and stand-alone program deepFPlearn that predicts the association between chemical structures and effects on the gene/pathway level using a combined deep learning approach. deepFPlearn uses a deep autoencoder for feature reduction before training a deep feed-forward neural network to predict the target association. We received good prediction qualities and showed that our feature compression preserves relevant chemical structural information. Using a vast chemical inventory (unlabeled data) as input for the autoencoder did not reduce our prediction quality but allowed capturing a much more comprehensive range of chemical structures. We predict meaningful—experimentally verified—associations of chemicals and effects on unseen data. deepFPlearn classifies hundreds of thousands of chemicals in seconds. We provide deepFPlearn as an open-source and flexible tool that can be easily retrained and customized to different application settings at https://github.com/yigbt/deepFPlearn.

Predicting the binding of small molecules to nuclear receptors using machine learning

Nuclear receptors (NRs) are important biological targets of endocrine-disrupting chemicals (EDCs). Identifying chemicals that can act as EDCs and modulate the function of NRs is difficult because of the time and cost of in vitro and in vivo screening to determine the potential hazards of the 100 000s of chemicals that humans are exposed to. Hence, there is a need for computational approaches to prioritize chemicals for biological testing. Machine learning (ML) techniques are alternative methods that can quickly screen millions of chemicals and identify those that may be an EDC. Computational models of chemical binding to multiple NRs have begun to emerge. Recently, a Nuclear Receptor Activity (NuRA) dataset, describing experimentally derived small-molecule activity against various NRs has been created. We have used the NuRA dataset to develop an ensemble of ML-based models to predict the agonism, antagonism, binding and effector binding of small molecules to nine different human NRs. We defined the applicability domain of the ML models as a measure of Tanimoto similarity to the molecules in the training set, which enhanced the performance of the developed classifiers. We further developed a user-friendly web server named 'NR-ToxPred' to predict the binding of chemicals to the nine NRs using the best-performing models for each receptor. This web server is freely accessible at http://nr-toxpred.cchem.berkeley.edu. Users can upload individual chemicals using Simplified Molecular-Input Line-Entry System, CAS numbers or sketch the molecule in the provided space to predict the compound's activity against the different NRs and predict the binding mode for each.

Review of in silico studies dedicated to the nuclear receptor family: Therapeutic prospects and toxicological concerns

Being in the center of both therapeutic and toxicological concerns, NRs are widely studied for drug discovery application but also to unravel the potential toxicity of environmental compounds such as pesticides, cosmetics or additives. High throughput screening campaigns (HTS) are largely used to detect compounds able to interact with this protein family for both therapeutic and toxicological purposes. These methods lead to a large amount of data requiring the use of computational approaches for a robust and correct analysis and interpretation. The output data can be used to build predictive models to forecast the behavior of new chemicals based on their in vitro activities. This atrticle is a review of the studies published in the last decade and dedicated to NR ligands in silico prediction for both therapeutic and toxicological purposes. Over 100 articles concerning 14 NR subfamilies were carefully read and analyzed in order to retrieve the most commonly used computational methods to develop predictive models, to retrieve the databases deployed in the model building process and to pinpoint some of the limitations they faced.

Towards machine learning discovery of dual antibacterial drug-nanoparticle systems

Artificial Intelligence/Machine Learning (AI/ML) algorithms may speed up the design of DADNP systems formed by Antibacterial Drugs (AD) and Nanoparticles (NP). In this work, we used IFPTML = Information Fusion (IF) + Perturbation-Theory (PT) + Machine Learning (ML) algorithm for the first time to study of a large dataset of putative DADNP systems composed by >165 000 ChEMBL AD assays and 300 NP assays vs. multiple bacteria species. We trained alternative models with Linear Discriminant Analysis (LDA), Artificial Neural Networks (ANN), Bayesian Networks (BNN), K-Nearest Neighbour (KNN) and other algorithms. IFPTML-LDA model was simpler with values of Sp ≈ 90% and Sn ≈ 74% in both training (>124 K cases) and validation (>41 K cases) series. IFPTML-ANN and KNN models are notably more complicated even when they are more balanced Sn ≈ Sp ≈ 88.5%-99.0% and AUROC ≈ 0.94-0.99 in both series. We also carried out a simulation (>1900 calculations) of the expected behavior for putative DADNPs in 72 different biological assays. The putative DADNPs studied are formed by 27 different drugs with multiple classes of NP and types of coats. In addition, we tested the validity of our additive model with 80 DADNP complexes experimentally synthetized and biologically tested (reported in >45 papers). All these DADNPs show values of MIC < 50 μg mL^-1 (cutoff used) better that MIC of AD and NP alone (synergistic or additive effect). The assays involve DADNP complexes with 10 types of NP, 6 coating materials, NP size range 5-100 nm vs. 15 different antibiotics, and 12 bacteria species. The IFPTML-LDA model classified correctly 100% (80 out of 80) DADNP complexes as biologically active. IFPMTL additive strategy may become a useful tool to assist the design of DADNP systems for antibacterial therapy taking into consideration only information about AD and NP components by separate.

SepPCNET: Deeping Learning on a 3D Surface Electrostatic Potential Point Cloud for Enhanced Toxicity Classification and Its Application to Suspected Environmental Estrogens

Deep learning (DL) offers an unprecedented opportunity to revolutionize the landscape of toxicity prediction based on quantitative structure-activity relationship (QSAR) studies in the big data era. However, the structural description in the reported DL-QSAR models is still restricted to the two-dimensional level. Inspired by point clouds, a type of geometric data structure, a novel three-dimensional (3D) molecular surface point cloud with electrostatic potential (SepPC) was proposed to describe chemical structures. Each surface point of a chemical is assigned its 3D coordinate and molecular electrostatic potential. A novel DL architecture SepPCNET was then introduced to directly consume unordered SepPC data for toxicity classification. The SepPCNET model was trained on 1317 chemicals tested in a battery of 18 estrogen receptor-related assays of the ToxCast program. The obtained model recognized the active and inactive chemicals at accuracies of 82.8 and 88.9%, respectively, with a total accuracy of 88.3% on the internal test set and 92.5% on the external test set, which outperformed other up-to-date machine learning models and succeeded in recognizing the difference in the activity of isomers. Additional insights into the toxicity mechanism were also gained by visualizing critical points and extracting data-driven point features of active chemicals.

Predicting Potential Endocrine Disrupting Chemicals Binding to Estrogen Receptor α (ERα) Using a Pipeline Combining Structure-Based and Ligand-Based in Silico Methods

The estrogen receptors α (ERα) are transcription factors involved in several physiological processes belonging to the nuclear receptors (NRs) protein family. Besides the endogenous ligands, several other chemicals are able to bind to those receptors. Among them are endocrine disrupting chemicals (EDCs) that can trigger toxicological pathways. Many studies have focused on predicting EDCs based on their ability to bind NRs; mainly, estrogen receptors (ER), thyroid hormones receptors (TR), androgen receptors (AR), glucocorticoid receptors (GR), and peroxisome proliferator-activated receptors gamma (PPARγ). In this work, we suggest a pipeline designed for the prediction of ERα binding activity. The flagged compounds can be further explored using experimental techniques to assess their potential to be EDCs. The pipeline is a combination of structure based (docking and pharmacophore models) and ligand based (pharmacophore models) methods. The models have been constructed using the Environmental Protection Agency (EPA) data encompassing a large number of structurally diverse compounds. A validation step was then achieved using two external databases: the NR-DBIND (Nuclear Receptors DataBase Including Negative Data) and the EADB (Estrogenic Activity DataBase). Different combination protocols were explored. Results showed that the combination of models performed better than each model taken individually. The consensus protocol that reached values of 0.81 and 0.54 for sensitivity and specificity, respectively, was the best suited for our toxicological study. Insights and recommendations were drawn to alleviate the screening quality of other projects focusing on ERα binding predictions.

Machine Learning in Predictive Toxicology: Recent Applications and Future Directions for Classification Models

In recent times, machine learning has become increasingly prominent in predictive toxicology as it has shifted from in vivo studies toward in silico studies. Currently, in vitro methods together with other computational methods such as quantitative structure-activity relationship modeling and absorption, distribution, metabolism, and excretion calculations are being used. An overview of machine learning and its applications in predictive toxicology is presented here, including support vector machines (SVMs), random forest (RF) and decision trees (DTs), neural networks, regression models, naïve Bayes, k-nearest neighbors, and ensemble learning. The recent successes of these machine learning methods in predictive toxicology are summarized, and a comparison of some models used in predictive toxicology is presented. In predictive toxicology, SVMs, RF, and DTs are the dominant machine learning methods due to the characteristics of the data available. Lastly, this review describes the current challenges facing the use of machine learning in predictive toxicology and offers insights into the possible areas of improvement in the field.

Comparing Machine Learning Models for Aromatase (P450 19A1)

Aromatase, or cytochrome P450 19A1, catalyzes the aromatization of androgens to estrogens within the body. Changes in the activity of this enzyme can produce hormonal imbalances that can be detrimental to sexual and skeletal development. Inhibition of this enzyme can occur with drugs and natural products as well as environmental chemicals. Therefore, predicting potential endocrine disruption via exogenous chemicals requires that aromatase inhibition be considered in addition to androgen and estrogen pathway interference. Bayesian machine learning methods can be used for prospective prediction from the molecular structure without the need for experimental data. Herein, the generation and evaluation of multiple machine learning models utilizing different sources of aromatase inhibition data are described. These models are applied to two test sets for external validation with molecules relevant to drug discovery from the public domain. In addition, the performance of multiple machine learning algorithms was evaluated by comparing internal five-fold cross-validation statistics of the training data. These methods to predict aromatase inhibition from molecular structure, when used in concert with estrogen and androgen machine learning models, allow for a more holistic assessment of endocrine-disrupting potential of chemicals with limited empirical data and enable the reduction of the use of hazardous substances.

Degradation mechanism for Zearalenone ring-cleavage by Zearalenone hydrolase RmZHD: A QM/MM study

The danger of zearalenone (ZEN) as an endocrine disruptor to humans and the environment has aroused increasing attention. In this study, we implemented the quantum mechanics/molecular mechanics (QM/MM) method to investigate the degradation mechanism of ZEN hydrolase (RmZHD) toward ZEN at the atomic level. The degradation process involves two concerted reaction pathways, where the active site contains a Ser-His-Glu triplet as a proton donor. With the Boltzmann-weighted average potential barriers of 18.1 and 21.5 kcal/mol, the process undergoes proton transfer and nucleophilic-substituted ring opening to form a hydroxyl product. Non-covalent interaction analyses elucidated hydrogen bonding between key amino acids with ZEN. The electrostatic influence analysis of 16 amino acids proposes residues Asp34 and His128 as the possible mutation target for future mutation design of enzyme RmZHD. An in-depth investigation of the protein environment of RmZHD can improve the bioremediation efficiency of endocrine disrupting chemicals.

Binding interactions of halo-benzoic acids, halo-benzenesulfonic acids and halo-phenylboronic acids with human transthyretin

The anionic form-dependent binding interaction of halo-phenolic substances with human transthyretin (hTTR) has been observed previously. This indicates that ionizable compounds should be the primary focus in screening potential hTTR disruptors. Here, the potential binding potency of halo-benzoic acids, halo-benzenesulfonic acids/sulfates and halo-phenylboronic acids with hTTR was determined and analyzed by competitive fluorescence displacement assay integrated with computational methods. The laboratorial results indicated that the three test groups of model compounds exhibited a distinct binding affinity to hTTR. All the tested halo-phenylboronic acids, some of the tested halo-benzoic acids and halo-benzenesulfonic acids/sulfates were shown to be inactive with hTTR. Other halo-benzoic acids and halo-benzenesulfonic acids/sulfates were moderate and/or weak hTTR binders. The binding affinity of halo-benzoic acids and halo-benzenesulfonic acids/sulfates with hTTR was similar. The low distribution ability of the model compounds from water to hTTR may be the reason why they exhibited the binding potency observed with hTTR. By introducing other highly hydrophobic compounds, we observed that the binding affinity between compounds and hTTR increased with increasing molecular hydrophobicity. Those results indicated that the highly hydrophobic halo-benzoic acids and halo-benzenesulfonic acids/sulfates may be high-priority hTTR disruptors. Finally, a binary classification model was constructed employing three predictive variables. The sensitivity (S_n), specificity (S_p), predictive accuracy (Q) values of the training set and validation set were >0.83, indicating that the model had good classification performance. Thus, the binary classification model developed here could be used to distinguish whether a given ionizable compound is a potential hTTR binder or not.

Assessment of the endocrine-disrupting effects of organophosphorus pesticide triazophos and its metabolites on endocrine hormones biosynthesis, transport and receptor binding in silico

Triazophos (TAP) was a widely used organophosphorus insecticide in developing countries. TAP could produce specific metabolites triazophos-oxon (TAPO) and 1-phenyl-3-hydroxy-1,2,4-triazole (PHT) and non-specific metabolites diethylthiophosphate (DETP) and diethylphosphate (DEP). The objective of this study involved computational approaches to discover potential mechanisms of molecular interaction of TAP and its major metabolites with endocrine hormone-related proteins using molecular docking in silico. We found that TAP, TAPO and DEP showed high binding affinity with more proteins and enzymes than PHT and DETP. TAP might interfere with the endocrine function of the adrenal gland, and TAP might also bind strongly with glucocorticoid receptors and thyroid hormone receptors. TAPO might disrupt the normal binding of androgen receptor, estrogen receptor, progesterone receptor and adrenergic receptor to their natural hormone ligands. DEP might affect biosynthesis of steroid hormones and thyroid hormones. Meanwhile, DEP might disrupt the binding and transport of thyroid hormones in the blood and the normal binding of thyroid hormones to their receptors. These results suggested that TAP and DEP might have endocrine disrupting activities and were potential endocrine disrupting chemicals. Our results provided further reference for the comprehensive evaluation of toxicity of organophosphorus chemicals and their metabolites.