{"title":"3CL<sup>pro</sup> of SARS-CoV-2 as a new target for bufadienolides: in silico and in vitro study.","authors":"Muzaffar Kayumov, Jamoliddin Razzokov, Mukhriddin Makhkamov, Murodjon Radjabov, Nurkhodja Mukhamedov, Makhmudjon Khakimov, Akmal M Asrorov, Okhunjon Khasanov, Ansor Yashinov, Mugrajitdin Tashmukhamedov, Ahmidin Wali, Abulimiti Yili, Sharafitdin Mirzaakhmedov","doi":"10.1007/s10822-025-00623-2","DOIUrl":null,"url":null,"abstract":"<p><p>The rapid spread of SARS-CoV-2 and its widespread public health implications have highlighted the urgent need for effective antiviral therapies. A promising strategy is to investigate natural compounds that may inhibit the key viral targets. In this work, we demonstrated the anti-SARS-CoV-2 activity of six bufadienolides, including bufalin (A), arenobufagin (B), gamabufotalin (C), telocinobufagin (D), marinobufagin (E), and bufarenogin (F) from the venom of the Central Asian green toad, Bufo viridis. Molecular docking assays revealed that all A-F bufadienolides bind to key residues (Thr26, His41, Met49, Met161, and Gln189) in the catalytic pocket of 3-chymotrypsin-like cysteine protease (3CL<sup>pro</sup>), an essential enzyme for viral replication and polyprotein processing. The stability of the protein-ligand complexes was then tested using molecular dynamics (MD) simulations, while the binding free energies were estimated using the umbrella sampling (US) technique. Compounds A (-49.8 ± 1.0 kJ/mol), C (-45.9 ± 1.8 kJ/mol), E (-45.6 ± 1.1 kJ/mol), and F (-45.8 ± 1.9 kJ/mol) had significantly higher binding affinities than compounds B (-14.6 ± 1.1 kJ/mol) and D (-10.3 ± 1.9 kJ/mol). In vitro enzymatic assays also confirmed these results, demonstrating that A and C exhibited potent inhibitory activity against 3CL<sup>pro</sup> with IC<sub>50</sub> values of 1.37 µM and 2 µM, respectively, compared to the other bufadienolides; however, they were less active than the positive control GC376 (IC<sub>50</sub> = 0.27 µM). The experimental results are consistent with the computational observations. In silico ADME profiling also revealed good pharmacokinetic properties, indicating that bufadienolides A-F are lead compounds for further antiviral drug development. Taken together, our results support the hypothesis that bufadienolides are SARS-CoV-2 3CL<sup>pro</sup> inhibitors and elucidate their mechanism of action, thereby laying the foundation for potential therapeutic advances against COVID-19.</p>","PeriodicalId":621,"journal":{"name":"Journal of Computer-Aided Molecular Design","volume":"39 1","pages":"41"},"PeriodicalIF":3.0000,"publicationDate":"2025-07-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":"0","resultStr":null,"platform":"Semanticscholar","paperid":null,"PeriodicalName":"Journal of Computer-Aided Molecular Design","FirstCategoryId":"99","ListUrlMain":"https://doi.org/10.1007/s10822-025-00623-2","RegionNum":3,"RegionCategory":"生物学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":null,"EPubDate":"","PubModel":"","JCR":"Q3","JCRName":"BIOCHEMISTRY & MOLECULAR BIOLOGY","Score":null,"Total":0}
引用次数: 0
Abstract
The rapid spread of SARS-CoV-2 and its widespread public health implications have highlighted the urgent need for effective antiviral therapies. A promising strategy is to investigate natural compounds that may inhibit the key viral targets. In this work, we demonstrated the anti-SARS-CoV-2 activity of six bufadienolides, including bufalin (A), arenobufagin (B), gamabufotalin (C), telocinobufagin (D), marinobufagin (E), and bufarenogin (F) from the venom of the Central Asian green toad, Bufo viridis. Molecular docking assays revealed that all A-F bufadienolides bind to key residues (Thr26, His41, Met49, Met161, and Gln189) in the catalytic pocket of 3-chymotrypsin-like cysteine protease (3CLpro), an essential enzyme for viral replication and polyprotein processing. The stability of the protein-ligand complexes was then tested using molecular dynamics (MD) simulations, while the binding free energies were estimated using the umbrella sampling (US) technique. Compounds A (-49.8 ± 1.0 kJ/mol), C (-45.9 ± 1.8 kJ/mol), E (-45.6 ± 1.1 kJ/mol), and F (-45.8 ± 1.9 kJ/mol) had significantly higher binding affinities than compounds B (-14.6 ± 1.1 kJ/mol) and D (-10.3 ± 1.9 kJ/mol). In vitro enzymatic assays also confirmed these results, demonstrating that A and C exhibited potent inhibitory activity against 3CLpro with IC50 values of 1.37 µM and 2 µM, respectively, compared to the other bufadienolides; however, they were less active than the positive control GC376 (IC50 = 0.27 µM). The experimental results are consistent with the computational observations. In silico ADME profiling also revealed good pharmacokinetic properties, indicating that bufadienolides A-F are lead compounds for further antiviral drug development. Taken together, our results support the hypothesis that bufadienolides are SARS-CoV-2 3CLpro inhibitors and elucidate their mechanism of action, thereby laying the foundation for potential therapeutic advances against COVID-19.
期刊介绍:
The Journal of Computer-Aided Molecular Design provides a form for disseminating information on both the theory and the application of computer-based methods in the analysis and design of molecules. The scope of the journal encompasses papers which report new and original research and applications in the following areas:
- theoretical chemistry;
- computational chemistry;
- computer and molecular graphics;
- molecular modeling;
- protein engineering;
- drug design;
- expert systems;
- general structure-property relationships;
- molecular dynamics;
- chemical database development and usage.