The high morbidity and mortality rate of Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2) infection arises majorly from the Acute Respiratory Distress Syndrome and “cytokine storm” syndrome, which is sustained by an aberrant systemic inflammatory response and elevated pro-inflammatory cytokines. Thus, phytocompounds with broad-spectrum anti-inflammatory activity that target multiple SARS-CoV-2 proteins will enhance the development of effective drugs against the disease. In this study, an in-house library of 117 steroidal plant-derived pregnanes (PDPs) was docked in the active regions of human glucocorticoid receptors (hGRs) in a comparative molecular docking analysis. Based on the minimal binding energy and a comparative dexamethasone binding mode analysis, a list of top twenty ranked PDPs docked in the agonist conformation of hGR, with binding energies ranging between −9.8 and −11.2 kcal/mol, was obtained and analyzed for possible interactions with the human Janus kinases 1 and Interleukins-6 and SARS-CoV-2 3-chymotrypsin-like protease, Papain-like protease and RNA-dependent RNA polymerase. For each target protein, the top three ranked PDPs were selected. Eight PDPs (bregenin, hirundigenin, anhydroholantogenin, atratogenin A, atratogenin B, glaucogenin A, glaucogenin C and glaucogenin D) with high binding tendencies to the catalytic residues of multiple targets were identified. A high degree of structural stability was observed from the 100 ns molecular dynamics simulation analyses of glaucogenin C and hirundigenin complexes of hGR. The selected top-eight ranked PDPs demonstrated high druggable potentials and favourable in silico ADMET properties. Thus, the therapeutic potentials of glaucogenin C and hirundigenin can be explored for further in vitro and in vivo studies.
Dual targeting of cytokine storm and viral replication in COVID-19 by plant-derived steroidal pregnanes: An in silico perspective