Structure-based design approach of potential BCL-2 inhibitors for cancer chemotherapy

Cell death pathways: molecular mechanisms and therapeutic targets for cancer

Cell death regulation is essential for tissue homeostasis and its dysregulation often underlies cancer development. Understanding the different pathways of cell death can provide novel therapeutic strategies for battling cancer. This review explores several key cell death mechanisms of apoptosis, necroptosis, autophagic cell death, ferroptosis, and pyroptosis. The research gap addressed involves a thorough analysis of how these cell death pathways can be precisely targeted for cancer therapy, considering tumor heterogeneity and adaptation. It delves into genetic and epigenetic factors and signaling cascades like the phosphatidylinositol 3-kinase/protein kinase B/mammalian target of rapamycin (PI3K/AKT/mTOR) and mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) pathways, which are critical for the regulation of cell death. Additionally, the interaction of the microenvironment with tumor cells, and particularly the influence of hypoxia, nutrient deprivation, and immune cellular interactions, are explored. Emphasizing therapeutic strategies, this review highlights emerging modulators and inducers such as B cell lymphoma 2 (BCL2) homology domain 3 (BH3) mimetics, tumour necrosis factor-related apoptosis-inducing ligand (TRAIL), chloroquine, and innovative approaches to induce ferroptosis and pyroptosis. This review provides insights into cancer therapy's future direction, focusing on multifaceted approaches to influence cell death pathways and circumvent drug resistance. This examination of evolving strategies underlines the considerable clinical potential and the continuous necessity for in-depth exploration within this scientific domain.

Pro-apoptotic gene BAX is a pan-cancer predictive biomarker for prognosis and immunotherapy efficacy

BACKGROUND

Apoptosis Regulator BCL2 Associated X (BAX) is a pro-apoptotic gene. Apoptosis is one of the important components of immune response and immune regulation. However, there is no systematic pan-cancer analysis of BAX.

METHODS

Original data of this study were downloaded from TCGA databases and GTEX databases. We conducted the gene expression analysis and survival analysis of BAX in 33 types of cancer via Gene Expression Profiling Interactive Analysis (GEPIA) database. Real-time PCR and immunohistochemistry (IHC) were further performed to examine the BAX expression in cancer cells and tissues. Moreover, the relationship between BAX and immune infiltration and gene alteration was studied by the Tumor Immune Estimation Resource (TIMER) and cBioPortal tools. Protein-protein interaction analysis was performed in the STRING database. Finally, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were utilized to evaluate the enrichment analysis.

RESULTS

BAX was highly expressed in most cancers and was associated with poor prognosis in nine cancer types. In addition, BAX showed significant clinical relevance, and the mRNA expression of BAX was also strongly associated with drug sensitivity of many drugs. Furthermore, BAX may participate in proliferation and metastasis of many cancers and was associated with methylation. Importantly, BAX expression was positively correlated with most immune infiltrating cells.

CONCLUSION

Our findings suggested that BAX can function as an oncogene and may be used as a potential predictive biomarker for prognosis and immunotherapy efficacy of human cancer, which could provide a new approach for cancer therapy.

Investigation of the insecticidal potential of curcumin derivatives that target the Helicoverpa armigera sterol carrier protein-2

Cotton bollworm (Helicoverpa armigera) is a highly polyphagous, widely prevalent, and persistent Old World insect pest that affects numerous important crops that are directly consumed by people, including tomato, cotton, pigeon pea, chickpea, rice, sorghum, and cowpea. Insects do not synthesize steroids but obtain them from their diet. Inhibition of dietary uptake of steroids by insects is a potentially effective insecticidal mechanism that should not be toxic to humans and other mammals, who synthesize their steroids. Ten curcumin derivatives were tested against H. armigera sterol carrier protein-2 (HaSCP-2) for their potential as insecticidal agents. Curcumin derivatives were initially docked at the binding site of HaSCP-2 to determine their binding affinities and plausible binding modes. The binding modes predominantly show hydrophobic interactions of derivatives with Phe53, Phe110, and Phe89 as core interacting residues in the active site. Validation of in silico results was carried out by performing a fluorescence binding and displacement assay to determine the binding affinities of curcumin derivatives. Among a collection of curcumin derivatives tested, Cur10 showed the lowest IC50 value of 9.64 μM, while Cur07 was 19.86 μM, and Cur06 was 20.79 μM. There was a significant negative correlation between the ability of the curcumin derivatives tested to displace the fluorescent probe from the sterol binding site of HaSCP-2 and to inhibit Sf9 insect cell growth in culture, which is consistent with the curcumin derivatives acting by the novel mechanism of blocking sterol uptake. Then molecular dynamics simulation studies validated the predicted binding modes and the interactions of curcumin derivatives with HaSCP-2 protein. In conclusion, these studies support the potential use of curcumin derivatives as insecticidal agents.

Recent trends and advances in novel formulations as an armament in Bcl-2/Bax targeted breast cancer

Breast cancer (BC) remains a significant health burden worldwide, necessitating the development of innovative therapeutic strategies. The B-cell lymphoma 2 (Bcl-2) family proteins, Bcl-2 and Bax, play a crucial role in regulating apoptosis and thus are promising targets for BC therapy. We focus on the recent advancements in novel formulations that specifically target Bcl-2/Bax pathway to combat BC. It provides an overview on biological functions of Bcl-2/Bax in apoptosis regulation, emphasizing their significance in pathogenesis and progression of the disease while covering the numerous therapeutic approaches aimed at modulating the Bcl-2/Bax pathway, including small-molecule inhibitors, peptides, gene-based therapies and other repurposed drugs harboured onto cutting-edge technologies and nanocarrier systems employed to enhance the targeted delivery of Bcl-2/Bax inhibitors tumor cells. These advanced formulations aim to improve therapeutic efficacy, minimize off-target effects, and overcome drug resistance, offering promising prospects in its treatment. In conclusion, it illuminates the diverse and evolving landscape of novel formulations as an essential armament in targeting these proteins while bridging and unravelling the obscurity of Bcl-2/Bax pathway-targeted drug delivery systems which are presently in their nascent stages of exploration for BC therapy which can benefit researchers, clinicians, and pharmaceutical scientists.

Porcupine quills keratin peptides induces G0/G1 cell cycle arrest and apoptosis via p53/p21 pathway and caspase cascade reaction in MCF-7 breast cancer cells

BACKGROUND

Porcupine quills, a by-product of porcupine pork, are rich in keratin, which is an excellent source of bioactive peptides. The objective of this study was to investigate the underlying mechanism of anti-proliferation effect of porcupine quills keratin peptides (PQKPs) on MCF-7 cells.

RESULTS

Results showed that PQKPs induced MCF-7 cells apoptosis by significantly decreasing the secretion level of anti-apoptosis protein Bcl-2 and increasing the secretion levels of pro-apoptosis proteins Bax, cytochrome c, caspase 9, caspase 3 and PARP. PQKPs also arrested the cell cycle at G0/G1 phase via remarkably reducing the protein levels of CDK4 and enhancing the protein levels of p53 and p21. High-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) analysis identified nine peptides with molecular weights less than 1000 Da in PQKPs. Molecular docking results showed that TPGPPT and KGPAC identified from PQKPs could bind with p53 mutant and Bcl-2 protein by conventional hydrogen bonds, carbon hydrogen bonds and van der Waals force. Furthermore, the anti-proliferation impact of synthesized peptides (TPGPPT and KGPAC) was shown in MCF-7 cells.

CONCLUSION

These findings indicated that PQKPs suppressed the proliferation of MCF-7 breast cancer cells by triggering apoptosis and G0/G1 cell cycle arrest. Moreover, the outcome of this study will bring fresh insights into the production and application of animal byproducts. © 2023 Society of Chemical Industry.

Ehretia Species Phytoconstituents as Potential Lead Compounds against Klebsiella pneumoniae Carbapenemase: A Computational Approach

The evolution of antibiotic-resistant carbapenemase has negatively impacted the management of critical healthcare-associated infections. K. pneumoniae carbapenemase-2- (KPC-2-) expressing bacteria have developed resistance to conventional therapeutic options, including those used as a last resort for life-threatening diseases. In this study, Ehretia species phytoconstituents were screened for their potential to inhibit KPC-2 protein using in silico approaches. Molecular docking was used to identify strong KPC-2 protein binding phytoconstituents retrieved from the literature. The best-docked conformation of the ligands was selected based on their glide energy and binding interactions. To determine their binding free energies, these hit compounds were subjected to molecular mechanics with generalized born and surface area (MM-GBSA) in the PRIME module. Pharmacological assessments of the ligands were performed to evaluate their drug-likeness. Molecular dynamic (MD) simulations were used to analyze the conformational stability of the selected druglike compounds within the active site of the KPC-2 protein. Overall, a total of 69 phytoconstituents were compiled from the literature. Fourteen of these compounds exhibited a stronger binding affinity for the protein target than the reference drugs. Four of these top hit compounds, DB09, DB12, DB28, and DB66, revealed the highest efficacy in terms of drug-likeness properties. The MD simulation established that among the druglike compounds, DB66 attained stable conformations after 150 ns simulation in the active site of the protein. We concluded that DB66 from Ehretia species could play a significant role in therapeutic efforts against KPC-2-expressing bacteria.

Exploiting the co-crystal ligands shape, features and structure-based approaches for identification of SARS-CoV-2 Mpro inhibitors

SARS-CoV-2 enters the host cell through the ACE2 receptor and replicates its genome using an RNA-Dependent RNA Polymerase (RDRP). The functional RDRP is released from pro-protein pp1ab by the proteolytic activity of Main protease (Mpro) which is encoded within the viral genome. Due to its vital role in proteolysis of viral polyprotein chains, it has become an attractive potential drug target. We employed a hierarchical virtual screening approach to identify small synthetic protease inhibitors. Statistically optimized molecular shape and color-based features (various functional groups) from co-crystal ligands were used to screen different databases through various scoring schemes. Then, the electrostatic complementarity of screened compounds was matched with the most active molecule to further reduce the hit molecules' size. Finally, five hundred eighty-seven molecules were docked in Mpro catalytic binding site, out of which 29 common best hits were selected based on Glide and FRED scores. Five best-fitting compounds in complex with Mpro were subjected to MD simulations to analyze their structural stability and binding affinities with Mpro using MM/GB(PB)SA models. Modeling results suggest that identified hits can act as the lead compounds for designing better active Mpro inhibitors to enhance the chemical space to combat COVID-19.Communicated by Ramaswamy H. Sarma.

In Silico Computational Studies of Bioactive Secondary Metabolites from Wedelia trilobata against Anti-Apoptotic B-Cell Lymphoma-2 (Bcl-2) Protein Associated with Cancer Cell Survival and Resistance

In the present study, the binding affinity of 52 bioactive secondary metabolites from Wedelia trilobata towards the anti-apoptotic B-cell lymphoma-2 (Bcl-2) protein (PDB: 2W3L) structure was identified by using in silico molecular docking and molecular dynamics simulation. The molecular docking results demonstrated that the binding energies of docked compounds with Bcl-2 protein ranged from -5.3 kcal/mol to -10.1 kcal/mol. However, the lowest binding energy (-10.1 kcal/mol) was offered by Friedelin against Bcl-2 protein when compared to other metabolites and the standard drug Obatoclax (-8.4 kcal/mol). The molecular dynamics simulations revealed that the Friedelin-Bcl-2 protein complex was found to be stable throughout the simulation period of 100 ns. Overall, the predicted Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) properties of Friedelin are relatively better than Obatoclax, with the most noticeable differences in many parameters where Friedelin has no AMES toxicity, hepatotoxicity, and skin sensitization. The ADMET profiling of selected compounds supported their in silico drug-likeness properties. Based on the computational analyses, the present study concluded that Friedelin of W. trilobata was found to be the potential inhibitor of the Bcl-2 protein, which merits attention for further in vitro and in vivo studies before clinical trials.

Computational design of MmpL3 inhibitors for tuberculosis therapy

Tuberculosis is a chronic communicable disease caused by Mycobacterium tuberculosis (Mtb) and spreads from lungs to lymphatic system. The cell wall of mycobacterium plays a prominent role in maintaining the virulence and pathogenicity and also acts as prime target for drug discovery. Hence, this study has put into emphasis with target MmpLs (Mycobacterial membrane proteins Large) which are significant for the growth and survival of Mycobacterium tuberculosis. MmpLs belongs to the resistance, nodulation and division (RND) protein superfamily. MmpL3 is the only MmpL deemed essential for the replication and viability of mycobacterial cells. For the study, we have selected SQ109 derivatives as Mmpl3 inhibitor, which holds non-covalent property. Structure-based pharmacophore model of MmpL3 target protein with SQ109 as co-crystallized ligand (PDB: 6AJG) was generated to screen the ligand database. Compounds with decent fitness score and pharmacophoric features were compared with standard drug and taken for molecular docking studies. Further prime molecular mechanics-Poisson-Boltzmann surface area (MM-GBSA) and induced fit calculations identified potential molecules for further drug-likeness screening. Overall computational calculations identified ZINC000000016638 and ZINC000000003594 as potential in silico MmpL3 inhibitors. Molecular dynamics simulations integrated with MM-PBSA free energy calculations identified that MmpL3-ZINC000000016638 complex was more stable. Study can be further extended for synthesis and biological evaluation, derivatization of active compound to identify potential and safe lead compounds for effective tuberculosis therapy.

Recent Advances in the Applications of Small Molecules in the Treatment of Multiple Myeloma

Therapy for multiple myeloma (MM), a hematologic neoplasm of plasma cells, has undergone remarkable changes over the past 25 years. Small molecules (molecular weight of less than one kDa), together with newer immunotherapies that include monoclonal antibodies, antibody-drug conjugates, and most recently, chimeric antigen receptor (CAR) T-cells, have combined to double the disease's five-year survival rate to over 50% during the past few decades. Despite these advances, the disease is still considered incurable, and its treatment continues to pose substantial challenges, since therapeutic refractoriness and patient relapse are exceedingly common. This review focuses on the current pipeline, along with the contemporary roles and future prospects for small molecules in MM therapy. While small molecules offer prospective benefits in terms of oral bioavailability, cellular penetration, simplicity of preparation, and improved cost-benefit considerations, they also pose problems of toxicity due to off-target effects. Highlighted in the discussion are recent developments in the applications of alkylating agents, immunomodulators, proteasome inhibitors, apoptosis inducers, kinesin spindle protein inhibitors, blockers of nuclear transport, and drugs that affect various kinases involved in intracellular signaling pathways. Molecular and cellular targets are described for each class of agents in relation to their roles as drivers of MM.

The Antiproliferative and Apoptotic Effects of a Novel Quinazoline Carrying Substituted-Sulfonamides: In Vitro and Molecular Docking Study

In order to investigate for a new effective and safe anticancer drug, we synthesized a novel series of quinazoline containing biologically active substituted-sulfonamide moiety at 3- position 4a-n. The structure of the newly prepared compounds was proved by microanalysis, IR, ¹H-NMR, ¹³C-NMR and mass spectral data. All the synthesized compounds were evaluated for their in vitro cytotoxic activity in numerous cancer cell lines including A549, HepG-2, LoVo and MCF-7 and normal HUVEC cell line. The two most active compounds 4d and 4f were then tested for their apoptosis induction using DNA content and Annexin V-FITC/PI staining. Moreover, apoptosis initiation was also confirmed using RT-PCR and Western blot. To further understand the binding preferences of quinazoline sulfonamides, docking simulations were used. Among the fourteen new synthesized compounds, we found that compounds 4d and 4f exerted the strongest cytotoxicity against MCF-7 cells with an IC₅₀ value of 2.5 and 5 μM, respectively. Flow cytometry data revealed the ability of compounds 4d and 4f to mediate apoptosis and arrest cell cycle growth at G1 phase. Furthermore, RT-PCR and Western blot results suggested that both 4d and 4f activates apoptotic cell death pathway in MCF-7 cells. Molecular docking assessments indicated that compounds 4d and 4f fit perfectly into Bcl2's active site. Based on the biological properties, we conclude that both compounds 4d and 4f could be used as a new type of anticancer agent, which provides a scientific basis for further research into the treatment of cancer.

Computational screening of natural compounds from Salvia plebeia R. Br. for inhibition of SARS-CoV-2 main protease

The novel Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV-2) has emerged to be the reason behind the COVID-19 pandemic. It was discovered in Wuhan, China and then began spreading around the world, impacting the health of millions. Efforts for treatment have been hampered as there are no antiviral drugs that are effective against this virus. In the present study, we have explored the phytochemical constituents of Salvia plebeia R. Br., in terms of its binding affinity by targeting COVID-19 main protease (M^pro) using computational analysis. Molecular docking analysis was performed using PyRx software. The ADMET and drug-likeness properties of the top 10 compounds showing binding affinity greater than or equal to - 8.0 kcal/mol were analysed using pkCSM and DruLiTo, respectively. Based on the docking studies, it was confirmed that Rutin and Plebeiosides B were the most potent inhibitors of the main protease of SARS-CoV-2 with the best binding affinities of - 9.1 kcal/mol and - 8.9 kcal/mol, respectively. Further, the two compounds were analysed by studying their biological activity using the PASS webserver. Molecular dynamics simulation analysis was performed for the selected protein-ligand complexes to confirm their stability at 300 ns. MM-PBSA provided the basis for analyzing the affinity of the phytochemicals towards M^pro by calculating the binding energy, and secondary structure analysis indicated the stability of protease structure when it is bound to Rutin and Plebeiosides B. Altogether, the study identifies Rutin and Plebeiosides B to be potent M^pro inhibitors of SARS-CoV-2.

Graphic abstract

Supplementary Information

The online version contains supplementary material available at 10.1007/s42535-021-00304-z.

Studies on ligand-based pharmacophore modeling approach in identifying potent future EGFR inhibitors

Epidermal growth factor receptor (EGFR) is a validated drug target for cancer chemotherapy. Mutations in EGFR are directly linked with the development of drug resistance and this has led for the development of newer drugs in quest for more efficacious inhibitors. The current research is focused on identifying potential and safe molecules as EGFR inhibitors by using both structure and ligand based computational approaches. In quest for finding newer moieties, we have developed a pharmacophore model utilizing drugs like lazertinib, osimertinib, nazartinib, avitinib, afatininb, and talazoparib that are known to inhibit EGFR along with their downstream signaling. Ligand-based pharmacophore model have been developed to screen the ZINC database through ZINCPharmer webserver. The server has identified 9482 best possible ligands with high pharmacophoric similarity i.e., RMSD value less than 0.2 Å. The top 10 ligands with the criteria of dock score(s) and interactions were further subjected to in silico ADMET studies giving two plausible ligands that were further subjected to Molecular Dynamics and MM/PBSA free energy calculations to ensure stability to the target site. Results deduced by in silico work in the current study may be corroborated biologically in the future. The current work, therefore, provides ample opportunity for computational and medicinal chemists to work in allied areas to facilitate the design and development of novel and more efficacious EGFR inhibitors for future experimental studies.

BH3 Mimetics in Hematologic Malignancies

Hematologic malignancies (HM) comprise diverse cancers of lymphoid and myeloid origin, including lymphomas (approx. 40%), chronic lymphocytic leukemia (CLL, approx. 15%), multiple myeloma (MM, approx. 15%), acute myeloid leukemia (AML, approx. 10%), and many other diseases. Despite considerable improvement in treatment options and survival parameters in the new millennium, many patients with HM still develop chemotherapy-refractory diseases and require re-treatment. Because frontline therapies for the majority of HM (except for CLL) are still largely based on classical cytostatics, the relapses are often associated with defects in DNA damage response (DDR) pathways and anti-apoptotic blocks exemplified, respectively, by mutations or deletion of the TP53 tumor suppressor, and overexpression of anti-apoptotic proteins of the B-cell lymphoma 2 (BCL2) family. BCL2 homology 3 (BH3) mimetics represent a novel class of pro-apoptotic anti-cancer agents with a unique mode of action—direct targeting of mitochondria independently of TP53 gene aberrations. Consequently, BH3 mimetics can effectively eliminate even non-dividing malignant cells with adverse molecular cytogenetic alterations. Venetoclax, the nanomolar inhibitor of BCL2 anti-apoptotic protein has been approved for the therapy of CLL and AML. Numerous venetoclax-based combinatorial treatment regimens, next-generation BCL2 inhibitors, and myeloid cell leukemia 1 (MCL1) protein inhibitors, which are another class of BH3 mimetics with promising preclinical results, are currently being tested in several clinical trials in patients with diverse HM. These pivotal trials will soon answer critical questions and concerns about these innovative agents regarding not only their anti-tumor efficacy but also potential side effects, recommended dosages, and the optimal length of therapy as well as identification of reliable biomarkers of sensitivity or resistance. Effective harnessing of the full therapeutic potential of BH3 mimetics is a critical mission as it may directly translate into better management of the aggressive forms of HM and could lead to significantly improved survival parameters and quality of life in patients with urgent medical needs.

Ligand-based design approach of potential Bcl-2 inhibitors for cancer chemotherapy

BACKGROUND AND OBJECTIVES

Overexpression of prosurvival Bcl-2 family members make tumor cells resistant to conventional cancer therapeutic agents. It is commonly observed feature in many different types of human tumors. Hence, small-molecules as Bcl-2 inhibitors may have a promising therapeutic potential for the treatment of human cancer. The given study focusses on development of novel and small Bcl-2 inhibitors using ligand-based drug design approach.

METHODS

Ligand based pharmacophore was generated using the PHASE tool of Schrödinger and screened ZINC database through ZINCPharmer webserver to identify compounds with similar features. Compounds having good fitness score were selected for molecular docking and binding interactions were compared with drugs in market as well as trials. QSAR model was generated using advanced AutoQSAR tool and validated for prediction of unknown compounds. QSAR prediction of in silico active identified three potential compounds and were subjected to investigate stability by molecular dynamics simulations and MM-PBSA binding energy calculations.

RESULTS

Study identified three in silico potential molecules with good stability and binding affinity. Further substructure search and pIC₅₀ value prediction has identified six more molecules. Total nine molecules have demonstrated good drug likeness features.

CONCLUSION

Final oral rat LD50 calculation of nine molecules has identified three hit molecules i.e., ZINC76760927, ZINC76768675 and ZINC52767796 for further in vitro and in vivo testing as safe and potential Bcl-2 inhibitors.

Screening of natural compounds from Cyperus rotundus Linn against SARS-CoV-2 main protease (Mpro): An integrated computational approach

Coronavirus disease 2019 (COVID-19) is a viral respiratory disease that has been spreading across the globe. The World Health Organization (WHO) declared it as a public health emergency. The treatment of COVID-19 has been hampered due to the lack of effective therapeutic efforts. Main Protease (M^pro) is a key enzyme in the viral replication cycle and its non-specificity to human protease makes it a potential drug target. Cyperus rotundus Linn, which belongs to the Cyperaceae family, is a traditional herbal medicine that has been widely studied for its antiviral properties. In this study, a computational approach was used to screen natural compounds from C. rotundus Linn using BIOVIA Discovery Suite and novel potential molecules against M^pro of SARS-CoV-2 were predicted. Molecular docking was performed using LibDock protocol and selected ligands were further subjected to docking analysis by CDOCKER. The docking scores of the selected ligands were compared with standard antiretroviral drugs such as lopinavir and ritonavir to assess their binding potentials. Interaction pharmacophore analysis was then performed for the compounds exhibiting good binding scores to evaluate their protein–ligand interactions. The selected protein–ligand complexes were subjected to molecular dynamics simulation for 50 ns. Results of binding free energy analysis revealed that two compounds—β-amyrin and stigmasta-5,22-dien-3-ol—exhibited the best binding interactions and stability. Finally, absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies were performed to understand the pharmacokinetic properties and safety profile of the compounds. The overall results indicate that the phytochemicals from Cyperus rotundus Linn, namely β-amyrin and stigmasta-5,22-dien-3-ol, can be screened as potential inhibitors of SARS-CoV-2 M^pro.

In silico analysis of Phyllanthus amarus phytochemicals as potent drugs against SARS-CoV-2 main protease

Phyllanthus amarus, also known as Bhui Korma in India, is well known for its medicinal properties and is used to treat several diseases worldwide. This study aims to identify phytochemicals from P. amarus and assess their anti-viral activity through in silico methods against the main protease (3CL^Pro/M^Pro) enzyme of the novel coronavirus. 190 compounds were obtained from literature and docked against 3CLPro and 16 compounds showed higher binding affinity with 3CL^Pro with their values lying between -8.9 ​kcal/mol to -9.6 ​kcal/mol. The top two compounds, Myricitrin (CID: 5352000) and Quercetin-3-O-glucuronide (CID: 12004528) gave high binding affinity values of -9.6 ​kcal/mol and -9.4 ​kcal/mol respectively and also display favourable binding interactions with the 3CL^Pro. Both the compounds were further subjected to molecular dynamics simulation and MM-PBSA based binding free energy calculations. ADMET and drug-likeness properties were studied to assess the pharmacokinetic properties of the compounds. Favourable pharmacokinetic results reinforced the applicability of the compounds assessed. Along with continuous studies being carried out with chemical compounds, research needs to expand into all areas, including the use of natural compounds as drug compounds. The identified hits from this study can be taken further for in vitro and in vivo studies to examine their efficacy against COVID-19.