Traditional herbal compounds as candidates to inhibit the SARS-CoV-2 main protease: an in silico study
Osmair Vital de Oliveira,Maria Cristina Andreazza Costa,Ricardo Marques da Costa,Rafael Giordano Viegas,Andrew S. Paluch,Márcia Miguel Castro Ferreira
DOI: https://doi.org/10.1080/07391102.2021.2023646
2022-01-04
Journal of Biomolecular Structure and Dynamics
Abstract:COVID-19, a disease caused by the SARS-CoV-2 virus, is responsible for a pandemic since March 2020 and it has no cure. Therefore, herein, different theoretical methods were used to obtain potential candidates from herbal compounds to inhibit the SARS-CoV-2 main protease (M<sup>pro</sup>). Initially, the 16 best-scored compounds were selected from a library containing 4066 ligands using virtual screening by molecular docking. Among them, six molecules (physalin B 5,6-epoxide (PHY), methyl amentoflavone (MAM), withaphysalin C (WPC), daphnoline or trilobamine (TRI), cepharanoline (CEP) and tetrandrine (TET)) were selected based on Lipinski's rule and ADMET analysis as criteria. These compounds complexed with the M<sup>pro</sup> were submitted to triplicate 100<b> </b>ns molecular dynamics simulations. RMSD, RMSF, and radius of gyration results show that the overall protein structure is preserved along the simulation time. The average ΔG<sub>binding</sub> values, calculated by the MM/PBSA method, were -41.7, -55.8, -45.2, -38.7, -49.3, and -57.9<b> </b>kcal/mol for the PHY-M<sup>pro</sup>, MAM-M<sup>pro</sup>, WPC-M<sup>pro</sup>, CEP-M<sup>pro</sup>, TRI-M<sup>pro</sup>, and TET-M<sup>pro</sup> complexes, respectively. Pairwise decomposition analyses revealed that the binding pocket is formed by His41-Val42, Met165-Glu166-Leu167, Asp187, and Gln189. The PLS regression model generated by QSPR analysis indicated that non-polar and polar groups with the presence of hydrogen bond acceptors play an important role in the herbal compounds-M<sup>pro</sup> interactions. Overall, we found six potential candidates to inhibit the SARS-CoV-2 M<sup>pro</sup> and highlighted key residues from the binding pocket that can be used for future drug design. Communicated by Ramaswamy H. Sarma.