MOLECULAR DOCKING ANALYSIS OF HETEROCYCLIC LIGANDS TARGETING NUDT5 PROTEIN IN BREAST CANCER
HTML Full TextMOLECULAR DOCKING ANALYSIS OF HETEROCYCLIC LIGANDS TARGETING NUDT5 PROTEIN IN BREAST CANCER
R. Govinda Rajan, S. Udhaya Vani, A. P. S. M. Krishna *, Ch. Dhanalakshmi Manimadhuri, G. Harika and S. K. Nujhat Aara
Department of Pharmaceutical Chemistry, Hindu College of Pharmacy, Guntur, Andhra Pradesh, India.
ABSTRACT: Breast cancer is the most common malignancy affecting women worldwide. Molecular docking is an important computational technique in drug discovery interaction between ligands and target proteins. In present study docking analysis binding affinity of 51 heterocyclic ligands against the breast cancer target protein NUDT5 (PDB ID:5NWH) are nucleotide- metabolizing enzymes (NUDIX) hydrolases. It plays a major role in hormone- dependent gene regulation and proliferation in breast cancer cells. The heterocyclic ligands docking analysis were analyzed by Autodock 4.0.1. The docking was validated by re docking co-crystallized inhibitor 9CH into active binding site of NUDT5 protein and evaluated binding energy, RMSD values and molecular interactions. The results were found ligand “GR39 exhibited the lowest binding energy” (-14.12kcal/mol), when compared with reference inhibitor 9CH [-10.59kcal/mol] and standard drugs like STD1 Tamoxifen (-9.65), STD2 Cyclophosphamide (-5.83), STD3 Doxorubicin (-11.74). However the present study is limited to molecular docking analyis. Further studies include in-vitro, in-vivo, ADME and cytotoxicity studies.
Keywords: Breast cancer, Molecular docking and NUDT5 protein
INTRODUCTION: Breast cancer is a malignant tumor that develops in the breast tissues, primarily in the ducts or lobules. It is the most frequently diagnosed cancer among women globally and accounts for a significant proportion of cancer-related deaths 1. It is a heterogenous disease characterized by the molecular subtypes and multiple histological include Ductal carcinoma in situ and Invasive ductal carcinoma most common (80%) 2.
At the molecular level, breast cancer is classified based on receptor status into Hormone receptor- positive and Triple-negative breast cancer. Breast cancer progression involves complex biological process such as Metastasis is primary cause of mortality, occurring years after initial treatment due to residual cancer stem cells and micro-environmental interactions 3.
Molecular docking is widely used in structure-based drug design to screen compounds and optimize drug candidates before experimental testing. The protein NUDT5 (NUDIX hydrolase 5) PDB ID:5NWH is a member of Nudix family of nucleotide-metabolizing enzymes. It functions primarily as an ADP-ribose pyrophosphatase, catalyzing the conversion of ADP-ribose into AMP OR ATP depending on cellular conditions. ATP plays Chromatin remodeling, Gene transcription and Hormone-dependent signaling pathways Fig. 1. This ATP production is particularly important in progesterone and estrogen signaling in breast cancer cells, linking metabolism directly to gene regulation 4. To overcome the problem and side effects of cancer, the present marketed drugs showing various adverse reactions, during treatments. In view of that, we want to develop new entities for the treatment of breast cancer cells 5.
FIG. 1: DEVELOPMENT OF POTENT NUDT5 INHIBITORS FROM A HIGH-THROUGHPUT SCREENING CAMPAIGN. A REPRESENTATION OF THE ENZYME-COUPLED MALACHITE GREEN ASSAY (MG ASSAY). ADPR IS HYDROLYZED TO AMP AND R5P BY NUDT5, THEN R5P IS CONVERTED TO FREE INORGANIC PHOSPHATE DETECTED BY THE MALACHITE GREEN REAGENT 5
MATERIALS AND METHODS:
Softwares: Chemdraw 8.0, Chemspider, Protein Data Bank, Pub Meb, Autodock 4.0.1, Discovery Studio V24, Crystal Protein-PDB CODE:5NWH Fig. 2 6.
FIG. 2: CRYSTAL PROTEIN STRUCTURE OF 5NWH
Ligands:
Standard Ligand: Inhibitor (9CH), STDI (Tamoxifen), STD2 (cyclophosphamide), STD3 (Doxorubicin) 7.
Test Ligand: 51 selected ligands of various heterocyclic derivatives were collected from thesis 8.
The 51 ligands were collected from various literatures, based on antioxidant activty, the selected ligands pharmacophores nucleus were presence in selected ligands such as chromen 2-one, Tetrahydropyrimide,-1H-pyrrole, 1,5 benzodiazepene, 1,4 benzodiazepene, pyrazinyl, pyrrolidine-1-yl, biphenyl, 1,3,5 triazine, 2H-chromen-3-yl, quinoline, phenoxyphenyl, 1-2—hydroxyphenyl, pyrimidin, tiazine, and imadazo [2,1-b] [1,3,4]thiadiazole 9.
Protein Data Bank Information of 5NWH:
Classification: Hydrolase Organism(s): Homo sapiens Expression System: Escherichia coli 10.
Mutation(s): NO
Study of Molecular Docking:
Protein Preparation: Macromolecule was downloaded from protein data bank website (RCSB PDB) 10 Download the macromolecule in PDB formate 11 open discovery studio for protein cleaning, ligand molecules, water molecules and other ions were removed using discovery studio. Open autodock 4.0.1 set preference and make default. Edit option and delete water, under edit select hydrogen and add polar hydrogen only and add kollman charges then check total residues and save molecule in PDBQT format.
Ligand Preparation: Structures of 51 selected ligand molecules, standard and reference inhibitor were drawn by using chemdraw 8.0 software 12. Save the structures in mol format. Covert the mol format structures to pdb format by using BIOVIA 13. Retrieve the standard ligand molecules from chemspider website (http://www.chemspider.com/) and download the 3d structure of inhibitor from the protein data bank 14.
Molecular Docking: The ligand prepared and docked against protein using autodock 4.0.1. Output -lamarckian GA -save the file in DPF format. Launch the auto dock. Open the DLG file in the folder and check the lowest binding energy, runs, RMSD and estimated inhibition constant 15.
Visualization: Open the PDBQT saved file in discovery studioand draw the 2D and 3D diagram 16.
Docking Validation: Docking validation was performed by re-docking the co-crystallized inhibitor 9CH into the active site of the NUDT5 protein (PDB ID: 5NWH) using AutoDock 4.0.1. The reference inhibitor showed a binding energy of -10.59 kcal/mol with an RMSD value of 1.557 Å and inhibition constant of 17.30 nM, indicating reliable docking accuracy. The ligand exhibited significant interactions with key active-site residues including TRP A28, TRP B46, ARG A51, GLU A166, and GLY B135. Since the RMSD value was below 2 Å, the docking protocol was considered validated and suitable for further docking studies 15, 16.
RESULTS AND DISCUSSION: It is of interest to design inhibitors for the breast cancer target NUDT5 using molecular docking based virtual screening followed by molecular docking 17.
STD1 (Tamoxifen) was shown more binding energy (-9.65), when compared with selected ligands of GR01, GR02, GR03, GR04, GR05, GR06, GR07, GR08, GR09, GR10, GR11, GR12, GR13, and GR1. Reference inhibitor (9CH) was shown more binding energy (-10.59), when compared with selected ligands of GR03, GR04, GR05, GR06, GR08, GR12, GR16, GR21, GR25, GR26, GR28, GR37, GR38, GR39, GR46, and GR47 respectively.
For docking results, amino acid residues linked with Reference inhibitor (9CH) binding position in the receptor at TRP A 28 (Pi-Pi Stacked Bond), TRP B 46 (Pi-Pi stacked bond), THR B 45 (Pi-Pi Stacked Bond), GLU B 47 (Carbon Hydrogen Bond), ARG A 51 (Conventional Hydrogen bond), GLU A 166 (Carbon Hydrogen Bond), GLY B 135 (Carbon Hydrogen Bond), ALA A 96 (Pi Alkyl Bond), ILE A 141 (Alkyl Bond). The Pose Run, RMSD value and the Inhibition Constant for Reference inhibitor (9CH) was found 41, 1.557 and 17.30 nM respectively.
GR16, GR17, GR18, GR19, GR20, GR21, GR22, GR23, GR25, GR26, GR28, GR30, GR37, GR38, GR39, GR41, GR44, GR46, GR47, GR49 Respectively: For docking results, amino acid residues with STD1 (Tamoxifen) binding position in the receptor at TRP A 28 (Pi-Pi T shaped bond), TRP B 46(Pi-Pi Stacked Bond), ARG A 51(Pi Alkyl bond), LEU A 98(Carbon Hydrogen Bond), GLU A 166 (Carbon Hydrogen Bond). The Pose Run, RMSD value and the Inhibition Constant for STD1 were found 6, 20.390 and 85.01nM respectively 18.
STD2 (cyclophosphamide) was shown more binding energy (-5.83), when compared with All the selected ligands. For docking results, amino acid residues with STD2 (cyclophosphamide) binding position in the receptor at TRP B 46 (Pi cation bond), ARG A 51 (Conventional Hydrogen Bond), GLY B 135 (Pi Cation Bond), TRP A 28 (Pi Alkyl bond). The Pose Run, RMSD value and the Inhibition Constant for STD2 were found 34, 19.875 and 53,670 NM respectively 19.
STD3 (Doxorubicin) was shown more binding energy (-11.74), when compared with selected ligands of GR25, GR37, GR38, GR39 respectively. For docking results, amino acid residues with STD3 (doxorubicin) binding position in the receptor at ASP B 164 (Conventional Hydrogen Bond), GLU B 166(Carbon Hydrogen Bond), GLY B 165 (Conventional Hydrogen Bond), TYR A 36 (Conventional Hydrogen Bond), TRP A 46(Pi-Pi T shaped bond), PHE B 165 (Carbon Hydrogen Bond), VAL B 168 (Alkyl bond), ARG B 51(Conventional Hydrogen Bond). The Pose Run, RMSD value and the Inhibition Constant for STD3 (doxorubicin) were found 47, 20.366 and 2.50nM respectively 20.
Grid values of co-ordinates axis X, Y, Z are 2.086, -9.948, -17.613
Spacing of atoms 0.375
Angle of Dimensions X, Y, Z – 126,126,126
Docking Studies: The selected ligands, standards and inhibitors interaction with protein were performed by using Autodock 4.0.1 and the output of docking results was calculated from dlg file Table 2. In comparison binding energy with Reference Inhibitor (-10.59), STD1 (-9.65), STD2 (-5.83), STD3 (-11.74); the selected ligand GR39 was showed lowest energy (-14.12). For docking results, amino acid residues with selected ligand GR39 binding position in the receptor at TRP B 46 (Napthyl ring, Pi-Pi stacked bond), ARG A 51(phenyl ethyl ring, Pi Cation bond), ARG A 84(Pi Cation bond), GLN B 15(Amino phenyl ring, Conventional Hydrogen Bond ), TYR B 36(Phenyl ring, Pi Donar Hydrogen Bond). The Pose Run, RMSD value and the Inhibition Constant for selected ligand GR39 were found 47, 30.726 and 44.59nM respectively. Fig. 3 showing the 2D and 3D structures of GR39 selected ligand interaction with the crystal protein 5NWH 21-22.
FIG. 3: 3D AND 2D STRUCTURES OF GR39 (E) (-1-(4-(4-(4-AMINOPHENYLAMINO)-6-(NAPHTHALEN-1-YLOXY)-1,3,5-TRIAZINE-2-YLAMINO) PHENYL)-3-(4-ETHYLPHENYL) PROP-2-EN-1-ONE)
FIG. 4: 3D AND 2D STRUCTURE OF INHIBITOR (9CH)
FIG. 5: 3D AND 2D STRUCTURE OF STD1 (TAMOXIFEN)
FIG. 6: 3D AND 2D STRUCTURE OF STD2 (CYCLOPHOSPHAMIDE)
FIG. 7: 3D AND 2D STRUCTURE OF STD3 (DOXORUBICIN)
TABLE 1: CHEMICAL PROPERTIES OF SELECTED LIGANDS, STANDARD DRUGS AND INHIBITOR’S
TABLE 2: DOCKING VALIDATION RESULTS OF SELECTED LIGANDS, STANDARD DRUGS AND INHIBITOR
| S. no. | Code | Lowest Binding Energy | Run | RMSD | Inhibition Constant (Nm) |
| 1 | GR01 | -10.11 | 10 | 20.842 | 39.01 |
| 2 | GR02 | -9.73 | 44 | 19.766 | 73.89 |
| 3 | GR03 | -10.96 | 42 | 17.240 | 9.23 |
| 4 | GR04 | -10.67 | 15 | 20.058 | 15.20 |
| 5 | GR05 | -11.14 | 25 | 30.720 | 6.81 |
| 6 | GR06 | -11.27 | 15 | 14.41 | 5.45 |
| 7 | GR07 | -9.88 | 2 | 17.564 | 57.58 |
| 8 | GR08 | -10.69 | 21 | 16.289 | 14.62 |
| 9 | GR09 | -9.73 | 28 | 15.760 | 73.67 |
| 10 | GR10 | -9.76 | 39 | 25.877 | 70.04 |
| 11 | GR11 | -10.13 | 35 | 16.210 | 37.83 |
| 12 | GR12 | -11.30 | 17 | 16.419 | 5.18 |
| 13 | GR13 | -10.03 | 9 | 26.957 | 44.16 |
| 14 | GR14 | -9.87 | 30 | 26.729 | 58.58 |
| 15 | GR15 | -8.44 | 42 | 27.418 | 651.27 |
| 16 | GR16 | -10.84 | 22 | 21.987 | 11.38 |
| 17 | GR17 | -10.37 | 4 | 21.906 | 25.18 |
| 18 | GR18 | -10.10 | 24 | 22.229 | 39.41 |
| 19 | GR19 | -9.84 | 6 | 16.500 | 61.08 |
| 20 | GR20 | -9.74 | 49 | 21.782 | 72.82 |
| 21 | GR21 | -10.63 | 36 | 28.176 | 16.10 |
| 22 | GR22 | -9.80 | 41 | 22.173 | 65.41 |
| 23 | GR23 | -10.59 | 30 | 29.195 | 17.24 |
| 24 | GR24 | -8.84 | 29 | 20.690 | 332.21 |
| 25 | GR25 | -13.14 | 33 | 24.165 | 232.57 |
| 26 | GR26 | -11.70 | 47 | 22.571 | 2.67 |
| 27 | GR27 | -9.56 | 18 | 17.957 | 98.70 |
| 28 | GR28 | -10.61 | 27 | 24.557 | 16.75 |
| 29 | GR29 | -9.56 | 19 | 17.767 | 98.93 |
| 30 | GR30 | -9.86 | 24 | 17.925 | 58.82 |
| 31 | GR31 | -9.28 | 45 | 25.546 | 158.00 |
| 32 | GR32 | -9.62 | 48 | 27.210 | 88.47 |
| 33 | GR33 | -9.21 | 12 | 17.264 | 177.49 |
| 34 | GR34 | -9.19 | 7 | 25.796 | 184.50 |
| 35 | GR35 | -9.28 | 20 | 26.988 | 158.16 |
| 36 | GR36 | -9.59 | 8 | 27.674 | 93.20 |
| 37 | GR37 | -13.33 | 24 | 19.650 | 170.11 |
| 38 | GR38 | -13.08 | 5 | 32.374 | 257.38 |
| 39 | GR39 | -14.12 | 47 | 30.726 | 44.59 |
| 40 | GR40 | -8.51 | 49 | 25.623 | 581.19 |
| 41 | GR41 | -9.75 | 34 | 20.465 | 71.56 |
| 42 | GR42 | -9.46 | 15 | 16.978 | 117.06 |
| 43 | GR43 | -8.24 | 26 | 25.157 | 917.00 |
| 44 | GR44 | -9.74 | 39 | 27.274 | 72.83 |
| 45 | GR45 | -8.12 | 34 | 18.677 | 1,130 |
| 46 | GR46 | -11.46 | 4 | 30.091 | 3.95 |
| 47 | GR47 | -11.35 | 15 | 27.740 | 4.80 |
| 48 | GR48 | -9.63 | 11 | 32.194 | 88.05 |
| 49 | GR49 | -8.76 | 25 | 20.865 | 379.97 |
| 50 | GR50 | -9.56 | 20 | 19.622 | 98.42 |
| 51 | GR51 | -9.17 | 18 | 20.684 | 188.56 |
| 52 | Inhibitor 9CH | -10.59 | 41 | 1.557 | 17.30 |
| 53 | STD1 Tamoxifen | -9.65 | 6 | 20.390 | 85.01 |
| 54 | STD2 Cyclophosphamide | -5.83 | 34 | 19.875 | 53,670 |
| 55 | STD3 Doxorubicin | -11.74 | 47 | 20.366 | 2.50 |
CONCLUSION: Breast cancer is one of the most common types of malignancies in women in world wide. The risk of breast increases with age Between 45 to 50 in females. So we conducted docking studies on the marketed drugs, Reference Inhibitor with selected ligands on breast cancer crystal protein (PDB: 5NWH).
In this study the crystal structure of 5NWH receptor bound to 9CH was replaced with 51 selected ligands as potential target for breast Cancer. So among 51 selected ligands GR39 (-14.12) showed lowest binding energy compared with Reference Inhibitor (-10.59), STD1 (-9.65), STD2 (-5.83), STD3 (-11.74). In future, studies will be performed to establish the in-vitro, in-vivo on breast cancer cell lines.
ACKNOWLEDGEMENT: Nil
CONFLICTS OF INTEREST: Nil
REFERENCES:
- Siegel RL, Miller KD, Fuchs HE and Jemal A: Cancer statistics, 2026. CA Cancer J Clin 2026; 76(1): 17-48.
- Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I and Jemal A: Global cancer statistics 2020. CA Cancer J Clin 2021; 71(3): 209-49.
- Waks AG and Winer EP: Breast cancer treatment: a review. JAMA 2019; 321(3): 288-300.
- Wright RHG and Beato M: NUDT5: a metabolic oncoprotein in hormone signaling and cancer. Trends Cancer 2018; 4(4): 268-72.
- Page BDG, Valerie NCK, Wright RHG, Wallner O, Isaksson R and Carter M: Targeted NUDT5 inhibitors block hormone signaling in breast cancer cells. Nat Commun 2018; 9: 250.
- Trott O and Olson AJ: AutoDock Vina: improving the speed and accuracy of docking. J Comput Chem 2010; 31(2): 455-61.
- Jordan VC: Tamoxifen: a most unlikely pioneering medicine. Nat Rev Drug Discov 2003; 2(3): 205-13.
- Wishart DS, Burley SK, Berman HM, Christie C, Duarte JM, Feng Z and Westbrook J: RCSB PDB: biological macromolecular structures. Nucleic Acids Res 2019; 47(1): 464-74.
- Patrick GL: An introduction to medicinal chemistry. Ed 7th Oxford: Oxford University Press 2021.
- Burley SK, Berman HM, Christie C, Duarte JM, Feng Z and Westbrook J: RCSB PDB: biological macromolecular structures. Nucleic Acids Res 2019; 47(1): 464-74.
- Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN and Weissig H: The Protein Data Bank. Nucleic Acids Res 2000; 28(1): 235-42.
- RCSB Protein Data Bank. Structure of human NUDT5 (PDB ID: 5NWH). Available from: https://www.rcsb.org/
- Morris GM, Huey R, Lindstrom W, Sanner MF, Belew RK and Goodsell DS: AutoDock4 and AutoDockTools4. J Comput Chem 2009; 30(16): 2785-91
- (http://www.chemspider.com/
- https://autodock.scripps.edu/
- Biovia: Discovery Studio Visualizer user guide. San Diego (CA): Dassault Systems 2024.
- Jain AN: Scoring functions for protein-ligand docking. Curr Protein Pept Sci 2006; 7(5): 407-20.
- Jordan VC and Brodie AMH: Development and evolution of therapies targeting estrogen receptor. Steroids 2007; 72(1): 7-25.
- Emadi A, Jones RJ and Brodsky RA: Cyclophosphamide and cancer. Nat Rev Clin Oncol 2009; 6(11): 638-47.
- Minotti G, Menna P, Salvatorelli E, Cairo G and Gianni L: Anthracyclines: doxorubicin. Pharmacol Rev 2004; 56(2): 185-229.
- Morris GM, Huey R and Lindstrom W: Autodock4 and Autodock tools; automated docking with selective receptor flexibility. J Comput Chem 2009; 30(16): 2785-2791
- Biovia DS: Discovery studio visualizer, version 2021.San diego: Dassault systems BIOVIA 2021.
How to cite this article:
Rajan RG, Vani SU, Krishna APSM, Manimadhuri CD, Harika G and Aara SKN: Molecular docking analysis of heterocyclic ligands targeting NUDT5 protein in breast cancer. Int J Pharm Sci & Res 2026; 17(9): 2660-75. doi: 10.13040/IJPSR.0975-8232.17(9).2660-75.
All © 2026 are reserved by International Journal of Pharmaceutical Sciences and Research. This Journal licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 3.0 Unported License.
Article Information
10
2660-2675
11121 KB
6
English
IJPSR
R. Govinda Rajan, S. Udhaya Vani, A. P. S. M. Krishna *, Ch. Dhanalakshmi Manimadhuri, G. Harika and S. K. Nujhat Aara
Department of Pharmaceutical Chemistry, Hindu College of Pharmacy, Guntur, Andhra Pradesh, India.
annavarapukrishna37@gmail.com
16 May 2026
27 May 2026
19 June 2026
10.13040/IJPSR.0975-8232.17(9).2660-75
01 September 2026



































































