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

Development of a Robust Read-Across Model for the Prediction of Biological Potency of Novel Peroxisome Proliferator-Activated Receptor Delta Agonists

A robust predictive model was developed using 136 novel peroxisome proliferator-activated receptor delta (PPARδ) agonists, a distinct subtype of lipid-activated transcription factors of the nuclear receptor superfamily that regulate target genes by binding to characteristic sequences of DNA bases. The model employs various structural descriptors and docking calculations and provides predictions of the biological activity of PPARδ agonists, following the criteria of the Organization for Economic Co-operation and Development (OECD) for the development and validation of quantitative structure-activity relationship (QSAR) models. Specifically focused on small molecules, the model facilitates the identification of highly potent and selective PPARδ agonists and offers a read-across concept by providing the chemical neighbours of the compound under study. The model development process was conducted on Isalos Analytics Software (v. 0.1.17) which provides an intuitive environment for machine-learning applications. The final model was released as a user-friendly web tool and can be accessed through the Enalos Cloud platform's graphical user interface (GUI).

Hepatotoxicity and lipid metabolism disorders of 8:2 polyfluoroalkyl phosphate diester in zebrafish: In vivo and in silico evidence

8:2 polyfluoroalkyl phosphate diester (8:2 diPAP) has been shown to accumulate in the liver, but whether it induces hepatotoxicity and lipid metabolism disorders remains largely unknown. In this study, zebrafish embryos were exposed to 8:2 diPAP for 7 d. Hepatocellular hypertrophy and karyolysis were noted after exposure to 0.5 ng/L 8:2 diPAP, suggesting suppressed liver development. Compared to the water control, 8:2 diPAP led to significantly higher triglyceride and total cholesterol levels, but markedly lower levels of low-density lipoprotein, implying disturbed lipid homeostasis. The levels of two peroxisome proliferator activated receptor (PPAR) subtypes (pparα and pparγ) involved in hepatotoxicity and lipid metabolism were significantly upregulated by 8:2 diPAP, consistent with their overexpression as determined by immunohistochemistry. In silico results showed that 8:2 diPAP formed hydrogen bonds with PPARα and PPARγ. Among seven machine learning models, Adaptive Boosting performed the best in predicting the binding affinities of PPARα and PPARγ on the test set. The predicted binding affinity of 8:2 diPAP to PPARα (7.12) was higher than that to PPARγ (6.97) by Adaptive Boosting, which matched well with the experimental results. Our results revealed PPAR - mediated adverse effects of 8:2 diPAP on the liver and lipid metabolism of zebrafish larvae.

Characterization of flavonoids with potent and subtype-selective actions on estrogen receptors alpha and beta

The initial step in estrogen-regulated transcription is the binding of a ligand to its cognate receptors, named estrogen receptors (ERα and ERβ). Phytochemicals present in foods and environment can compete with endogenous hormones to alter physiological responses. We screened 224 flavonoids in our engineered biosensor ERα and ERβ PRL-array cell lines to characterize their activity on several steps of the estrogen signaling pathway. We identified 83 and 96 flavonoids that can activate ERα or ERβ, respectively. While most act on both receptors, many appear to be subtype-selective, including potent flavonoids that activate ER at sub-micromolar concentrations. We employed an orthogonal assay using a transgenic zebrafish in vivo model that validated the estrogenic potential of these compounds. To our knowledge, this is the largest study thus far on flavonoids and the ER pathway, facilitating the identification of a new set of potential endocrine disruptors acting on both ERα and ERβ.

Prediction of the Interactions of a Large Number of Per- and Poly-Fluoroalkyl Substances with Ten Nuclear Receptors

Per- and poly-fluoroalkyl substances (PFASs) are persistent, toxic chemicals that pose significant hazards to human health and the environment. Screening large numbers of chemicals for their ability to act as endocrine disruptors by modulating the activity of nuclear receptors (NRs) is challenging because of the time and cost of in vitro and in vivo experiments. For this reason, we need computational approaches to screen these chemicals and quickly prioritize them for further testing. Here, we utilized molecular modeling and machine-learning predictions to identify potential interactions between 4545 PFASs with ten different NRs. The results show that some PFASs can bind strongly to several receptors. Further, PFASs that bind to different receptors can have very different structures spread throughout the chemical space. Biological validation of these in silico findings should be a high priority.

Screening of the Antagonistic Activity of Potential Bisphenol A Alternatives toward the Androgen Receptor Using Machine Learning and Molecular Dynamics Simulation

Over the past few decades, extensive research has indicated that exposure to bisphenol A (BPA) increases the health risks in humans. Toxicological studies have demonstrated that BPA can bind to the androgen receptor (AR), resulting in endocrine-disrupting effects. In recent investigations, many alternatives to BPA have been detected in various environmental media as major pollutants. However, related experimental evaluations of BPA alternatives have not been systematically implemented for the assessment of chemical safety and the effects of structural characteristics on the antagonistic activity of the AR. To promote the green development of BPA alternatives, high-throughput toxicological screening is fundamental for prioritizing chemical tests. Therefore, we proposed a hybrid deep learning architecture that combines molecular descriptors and molecular graphs to predict AR antagonistic activity. Compared to previous models, this hybrid architecture can extract substantial chemical information from various molecular representations to improve the model's generalization ability for BPA alternatives. Our predictions suggest that lignin-derivable bisguaiacols, as alternatives to BPA, are likely to be nonantagonist for AR compared to bisphenol analogues. Additionally, molecular dynamics (MD) simulations identified the dihydrotestosterone-bound pocket, rather than the surface, as the major binding site of bisphenol analogues. The conformational changes of key helix H12 from an agonistic to an antagonistic conformation can be evaluated qualitatively by accelerated MD simulations to explain the underlying mechanism. Overall, our computational study is helpful for toxicological screening of BPA alternatives and the design of environmentally friendly BPA alternatives.

Determination of Molecule Category of Ligands Targeting the Ligand-Binding Pocket of Nuclear Receptors with Structural Elucidation and Machine Learning

The mechanism of transcriptional activation/repression of the nuclear receptors (NRs) involves two main conformations of the NR protein, namely, the active (agonistic) and inactive (antagonistic) conformations. Binding of agonists or antagonists to the ligand-binding pocket (LBP) of NRs can regulate the downstream signaling pathways with different physiological effects. However, it is still hard to determine the molecular type of a LBP-bound ligand because both the agonists and antagonists bind to the same position of the protein. Therefore, it is necessary to develop precise and efficient methods to facilitate the discrimination of agonists and antagonists targeting the LBP of NRs. Here, combining structural and energetic analyses with machine-learning (ML) algorithms, we constructed a series of structure-based ML models to determine the molecular category of the LBP-bound ligands. We show that the proposed models work robustly and with high accuracy (ACC > 0.9) for determining the category of molecules derived from docking-based and crystallized poses. Furthermore, the models are also capable of determining the molecular category of ligands with dual opposite functions on different NRs (i.e., working as an agonist in one NR target, whereas functioning as an antagonist in another) with reasonable accuracy. The proposed method is expected to facilitate the determination of the molecular properties of ligands targeting the LBP of NRs with structural interpretation.

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.