Main Article Content
Abstract
Triple-negative breast cancer (TNBC) is an aggressive subtype with limited targeted therapies, emphasizing the urgent need for novel molecular inhibitors. Epidermal growth factor receptor (EGFR) is overexpressed in TNBC, but current inhibitors show limited clinical efficacy, indicating a significant therapeutic gap. This study aims to evaluate cardiac glycosides from Vernonia amygdalina as EGFR inhibitors using an integrated in silico approach. Docking results demonstrated that vernoniosides A3 exhibited the lowest binding energy of -12.2 kcal/mol, outperforming the native ligand (-9.0 kcal/mol) and gefitinib (-8.5 kcal/mol). Molecular dynamics (MD) simulations revealed that EGFR-vernoniosides A3 complex remained stable, with RMSD values within stable range (0.7-3.8 Å), RMSF values (0.4-8.3 Å), radius of gyration (Rg) (6.0-7.0 Å), solvent-accessible surface area (SASA) (830-871 Ų), while forming up to 15 hydrogen bonds. However, MM-PBSA analysis showed that binding free energy of vernoniosides A3 (321.60 ± 5084.43 kcal/mol) was lower than the native ligand (569.25 ± 8997.10 kcal/mol) and gefitinib (623.87 ± 9862.73 kcal/mol), indicating weaker thermodynamic binding strength. This discrepancy arises as docking measures potential binding in a static state, while MM-PBSA reflects dynamic, realistic binding strength. These results suggest that vernoniosides A3 shows promising binding affinity but requires further optimization and validation for TNBC therapy.
Keywords
Article Details

This work is licensed under a Creative Commons Attribution 4.0 International License.
References
- [1] Bray, F., Laversanne, M., Sung, H., Ferlay, J., Siegel, R.L., Soerjomataram, I., et al. (2024). Global cancer statistics 2022: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin, 74(3), 229-263.
- [2] Obidiro, O., Battogtokh, G., Akala, E. O. (2023). Triple Negative Breast Cancer Treatment Options and Limitations: Future Outlook. Pharmaceutics, 15(7). 1-26.
- [3] Li, X., Yang, J., Peng, L., Sahin, A. A., Huo, L., Ward, K. C., et al. (2016). Triple-negative breast cancer has worse overall survival and cause-specific survival than non-triple-negative breast cancer. Breast Cancer Res Treat, 161(2), 279-287.
- [4] Yamaoka, T., Ohba, M., Ohmori, T. (2017). Molecular-targeted therapies for epidermal growth factor receptor and its resistance mechanismse. Int J Mol Sci, 18(11), 1-22.
- [5] Sharma, B., Singh, V. J., Chawla, P. A. (2021). Epidermal growth factor receptor inhibitors as potential anticancer agents: An update of recent progress. Bioorg Chem, 116, 1-29.
- [6] Hashmi, A. A., Naz, S., Hashmi, S. K, Irfan, M., Hussain, Z. F., Khan, E. Y., et al. (2019). Epidermal growth factor receptor (EGFR) overexpression in triple-negative breast cancer: association with clinicopathologic features and prognostic parameters. Surgical and Experimental Pathology, 2(6), 1-7.
- [7] Gluz, O., Liedtke, C., Gottschalk, N., Pusztai, L., Nitz, U., Harbeck, N. (2009). Triple-negative breast cancer - Current status and future directions. Annals of Oncology, 20(12), 1913-1927.
- [8] Cohen, M. H., Williams, G. A., Sridhara, R., Chen, G., Pazdur, R. (2003). FDA Drug Approval Summary: Gefitinib (ZD1839) (Iressa ® ) Tablets. Oncologist, 8(4), 303-306.
- [9] Cohen, M. H., Johnson, J. R., Chen, Y. F., Sridhara, R., Pazdur, R. (2005). Regulatory Issues-FDA FDA Drug Approval Summary: Erlotinib (Tarceva ® ) Tablets. Oncologist, 10(7), 461-466.
- [10] Sarijowan, T. P. D., Bodhi, W., Lebang, J. S. (2022). Uji Aktivitas Antibakteri Daun Afrika Terhadap Pertumbuhan Bakteri Staphylococcus aureus dan Pseudomonas aeruginosa. PHARMACON, 11(4), 1678-1685.
- [11] Yani, N dan Thristy, I. (2021). Perbandingan Efektifitas Ekstrak Etanol Daun Afrika (Vernonia amygdalina Del) Dengan Simvastatin Terhadap Kadar Trigliserida Tikus Jantan Galur Wistar Yang Diinduksi Kuning Telur. JURNAL PANDU HUSADA, 2(1), 1-7.
- [12] Wang, W. T., Liao, S. F., Wu, Z. L., Chang, C. W., Wu, J. Y. (2020). Simultaneous study of antioxidant activity, DNA protection and anti-inflammatory effect of Vernonia amygdalina leaves extracts. PLoS One, 15(7), 1-17.
- [13] Pham, E. C., Doan, V. V., Thi, T. V. L., Ngo, C. V., Van, L. V. (2024). In vivo and in silico antihypertensive, anti-inflammatory, and analgesic activities of Vernonia amygdalina Del. leaf extracts. Heliyon, 10(19), 1-26.
- [14] Lifiani, R., Harahap, U., Hasibuan, P. A. Z., Satria, D. (2018). Anticancer effect of african leaves (Vernonia amygdalina del.) to T47D cell resistant. Asian Journal of Pharmaceutical and Clinical Research, 11(1), 4-7.
- [15] Hasibuan, P. A. Z., Harahap, U., Sitorus, P., Satria, D. (2020). The anticancer activities of Vernonia amygdalina Delile. Leaves on 4T1 breast cancer cells through phosphoinositide 3-kinase (PI3K) pathway. Heliyon, 6(7), 1-5.
- [16] Satria, D., Hasibuan, P. A. Z., Muhammad, M., Waruwu, S. B., Utomo, R. Y., Ghoran, S. H. (2024). Cytotoxic and apoptotic effect of Vernonia amygdalina Delile. fractions against Hs578t triple-negative breast cancer cell lines. Phytomedicine Plus, 4(4), 1-8.
- [17] Hasibuan, P. A. Z., Sitorus, R. K. U. A. B., Hermawan, A., Huda, F., Waruwu, S. B., Satria, D. (2024). Anticancer activity of the ethylacetate fraction of Vernonia amygdalina Delile towards overexpression of HER-2 breast cancer cell lines. Pharmacia, 71, 1-8.
- [18] Tafrihani, A. S., Hanif, N., Yoga, I. M. B. K., Irmasari, I., Fakih, T. M., Novitasari, D., et al. (2025). A computational study of cardiac glycosides from Vernonia amygdalina as PI3K inhibitors for targeting HER2 positive breast cancer. J Comput Aided Mol Des, 39(1), 1-13.
- [19] Hermawan, A., Satria, D., Hasibuan, P. A. Z., Huda, F., Tafrihan, A. S., Fatimah, N., et al. (2024). Identification of potential target genes of cardiac glycosides from Vernonia amygdalina Delile in HER2+ breast cancer cells. South African Journal of Botany, 164, 401-418.
- [20] Libera, J. L. V., Verdugo, F. D., Jiménez, A. V., Vivanco, G. N., Caballero, J. (2020). LigRMSD: A web server for automatic structure matching and RMSD calculations among identical and similar compounds in protein-ligand docking. Bioinformatics, 36(9), 2912-2914.
- [21] Fadlan, A., Nusantoro, Y. R. (2021). The Effect of Energy Minimization on The Molecular Docking of Acetone-Based Oxindole Derivatives. JKPK (Jurnal Kimia Dan Pendidikan Kimia), 6(1), 69-77.
- [22] Planken, S., Behenna, D. C., Nair, S. K., Johnson, T. O., Nagata, A., Almaden, C., et al. (2017). Discovery of N-((3R,4R)-4-Fluoro-1-(6-((3-methoxy-1-methyl-1H-pyrazol-4-yl)amino)-9-methyl-9H-purin-2-yl)pyrrolidine-3-yl)acrylamide (PF-06747775) through Structure-Based Drug Design: A High Affinity Irreversible Inhibitor Targeting Oncogenic EGFR Mutants with Selectivity over Wild-Type EGFR. J Med Chem, 60, 3002-3019.
- [23] Zhao, Z., Xie, L., Bourne, P. E. (2019). Structural Insights into Characterizing Binding Sites in Epidermal Growth Factor Receptor Kinase Mutants. J Chem Inf Model, 59(1), 453-462.
- [24] Purba, E. R., Saita, E. I., Maruyama, I. N. (2017). Activation of the EGF receptor by ligand binding and oncogenic mutations: The “rotation model”. Cells, 6(2), 1-19.
- [25] Alzain, A. A., Almogaddam, M. A., Yousif, R., Alqarni, M. H., Foudah, A. I., Osman, W., et al. (2025). Molecular Docking, Molecular Dynamics Simulation, and Pharmacophore-Based Virtual Screening Unveil Natural Compounds with TIM-3 Inhibitory Activity. J Pharm Bioallied Sci, 17(2), S1882-S1887.
- [26] Alsaady, I. M., Gattan, H. S., Aljahdali, S. M., Alruhaili, M. H., Dwived, V. D., Azhar, E. I. (2026). Conformational dynamics and binding free energy analyses unveil a stable flavonoid inhibitor of dengue virus NS5 polymerase. Scientific Reports, 16(1), 7761.
- [27] Nurlelasari, N., Firdaus, A. R. R., Harneti, D., Indrayati, N., Baroroh, U., Yusuf, M. (2020). Computational study of the potential molecular target for antibreast cancer activity of limonoid derivatives from chisocheton SP. J Appl Pharm Sci, 10(2), 7-12.
- [28] Isa, M. A, Kappo, A. P. (2025). Exploring phytoconstituents through molecular dynamics simulation: uncovering potential inhibitors for multiple targeted pathways in breast cancer. J Proteins Proteom, 16, 125-140.
- [29] Sharanya, C. S., Sasikala, W. D., Sathi, S. N., Natarajan, K. (2024). Computational screening combined with well tempered metadynamics simulations identifies potential TMPRSS2 inhibitors. Scientific Reports, 14(1), 16197.
- [30] Dauda, G., Olorukooba, A. B., Adebisi, I. M., Ogundele, O. I., Jibril, M. M., Igharo, A. K., et al. (2026). Molecular docking and dynamics identify novel high-affinity plasmepsin II inhibitors from neem phytochemicals for antimalarial drug development. Discover Chemistry, 3(6), 1-26.
- [31] Sahihi, M., Faraudo, J. (2022). Molecular Dynamics Simulations of Adsorption of SARS-CoV-2 Spike Protein on Polystyrene Surface. J Chem Inf Model, 62(16), 3814-3824.
- [32] Sinha, P., Yadav, A. K. (2023). In silico identification and molecular dynamic simulations of derivatives of 6,6-dimethyl-3-azabicyclo[3.1.0]hexane-2-carboxamide against main protease 3CLpro of SARS-CoV-2 viral infection. J Mol Model, 29(130), 1-12.
- [33] Alisaac, A. (2025). In silico analysis of quorum sensing modulators: Insights into molecular docking and dynamics and potential therapeutic applications. PLOS One, 20(6), 1-27.
- [34] Mahfud, R. (2024). Molecular Dynamics Computational Study of Sustainable Green Surfactant for Application in Chemical Enhanced Oil Recovery. ACS Omega, 9(25), 27177-27191.
- [35] Frimayanti, N., Nasution, M. R., Etavianti, E. (2021). Molecular docking and molecular dynamic simulation of 1,5-benzothiazepine chalcone derivative compounds as potential inhibitors for Zika virus helicase. Jurnal Riset Kimia, 12(1), 44–52.
- [36] Singh, M. B., Jain, P., Tomar, J., Kumar, V., Bahadur, I., Arya, D. K., et al. (2022). An In Silico investigation for acyclovir and its derivatives to fight the COVID-19: Molecular docking, DFT calculations, ADME and td-Molecular dynamics simulations. Journal of the Indian Chemical Society, 99(5), 100433
- [37] Nguyen, V. K. T., Camproux, A. C. (2025). Computational modeling of protein – ligand interactions: From binding site identification to pose prediction and beyond. Current Opinion in Structural Biology, 95, 1-3
- [38] Irsal, R. A. P., Gholam, G. M., Firdaus, D. A., Liwanda, N., Chairunisa, F. (2024). Molecular Docking and Dynamics of Xylocarpus granatum as A Potential Parkinson’s Drug Targeting Multiple Enzymes. Borneo Journal of Pharmacy, 7(2), 161–71.
