journal article·2025-02-01Open access

Revitalizing Recovery: Unveiling the Potential of Apigenin and Related Flavonoids in Long COVID-19 Therapy Through Molecular Dynamics Simulation

Muchtaridi Muchtaridi, Riska Prasetiawati, Siti Ajah Alawiah, Shela Salsabila, Taufik Muhammad Fakih, Rina Fajri Nuwarda, Nur Kusaira Khairul Ikram

Applied Sciences· Volume 15· Issue 3

Abstract

Long COVID-19, also known as post-acute sequelae of SARS-CoV-2 infection (PASC), involves symptoms or effects that persist for more than 4 weeks after the initial SARS-CoV-2 infection. One contributing factor to this condition is the disruption in the expression of the antioxidant enzyme Nuclear Factor Erythroid-2 (Nrf2) induced by the COVID-19 infection. Apigenin and related flavonoids, known for their diverse pharmacological activities, including potent antioxidant properties, have emerged as promising candidates for Long COVID-19 therapy. These compounds, particularly apigenin, are recognized for their ability to modulate oxidative stress and inflammation, making them potential activators of the Nrf2 pathway. This study aims to predict the activity of apigenin and its related flavonoids as Nrf2 activators using molecular modeling and molecular dynamics (MD) techniques, providing insights into their therapeutic potential in managing Long COVID-19. The findings from the molecular modeling analysis indicate that apigenin has a favorable affinity, with a free energy value (ΔG) of −6.40 kcal/mol. Additionally, MD simulation results demonstrate the strong stability of the Keap1-apigenin complex, with an average Root Mean Square Deviation (RMSD) value below 0.20 nm and the lowest average Root Mean Square Fluctuation (RMSF) value of 0.86 nm. Using the Molecular Mechanics/Generalized Born Surface Area (MM-GBSA) calculation method, the binding affinity of the Keap1-apigenin complex yields a lower free energy value (ΔG) of −67.039 kJ/mol, consistent with the molecular modeling results. Apigenin also exhibits the ability to inhibit the binding of Nrf2 to the hydrophobic surface of Keap1, with a total energy value of 993.266 kcal/mol and binding affinity value of −1.162 kJ/mol through peptide−receptor docking. In conclusion, the comprehensive results suggest that apigenin has the potential to be a lead compound for developing Nrf2 activators specifically designed for Long COVID-19 therapy.
Source: crossref · Quality: baseline

Research context

OpenAlex topics: Flavonoids in Medical Research · Computational Drug Discovery Methods · Phytochemistry and Biological Activities

SciVX discovery tags are terminology-based. OpenAlex topics are source-provided classifications.

Affiliations

Department of Pharmaceutical Analysis and Medicinal Chemistry, Faculty of Pharmacy, Universitas Padjadjaran, Jalan Raya Bandung Sumedang Km 21, Sumedang 45363, Indonesia

Research Collaboration Centre for Radiopharmaceuticals Theranostic, National Research and Innovation Agency (BRIN), Jalan Raya Bandung Sumedang Km 21, Sumedang 45363, Indonesia

Department of Pharmacy, Faculty of Mathematics and Natural Sciences, Universitas Garut, Jalan Jati 42B, Garut 44151, Indonesia

Department of Pharmacy, Faculty of Mathematics and Natural Sciences, Universitas Islam Bandung, Jalan Rangga Gading 8, Bandung 40116, Indonesia

Institute of Biological Sciences, Faculty of Science, Universiti Malaya, Kuala Lumpur 50603, Malaysia

Funding

Universitas Padjadjaran, Indonesia

References (82)

Show source-provided references
1. Mboussou, F., Farham, B., Nsasiirwe, S., Atagbaza, A., Oyaole, D., Atuhebwe, P.L., Alegana, V., Osei-Sarpong, F., Bwaka, A., and Paluku, G. (2023). COVID-19 Vaccination in the WHO African Region: Progress Made in 2022 and Factors Associated. Vaccines, 11.10.3390/vaccines11051010 ↗
2. Real-time allelic assays of SARS-CoV-2 variants to enhance sewage surveillance10.1016/j.watres.2022.118686 ↗
3. Determine the most common clinical symptoms in COVID-19 patients: A systematic review and meta-analysis
4. Sex differences in sequelae from COVID-19 infection and in long COVID syndrome: A review10.1080/03007995.2022.2081454 ↗
5. Transcription factor Nrf2 as a potential therapeutic target for COVID-1910.1007/s12192-022-01296-8 ↗
6. He, F., Ru, X., and Wen, T. (2020). NRF2, a Transcription Factor for Stress Response and Beyond. Int. J. Mol. Sci., 21.10.3390/ijms21134777 ↗
7. Advances and challenges in therapeutic targeting of NRF210.1016/j.tips.2022.12.003 ↗
8. Antioxidant response elements: Discovery, classes, regulation and potential applications10.1016/j.redox.2018.05.002 ↗
9. The Association Between Antioxidants and COVID-19 Outcomes: A Systematic Review on Observational Studies10.1007/s12011-023-03588-1 ↗
10. Antioxidants and COVID-19
11. Hossain, M.A., Rahman, M.H., Sultana, H., Ahsan, A., Rayhan, S.I., Hasan, M.I., Sohel, M., Somadder, P.D., and Moni, M.A. (2023). An integrated in-silico Pharmaco-BioInformatics approaches to identify synergistic effects of COVID-19 to HIV patients. Comput. Biol. Med., 155.10.1016/j.compbiomed.2023.106656 ↗
12. Literatur Riview: Flavonoid Activity as an Anti-Cancer Compound
13. Molecular Docking Analysis of Phytochemicals from Tinospora Cordifolia as Potential Inhibitor against SARS-CoV-2
14. Drug repurposing for identification of potential spike inhibitors for SARS-CoV-2 using molecular docking and molecular dynamics simulations10.1016/j.ymeth.2022.02.004 ↗
15. The upshot of Polyphenolic compounds on immunity amid COVID-19 pandemic and other emerging communicable diseases: An appraisal10.1007/s13659-020-00271-z ↗
16. Molecular Dynamics Simulation of Cellulose Synthase Subunit D Octamer with Cellulose Chains from Acetic Acid Bacteria: Insight into Dynamic Behaviors and Thermodynamics on Substrate Recognition10.1021/acs.jctc.0c01027 ↗
17. Mechanistic insight of α-mangostin encapsulation in 2-hydroxypropyl-β-cyclodextrin for solubility enhancement10.1080/07391102.2023.2214237 ↗
18. Pharmacophore-guided virtual screening and dynamic simulation of Kallikrein-5 inhibitor: Discovery of potential molecules for rosacea therapy10.1016/j.imu.2022.100844 ↗
19. Molecular Docking and Interaction Analysis of Propolis Compounds Against SARS-CoV-2 Receptor10.11594/jtls.12.02.08 ↗
20. Computational protein–ligand docking and virtual drug screening with the AutoDock suite10.1038/nprot.2016.051 ↗
21. Pantsar, T., and Poso, A. (2018). Binding affinity via docking: Fact and fiction. Molecules, 23.10.3390/molecules23081899 ↗
22. Varma, A.K., Patil, R., Das, S., Stanley, A., Yadav, L., and Sudhakar, A. (2010). Optimized hydrophobic interactions and hydrogen bonding at the target-ligand interface leads the pathways of Drug-Designing. PLoS ONE, 5.10.1371/journal.pone.0012029 ↗
23. Sang, P., Chen, Y.Q., Liu, M.T., Wang, Y.T., Yue, T., Li, Y., Yin, Y.R., and Yang, L.Q. (2022). Electrostatic Interactions Are the Primary Determinant of the Binding Affinity of SARS-CoV-2 Spike RBD to ACE2: A Computational Case Study of Omicron Variants. Int. J. Mol. Sci., 23.10.3390/ijms232314796 ↗
24. Molecular docking of potential SARS-CoV-2 papain-like protease inhibitors10.1016/j.bbrc.2020.11.083 ↗
25. Discovery of Potential Prolyl-tRNA Synthetase Allosteric Inhibitor Through Virtual Screening and In Vitro Assay against Plasmodium falciparum10.35516/jjps.v16i4.1027 ↗
26. In-silico study for the identification of potential destabilizers between the spike protein of SARS-CoV-2 and human ACE-210.1016/j.imu.2023.101278 ↗
27. Structure of the Keap1:Nrf2 interface provides mechanistic insight into Nrf2 signaling10.1038/sj.emboj.7601243 ↗
28. Guedes, I.A., Costa, L.S.C., dos Santos, K.B., Karl, A.L.M., Rocha, G.K., Teixeira, I.M., Galheigo, M.M., Medeiros, V., Krempser, E., and Cust?dio, F.L. (2021). Drug design and repurposing with DockThor-VS web server focusing on SARS-CoV-2 therapeutic targets and their non-synonym variants. Sci. Rep., 11.10.1038/s41598-021-84700-0 ↗
29. Huang, Y., Yang, W., Yang, L., Wang, T., Li, C., Yu, J., Zhang, P., Yin, Y., Li, R., and Tao, K. (2023). Nrf2 inhibition increases sensitivity to chemotherapy of colorectal cancer by promoting ferroptosis and pyroptosis. Sci. Rep., 13.10.1038/s41598-023-41490-x ↗
30. Molecular docking study of lamellarin analogues and identification of potential inhibitors of HIV-1 integrase strand transfer complex by virtual screening10.1016/j.heliyon.2019.e02811 ↗
31. A BOILED-Egg To Predict Gastrointestinal Absorption and Brain Penetration of Small Molecules10.1002/cmdc.201600182 ↗
32. pkCSM: Predicting Small-Molecule Pharmacokinetic and Toxicity Properties Using Graph-Based Signatures10.1021/acs.jmedchem.5b00104 ↗
33. Kalasariya, H.S., Patel, N.B., Gacem, A., Alsufyani, T., Reece, L.M., Yadav, V.K., Awwad, N.S., Ibrahium, H.A., Ahn, Y., and Yadav, K.K. (2022). Marine Alga Ulva fasciata-Derived Molecules for the Potential Treatment of SARS-CoV-2: An In Silico Approach. Mar. Drugs, 20.10.3390/md20090586 ↗
34. in silico Analysis of 4-((1-(3-Nitrophenyl)-5-oxo-3-phenyl-4,5-dihydro-1H-pyrazol-4-yl)methyl)benzoic Acid: An Emerging 3-CLpro Non-peptidic Inhibitors for COVID-1910.14233/ajomc.2022.AJOMC-P376 ↗
35. The Prediction of Pharmacokinetic Properties of Compounds in Hemigraphis alternata (Burm. F.) T. Ander Leaves Using pkCSM10.22146/ijc.73117 ↗
36. Allemailem, K.S., Almatroudi, A., Alharbi, H.O.A., AlSuhaymi, N., Alsugoor, M.H., Aldakheel, F.M., Khan, A.A., and Rahmani, A.H. (2024). Apigenin: A Bioflavonoid with a Promising Role in Disease Prevention and Treatment. Biomedicines, 12.10.3390/biomedicines12061353 ↗
37. Tsiailanis, A.D., Tellis, C.C., Papakyriakopoulou, P., Kostagianni, A.D., Gkalpinos, V., Chatzigiannis, C.M., Kostomitsopoulos, N., Valsami, G., Tselepis, A.D., and Tzakos, A.G. (2023). Development of a Novel Apigenin Dosage form as a Substitute for the Modern Triple Antithrombotic Regimen. Molecules, 28.10.3390/molecules28052311 ↗
38. Prediction of drug candidate from Rosmarinus officinalis L to inhibit IL-6R, IL-1R1, and TNF-α: In silico study10.30574/wjarr.2024.21.2.0417 ↗
39. Strategically Placed Trifluoromethyl Substituent in the Realm of Antitubercular Drug Design10.2174/1574885513666180906101732 ↗
40. Tripathi, P., Ghosh, S., and Talapatra, S.N. (2024, August 22). Bioavailability Prediction of Phytochemicals Present in Calotropis Procera (Aiton) R. Br. by Using Swiss-ADME Tool. Available online: http://www.swissadme.ch/.
41. Molecular Dynamic Simulation and 3d-pharmacophore Modeling of Alpha Mangostin and Its Derivatives against Estrogen Alpha Receptor10.2174/1570180820666230220122600 ↗
42. Arab, F.L., Hoseinzadeh, A., Mohammadi, F.S., Rajabian, A., Faridzadeh, A., and Mahmoudi, M. (2024). Immunoregulatory effects of nanocurcumin in inflammatory milieu: Focus on COVID-19. Biomed. Pharmacother., 171.10.1016/j.biopha.2024.116131 ↗
43. Describing Chemical Reactivity with Frontier Molecular Orbitalets10.1021/jacsau.2c00085 ↗
44. Frontier molecular orbital, molecular structure and Thermal properties of 2,4,6,8-tetramethyl-2,3,6,7-tetrahydro-s-indacene-1,5-dione using DFT calculation10.1016/j.matpr.2022.03.215 ↗
45. Chemical Reactivity as Described by Quantum Chemical Methods10.3390/i3040276 ↗
46. Electrophilicity and nucleophilicity scales at different DFT computational levels10.1002/poc.4503 ↗
47. Domingo, L.R., Ríos-Gutiérrez, M., and Pérez, P. (2016). Applications of the Conceptual Density Functional Theory Indices to Organic Chemistry Reactivity. Molecules, 21.10.3390/molecules21060748 ↗
48. Targeting of protein ubiquitination by BTB-Cullin 3-Roc1 ubiquitin ligases10.1038/ncb1056 ↗
49. Nrf2-Keap1 signaling in oxidative and reductive stress10.1016/j.bbamcr.2018.02.010 ↗
50. Guan, T., Bian, C., and Ma, Z. (2023). In vitro and in silico perspectives on the activation of antioxidant responsive element by citrus-derived flavonoids. Front. Nutr., 10.10.3389/fnut.2023.1257172 ↗
51. Huang, W., Zhong, Y., Gao, B., Zheng, B., and Liu, Y. (2023). Nrf2-mediated therapeutic effects of dietary flavones in different diseases. Front. Pharmacol., 14.10.3389/fphar.2023.1240433 ↗
52. Mendonca, P., and Soliman, K.F.A. (2020). Flavonoids Activation of the Transcription Factor Nrf2 as a Hypothesis Approach for the Prevention and Modulation of SARS-CoV-2 Infection Severity. Antioxidants, 9.10.3390/antiox9080659 ↗
53. Therapeutic and prophylactic effect of flavonoids in post-COVID-19 therapy10.1002/ptr.7436 ↗
54. Antiviral activities of olive oil apigenin and taxifolin against SARS-CoV-2 RNA-dependent RNA polymerase (RdRP): In silico, pharmacokinetic, ADMET, and in-vitro approaches10.1080/23311932.2023.2236828 ↗
55. KEAP1-NRF2 signalling and autophagy in protection against oxidative and reductive proteotoxicity10.1042/BJ20150568 ↗
56. Apigenin prevents metabolic syndrome in high-fructose diet-fed mice by Keap1-Nrf2 pathway10.1016/j.biopha.2018.06.108 ↗
57. Pharmacological stimulation of nuclear factor (erythroidderived 2)-like 2 translation activates antioxidant responses10.1074/jbc.M116.770925 ↗
58. Golonko, A., Olichwier, A.J., Swislocka, R., Szczerbinski, L., and Lewandowski, W. (2023). Why Do Dietary Flavonoids Have a Promising Effect as Enhancers of Anthracyclines? Hydroxyl Substituents, Bioavailability and Biological Activity. Int. J. Mol. Sci., 24.10.3390/ijms24010391 ↗
59. Auti, P.S., Jagetiya, S., and Paul, A.T. (2023). Chromone Containing Hybrid Analogs: Synthesis and Applications in Medicinal Chemistry. Chem. Biodivers., 20.10.1002/cbdv.202300587 ↗
60. Therapeutical properties of apigenin: A review on the experimental evidence and basic mechanisms10.1080/10942912.2023.2236329 ↗
61. Liu, W.-N., Shi, J., Fu, Y., and Zhao, X.-H. (2019). The Stability and Activity Changes of Apigenin and Luteolin in Human Cervical Cancer Hela Cells in Response to Heat Treatment and Fe2+/Cu2+ Addition. Foods, 8.10.3390/foods8080346 ↗
62. Structure, Conformation, and Electronic Properties of Apigenin, Luteolin, and Taxifolin Antioxidants. A First Principle Theoretical Study10.1021/jp035901j ↗
63. PubChem 2023 update10.1093/nar/gkac956 ↗
64. Quantum ESPRESSO toward the exascale10.1063/5.0005082 ↗
65. Small molecules inhibit the interaction of Nrf2 and the Keap1 Kelch domain through a non-covalent mechanism10.1016/j.bmc.2013.04.019 ↗
66. Understanding the binding interaction between phenyl boronic acid P1 and sugars: Determination of association and dissociation constants using S-V plots, steady-state spectroscopic methods and molecular docking10.1002/bio.3931 ↗
67. Forli, W., Halliday, S., Belew, R., and Olson, A. (2024, June 18). AutoDock, Available online: https://autodock.scripps.edu/.
68. OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules10.1021/acs.jctc.8b01026 ↗
69. How to improve docking accuracy of AutoDock4.2: A case study using different electrostatic potentials10.1021/ci300417y ↗
70. Biovia (2017). Dassault Systemes BIOVIA, Discovery Studio Modeling Environment, Release 2017, Dassault Systems.
71. α-Mangostin and its derivatives against estrogen receptor alpha10.1080/07391102.2020.1841031 ↗
72. Budiman, A., Fakih, T.M., Megantara, S., Muchtaridi, M., and Aulifa, D.L. (2025). Leveraging Non-Covalent Adaptable Networks to Stabilize Drug-Polymer Systems in Supersaturated Solutions: A Computational and Experimental Approach. Applied Sciences, 15.10.3390/app15010307 ↗
73. van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD10.1021/acs.jctc.3c01164 ↗
74. Peptide bond planarity constrains hydrogen bond geometry and influences secondary structure conformations10.1016/j.crstbi.2020.11.002 ↗
75. Eastman, P., Swails, J., Chodera, J.D., McGibbon, R.T., Zhao, Y., Beauchamp, K.A., Wang, L.P., Simmonett, A.C., Harrigan, M.P., and Stern, C.D. (2017). OpenMM 7: Rapid development of high performance algorithms for molecular dynamics. PLoS Comput. Biol., 13.10.1371/journal.pcbi.1005659 ↗
76. Alzain, A.A., Mukhtar, R.M., Abdelmoniem, N., Shoaib, T.H., Osman, W., Alsulaimany, M., Aljohani, A.K.B., Almadani, S.A., Alsaadi, B.H., and Althubyani, M.M. (2023). Modulation of NRF2/KEAP1-Mediated Oxidative Stress for Cancer Treatment by Natural Products Using Pharmacophore-Based Screening, Molecular Docking, and Molecular Dynamics Studies. Molecules, 28.10.3390/molecules28166003 ↗
77. Automated Topology Builder Version 3.0: Prediction of Solvation Free Enthalpies in Water and Hexane10.1021/acs.jctc.8b00768 ↗
78. Shiga, M., and Sorai, M. (2023, January 3). Salinity Effect on the Interfacial Tension of CO2-Brine: Estimation and Investigation by Molecular Dynamics Simulations. Proceedings of the International Petroleum Technology Conference, Bangkok, Thailand.10.2523/IPTC-23016-MS ↗
79. Multiple Staggered Mesh Ewald: Boosting the Accuracy of the Smooth Particle Mesh Ewald Method10.1021/acs.jctc.6b00701 ↗
80. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities10.1517/17460441.2015.1032936 ↗
81. CHARMM-GUI Implicit Solvent Modeler for Various Generalized Born Models in Different Simulation Programs10.1021/acs.jpcb.2c05294 ↗
82. Phytoconstituents from ten natural herbs as potent inhibitors of main protease enzyme of SARS-COV-2: In silico study10.1016/j.phyplu.2021.100083 ↗