In silico identification of VDR as a new potential target of Thymoquinone

Background and Aim: in recent decades, Nigella sativa and its main compound, thymoquinone, have been attracted much attention among civilizations and medical researchers with their therapeutic effects on many diseases including asthma, fatty liver, hypertension, diabetes, inflammation, cough, bronchitis, Headaches, eczema, fever, dizziness and flu and a variety of cancer types, but many of thymoquinone mechanisms remain unclear. In the present study, computational chemical and biological approaches have been used to identify the potential targets of thymoquinone. in silico target fishing, gene expression data analysis, Quantitative Structure–Activity Relationship (QSAR) and docking studies were the methods used in this study. The results showed that Thymoquinone is a potent antagonist of vitamin D۳ receptor (VDR).Methods: ۱. Gene Expression Data Analysis: in this section, the gene expression data information for TQ-treated HCT۱۱۶ colorectal cancer cells, with the name of GSE۱۲۲۸۶۰_Fold-change_data.txt.gz was downloaded from the Gene Expression Omnibus (GEO) database available at https://www.ncbi.nlm.nih.gov/geo/ with the accession number of GSE۱۲۲۸۶۰ and the data were adjusted based on P-value (less than ۰.۰۵) for TQ treated cells and then the symbols of negatively regulated gene (fold change < ۰) were submitted to EnrichR available at http://amp.pharm.mssm.edu/Enrichr for enrichment analysis. ۲. Target Fishing by PharmaMapper: The PharmMapper online tool (available at http://www.lilab-ecust.cn/pharmmapper/) is a web server for potential drug target identification by reversed pharmacophore matching the query compound against an in-house pharmacophore model database (Wang et al., ۲۰۱۷). Pharmacophore modeling is a broadly used ligand-based method in drug target identification. It refers to a protein−ligand interaction pattern corresponding to a desired pharmacological effect and can be considered as the largest common denominator shared by a set of active molecules (Wang et al., ۲۰۱۷). The SDF file of TQ was downloaded from PubChem database (CID:۱۰۲۸۱) and uploaded to the PharmMapper server. The submission started by selecting “Human Protein Targets Only (v۲۰۱۰, ۲۲۴۱)” as targets set and other parameters were not changed. ۳. Analyzing Protein-Ligand Interaction Using Autodock۴: AutoDock is a suite of C programs used to predict the bound conformations of a small, flexible ligand to a macromolecular target of known structure. The technique combines simulated annealing for conformation searching with a rapid grid‐ based method of energy evaluation (Goodsell et al., ۱۹۹۶). in the beginning the SDF file of TQ was converted to PDB file using OpenBabel and the crystal structure of vitamin D۳ receptor (VDR) was downloaded from protein data bank (PDB) database (PDB ID: ۱DB۱). The structure was adjusted and all non-standard elements were deleted using ViewerLite software. Autodock۴ and Autodocktools were employed to perform docking with flexible ligand and rigid protein (Huey et al., ۲۰۱۱). The results were analyzed with Ligplot+ and Notepad++ after getting the ligand-protein complex written by Autodocktools and converted to PDB using Discovery Studio ۴.۵ software. ۴. Prediction of IC۵۰ by QSAR Modeling: Quantitative Structure-Activity Relationship (QSAR) describes how a known biological activity can be considered as a function of the molecular descriptors derived from the chemical structure of a set of molecules (Bajorath, ۲۰۱۱). Many physiological activities of a molecule can be related to their chemical structure. Molecular descriptors as numerical values associated with the chemical constitution for correlation of chemical structure with various physical properties, chemical reactivity, or biological activity are termed as descriptors for a QSAR model (Roy et al., ۲۰۱۵). Half Maximum Inhibitory Concentration (IC۵۰) is the concentration of an inhibitor that is required to inhibit ۵۰% of the enzyme activity in vitro. In this section, the structures and Information of compounds with known IC۵۰ values against VDR were downloaded from Binding database (available at https://www.bindingdb.org/bind/index.jsp) in SDF and TSV formats, respectively. PaDEL-Descriptor software was used to define molecular descriptors of downloaded structures and thymoquinone SDF file, and then the data were adjusted based on target source organisms (homo sapiens) by Excel. S-MLR software was employed to select descriptors for QSAR equations and finally, building the QSAR model and IC۵۰ value prediction for TQ were performed by Chemoface software using “Partial Least Square” as regression method.Results: ۱. Gene Expression Data Analysis: based on TRRUST Transcription Factors ۲۰۱۹, Human VDR was identified as the transcription factor for the four down regulated genes DKK۴, ERBB۲, DDIT۴ and EGFR with the P-value, Adjusted p-value, Odd Ratio and Combined score of ۰.۰۰۰۱۴۴۰, ۰.۰۰۵۴۸۱, ۱۵.۰۰ and ۱۳۲.۶۷ respectively. ۲. Target Fishing by PharmaMapper: VDR and the pharmacophore model of ۱DB۱ (VDR complexed to ۱,۲۵ DIHYDROXY VITAMIN D۳ structure at PDB database) were predicted by PharmMapper as potential target of TQ with ۱۰ similar pharmacophoric features and fit score of ۳.۸۱۴. ۱,۲۵ DIHYDROXY VITAMIN D۳ or Calcitriol has been identified as VDR antagonist (Reinhart, ۲۰۰۴). ۳. Analyzing Protein-Ligand Interaction Using Autodock۴: the docking results showed that VDR was targeted by TQ by making hydrogen bonds at tyrA:۱۴۳ and serA:۲۳۷ positions with respectively the estimated free energy of binding, final intermolecular energy and final total internal energy of -۵.۳۱, -۵.۶۱ and -۰.۲۳ kcal/mol. ۴. Prediction of IC۵۰ by QSAR Modeling: the predicted IC۵۰ of TQ against was ۲۶۷.۹۹۱۵ nM, while the reported IC۵۰ values for Calcitriol against VDR on binding database range from ۰.۰۶۰۰ nM to ۱۳ nM.Conclusion: In recent decades, computational chemogenomics has emerged as an interdisciplinary field of biomedical and pharmaceutical research aimed at identifying ligands or potential regulators for a wide range of gene products. This approach has accelerated studies of the function of gene products in living things and pharmaceutical discoveries (Brown and Bajorath, ۲۰۱۹). The basic assumption that guides all computational approaches in chemogenomics is that similar compounds bind to similar targets and therefore show a similar binding profile and conversely targets that bind similar ligands have similar binding sites (Jacoby, ۲۰۱۳). A wide range of vegetables, fruits, spices and herbs from long time ago have been known for a number of therapeutic and preventative uses and contain various bioactive compounds. Plants are a natural source of many available drugs that have found their place among modern treatments (Kaur et al., ۲۰۱۸). In recent years, using bioinformatics and combining computational simulations and predictions with laboratory research results, has made it possible to more effectively explore pharmaceutical compounds and predict their role in activity regulation of important therapeutic targets (Zheng et al., ۲۰۱۸). Synergistic effect of both Vitamin D and TQ as anticancer agents has been reported in both in vitro and in vivo studies (Kaseb et al., ۲۰۰۷, Yi et al., ۲۰۰۸, Ordonez-Moran et al., ۲۰۰۵, Hummel et al., ۲۰۱۳, Mohamed et al., ۲۰۱۷). In this study we identified VDR as the potential target of TQ with similar binding profile to Calcitriol as a synthetic physiologically active analog of vitamin D, specifically the vitamin D۳ and VDR complex. These results help us to use chemical similarity as a probe of protein function, predict off-target effects of known drug Calcitriol and predict new biological targets of the both compounds.

کلیدواژه ها:

Thymoquinone ، VDR ، Calcitriol ، Computational Chemogenomics ، Nigella sativa