Transformer-Based Molecular Generative Model for Antiviral Drug Design.

Since the Simplified Molecular Input Line Entry System (SMILES) is oriented to the atomic-level representation of molecules and is not friendly in terms of human readability and editable, however, IUPAC is the closest to natural language and is very friendly in terms of human-oriented readability and performing molecular editing, we can manipulate IUPAC to generate corresponding new molecules and produce programming-friendly molecular forms of SMILES. In addition, antiviral drug design, especially analogue-based drug design, is also more appropriate to edit and design directly from the functional group level of IUPAC than from the atomic level of SMILES, since designing analogues involves altering the R group only, which is closer to the knowledge-based molecular design of a chemist. Herein, we present a novel data-driven self-supervised pretraining generative model called "TransAntivirus" to make select-and-replace edits and convert organic molecules into the desired properties for design of antiviral candidate analogues. The results indicated that TransAntivirus is significantly superior to the control models in terms of novelty, validity, uniqueness, and diversity. TransAntivirus showed excellent performance in the design and optimization of nucleoside and non-nucleoside analogues by chemical space analysis and property prediction analysis. Furthermore, to validate the applicability of TransAntivirus in the design of antiviral drugs, we conducted two case studies on the design of nucleoside analogues and non-nucleoside analogues and screened four candidate lead compounds against anticoronavirus disease (COVID-19). Finally, we recommend this framework for accelerating antiviral drug discovery.

Structure based virtual screening and molecular simulation study of FDA-approved drugs to inhibit human HDAC6 and VISTA as dual cancer immunotherapy.

Cancer immunotherapy has significantly contributed to the treatment of various types of cancers mainly by targeting immune checkpoint inhibitors (ICI). Among them, V-domain immunoglobulin suppressor of T cell activation (VISTA) has been explored as a promising therapeutic target. Besides, histone deacetylase 6 (HDAC6) has been demonstrated to be efficacious target for several cancers. The current theoretical work was performed to explore the virtual repurposing of the FDA-approved drugs as inhibitors against these two (VISTA and HDAC6) cancers therapeutic targets. The crystal structure of the two proteins were downloaded from PDB and subjected to virtual screening by DrugRep webserver while using FDA-approved drugs library as ligands database. Our study revealed that Oxymorphone and Bexarotene are the top-ranked inhibitors of VISTA and HDAC6, respectively. The docking score of Bexarotene was predicted as - 10 kcal/mol while the docking score of Oxymorphone was predicted as - 6.2 kcal/mol. Furthermore, a total of 100 ns MD simulation revealed that the two drugs Oxymorphone and Bexarotene formed stable complexes with VISTA and HDAC6 drug targets. As compared to the standard drug the two drugs Oxymorphone and Bexarotene revealed great stability during the whole 100 ns MD simulation. The binding free energy calculation further supported the Root Mean Square Deviation (RMSD) result which stated that as compared to the ref/HDAC6 (- 18.0253 ± 2.6218) the binding free energy score of the Bexarotene/HDAC6 was good (- 51.9698 ± 3.1572 kcal/mol). The binding free energy score of Oxymorphone/VISTA and Ref/VISTA were calculated as - 36.8323 ± 3.4565, and - 21.5611 ± 4.8581 respectively. In conclusion, the two drugs deserve further consideration as cancer treatment option.

Corrigendum: Exploring the molecular mechanism of hepatitis virus inducing hepatocellular carcinoma by microarray data and immune infiltrates analysis.

[This corrects the article DOI: 10.3389/fimmu.2022.1032819.].

Application of a deep generative model produces novel and diverse functional peptides against microbial resistance.

Antimicrobial resistance could threaten millions of lives in the immediate future. Antimicrobial peptides (AMPs) are an alternative to conventional antibiotics practice against infectious diseases. Despite the potential contribution of AMPs to the antibiotic's world, their development and optimization have encountered serious challenges. Cutting-edge methods with novel and improved selectivity toward resistant targets must be established to create AMPs-driven treatments. Here, we present AMPTrans-lstm, a deep generative network-based approach for the rational design of AMPs. The AMPTrans-lstm pipeline involves pre-training, transfer learning, and module identification. The AMPTrans-lstm model has two sub-models, namely, (long short-term memory) LSTM sampler and Transformer converter, which can be connected in series to make full use of the stability of LSTM and the novelty of Transformer model. These elements could generate AMPs candidates, which can then be tailored for specific applications. By analyzing the generated sequence and trained AMPs, we prove that AMPTrans-lstm can expand the design space of the trained AMPs and produce reasonable and brand-new AMPs sequences. AMPTrans-lstm can generate functional peptides for antimicrobial resistance with good novelty and diversity, so it is an efficient AMPs design tool.

Exploring the molecular mechanism of hepatitis virus inducing hepatocellular carcinoma by microarray data and immune infiltrates analysis.

The number of new cases of hepatocellular carcinoma (HCC) worldwide reached 910,000, ranking the sixth, 80% HCC is associated with viruses, so exploring the molecular mechanism of viral carcinogenicity is imperative. The study showed that both HBV and HCV associated HCC and non-viral HCC have the same molecular phenotype (low gene expression and inhibition of immune pathways), but in the tumor immune micro-environment, there is excessive M2-type macrophage polarization in virus-associated hepatocellular carcinoma. To address this phenomenon, the data sets were analyzed and identified five hub genes (POLR2A, POLR2B, RPL5, RPS6, RPL23A) involved in viral gene expression and associated with PI3K-Akt-mTOR pathway activation by six algorithms. In addition, numerous studies have reported that M2-type macrophages participate in the hepatic fibro-pathological process of the development of HCC and are regulated by the PI3K-Akt-mTOR pathway. On this basis, the study showed that hepatitis virus causes abnormal expression of hub genes, leading to the activation of the pathway, which in turn promote the differentiation of M2-type macrophages and eventually promote the formation of liver fibrosis, leading to the occurrence of HCC. In addition, these hub genes are regulated by transcription factors and m6A enzyme, and have good prognosis and diagnostic value. With regard to drug reuse, the results suggest that patients with virus-related HCC for whom Cytidine triphosphate disodium salt and Guanosine-5'-Triphosphate are used as supplementary therapy, and may have a better prognosis. In conclusion, the study has identified novel molecules that are carcinogenic to hepatitis viruses and are expected to serve as molecular markers and targets for diagnosis and treatment.

Bergamottin and PAP-1 Induced ACE2 Degradation to Alleviate Infection of SARS-CoV-2.

Angiotensin-converting enzyme 2 (ACE2), a functional receptor for SARS-CoV, now appears likely to mediate 2019-nCoV entry into human cells. However, inhibitors such as PAP-1 and bergamottin have been discovered; both of them can preferentially bind to ACE2, prevent RBD Spike S protein from binding to ACE2, and reduce the binding sites for RBD Spike S protein. In addition, we investigated the binding energy of PAP-1 and bergamottin with ACE2 through molecular docking with bio-layer interferometry (BLI) and found relatively high binding affinity (KD = 48.5 nM, 53.1 nM) between the PAP-1 and bergamottin groups. In addition, the nanomolar fraction had no effect on growth of the AT-II cell, but 150 µM PAP-1 and 75 µM bergamottin inhibited the proliferation of AT-II cells in vitro by 75% and 68%, respectively. Meanwhile, they significantly reduced ACE2 mRNA and proteins by 67%, 58% and 55%, 41%, respectively. These results indicate that psoralen compounds PAP-1 and bergamottin binding to ACE2 protein could be further developed in the fight against COVID-19 infection during the current pandemic. However, attention should be paid to the damage to human alveolar type II epithelial cells.

Development of an innovative data-driven system to generate descriptive prediction equation of dielectric constant on small sample sets.

Dielectric constant (DC, ε) is a fundamental parameter in material sciences to measure polarizability of the system. In industrial processes, its value is an imperative indicator, which demonstrates the dielectric property of material and compiles information including separation information, chemical equilibrium, chemical reactivity analysis, and solubility modeling. Since, the available ε-prediction models are fairly primitive and frequently suffer from serious failures especially when deals with strong polar compounds. Therefore, we have developed a novel data-driven system to improve the efficiency and wide-range applicability of ε using in material sciences. This innovative scheme adopts the correlation distance and genetic algorithm to discriminate features' combination and avoid overfitting. Herein, the prediction output of the single ML model as a coding to estimate the target value by simulating the layer-by-layer extraction in deep learning, and enabling instant search for the optimal combination of features is recruited. Our model established an improved correlation value of 0.956 with target as compared to the previously available best traditional ML result of 0.877. Our framework established a profound improvement, especially for material systems possessing ε value >50. In terms of interpretability, we have derived a conceptual computational equation from a minimum generating tree. Our innovative data-driven system is preferentially superior over other methods due to its application for the prediction of dielectric constants as well as for the prediction of overall micro and macro-properties of any multi-components complex.

3D-QSAR-Based Pharmacophore Modeling, Virtual Screening, and Molecular Dynamics Simulations for the Identification of Spleen Tyrosine Kinase Inhibitors.

Spleen tyrosine kinase (SYK) is an essential mediator of immune cell signaling and has been anticipated as a therapeutic target for autoimmune diseases, notably rheumatoid arthritis, allergic rhinitis, asthma, and cancers. Significant attempts have been undertaken in recent years to develop SYK inhibitors; however, limited success has been achieved due to poor pharmacokinetics and adverse effects of inhibitors. The primary goal of this research was to identify potential inhibitors having high affinity, selectivity based on key molecular interactions, and good drug-like properties than the available inhibitor, fostamatinib. In this study, a 3D-QSAR model was built for SYK based on known inhibitor IC50 values. The best pharmacophore model was then used as a 3D query to screen a drug-like database to retrieve hits with novel chemical scaffolds. The obtained compounds were subjected to binding affinity prediction using the molecular docking approach, and the results were subsequently validated using molecular dynamics (MD) simulations. The simulated compounds were ranked according to binding free energy (ΔG), and the binding affinity was compared with fostamatinib. The binding mode analysis of selected compounds revealed that the hit compounds form hydrogen bond interactions with hinge region residue Ala451, glycine-rich loop residue Lys375, Ser379, and DFG motif Asp512. Identified hits were also observed to form a desirable interaction with Pro455 and Asn457, the rare feature observed in SYK inhibitors. Therefore, we argue that identified hit compounds ZINC98363745, ZINC98365358, ZINC98364133, and ZINC08789982 may help in drug design against SYK.

De novo molecular design with deep molecular generative models for PPI inhibitors.

We construct a protein-protein interaction (PPI) targeted drug-likeness dataset and propose a deep molecular generative framework to generate novel drug-likeness molecules from the features of the seed compounds. This framework gains inspiration from published molecular generative models, uses the key features associated with PPI inhibitors as input and develops deep molecular generative models for de novo molecular design of PPI inhibitors. For the first time, quantitative estimation index for compounds targeting PPI was applied to the evaluation of the molecular generation model for de novo design of PPI-targeted compounds. Our results estimated that the generated molecules had better PPI-targeted drug-likeness and drug-likeness. Additionally, our model also exhibits comparable performance to other several state-of-the-art molecule generation models. The generated molecules share chemical space with iPPI-DB inhibitors as demonstrated by chemical space analysis. The peptide characterization-oriented design of PPI inhibitors and the ligand-based design of PPI inhibitors are explored. Finally, we recommend that this framework will be an important step forward for the de novo design of PPI-targeted therapeutics.

Isolation and Characterization of an Acidic, Salt-Tolerant Endoglucanase Cel5A from a Bacterial Strain Martelella endophytica YC6887 Genome.

A Martelella endophytica (M. endophytica) strain YC6887 was previously isolated from the roots of a halophyte, Rosa rugosa, which was sequenced and characterized. The genomic and proteomic analysis showed a carbohydrate-degrading enzyme, endoglucanase Cel5A which was further characterized. The protein analysis revealed that this endoglucanase belongs to glycosidic hydrolase family 5 (GH5) with catalytic domain. This gene encodes 349-residue polypeptide and shows closest similarity with cellulases of other Martelella species. The protein was purified to homogeneity and shown that it was a 39 kDa protein. The purified recombinant Cel5A endoglucanase exhibited maximum activity at 50 °C and pH 4.5. The enzyme was salt tolerant and retained more than 50% residual activity up to 15% NaCl. The homology model structure of Cel5A displayed that it is stable and compact protein structure consisting of eleven α-helical structures and eight ß-sheets. According to the predicted ligand binding site after superimposition with Pseudomonas stutzeri endoglucanase Cel5A (PDB ID: 4LX4), it consisted of five amino acid Asn157, Tyr116, Glu158, Glu270 and Trp303 that may be the expected active site of Cel5A from YC6887. This presented that our strain M. endophytica YC6887 that produces cellulase partially degrade the insoluble polysaccharides into reducing sugars.

Synthesis and Biological Evaluation of Amoxicillin Loaded Hybrid Material Composite Spheres Against Methicillin-Resistant Staphylococcus aureus.

BACKGROUND: Incoherent use of antibiotics has led toward resistance in MRSA, becoming multidrug-resistant with a high rate of virulence in the community and hospital settings. OBJECTIVE: Synergistic anti-MRSA activity was investigated in this study for hybrid material composite spheres of amoxicillin, Ag nanoparticles, and chitosan, which were prepared by one-step synthesis method, and various characterizations were performed. METHODS: Antimicrobial-susceptibility assay on MRSA was achieved by disc diffusion and agar dilution techniques, while agar well diffusion was used for hybrid composite spheres. The in vitro and cytotoxicity studies were conducted on the skin abrasion mouse model and MTT assay on RD cell, respectively. RESULTS: All isolates showed resistance to the tested antibiotics except vancomycin. MIC against MRSA showed high resistance with amoxicillin from 4 to 128 mg L-1. The mean diameter of chitosan spheres and Ag nanoparticles was 02 mm and 277 nm, respectively. Morphology of spheres was uneven, varied, porous, and irregular in SEM, and Ag nanoparticles presence and formation was also seen in the micrograph. No substantial interface among drug, nanoparticles, and polymer was found in XRD, and IR showed characteristic peaks of all compounds in the formulation. The in vitro assay showed augmented anti-MRSA activity with amoxicillin loaded hybrid composite spheres (22-29 mm). A significant reduction in microbial burden (~6.5 log10 CFU mL-1) was seen in vivo with loaded hybrid composite spheres formulation. The MTT assay indicated no potential cytotoxicity with hybrid composite spheres. CONCLUSION: The synergistic effect of Amoxycillin in the current study predicts a promising hybrid formulation with enhanced anti-MRSA activity.

In Silico Study Identified Methotrexate Analog as Potential Inhibitor of Drug Resistant Human Dihydrofolate Reductase for Cancer Therapeutics.

Drug resistance is a core issue in cancer chemotherapy. A known folate antagonist, methotrexate (MTX) inhibits human dihydrofolate reductase (hDHFR), the enzyme responsible for the catalysis of 7,8-dihydrofolate reduction to 5,6,7,8-tetrahydrofolate, in biosynthesis and cell proliferation. Structural change in the DHFR enzyme is a significant cause of resistance and the subsequent loss of MTX. In the current study, wild type hDHFR and double mutant (engineered variant) F31R/Q35E (PDB ID: 3EIG) were subject to computational study. Structure-based pharmacophore modeling was carried out for wild type (WT) and mutant (MT) (variant F31R/Q35E) hDHFR structures by generating ten models for each. Two pharmacophore models, WT-pharma and MT-pharma, were selected for further computations, and showed excellent ROC curve quality. Additionally, the selected pharmacophore models were validated by the Guner-Henry decoy test method, which yielded high goodness of fit for WT-hDHFR and MT-hDHFR. Using a SMILES string of MTX in ZINC15 with the selections of 'clean', in vitro and in vivo options, 32 MTX-analogs were obtained. Eight analogs were filtered out due to their drug-like properties by applying absorption, distribution, metabolism, excretion, and toxicity (ADMET) assessment tests and Lipinski's Rule of five. WT-pharma and MT-pharma were further employed as a 3D query in virtual screening with drug-like MTX analogs. Subsequently, seven screening hits along with a reference compound (MTX) were subjected to molecular docking in the active site of WT- and MT-hDHFR. Through a clustering analysis and examination of protein-ligand interactions, one compound was found with a ChemPLP fitness score greater than that of MTX (reference compound). Finally, a simulation of molecular dynamics (MD) identified an MTX analog which exhibited strong affinity for WT- and MT-hDHFR, with stable RMSD, hydrogen bonds (H-bonds) in the binding site and the lowest MM/PBSA binding free energy. In conclusion, we report on an MTX analog which is capable of inhibiting hDHFR in wild type form, as well as in cases where the enzyme acquires resistance to drugs during chemotherapy treatment.

Effectiveness of occupational therapy in improving activities of daily living performance in complete cervical tetraplegic patients; A quasi experimental study.

OBJECTIVE: To determine the effectiveness of Occupational therapy (OT) in improving activities of daily living performance in complete cervical tetraplegic patients. METHODS: A quasi experimental study was conducted at Paraplegic Center Peshawar from May 2018 to March 2019. Seventy two spinal cord injury patients with complete cervical tetraplegia having age 18-60 years were included in the study using non probability convenience sampling technique. A trained Occupational therapist conducted two occupational therapy sessions per day, three times a week on alternative days for a period of six weeks. Self-care portion of Spinal cord injury independence measure (SCIM) was used to collect pre and post data. Data was analyzed using SPSS version 20. RESULTS: Mean age of the participants was 30.21±13.52.Majority of the participants were (n= 61, 84.7%) male while remaining 11 (15.3%) participants were female. Pre self-care total score was 0.39±0.987 and post self-care total score was 7.17±5.536. There was significant differences (P value <0.05) between pre and post scores of feeding, upper & lower body bathing, upper & lower body dressing, grooming and total self-care scores. CONCLUSION: Occupational therapy significantly improves activities of daily living performance in complete cervical tetraplegic spinal cord injury patients.

Correction to: Exploration for novel inhibitors showing back-to-front approach against VEGFR-2 kinase domain (4AG8) employing molecular docking mechanism and molecular dynamics simulations.

Following publication of the original article [1], the authors reported errors in Figure 3, Figure 14a, Figure 18, Figure 19b, Additional file 3 and Additional file 7.

Biochemical and Structural Insights Concerning Triclosan Resistance in a Novel YX7K Type Enoyl-Acyl Carrier Protein Reductase from Soil Metagenome.

Enoyl-acyl carrier protein reductase (ENR) catalyzes the last reduction step in the bacterial type II fatty acid biosynthesis cycle. ENRs include FabI, FabL, FabL2, FabK, and FabV. Previously, we reported a unique triclosan (TCL) resistant ENR homolog that was predominant in obligate intracellular pathogenic bacteria and Apicomplexa. Herein, we report the biochemical and structural basis of TCL resistance in this novel ENR. The purified protein revealed NADH-dependent ENR activity and shared similarity to prototypic FabI. Thus, this metagenome-derived ENR was designated FabI2. Unlike other prototypic bacterial ENRs with the YX6K type catalytic domain, FabI2 possessed a unique YX7K type catalytic domain. Computational modeling followed by site-directed mutagenesis revealed that mild resistance (20 µg/ml of minimum inhibitory concentration) of FabI2 to TCL was confined to the relatively less bulky side chain of A128. Substitution of A128 in FabI2 with bulky valine (V128) elevated TCL resistance. Phylogenetic analysis further suggested that the novel FabI2 and prototypical FabI evolved from a common short-chain dehydrogenase reductase family. To our best knowledge, FabI2 is the only known ENR shared by intracellular pathogenic prokaryotes, intracellular pathogenic lower eukaryotes, and a few higher eukaryotes. This suggests that the ENRs of prokaryotes and eukaryotes diverged from a common ancestral ENR of FabI2.

Employment status of person with disability in Government Departments in Khyber Pakhtunkhwa Pakistan.

OBJECTIVE: The objective of the study was to determine whether anticipated number of person with disabilities was employed in different government departments of Khyber Pakhtunkhwa (KPK), Pakistan. METHODS: Using across sectional survey, data from 16 different departments of provincial government of KPK was collected by volunteers of 'Friends of Paraplegics. RESULTS: Out of total 1, 71,137 Govt. employees, only 1151 (0.67%) were person with disabilities. None of the included departments fulfilled 2% allocated job quota for person with disabilities. Majority of the employees included in study were from Elementary & Secondary Education Department (n=140345) and Agriculture, Livestock & Cooperation Department (n=14315). The number of person with disabilities in these departments were 960 (0.68%) and 68 (0.48%) respectively. The highest percentage of person with disabilities were working in Higher Education Archives & Libraries Department (1.65%)followed by Law, Parliamentary Affairs & Human Rights Department (1.42%), Planning & Development Department (1.39%) and Administration & Establishment Department (1.16%). CONCLUSION: It is concluded that a small number of person with disabilities are employed in different Govt. Departments. Moreover, these departments has not ensured providing 2% job quota for person with disabilities.

Spinal cord injury due to fall from electricity poles after electrocution.

OBJECTIVE: The objective of this study was to report epidemiology, complications and rehabilitation outcomes of patients who sustained spinal cord injury (SCI) due to fall from electricity poles after electrocution. METHODS: A prospective observational study was conducted in which patients admitted to Paraplegic Centre Peshawar from July 2016 to July 2018 who sustained SCI due to fall from electricity poles after electrocution were included. Of total 852 patients, 39 (4.58%) sustained SCI due to fall from electricity poles after electrocution. Two patients were excluded and data of 37 patients was analyzed. RESULTS: The mean age of the participants was 35.03±13.47 years. Twenty-two (59.4%) patients had associated burns on different parts of body. Twenty-seven (72.9%) had pressure ulcers, 31 (83.8%) had spasticity, 18 (48.6%) had neuropathic pain and 2 (5.4%) had limb amputations due to injury. Mean Spinal cord injury independence measure score at the time of discharge was 53.4±5.7. CONCLUSION: SCI due to fall from electricity poles after electrocution is rare however combined effect of injury by electricity along with fall from electricity poles are associated with severe complications. Rehabilitation outcomes in these patients are also minimal.

Computational Simulations Identified Two Candidate Inhibitors of Cdk5/p25 to Abrogate Tau-associated Neurological Disorders.

Deregulation of Cdk5 is a hallmark in neurodegenerative diseases and its complex with p25 forms Cdk5/p25, thereby causes severe neuropathological insults. Cdk5/p25 abnormally phosphorylates tau protein, and induces tau-associated neurofibrillary tangles in neurological disorders. Therefore, the pharmacological inhibition of Cdk5/p25 alleviates tau-associated neurological disorders. Herein, computational simulations probed two candidate inhibitors of Cdk5/p25. Structure-based pharmacophore investigated the essential complementary chemical features of ATP-binding site of Cdk5 in complex with roscovitine. Resultant pharmacophore harbored polar interactions with Cys83 and Asp86 residues and non-polar interactions with Ile10, Phe80, and Lys133 residues of Cdk5. The chemical space of selected pharmacophore was comprised of two hydrogen bond donors, one hydrogen bond acceptor, and three hydrophobic features. Decoy test validation of pharmacophore obtained highest Guner-Henry score (0.88) and enrichment factor score (7.23). The screening of natural product drug-like databases by validated pharmacophore retrieved 1126 compounds as candidate inhibitors of Cdk5/p25. The docking of candidate inhibitors filtered 10 molecules with docking score >80.00 and established polar and non-polar interactions with the ATP-binding site residues of Cdk5/p25. Finally, molecular dynamics simulation and binding free energy analyses identified two candidate inhibitors of Cdk5/p25. During 30 ns simulation, the candidate inhibitors established <3.0 Šroot mean square deviation and stable hydrogen bond interactions with the ATP-binding site residues of Cdk5/p25. The final candidate inhibitors obtained lowest binding free energies of -122.18 kJ/mol and - 117.26 kJ/mol with Cdk5/p25. Overall, we recommend two natural product candidate inhibitors to target the pharmacological inhibition of Cdk5/p25 in tau-associated neurological disorders.

Pharmacotherapeutics and Molecular Mechanism of Phytochemicals in Alleviating Hormone-Responsive Breast Cancer.

Breast cancer (BC) is the leading cause of death among women worldwide devoid of effective treatment. It is therefore important to develop agents that can reverse, reduce, or slow the growth of BC. The use of natural products as chemopreventive agents provides enormous advantages. The aim of the current investigation is to determine the efficacy of the phytochemicals against BC along with the approved drugs to screen the most desirable and effective phytocompound. In the current study, 36 phytochemicals have been evaluated against aromatase to identify the potential candidate drug along with the approved drugs employing the Cdocker module accessible on the Discovery Studio (DS) v4.5 and thereafter analysing the stability of the protein ligand complex using GROningen MAchine for Chemical Simulations v5.0.6 (GROMACS). Additionally, these compounds were assessed for the inhibitory features employing the structure-based pharmacophore (SBP). The Cdocker protocol available with the DS has computed higher dock scores for the phytochemicals complemented by lower binding energies. The top-ranked compounds that have anchored with key residues located at the binding pocket of the protein were subjected to molecular dynamics (MD) simulations employing GROMACS. The resultant findings reveal the stability of the protein backbone and further guide to comprehend on the involvement of key residues Phe134, Val370, and Met374 that mechanistically inhibit BC. Among 36 compounds, curcumin, capsaicin, rosmarinic acid, and 6-shogaol have emerged as promising phytochemicals conferred with the highest Cdocker interaction energy, key residue interactions, stable MD results than reference drugs, and imbibing the key inhibitory features. Taken together, the current study illuminates the use of natural compounds as potential drugs against BC. Additionally, these compounds could also serve as scaffolds in designing and development of new drugs.