Coronaviruses are contagious pathogens primarily responsible for respiratory and intestinal infections. Research efforts to develop antiviral agents against coronavirus demonstrated the main protease (Mpro) protein may represent effective drug target. X-ray crystallographic structure of the SARS-CoV2 Mpro protein demonstrated the significance of Glu166, Cys141, and His41 residues involved in protein dimerization and its catalytic function. We performed in silico screening of compounds from Curcuma longa L. (Zingiberaceae family) against Mpro protein inhibition. Employing a combination of molecular docking, scoring functions, and molecular dynamics simulations, 267 compounds were screened by docking on Mpro crystallographic structure. Docking score and interaction profile analysis exhibited strong binding on the Mpro catalytic domain with compounds C1 (1E,6E)-1,2,6,7-tetrahydroxy-1,7-bis(4-hydroxy-3-methoxyphenyl)hepta-1,6-diene-3,5-dione) and C2 (4Z,6E)‐1,5‐dihydroxy‐1,7‐bis(4‐hydroxyphenyl)hepta‐4,6‐dien‐3‐one as lead agents. Compound C1 and C2 showed minimum binding score (-9.08 and -8.07 kcal/mole) against Mpro protein in comparison to shikonin and lopinavir (≈ -5.4 kcal/mole) a standard Mpro inhibitor. Furthermore, principal component analysis, free energy landscape and protein-ligand energy calculation studies revealed that these two compounds strongly bind to the catalytic core of the Mpro protein with higher efficacy than lopinavir, a standard antiretroviral of the protease inhibitor class. Taken together, this structure based optimization has provided lead on two natural Mpro inhibitors for further testing and development as therapeutics against human coronavirus.

The present study explores the SARS-CoV-2 drugable target inhibition efficacy of phytochemicals from Indian medicinal plants using molecular docking, molecular dynamics (MD) simulation, and MM-PBSA analysis. A total of 130 phytochemicals were screened against SARS-CoV-2 Spike (S)-protein, RNA-dependent RNA polymerase (RdRp), and Main protease (M^pro). Result of molecular docking showed that Isoquercetin potentially binds with the active site/protein binding site of the Spike, RdRP, and Mpro targets with a docking score of -8.22, -6.86, and -9.73 kcal/mole, respectively. Further, MS 3, 7-Hydroxyaloin B, 10-Hydroxyaloin A, showed -9.57, -7.07, -8.57 kcal/mole docking score against Spike, RdRP, and M^pro targets respectively. The MD simulation was performed to study the favorable confirmation and energetically stable complex formation ability of Isoquercetin and 10-Hydroxyaloin A phytochemicals in M^pro-unbound/ligand bound/standard inhibitor bound system. The parameters such as RMSD, RMSF, Rg, SASA, Hydrogen-bond formation, energy landscape, principal component analysis showed that the lead phytochemicals form stable and energetically stabilized complex with the target protein. Further, MM-PBSA analysis was performed to compare the Gibbs free energy of the M^pro-ligand bound and standard inhibitor bound complexes. The analysis revealed that the His-41, Cys145, Met49, and Leu27 amino acid residues were majorly responsible for the lower free energy of the complex. Drug likeness and physiochemical properties of the test compounds showed satisfactory results. Taken together, the study concludes that that the Isoquercetin and 10-Hydroxyaloin A phytochemical possess significant efficacy to bind SARS-Cov-2 M^pro active site. The study necessitates further in vitro and in vivo experimental validation of these lead phytochemicals to assess their anti-SARS-CoV-2 potential.