Drug Targets Research Articles

Adipose Triglyceride Lipase–Mediated Adipocyte Lipolysis Exacerbates Acute Pancreatitis Severity in Mouse Models and Patients

Dyslipolysis of adipocytes has played a critical role in various diseases. Adipose triglyceride lipase (ATGL) is a rate-limiting enzyme in adipocyte autonomous lipolysis. However, whether the degree of adipocyte lipolysis relates to the prognoses in acute pancreatitis (AP) and the role of ATGL mediated lipolysis in the pathogenesis of AP remains elusive. The visceral adipose tissue (VAT) consumption rate (VATR) in the acute stage was measured in both AP patients and mouse models. Lipolysis levels and ATGL expression were detected in caerulein (CER)-induced AP models. CL316,243, a lipolysis stimulator, and adipose tissue-specific ATGL knockout (AAKO) mice were used to further investigate the role of lipolysis in AP. The ATGL-specific inhibitor, Atglistatin, was performed in C57Bl/6N and ob/ob AP models. This study found that increased VATR in the acute phase was independently associated with adverse prognoses in AP patients, which was validated in mice AP models. Lipolysis of adipocytes was elevated in AP mice. Stimulation of lipolysis could aggravate AP. Genetic blockage of ATGL specifically in adipocytes was able to alleviate the damage to AP. The application of Atglistatin could effectively protect against AP in both lean and obese mice. These findings demonstrated that ATGL-mediated adipocyte lipolysis exacerbates AP and highlighted the therapeutic potential of ATGL as a drug target for AP.

Structural basis of human mpox viral DNA replication inhibition by brincidofovir and cidofovir

Mpox virus has wildly spread over 108 non-endemic regions in the world since May 2022. DNA replication of mpox is performed by DNA polymerase machinery F8-A22-E4, which is known as a great drug target. Brincidofovir and cidofovir are reported to have broad-spectrum antiviral activity against poxviruses, including mpox virus in animal models. However, the molecular mechanism is not understood. Here we report cryogenic electron microscopy structures of mpox viral F8-A22-E4 in complex with a DNA duplex, or dCTP and the DNA duplex, or cidofovir diphosphate and the DNA duplex at resolution of 3.22, 2.98 and 2.79 Å, respectively. Our structural work and DNA replication inhibition assays reveal that cidofovir diphosphate is located at the dCTP binding position with a different conformation to compete with dCTP to incorporate into the DNA and inhibit DNA synthesis. Conformation of both F8-A22-E4 and DNA is changed from the pre-dNTP binding state to DNA synthesizing state after dCTP or cidofovir diphosphate is bound, suggesting a coupling mechanism. This work provides the structural basis of DNA synthesis inhibition by brincidofovir and cidofovir, providing a rational strategy for new therapeutical development for mpox virus and other pox viruses.

Cryo-EM structure of I domain-containing integrin αEβ7

The integrin family is a transmembrane receptor that plays critical roles in the cell-cell and cell-extracellular matrix adhesion, signal transduction such as cell cycle regulation, organization of the intracellular cytoskeleton, and immune responses. Consequently, dysfunction of integrins is associated with a wide range of human diseases, including cancer and immune diseases, which makes integrins therapeutic targets for drug discovery. Here we report the cryo-EM structure of the α-I domain-containing full-length integrin αEβ7, which is expressed in the leukocytes of the immune system and a drug target for inflammatory bowel disease (IBD). The structure reveals the half-bent conformation, an intermediate between the close and the open conformation, while the α-I domain responsible for the ligand binding covers the headpiece domain by a unique spatial arrangement. Our results provide the structural information for the drug design targeting IBD.

Artemisinin derivatives as potential drug candidates against Mycobacterium tuberculosis: insights from molecular docking, MD simulations, PCA, MM/GBSA and ADMET analysis

ABSTRACT Tuberculosis, one of the most ancient and formidable infectious diseases, primarily arises from the pathogenic bacterium Mycobacterium tuberculosis (Mtb). Phosphoenolpyruvate carboxykinase (Pck), a potential drug target, is essential for the growth of Mtb that involves in the pathway of gluconeogenesis at the centre of phosphoenolpyruvate-pyruvate-oxaloacetate node. This study aims to understand drug-like properties of the compounds derived from artemisinin, and to investigate their inhibitory roles against Mtb. Molecular docking was performed on a set of 56 artemisinin compounds to identify their binding efficacy to the Mtb protein target. Subsequently, the selected top three complexes were subjected to molecular dynamics simulations. By the trajectory analysis, root mean square deviation and radius of gyration indicated the compactness of the systems without much fluctuation. The principal component analysis revealed that the complexes were less dynamic and energetically more favourable, and Gibbs energy landscape revealed the favourable energetic transitions between conformations. The artemisinin dimer primary alcohol holding good bioavailability scores appears highly stable in complex with Pck (MM/GBSA of −37.67 kcal/mol) that stands as a potential inhibitor of the Mtb target. However, further preclinical experiments and investigations are necessary to evaluate the intrinsic properties and to confirm effectiveness of the drug candidate.

Identifying Potential Alzheimer's Disease Therapeutics through GSK-3β Inhibition: A Molecular Docking and Dynamics Approach

Emerging as a promising drug target for Alzheimer's disease (AD) therapy, glycogen synthase kinase 3β (GSK-3β) has garnered attention. This study sought to rigorously scrutinize a compendium of natural compounds retrieved from the ZINC database through pharmacodynamic experiments, employing a 1H-indazole-3-carboxamide (INDZ) scaffold, to identify compounds capable of inhibiting the GSK-3β protein. Utilizing a multi-step approach, the study involved pharmacophore analysis, followed by molecular docking to select five promising ligands for further investigation. Subsequently, ESMACS simulations were employed to assess the stability of the ligand-protein interactions. Evaluation of the binding modes and free energy of the ligands revealed that five compounds (2a-6a) exhibited crucial interactions with the active site residues. Furthermore, various methodologies, including hydrogen bond and clustering analyses, were utilized to ascertain their inhibitory potential and elucidate the factors contributing to ligand binding in the protein's active site. The findings from MMPBSA/GBSA analysis indicated that these five selected small molecules closely approached the IC50 value of the reference ligand (OH8), yielding energy values of −34.85, −32.58, −31.71, and −30.39kcal/mol, respectively. Additionally, an assessment of the interactions using hydrogen bond and dynamic analyses delineated the effective binding of the ligands with the binding pockets in the protein. Through computational analysis, we obtained valuable insights into the molecular mechanisms of GSK-3β, aiding in the development of more potent inhibitors.

Effect of luteolin on glioblastoma's immune microenvironment and tumor growth suppression

BackgroundGlioblastoma is the most malignant and prevalent primary human brain tumor, and the immunological microenvironment controlled by glioma stem cells is one of the essential elements contributing to its malignancy. The use of medications to ameliorate the tumor microenvironment may give a new approach for glioma treatment. MethodsGlioma stem cells were separated from clinical patient-derived glioma samples for molecular research. Other studies, including CCK8, EdU, Transwell, and others, supported luteolin's ability to treat glioma progenitor cells. Network pharmacology and molecular docking models were used to study the drug target, and qRT-PCR, WB, and IF were used to evaluate the molecular mechanism. Intracranial xenografts were examined using HE and IHC, while macrophage polarization was examined using FC. ResultsWe originally discovered that luteolin inhibits glioma stem cells. IL6 released by glioma stem cells is blocked during medication action and inhibits glioma stem cell proliferation and invasion via the IL6/STAT3 signaling pathway. Additionally, luteolin inhibits the secretion of TGFβ1, affects the polarization function of macrophages in the microenvironment, inhibits the polarization of M2 macrophages in TAM, and further inhibits various functions of glioma stem cells by affecting the IL6/STAT3 signaling pathway, luteolin crosstalk TGFβ1/SMAD3 signaling pathway, and so on. ConclusionsThrough the suppression of the immunological microenvironment and inhibition of the IL6/STAT3 signaling pathway, our study determined the inhibitory effect of luteolin on glioma stem cells. This medication's dual inhibitory action, which has a significant negative impact on the glioma stem cells' malignant process, makes it both a viable anti-glioma medication and a candidate for targeted glioma microenvironment therapy.

Pentacyclic triterpenoids as potential ACL inhibitors from the rare medicinal plant Semiliquidambar cathayensis

An extensive phytochemical investigation on the rare medicinal plant Semiliquidambar cathayensis (family: Hamamelidaceae) led to the isolation of four new (1–4, named semiliquidacids A–D, respectively) and 25 related known pentacyclic triterpenoids. The new structures with absolute configurations were elucidated by spectroscopic methods, electronic circular dichroism (ECD) calculations, and single-crystal X-ray diffraction analysis. Compound 1 represents the first naturally occurring ursane-type triterpenoid featuring an uncommon C-25 formyl group. Compound 4 and oleanolic acid (13) exhibited remarkable inhibitory effects against the ATP-citrate lyase (ACL, an emerging drug target for hyperlipidemia and related metabolic disorders) with IC50 values of 6.5 and 11.9 μM, respectively. The molecular interaction and binding mode between the bioactive triterpenoids and ACL were elaborated by conducting a molecular docking study. Meanwhile, the chemotaxonomic significance of the isolated triterpenoids has been briefly discussed.

O-Mannosyltransferase CfPmt4 Regulates the Growth, Development and Pathogenicity of Colletotrichum fructicola

Camellia oleifera is a woody, edible-oil plant native to China. Anthracnose is the major disease of Ca. oleifera, and Colletotrichum fructicola is the main epidemic pathogen. Our previous research indicated that CfHac1 (homologous to ATF/CREB1) and CfGcn5 (general control nonderepressible 5, Gcn5) are integral to key cellular processes that govern fungal development and pathogenesis. Further transcriptomic analyses of the CfHac1 and CfGcn5 mutants, particularly under conditions of endoplasmic reticulum (ER) stress, hold the potential to unveil additional genes implicated in this critical cellular response. We identified all OST/PMT (oligosaccharyltransferase/Protein O-Mannosyltransferases) genes in C. fructicola and analyzed their expression levels. To elucidate novel glycosylation-related genes that may be important for the virulence of C. fructicola, we took an unbiased transcriptomic approach comparing wild-type and the ∆Cfhac1 mutant. Notably, all OST/PMT genes were induced by dithiothreitol and down-regulated in the ΔCfhac1 mutant, yet only the CfPMT4 (Protein O-Mannosyltransferases 4) gene (A04626) was unaffected in the ΔCfgcn5. The results of targeted gene deletion experiments indicate that CfPMT4 plays a crucial role in both vegetative growth and conidiation. Additionally, our investigation revealed that the ΔCfpmt4 exhibits deficiencies in appressorium formation, as well as in its response to cell wall integrity and endoplasmic reticulum stresses. Furthermore, the mutant displayed impaired glycogen metabolism, which may contribute to reduced penetration ability. Overall, CfPmt4, an O-mannosyltransferase, controls the growth, development, and pathogenicity of Colletotrichum fructicola. Understanding the function of the CfPMT4 homolog could provide a potential molecular target for controlling Ca. oleifera anthracnose.

Epicardial SMARCA4 deletion exacerbates cardiac injury in myocardial infarction and is related to the inhibition of epicardial epithelial-mesenchymal transition

The reactivated adult epicardium produces epicardium-derived cells (EPDCs) via epithelial-mesenchymal transition (EMT) to benefit the recovery of the heart after myocardial infarction (MI). SMARCA4 is the core catalytic subunit of the chromatin re-modeling complex, which has the potential to target some reactivated epicardial genes in MI. However, the effects of epicardial SMARCA4 on MI remain uncertain. This study found that SMARCA4 was activated over time in epicardial cells following MI, and some of activated cells belonged to downstream differentiation types of EPDCs. This study used tamoxifen to induce lineage tracing and SMARCA4 deletion from epicardial cells in Wt1-CreER;Smarca4fl/fl;Rosa26-RFP adult mice. Epicardial SMARCA4 deletion reduces the number of epicardial cells in adult mice, which was related to changes in the activation, proliferation, and apoptosis of epicardial cells. Epicardial SMARCA4 deletion reduced collagen deposition and angiogenesis in the infarcted area, exacerbated cardiac injury in MI. The exacerbation of cardiac injury was related to the inhibition of generation and differentiation of EPDCs. The alterations in EPDCs were associated with inhibited transition between E-CAD and N-CAD during the epicardial EMT, coupled with the down-regulation of WT1, SNAIL1, and PDGF signaling. In conclusion, this study suggests that Epicardial SMARCA4 plays a critical role in cardiac injury caused by MI, and its regulatory mechanism is related to epicardial EMT. Epicardial SMARCA4 holds potential as a novel molecular target for treating MI.

Revisiting the Role of B-RAF Kinase as a Therapeutic Target in Melanoma.

Malignant melanoma is the rarest but most aggressive and deadly skin cancer. Melanoma is the result of a malignant transformation of melanocytes, which leads to their uncontrolled proliferation. Mutations in the mitogen-activated protein kinase (MAPK) pathway, which are crucial for the control of cellular processes, such as apoptosis, division, growth, differentiation, and migration, are one of its most common causes. BRAF kinase, as one of the known targets of this pathway, has been known for many years as a prominent molecular target in melanoma therapy, and the following mini-review outlines the state-of-the-art knowledge regarding its structure, mutations and mechanisms.

DTNPD: A comprehensive database of drugs and targets for neurological and psychiatric disorders

In response to the shortcomings in data quality and coverage for neurological and psychiatric disorders (NPDs) in existing comprehensive databases, this paper introduces the DTNPD database, specifically designed for NPDs. DTNPD contains detailed information on 30 NPDs types, 1847 drugs, 514 drug targets, 64 drug combinations, and 61 potential target combinations, forming a network with 2389 drug-target associations. The database is user-friendly, offering open access and downloadable data, which is crucial for network pharmacology studies. The key strength of DTNPD lies in its robust networks of drug and target combinations, as well as drug-target networks, facilitating research and development in the field of NPDs. The development of the DTNPD database marks a significant milestone in understanding and treating NPDs. For accessing the DTNPD database, the primary URL is http://dtnpd.cnsdrug.com, complemented by a mirror site available at http://dtnpd.lyhbio.com.

E2F transcription factors promote tumorigenicity in pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is one of the most lethal cancers with limited treatment options, illustrating an urgent need to identify new drugable targets in PDACs. Using the similarities between tumor development and normal embryonic development, which is accompanied by rapid cell expansion, we aimed to identify and characterize embryonic signaling pathways that were reinitiated during tumor formation and expansion. Here, we report that the transcription factors E2F1 and E2F8 are potential key regulators in PDAC. E2F1 and E2F8 RNA expression is mainly localized in proliferating cells in the developing pancreas and in malignant ductal cells in PDAC. Silencing of E2F1 and E2F8 in PANC-1 pancreatic tumor cells inhibited cell proliferation and impaired cell spreading and migration. Moreover, loss of E2F1 also affected cell viability and apoptosis with E2F expression in PDAC tissues correlating with expression of apoptosis and mitosis pathway genes, suggesting that E2F factors promote cell cycle regulation and tumorigenesis in PDAC cells. Our findings illustrate that E2F1 and E2F8 transcription factors are expressed in pancreatic progenitor and PDAC cells, where they contribute to tumor cell expansion by regulation of cell proliferation, viability, and cell migration making these genes attractive therapeutic targets and potential prognostic markers for pancreatic cancer.

Drug-target Mendelian randomization analysis supports lowering plasma ANGPTL3, ANGPTL4, and APOC3 levels as strategies for reducing cardiovascular disease risk

Abstract Background and Aims APOC3, ANGPTL3, and ANGPTL4 are circulating proteins that are actively pursued as pharmacological targets to treat dyslipidemia and reduce the risk of atherosclerotic cardiovascular disease. Here, we used human genetic data to compare the predicted therapeutic and adverse effects of APOC3, ANGPTL3, and ANGPTL4 inactivation. Methods We conducted drug-target Mendelian randomization analyses using variants in proximity to the genes associated with circulating protein levels to compare APOC3, ANGPTL3, and ANGPTL4 as drug targets. We obtained exposure and outcome data from large-scale genome-wide association studies and used generalized least squares to correct for linkage disequilibrium-related correlation. We evaluated five primary cardiometabolic endpoints and screened for potential side effects across 694 disease-related endpoints, 43 clinical laboratory tests, and 11 internal organ MRI measurements. Results Genetically lowering circulating ANGPTL4 levels reduced the odds of coronary artery disease (CAD) (odds ratio, 0.57 per s.d. protein [95%CI,0.47–0.70]) and type 2 diabetes (T2D) (odds ratio, 0.73 per s.d. protein [95%CI,0.57–0.94]). Genetically lowering circulating APOC3 levels also reduced the odds of CAD (odds ratio, 0.90 per s.d. protein [95%CI,0.82–0.99]). Genetically lowered ANGPTL3 levels via common variants were not associated with CAD. However, meta-analysis of deleterious variants revealed that ANGPTL3 inactivation protected against CAD (odds ratio, 0.79 per allele [95%CI, 0.69–0.90]). Analysis of lowered ANGPTL3, ANGPTL4, and APOC3 levels did not identify important safety concerns. Conclusion Human genetic evidence suggests that therapies aimed at reducing circulating levels of ANGPTL3, ANGPTL4, and APOC3 reduce the risk of CAD. ANGPTL4 lowering may also reduce the risk of T2D.

Structure-Function Analysis of the Essential Mycobacterium tuberculosis P450 Drug Target, CYP121A1.

Cytochrome P450 CYP121A1 is a well-known drug target against Mycobacterium tuberculosis, the human pathogen that causes the deadly disease tuberculosis (TB). CYP121A1 is a unique P450 enzyme because it uses classical and non-classical P450 catalytic processes and has distinct structural features among P450s. However, a detailed investigation of CYP121A1 protein structures in terms of active site cavity dynamics and key amino acids interacting with bound ligands has yet to be undertaken. To address this research knowledge gap, 53 CYP121A1 crystal structures were investigated in this study. Critical amino acids required for CYP121A1's overall activity were identified and highlighted this enzyme's rigid architecture and substrate selectivity. The CYP121A1-fluconazole crystal structure revealed a novel azole drug-P450 binding mode in which azole heme coordination was facilitated by a water molecule. Fragment-based inhibitor approaches revealed that CYP121A1 can be inhibited by molecules that block the substrate channel or by directly interacting with the P450 heme. This study serves as a reference for the precise understanding of CYP121A1 interactions with different ligands and the structure-function analysis of P450 enzymes in general. Our findings provide critical information for the synthesis of more specific CYP121A1 inhibitors and their development as novel anti-TB drugs.

Biophysical profiling and in vitro study of apocarotenoids with calf thymus DNA to explore their efficacy in adjuvant cancer therapy

DNA is one of the important targets for therapeutic drugs. Investigating the interaction between these drug molecules and DNA yields valuable insights for the further advancement of medication formulations. In this study, we explore the interaction between calf thymus-DNA and temozolomide with apocarotenoids such as bixin and crocin, which are known for their therapeutic activities. The combinatory effect of the compounds on ct-DNA was studied through diverse spectroscopic methods including UV-visible, circular dichroism studies coupled with viscosity and cell cycle analysis. UV-visible spectroscopy evaluation indicates that apocarotenoids and the temozolomide combination display a non-intercalative binding mechanism. The circular dichroism study revealed that the interaction between DNA and apocarotenoids, as well as their combinations with temozolomide, did not present any significant alterations in the structure of DNA suggesting a possible groove binding interaction. The viscosity analyses verified that ct-DNA interacts with apocarotenoids along with their temozolomide combination through a groove-binding mechanism. This study offers a more comprehensive understanding of the potential mechanism of interaction between ct-DNA and natural compounds with the addition of temozolomide.

Drug Screening and Target Prediction Based on Machine Learning

Drug Screening and Target Prediction Based on Machine Learning

QUBO Problem Formulation of Fragment-Based Protein–Ligand Flexible Docking

Protein–ligand docking plays a significant role in structure-based drug discovery. This methodology aims to estimate the binding mode and binding free energy between the drug-targeted protein and candidate chemical compounds, utilizing protein tertiary structure information. Reformulation of this docking as a quadratic unconstrained binary optimization (QUBO) problem to obtain solutions via quantum annealing has been attempted. However, previous studies did not consider the internal degrees of freedom of the compound that is mandatory and essential. In this study, we formulated fragment-based protein–ligand flexible docking, considering the internal degrees of freedom of the compound by focusing on fragments (rigid chemical substructures of compounds) as a QUBO problem. We introduced four factors essential for fragment–based docking in the Hamiltonian: (1) interaction energy between the target protein and each fragment, (2) clashes between fragments, (3) covalent bonds between fragments, and (4) the constraint that each fragment of the compound is selected for a single placement. We also implemented a proof-of-concept system and conducted redocking for the protein–compound complex structure of Aldose reductase (a drug target protein) using SQBM+, which is a simulated quantum annealer. The predicted binding pose reconstructed from the best solution was near-native (RMSD=1.26 Å), which can be further improved (RMSD=0.27 Å) using conventional energy minimization. The results indicate the validity of our QUBO problem formulation.

Similar excitability through different sodium channels and implications for the analgesic efficacy of selective drugs

Nociceptive sensory neurons convey pain-related signals to the CNS using action potentials. Loss-of-function mutations in the voltage-gated sodium channel NaV1.7 cause insensitivity to pain (presumably by reducing nociceptor excitability) but clinical trials seeking to treat pain by inhibiting NaV1.7 pharmacologically have struggled. This may reflect the variable contribution of NaV1.7 to nociceptor excitability. Contrary to claims that NaV1.7 is necessary for nociceptors to initiate action potentials, we show that nociceptors can achieve similar excitability using different combinations of NaV1.3, NaV1.7, and NaV1.8. Selectively blocking one of those NaV subtypes reduces nociceptor excitability only if the other subtypes are weakly expressed. For example, excitability relies on NaV1.8 in acutely dissociated nociceptors but responsibility shifts to NaV1.7 and NaV1.3 by the fourth day in culture. A similar shift in NaV dependence occurs in vivo after inflammation, impacting ability of the NaV1.7-selective inhibitor PF-05089771 to reduce pain in behavioral tests. Flexible use of different NaV subtypes exemplifies degeneracy – achieving similar function using different components – and compromises reliable modulation of nociceptor excitability by subtype-selective inhibitors. Identifying the dominant NaV subtype to predict drug efficacy is not trivial. Degeneracy at the cellular level must be considered when choosing drug targets at the molecular level.

Exploring the Antiviral Potential of Natural Compounds against Influenza: A Combined Computational and Experimental Approach.

The influenza A virus nonstructural protein 1 (NS1), which is crucial for viral replication and immune evasion, has been identified as a significant drug target with substantial potential to contribute to the fight against influenza. The emergence of drug-resistant influenza A virus strains highlights the urgent need for novel therapeutics. This study proposes a combined theoretical criterion for the virtual screening of molecular libraries to identify candidate NS1 inhibitors. By applying the criterion to the ZINC Natural Product database, followed by ligand-based virtual screening and molecular docking, we proposed the most promising candidate as a potential NS1 inhibitor. Subsequently, the selected natural compound was experimentally evaluated, revealing measurable virus replication inhibition activity in cell culture. This approach offers a promising avenue for developing novel anti-influenza agents targeting the NS1 protein.

Drug repurposing: identification of SARS-CoV-2 potential inhibitors by virtual screening and pharmacokinetics strategies.

The coronavirus disease 2019 (COVID-19) pandemic caused global health, economic, and population loss. Variants of the coronavirus contributed to the severity of the disease and persistent rise in infections. This study aimed to identify potential drug candidates from fifteen approved antiviral drugs against SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike protein (6M0J) using virtual screening and pharmacokinetics to gain insights into COVID-19 therapeutics. We employed drug repurposing approach to analyze binding performance of fifteen clinically approved antiviral drugs against the main protease of SARS-CoV-2 (6LU7), SARS-CoV (5B6O), and SARS-CoV-2 spike proteins bound to ACE-2 receptor (6M0J), to provide an insight into the therapeutics of COVID-19. AutoDock Vina was used for docking studies. The binding affinities were calculated, and 2-3D structures of protein-ligand interactions were drawn. Rutin, hesperidin, and nelfinavir are clinically approved antiviral drugs with high binding affinity to proteins 6LU7, 5B6O, and 6M0J. These ligands have excellent pharmacokinetics, ensuring efficient absorption, metabolism, excretion, and digestibility. Hesperidin showed the most potent interaction with spike protein 6M0J, forming four H-bonds. Nelfinavir had a high human intestinal absorption (HIA) score of 0.93, indicating maximum absorption in the body and promising interactions with 6LU7. Our results indicated that rutin, hesperidin, and nelfinavir had the highest binding results against the proposed drug targets. The computational approach effectively identified SARS-CoV-2 inhibitors. COVID-19 is still a recurrent threat globally and predictive analysis using natural compounds might serve as a starting point for new drug development against SARS-CoV-2 and related viruses.