Volume 6, Issue 1 (2020)                   Pharm Biomed Res 2020, 6(1): 27-36 | Back to browse issues page


XML Print


Download citation:
BibTeX | RIS | EndNote | Medlars | ProCite | Reference Manager | RefWorks
Send citation to:

Elengoe A, Sundramoorthy N D. Molecular Docking of Curcumin With Breast Cancer Cell Line Proteins. Pharm Biomed Res 2020; 6 (1) :27-36
URL: http://pbr.mazums.ac.ir/article-1-254-en.html
1- Department of Biotechnology, Faculty of Sciences, Lincoln University College, Selangor, Malaysia.
2- Department of Sciences and Biotechnology, Faculty of Engineering and Life Sciences, University of Selangor, Selangor, Malaysia.
Full-Text [PDF 1207 kb]   (1538 Downloads)     |   Abstract (HTML)  (5265 Views)
Full-Text:   (1363 Views)
Introduction
Breast cancer is a known-recognized malignant tumor, causing the highest mortality rate among cancers in women [1]. The reason is the lack of specific signs and symptoms at the early stage of this cancer and its aggressive nature. Currently, breast cancer treatments such as chemotherapy, surgery, and radiotherapy have not been effective. Moreover, these treatments have side effects such as liver, kidney, heart failure and mutation to healthy cells [2].
Plant compounds can be candidates for finding new anticancer agents [3, 4, 5, 6]. Curcumin, which is an active compound from turmeric, has been accounted to have an anti-cancer and chemoprevention effect on breast cancer. According to previous studies, curcumin is a useful herb in cancer treatment. It affects cancer growth, development, and spread at the molecular level. Besides, it kills cancer cells and reduces the growth of new blood vessels in tumors (angiogenesis) and spread of cancer (metastasis). Curcumin exerts its anticancer impact through a complicated molecular signaling network, including multiplication, Estrogen Receptor (ER), and human epidermal growth factor receptor 2 (HER2) pathways [7].
Curcumin is an anti-inflammatory molecule with anticancer properties, similar to fish oil. According to previous studies, curcumin is a potential therapeutic agent in diseases such as inflammatory bowel syndrome, arthritis, and chronic anterior uveitis [8] and prevents heart disease, Alzheimer’s, and cancer. It also alleviates the symptoms of depression.
Computational biology tool helps design substrate-based drug and study the interaction between the phytocompounds and cancer cell proteins. This computer-based approach reduces cost and saves energy and time.
This study aims to generate Three-Dimensional (3D) structures of proteins of breast cancer cell lines, i.e. nuclear factor NF-kappa-B-p105 subunit, caspase 3, and Mucosal Addressin Cell Adhesion Molecule 1 (MADCAM1) and then evaluate their binding affinities with curcumin through local docking.

Materials and Methods
Target sequence
The complete amino acid sequences of nuclear factor NF-kappa-B-p105 subunit, caspase 3, and MADCAM1 were obtained from the National Centre of Biotechnology Information (NCBI) [9]. The template protein PDB file and amino acid sequence in the FASTA format were downloaded from the Protein Data Bank [10].

Homology modeling of protein structure
To date, the 3D models of nuclear factor NF-kappa-B-p105 subunit, caspase 3, and MADCAM1 are not available in Protein Data Bank. Hence, the protein models were generated using the Swiss Model and visualized by the PyMol software [11].

Physiochemical characterization of the protein models
The number of disulfide bonds was calculated using the Cys_Rec program [12]. Protein structural analysis was carried out using the Expasy’s ProtParam Proteomics server [13].

Secondary structure prediction of the protein structure
The secondary structure features were predicted with Self Optimized Prediction Method with Alignment (SOPMA) [14].

Validation of the protein models
The protein structures were validated with PROCHECK by Ramachandran plot analysis [15]. The protein models were further analyzed by ProQ [16], ERRAT [17], and Verify3D [18].

Identification of active sites
The protein models were submitted to active site prediction server-SCFBio to identify the binding sites of the nuclear factor NF-kappa-B-p105 subunits, caspase 3, and MADCAM1 [19].

Protein-ligand docking
The 3D protein model of nuclear kappa NF-kappa-B-p105 subunit, caspase 3, and MADCAM1 subunit were performed for molecular docking with the 3D structure of curcumin using BSP-Slim server [20]. The top-ranked conformation with the lowest binding energy was selected. The same docking simulation approach was performed with the other two protein molecules.

Results

Physiochemical characterization of the protein models
The calculated isoelectric point (pI) value for MADCAM1 and nuclear factor NF-kappa-B-p105 subunits were less than 7, which indicated the acidic characteristics while the pI for caspase 3 was more than 7. Furthermore, the molecular weight of caspase 3, MADCAM1, and nuclear factor NF-kappa-B-p105 were 31641.92, 40155.28, and 1053561.00 Da, respectively. The extinction coefficient was calculated using the extent of light being absorbed by the protein at a particular wavelength range of 27,960-60, and 740 M/cm. The total number of negatively charged residues are (ASP+GLU) is -R, where the total number of positively charged residues are (ARG+LYS) is +R. There were a total of 23 negatively charged and 27 positively charged residues for caspase 3, while 40 negatively charged and 26 positively charged residues for MADCAM1. There are 133  negatively charged and 93 positively charged residues in the amino acid sequence of nuclear factor NF-kappa-B-p105 subunit.
Based on the Expasy’s ProtParam instability index, the nuclear factor NF-kappa-B-p105 subunit was stable because its value was less than 40. Caspase 3 and MADCAM1 had an instability index value of 41.13 and 63.75, respectively. Table 1 presents the aliphatic indexes of all the mutant proteins. The very low grand average of hydropathicity (GRAVY) index (a negative value GRAVY) of all proteins demonstrated their hydrophilic nature.

The thermostability of a molecule increases in the presence of arginine in the protein structure by giving more electrostatic interactions through their guanidine group. Therefore, all proteins have the same stability. Furthermore, the Cys_Rec analysis demonstrated that the number of disulfide bonds was 11 in nuclear factor NF-kappa-B-p105 subunit contrasted with just 8 in caspase 3 and 7 in MADCAM1 (Table 2).

Secondary structure prediction of the protein structures

SOPMA view shows the presence of an alpha helix dominated among secondary structure elements followed by random coils, extended strand, and beta turns at various positions in all protein models (Table 3). Helix structure enhances the stability of the protein structure. The result from this analysis indicates that nuclear factor NF-kappa-B-p105 subunits, caspase 3, and MADCAM1 consist of 15, 18, and 19 α helices, respectively. Table 4 presents the shortest and longest α helix of the proteins.


Validation of the protein models
The 3D protein models were evaluated by Ramachandran plot calculations using PROCHECK software for stereochemical quality and geometry of protein. PROCHECK analysis shows that the quality of the MADCAM1 was excellent. The PROCHECK study confirmed that the residues of all the protein structures are in the most favorable area which was more than 80% (Table 5 and Figure 1).


ProQ was used to evaluate the quality of protein models using the Levitt-Gerstein (LG) score and maximum subarray (MaxSub). The results demonstrated the anticipated LG score (>1.5: good model) and predicted the MaxSub score (>0.1 good model) for all protein structures. The protein models were in the acceptable range to generate a good structure.
ERRAT analysis depends on the quality factor for nonbonded atomic interactions in the 3D protein structures. The higher scores suggest the better quality of the protein model. For a good model, the acceptable range is more than 50%. In this study, the ERRAT score for caspase 3 was the highest (92.5%) (Figure 2).

Furthermore, this confirmed that the caspase 3 structure has a good high resolution and quality contrasted with other protein models. The score for MADCAM1 was 86.928% while it was 80.137% for nuclear factor NF-kappa-B-p105 subunit which was the lowest score.
Verify 3D server showed 96.81%, 96.63%, and 95.69 % of the residues in the nuclear factor NF-kappa-B-p105 subunit, caspase 3, and MADCAM1, respectively. They had an average 3D-1D score of more than 0.2, showing that all protein structures had good quality (Figure 3). Thus, docking analysis with ligand will proceed.

Identification of active sites
The size, protein volume, and the residues of active sites in caspase 3, MADCAM1, and nuclear factor NF-kappa-B-p105 subunit were obtained using active site prediction server (Table 6).

Molecular docking

In this study, nuclear factor NF-kappa-B-p105 subunit, caspase 3, and MADCAM1 were successfully docked with curcumin using the BSP-slim server. The best docking orientation was chosen based on the lowest energy value (DGbind) and negative. Docking analysis showed a strong bond between nuclear factor NF-kappa-B-p105 subunit and ligand because it had the lowest energy value (3.165) compared with the other protein models (Table 7).

Discussion
According to Sohilait et al. study (2017), three analogues of curcumin, i.e. 1-(1,3-benzodioxol-5-yl)-5-(4-hydroxy-3-methoxyphenyl) penta-1, 4-dien-3-one, 1-(3,4 dimethoxyphenyl)-5-(4-nitrophenyl) penta-1, 4-dien-3-one, and 1-(4-hydroxy-3- methoxyphenyl)-5-(4-methoxyphenyl) penta-1, 4-dien-3-one could successfully dock into COX-2 active site. Their binding energies were -8.2, -7.6, and -7.5 kcal/mol, respectively. Curcumin analogues were analysed for COX-2 inhibitory as anti-inflammatory activities [21].
Also, proanthocyanidin (PAC) exhibited better binding affinity with BCL-XL and CDK2 proteins. Their binding energies were -5.23 and -5.17 kcal/mol. This study shows that the proanthocyanidin compound has anti-cancer activity. Therefore, the PAC compound could be a potential drug for colon cancer treatment [22].
Based on Pusphalatha et al. study (2017), the binding score results of curcumin with breast cancer proteins (HER2, estrogen receptor, ERBB2, tyrosine kinase, and HSP90) were -12.1, -9.92, -8.30, -8.22, and -8.05 kcal/mol, respectively. The highest docking score of curcumin with HER2 was -12.1 kcal/mol. The more negative docking score indicates stronger binding between the plant compound and breast cancer protein. This study revealed that curcumin had a strong bond with HER2 protein [23].
According to Iheamgham et al. study (2019) [24], seven phytochemical compounds (7-hydroxy-4′-methoxy-3,11-dehydrohomoisoflavanone, 4,4′-dihydroxy-2’-methoxy-chalcone, 7,4′-dihydroxy-3,11-dehydrohomoisoflavanone, luteolin, quercetin-3-methyl, kaempferol-3-O-β-d-xylopyranoside, and kaempferol-3-O-α-l-rhamnopyranosyl-(1 → 2)-β-D-xylopyranoside) from young twigs and leaves of C bonduc had the strong interactions with cancer target proteins such as Tyrosine Kinase (TK), Vascular Endothelial Growth Factor (VEGF) and Matrix Metalloproteinases (MMP). These phytocompounds can be designed as putative lead compounds for developing new anti-cancer drugs.
Based on Mishra et al. study (2019) [25], Flavonoid Compound (FRC-1) of the plant (Launaea procumbens) acts as a potential anticancer agent for the target protein, tyrosine-protein kinase, Fyn. Molecular docking results showed that the binding between the phytocompound and target protein was strong. FRC-1 has high possibilities of being an anticancer drug.
    In this study, the nuclear factor NF-kappa-B-p105 subunit strongly interacted with the curcumin due to its lowest energy value (3.165) compared with the other protein models.

Conclusion
Nuclear factor NF-kappa-B-p105 subunits, caspase 3, and MADCAM1 were docked successfully with the curcumin. Curcumin can be a potential medication for breast cancer. Further studies of this research can be conducted through lab work.

Ethical Considerations
Compliance with ethical guidelines

There was no ethical considerations to be considered in this research.

Funding
This work was supported by Department of Biotechnology, Faculty of Science, Lincoln University College.

Authors' contributions
Conceptualization, methodology, supervision, and funding acquisition: Asita Elengoe; Investigation, writing-review & editing, resources, and writing-original draft: Asita Elengoe and Nishalini Devi Sundramoorthy.

Conflict of interest
The authors declared no conflict of interest.


References
  1. World Health Organization (WHO) [Internet]. [Cited: 12 September 2018]. Available from: http://www.who.int/mediacentre/factsheets/fs297/en/.
  2. Hawkins LK, Hermiston T. Gene delivery from the E3 region of replicating human adenovirus: Evaluation of the E3B region. Gene Ther. 2001; 8:1142-8. [DOI:10.1038/sj.gt.3301507] [PMID]
  3. David M, Ain QU, Ahmad M, Zaman W, Jahan S. A biochemical and histological approach to study antifertility effects of methanol leaf extract of Asplenium dalhousiae Hook in adult male rats. Andrologia. 2019; 51(6):e13262. [DOI:10.1111/and.13262]
  4. Zaman W, Ahmad M, Zafar H, Amina, Lubna F, Ullah S, et al. The quest for some novel antifertility herbals used as male contraceptives in district Shangla, Pakistan. Acta Ecologica Sinica. 2020; 40(1):102-12. [DOI:10.1016/j.chnaes.2019.05.017]
  5. Zeb S, Ali A, Zaman W, Zeb S, Ali S, Ullah F, Shakoor A. Pharmacology, taxonomy, and phytochemistry of the genus artemisia specifically from Pakistan: A comprehensive review. Pharmaceut Biomed Res. 2017; 4(4):1-10. [DOI:10.18502/pbr.v4i4.543]
  6. Zaman W, Shah SN, Ullah F, Ayaz A, Ahmad M, Ali A. Systematic approach to the correct identification of Asplenium dalhousiae (Aspleniaceae) with their medicinal uses. Microsc Res Tech. 2019; 82(4):459-65. [DOI:10.1002/jemt.23189] [PMID]
  7. Somasundaram S, Edmund NA, Moore DT, Small GW, Shi YY, Orlowski RZ. Dietry curcumin inhibits chemotheraphy-induced apoptosis in models of human breast cancer. Cancer Res. 2002; 62(13):3868-75.
  8. Tathagata C, Suman P. Curcumin induces apotosis in human breast cancer cell through p53-dependent Bax induction. Febs Press. 2002; 512(13):334-40. [DOI:10.1016/S0014-5793(02)02292-5]
  9. National Center for Biotechnology Information Search database (NCBI) [Internet]. 2017 [Cited 2020 Apr 13]. Available from: https://www.ncbi.nlm.nih.gov/.%20[cited:%2010th%20July%202017].
  10. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H. The Protein Data Bank, et al. Nucleic Acids Res. 2000; 28(1):235-42. [DOI:10.1093/nar/28.1.235.] [PMID] [PMCID]
  11. Delano WL. PyMOL: An open-source molecular graphics tool. [Internet]. 2001 [Cited 2020 Apr1 3]. Available from: https://www.ccp4.ac.uk/newsletters/newsletter40/11_pymol.pdf
  12. Gasteiger E, Hoogland C, Gattiker A, Duvaud S, Wilkins MR, Appel RD, et al. Protein identification and analysis tools on the ExPASy server. The Proteomics Protocols Handbook. Totowa: Humana Press; 2005. [DOI:10.1385/1-59259-890-0:571]
  13. Roy S, Maheshwari N, Chauhan R, Sen NK, Sharma A. Structure prediction and functional characterization of secondary metabolite proteins of Ocimum. Bioinformation. 2011; 6(8):315-9. [DOI:10.6026/97320630006315] [PMID] [PMCID]
  14. Geourjon C, Deleage G. SOPMA: significant improvements in protein secondary structure prediction by consensus prediction from multiple alignments. Comput Appl Biosci. 1995; 11:681-4. [DOI:10.1093/bioinformatics/11.6.681] [PMID]
  15. Laskowski RA, MacArthur MW, Moss DS, Thornton JM.: PROCHECK: A program to check the stereo chemical quality of protein structures. J Appl Cryst. 1993; 26:283-91. [DOI:10.1107/S0021889892009944]
  16. Wallner B, Elofsson A. Can correct protein models be identified? Protein Sci. 2003; 12:1073-86. [DOI:10.1110/ps.0236803] [PMID] [PMCID]
  17. Colovos C, Yeates TO. Verification of protein structures: Patterns of non-bonded atomic interactions. Protein Sci. 1993; 2:1511-9. [DOI:10.1002/pro.5560020916] [PMID] [PMCID]
  18. Eisenberg D, Luthy R, Bowie JU. VERIFY3D: Assessment of protein models with three-dimensional profiles. Methods Enzymol. 1997; 277:396-404. [DOI:10.1016/S0076-6879(97)77022-8]
  19. Jayaram B. Active site prediction server [Internet]. 2014 [Cited 2020 Apr. 13]. Available from: http://www.scfbio-iitd.res.in/dock/ActiveSite.jsp
  20. Hui SL, Yang Z. BSP-SLIM: A blind low-resolution ligand-protein docking approach using theoretically predicted protein structures. Proteins. 2012; 80(1):93-110. [DOI:10.1002/prot.23165] [PMID] [PMCID]
  21. Sohilait MR, Pranowo HD, Haryadi, W. Molecular docking analysis of curcumin analogues with COX-2. Bioinformation. 2014; 13(11):356-9. [DOI:10.6026/97320630013356] [PMID] [PMCID]
  22. Déziel BA, Patel K, Neto C, Gottschall‐Pass K, Hurta RA. Proanthocyanidins from the American Cranberry (Vaccinium macrocarpon) inhibit matrix metalloproteinase‐2 and matrix metalloproteinase‐9 activity in human prostate cancer cells via alterations in multiple cellular signalling pathways. J Cell Biochem. 2010; 111(3):742-54. [DOI:10.1002/jcb.22761] [PMID]
  23. Pushpalatha R, Selvamuthukumar S, Kilimozhi D. Comparative in silico docking analysis of curcumin and resveratrol on breast cancer proteins and their synergistic effect on mcf-7 cell line. J Young Pharm. 2017; 9(4):480-5. [DOI:10.5530/jyp.2017.9.94]
  24. Iheagwam FN, Ogunlana OO, Ogunlana OE, Isewon I, Oyelade J. Potential anti-cancer flavonoids isolated from caesalpinia bonduc young twigs and leaves: Molecular docking and in silico studies. Bioinform Biol Insights. 2019; 13:1177932218821371. [DOI:10.1177/1177932218821371] [PMID] [PMCID]
  25. Mishra GJ, Prajapati C, Patel M, Reddy MN. Comparative molecular docking analysis of flavonoid component from launaea procumbens (Roxb.) against tyrosine-protein kinase fyn. Int J Pharm Sci Res. 2019; 10(2):537-45.
 
 
Type of Study: Original Research | Subject: Biotechnology

Add your comments about this article : Your username or Email:
CAPTCHA

Rights and permissions
Creative Commons License This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License.

© 2024 CC BY-NC 4.0 | Pharmaceutical and Biomedical Research

Designed & Developed by : Yektaweb