Isopetasin and S-isopetasin as novel P-glycoprotein inhibitors against multidrug-resistant cancer cells

BackgroundA major problem of cancer treatment is the development of multidrug resistance (MDR) to chemotherapy. MDR is caused by different mechanisms such as the expression of the ABC-transporters P-glycoprotein (P-gp, MDR1, ABCB1) and breast cancer resistance protein (BCRP, ABCG2). These transporters efflux xenobiotic toxins, including chemotherapeutics, and they were found to be overexpressed in different cancer types.PurposeIdentification of novel molecules that overcome MDR by targeting ABC-transporters.MethodsResazurin reduction assay was used for cytotoxicity test. AutoDock 4.2. was used for molecular docking. The function of P-gp and BCRP was tested using a doxorubicin uptake assay and an ATPase assay. ROS generation was detected using flow cytometry for the measurement of H₂DCFH-DA fluorescence. Annexin/PI staining was applied for the detection of apoptosis. Bioinformatic analyses were performed using LigandScout 3.12. software and DataWarrior software.ResultsIn our search for new molecules that selectively act against resistant phenotypes, we identified isopetasin and S-isopetasin, which are bioactive natural products from Petasites formosanus. They exerted collateral sensitivity towards leukemia cells with high P-gp expression in CEM/ADR5000 cells, compared to sensitive wild-type CCRF-CEM leukemia cells. Also, they revealed considerable activity towards breast cancer cells overexpressing breast cancer resistance protein, MDA-MB-231-BCRP clone 23. This motivated us to investigate whether the function of P-gp was inhibited. In-silico results showed the compounds bound with high affinity and interacted with key amino acid residues in P-gp . Then, we found that the two compounds increased doxorubicin accumulation in P-gp overexpressing CEM/ADR5000 by three-fold compared to cells without inhibitor. P-gp-mediated drug efflux was ATP-dependent. Isopetasin and S-isopetasin increased the ATPase activity of human P-gp in a comparable fashion as verapamil used as control P-gp inhibitor. As isopetasin and S-isopetasin exerted dual roles, first as cytotoxic compounds and then as P-gp inhibitors, we suggested that their P-gp inhibition is part of a larger complex of mechanisms to induce cell death in cancer patients. P-gp dysfunction induces mitochondrial stress to generate ATP. Upon continuing stress by P-gp inhibition, the mitochondria generate reactive oxygen species (ROS). Initially established for verapamil, this theory was validated in the present study for isopetasin and S-isopetasin, as treatment with the two candidates increased ROS levels in CEM/ADR5000 cells followed by apoptosis.ConclusionOur study highlights the importance of isopetasin and S-isopetasin as novel ROS-generating and apoptosis-inducing P-gp inhibitors.     

Identification of novel CDK2 inhibitors by a multistage virtual screening method based on SVM,pharmacophore and docking model

Abstract

Cyclin-dependent kinase 2 (CDK2) is the family of Ser/Thr protein kinases that has emerged as a highly selective with low toxic cancer therapy target. A multistage virtual screening method combined by SVM, protein-ligand interaction fingerprints (PLIF) pharmacophore and docking was utilised for screening the CDK2 inhibitors. The evaluation of the validation set indicated that this method can be used to screen large chemical databases because it has a high hit-rate and enrichment factor (80.1% and 332.83 respectively). Six compounds were screened out from NCI, Enamine and Pubchem database. After molecular dynamics and binding free energy calculation, two compounds had great potential as novel CDK2 inhibitors and they also showed selective inhibition against CDK2 in the kinase activity assay.     

基于处方挖掘与分子动力学模拟筛选严重急性呼吸综合征冠状病毒2潜在抑制剂分子的研究

目的 从中药筛选具有潜在抑制严重急性呼吸综合征冠状病毒2 （SARS-CoV-2） 活性的成分，进一步从原子水平揭 示其抑制SARS-CoV-2 表面刺突蛋白（S 蛋白） 受体结合域（RBD） 与血管紧张素转化酶2 （ACE2） 结合的内在机制。 方法 检索新型冠状病毒（简称“新冠肺炎”） 治疗中药处方，构建“新冠肺炎中药候选活性成分数据库”。用具有ACE2 抑制活性的小分子化合物构建HipHop药效团模型，并对“新冠肺炎中药候选活性成分数据库”中活性成分筛选。采用分子 对接和分子动力学模拟方法研究候选活性成分与ACE2 的结合方式及其对SARS-CoV-2 S 蛋白与ACE2 识别的影响。 结果 本文通过中药处方挖掘和分子动力学模拟，从143 个新冠肺炎治疗中药处方中筛选出10 种可与SARS-CoV-2 S 蛋白/ 人源ACE2 识别位点结合的中药成分。其中，枇杷叶主要活性成分23-trans-p-coumaryhormentic acid 与ACE2 具有最高的亲和 力，且23-trans-p-coumaryhormentic acid 的结合可有效阻断SARS-CoV-2 S蛋白与宿主细胞ACE2 的结合。结论 本文通过虚 拟筛选发现了SARS-CoV-2 潜在抑制剂分子23-trans-p-coumaryhormentic acid，同时从原子水平预测了其抑制SARS-CoV-2 S 蛋白与ACE2 结合的内在机制，这将为SARS-CoV-2 特异性抗病毒药物的研发提供理论依据。     

Novel inhibitors of human glucose-6-phosphate dehydrogenase (HsG6PD) affect the activity and stability of the protein

BackgroundThe pentose phosphate pathway (PPP) has received significant attention because of the role of NADPH and R-5-P in the maintenance of cancer cells, which are necessary for the synthesis of fatty acids and contribute to uncontrollable proliferation. The HsG6PD enzyme is the rate-limiting step in the oxidative branch of the PPP, leading to an increase in the expression levels in tumor cells; therefore, the protein has been proposed as a target for the development of new molecules for use in cancer.MethodsThrough in vitro studies, we assayed the effects of 55 chemical compounds against recombinant HsG6PD. Here, we present the kinetic characterization of four new HsG6PD inhibitors as well as their functional and structural effects on the protein. Furthermore, molecular docking was performed to determine the interaction of the best hits with HsG6PD.ResultsFour compounds, JMM-2, CCM-4, CNZ-3, and CNZ-7, were capable of reducing HsG6PD activity and showed noncompetitive and uncompetitive inhibition. Moreover, experiments using circular dichroism and fluorescence spectroscopy showed that the molecules affect the structure (secondary and tertiary) of the protein as well as its thermal stability. Computational docking analysis revealed that the interaction of the compounds with the protein does not occur at the active site.ConclusionsWe identified two new compounds (CNZ-3 and JMM-2) capable of inhibiting HsG6PD that, compared to other previously known HsG6PD inhibitors, showed different mechanisms of inhibition.General significanceScreening of new inhibitors for HsG6PD with a future pharmacological approach for the study and treatment of cancer.     

In silico hit optimization toward AKT inhibition: fragment-based approach,molecular docking and molecular dynamics study

Abstract

Protein kinase B also known as AKT is a cardinal node in different signaling pathways that regulates diverse cell processes. AKT has three isoforms that share high homology. Hyperactivation of each isoform is related with different types of cancer. This work describes the computational search for new inhibitors using a hit optimization process of the previously reported AKT pan inhibitor, a 2,4,6-trisubstituted pyridine. A database of new molecules was proposed using a variant of fragment-based docking methodology and previous reported considerations. Molecular docking followed by molecular dynamics studies were performed to select the best compounds and analyze their behavior. Protein–ligand complexes energy was calculated using molecular mechanics Poisson–Boltzmann surface area protocol. Further, proposed molecules were compared with the ChEMBL database of compounds assayed against AKT. Data analysis leads to determine the structural requirements necessary for a favorable interaction of the proposed ligands with the AKT pocket. Molecular dynamics data suggested that the pK_a of the ligands is important for the stability in the AKT pocket. Molecular similarity analysis shows that proposed ligands have not been previously reported. Thus, ligands with high docking scores and favorable behavior on molecular dynamics simulations are proposed as potential AKT inhibitors.     

Resveratrol-salicylate derivatives as selective DNMT3 inhibitors and anticancer agents

Resveratrol is a natural polyphenol with plethora of biological activities. Resveratrol has previously shown to decrease DNA-methyltransferase (DNMT) enzymes expression and to reactivate silenced tumor suppressor genes. Currently, it seems that no resveratrol analogs have been developed as DNMT inhibitors. Recently, we reported the synthesis of resveratrol-salicylate derivatives and by examining the chemical structure of these analogs, we proposed that these compounds could exhibit DNMT inhibition especially that they resembled NSC 14778, a compound we previously identified as a DNMT inhibitor by virtual screening. Indeed, using in vitro DNMT inhibition assay, some of the resveratrol-salicylate analogs we screened in this work that showed selective inhibition against DNMT3 enzymes which were greater than resveratrol. A molecular docking study revealed key binding interactions with DNMT3A and DNMT3B enzymes. In addition, the most active analog, 10 showed considerable cytotoxicity against three human cancer cells; HT-29, HepG2 and SK-BR-3, which was greater than resveratrol. Further studies are needed to understand the anticancer mechanisms of these derivatives.     

Design,synthesis and biological evaluation of trinary benzocoumarin-thiazoles-azomethines derivatives as effective and selective inhibitors of alkaline phosphatase

Design, synthesis and characterization of new trinary Benzocoumarin-Thiazoles-Azomethine derivatives having three bioactive scaffolds in a single structural unit were carried out. The newly synthesized molecules were investigated for the inhibitory activity on human tissue nonspecific alkaline phosphatase (h-TNAP) and human intestinal alkaline phosphatase (h-IAP) isozymes. All the tested compounds exhibited the potent inhibition profile on both isozymes of alkaline phosphatase i.e., h-TNAP and h-IAP. Molecular docking studies were performed to explore the putative binding mode of interactions of selective inhibitors. Moreover, the synthesized derivatives were evaluated against cervical cancer cell line, HeLa and a few compounds exhibited significant inhibition in the range of 21.0–69.7%. The derivatives can be potential and selective alkaline phosphatase inhibitors for future studies.     

Insight into the Enzyme-Inhibitor Interactions of the First Experimentally Determined Human Aromatase

Abstract

Aromatase is an important pharmacological target in the anti-cancer therapy as the intratumoral aromatase is the source of local estrogen production in breast cancer tissues. Suppression of estrogen biosynthesis by aromatase inhibition represents an effective approach for the treatment of hormone-sensitive breast cancer. Because of the membrane-bound character and heme-binding instability, no crystal structure of aromatase was reported for a long time, until recently when crystal structure of human placental aromatase cytochrome P450 in complex with androstenedione was deposited in PDB. The present study is towards understanding the structural and functional characteristics of aromatase to address unsolved mysteries about this enzyme and elucidate the exact mode of binding of aromatase inhibitors. We have performed molecular docking simulation with twelve different inhibitors (ligands), which includes four FDA approved drugs; two flavonoids; three herbal compounds and three compounds having biphenyl motif with known IC₅₀ values into the active site of the human aromatase enzyme. All ligands showed favorable interactions and most of them seemed to interact to hydrophobic amino acids Ile133, Phe134, Phe221, Trp224, Ala306, Val370, Val373, Met374 and Leu477 and hydrophilic Arg115 and neutral Thr310 residues. The elucidation of the actual structure-function relationship of aromatase and the exact binding mode described in this study will be of significant interest as its inhibitors have shown great promise in fighting breast cancer.     

Molecular docking, 3D-QSAR,molecular dynamics,synthesis and anticancer activity of tyrosine kinase 2 (TYK 2) inhibitors

Abstract

A therapeutic rationale is proposed by selectively targeting tyrosine kinase 2 (TYK 2) to obtain potent TYK 2 inhibitors by molecular modeling studies. In the present study, we have taken tyrosine kinase (TYK 2) inhibitors and carried out molecular docking, 3?D quantitative structure–activity relationship (3D-QSAR) analysis and molecular dynamics (MD). Based on the 3D-QSAR results thirteen new compounds (R-1 to R-13) were designed and synthesized in good yields. The synthesized molecules were evaluated for their in vitro anticancer activity against LnCap and A549 cell lines. The molecules R-1, R-3, R-5, R-7, and R-10 exhibited considerable anti cancer activity.     

Kinetic and structural studies on the interactions of Torpedo californica acetylcholinesterase with two donepezil-like rigid analogues

Acetylcholinesterase inhibitors were introduced for the symptomatic treatment of Alzheimer’s disease (AD). Among the currently approved inhibitors, donepezil (DNP) is one of the most preferred choices in AD therapy. The X-ray crystal structures of Torpedo californica AChE in complex with two novel rigid DNP-like analogs, compounds 1 and 2, have been determined. Kinetic studies indicated that compounds 1 and 2 show a mixed-type inhibition against TcAChE, with K_i values of 11.12?±?2.88 and 29.86?±?1.12?nM, respectively. The DNP rigidification results in a likely entropy-enthalpy compensation with solvation effects contributing primarily to AChE binding affinity. Molecular docking evidenced the molecular basis for the binding of compounds 1 and 2 to the active site of β-secretase-1. Overall, these simplified DNP derivatives may represent new structural templates for the design of lead compounds for a more effective therapeutic strategy against AD by foreseeing a dual AChE and BACE-1 inhibitory activity.     

Computational aided mechanistic understanding of Camellia sinensis bioactive compounds against co-chaperone p23 as potential anticancer agent

Co-chaperon p23 has been well established as molecular chaperon for the heat shock protein 90 (Hsp90) that further leads to immorality in cancer cells by providing defense against Hsp90 inhibitors, and as stimulating agent for generating overexpressed antiapoptotic proteins, that is, Hsp70 and Hsp27. The natural compounds such as catechins from Camellia sinensis (green tea) are also well known for inhibition activity against various cancer. However, molecular interaction profile and potential lead bioactive compounds against co-chaperon p23 from green tea are not yet reported. To this context, we study the various secondary metabolites of green tea against co-chaperon p23 using structure-based virtual screening from Traditional Chinese Medicine (TCM) database. Following 26 compounds were obtained from TCM database and further studied for extra precision molecular docking that showed binding score between −10.221 and −2.276 kcal/mol with co-chaperon p23. However, relative docking score to known inhibitors, that is, ailanthone (−4.54 kcal/mol) and gedunin ( 3.60 kcal/mol) along with ADME profile analysis concluded epicatechin (−7.013 kcal/mol) and cis-theaspirone (−4.495 kcal/mol) as potential lead inhibitors from green tea against co-chaperone p23. Furthermore, molecular dynamics simulation and molecular mechanics generalized born surface area calculations validated that epicatechin and cis-theaspirone have significantly occupied the active region of co-chaperone p23 by hydrogen and hydrophobic interactions with various residues including most substantial amino acids, that is, Thr90, Ala94, and Lys95. Hence, these results supported the fact that green tea contained potential compounds with an ability to inhibit the cancer by disrupting the co-chaperon p23 activity.     

Identifying the novel inhibitors of lysine-specific demethylase 1 (LSD1) combining pharmacophore-based and structure-based virtual screening

Abstract

Lysine-specific demethylase 1 (LSD1) has been reported to connect with a range of solid tumors. Thus, the exploration of LSD1 inhibitors has emerged as an effective strategy for cancer treatment. In this study, we constructed a pharmacophore model based on a series of flavin adenine dinucleotide (FAD)-competing inhibitors bearing triazole???dithiocarbamate scaffold combining docking, structure–activity relationship (SAR) study, and molecular dynamic (MD) simulation. Meanwhile, another pharmacophore model was also constructed manually, relying on several speculated substrate-competing inhibitors and reported putative vital interactions with LSD1. On the basis of the two pharmacophore models, multi-step virtual screenings (VSs) were performed against substrate-binding pocket and FAD-binding pocket, respectively, combining pharmacophore-based and structure-based strategy to exploit novel LSD1 inhibitors. After bioassay evaluation, four compounds among 21 hits with diverse and novel scaffolds exhibited inhibition activity at the range of 3.63–101.43?μM. Furthermore, substructure-based enrichment was performed, and four compounds with a more potent activity were identified. After that, the time-dependent assay proved that the most potent compound with IC₅₀ 2.21?μM inhibits LSD1 activity in a manner of time-independent. In addition, the compound exhibited a cellular inhibitory effect against LSD1 in MGC-803 cells and may inhibit cell migration and invasion by reversing EMT in cultured gastric cancer cells. Considering the binding mode and SAR of the series of compounds, we could roughly deem that these compounds containing 3-methylxanthine scaffold act through occupying substrate-binding pocket competitively. This study presented a new starting point to develop novel LSD1 inhibitors.     

Design,Synthesis and Molecular Docking Studies of Some Tetrahydropyrimidine Derivatives as Possible Fascin Inhibitors

Eight derivatives of tetrahydropyrimidine scaffold were designed and prepared as hybrid compounds possessing the structural features of both monastrol as an anticancer drug and nifedipine as a fascin blocking agent. All of the compounds were evaluated for their cytotoxic potency and the ability to inhibit 4T1 breast cancer cells migration. Then, they were investigated in silico for their ability to inhibit the fascin protein using molecular docking simulation. The most potent compound was 4d and the weakest one was 4a according to the in vitro cytotoxicity assay. The corresponding IC₅₀ values were 193.70 and 248.75 μm , respectively. The least cytotoxic compound ( 4a ) was one of the strongest ones in binding to the fascin binding site according to the molecular docking results. 4a and 4e inhibited the 4T1 cells migration better than other compounds. They were more potent than nifedipine in inhibiting the migration process. In silico studies proved 4h to be the most potent fascin inhibitor in terms of ΔG_bind although it was not inhibiting migration. The controversy between the in vitro and in silico results may cancel the theory of the involvement of the fascin inhibition in the migration inhibition. However, the considerable antimigratory effects of some of the synthesized compounds encourage performing further in vivo experiments to introduce novel tumor metastasis inhibitors.     

Novel inhibitors of Rho-kinase mediated neuroinflammatory pathways and their potential application in recovery of injured spinal cord

Abstract

Spinal cord injury (SCI) involves damage to any part of the spinal cord which results in temporary or permanent changes in its function. Spinal cord secondary injury activates Rho-associated protein kinase 2 (ROCK2), which is involved in neuroinflammation and cell death by mediating secretion of inflammatory cytokines such as tumor necrosis factor-α (TNF-α), interleukin-1 beta (IL-1β), interleukin-2 (IL-2), and CXC chemokines. Here we evaluated potential inhibitors of ROCK2, Caspase-1, and TNF-α from Cissus quadrangularis derived natural compounds and compared them with structural analogues of quadrangularin by molecular docking, followed by correlation using molecular dynamic simulations studies. The results clearly demonstrate that the naturally derived compounds, quadrangularin and luteolin potentially inhibit ROCK2 and Caspase-1 with high binding affinity, and showed stable conformation throughout simulation trajectory period. Interestingly, quadrangularin and its structural analogues demonstrate effective binding affinity against ROCK2, caspase-1, and TNF-α when compared to their respective known inhibitors. From our studies, we can infer that natural compounds derived from C. quadrangularis are potentially capable of inhibitory activity against ROCK2, Caspase-1, and TNF-α. These findings could help in identifying novel therapeutic drugs targeting SCI.

Communicated by Ramaswamy H. Sarma     

Discovery of novel anti-breast cancer agents derived from deguelin as inhibitors of heat shock protein 90 (HSP90)

A series of O-substituted analogues of the B,C-ring truncated scaffold of deguelin were designed as C-terminal inhibitors of heat shock protein 90 (HSP90) and investigated as novel antiproliferative agents against HER2-positive breast cancer. Among the synthesized compounds, compound 80 exhibited significant inhibition in both trastuzumab-sensitive and trastuzumab-resistant breast cancer cells, whereas compound 80 did not show any cytotoxicity in normal cells. Compound 80 markedly downregulated the expression of the major client proteins of HSP90 in both cell types, indicating that the cytotoxicity of 80 in breast cancer cells is attributed to the destabilization and inactivation of HSP90 client proteins and that HSP90 inhibition represents a promising strategy to overcome trastuzumab resistance. A molecular docking study of 80 with the homology model of a HSP90 homodimer showed that 80 fit nicely in the C-terminal domain with a higher electrostatic complementary score than that of ATP.     

Design,synthesis and bioevaluation of novel umbelliferone analogues as potential mushroom tyrosinase inhibitors

Abstract

A series of umbelliferone analogues were synthesized and their inhibitory effects on the DPPH and mushroom tyrosinase were evaluated. The results showed that some of the synthesized compounds exhibited significant mushroom tyrosinase inhibitory activities. Especially, 2-oxo-2-[(2-oxo-2H-chromen-7-yl)oxy]ethyl-2,4-dihydroxybenzoate (4e) bearing 2,4-dihydroxy substituted phenyl ring exhibited the most potent tyrosinase inhibitory activity with IC₅₀ value 8.96?µM and IC₅₀ value of kojic acid is 16.69. The inhibition mechanism analyzed by Lineweaver–Burk plots revealed that the type of inhibition of compound 4e on tyrosinase was non-competitive. The docking study against tyrosinase enzyme was also performed to determine the binding affinity of the compounds. The compounds 4c and 4e showed the highest binding affinity with active binding site of tyrosinase. The initial structure activity relationships (SARs) analysis suggested that further development of such compounds might be of interest. The statistics of our results endorses that compounds 4c and 4e may serve as a structural template for the design and development of novel tyrosinase inhibitors.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.     

Novel triazolothiadiazines act as potent anticancer agents in liver cancer cells through Akt and ASK-1 proteins

Newly designed triazolothiadiazines incorporating with structural motifs of nonsteroidal analgesic anti-inflammatory drugs were synthesized and screened for their bioactivity against epithelial cancer cells. Compounds with bioactivities less then ∼5 μM (IC₅₀) were further analyzed and showed to induce apoptotic cell death and SubG₁ cell cycle arrest in liver cancer cells. Among this group, two compounds (1g and 1h) were then studied to identify the mechanism of action. These molecules triggered oxidative stress induced apoptosis through ASK-1 protein activation and Akt protein inhibition as demonstrated by downstream targets such as GSK3β, β-catenin and cyclin D1. QSAR and molecular docking models provide insight into the mechanism of inhibition and indicate the optimal direction of future synthetic efforts. Furthermore, molecular docking results were confirmed with in vitro COX bioactivity studies. This study demonstrates that the novel triazolothiadiazine derivatives are promising drug candidates for epithelial cancers, especially liver cancer.     

Virtual screening to identify novel potential inhibitors for Glutamine synthetase of Mycobacterium tuberculosis

Abstract

Glutamine synthetase (GS) of Mycobacterium tuberculosis (Mtb) is an essential enzyme which is involved in nitrogen metabolism and cell wall synthesis. It is involved in the inhibition of phagosome-lysosome fusion by preventing acidification. Targeting GS can be helpful to control the infection of Mtb. In order to identify potential inhibitors, we screened chemical libraries (56,400 compounds of ChEMBL anti-mycobacterial, 1596 FDA approved drugs, 419 Natural product and 916 phytochemical) against this target using the virtual screening approach. Screening by molecular docking identified ten top-ranked compounds as GSMtb inhibitors and they were compared with known inhibitors (as control). Since GS enzyme (GSHs) is also present in human. We have compared the protein sequence of GS from Mtb and human using the P-BLAST in NCBI. We found ～27% identity in between these two sequences, so we also compared the binding affinity of inhibitor between Mtb and human. Finally, we identified top two compounds namely CHEMBL387509, CHEMBL226198 from ChEMBL anti-mycobacterial dataset, and Eriocitrin and Malvidin from phytochemical dataset which showed lees binding affinity towards GSHs whereas Pamidronate, and Phentermine from FDA approved drugs and (-)-Quinic Acid, Hexopyranuronic acid, Quebrachit, and Castanospermine from natural product showed protein-ligand interaction with Mtb protein while no interaction with GSHs. The top two docked complexes were subjected to molecular dynamic simulation to understand the stability of the molecule. Further, we calculated the binding free energy of the docked complex and analyzed hydrogen bond, salt bridge, pie stacking, and hydrophobic interaction in the docking region. These ligands exhibited very good binding affinity GSMtb enzymes. Therefore, these ligands are novel and drug-likeness compounds, and they may be potential inhibitors of M tuberculosis.

Communicated by Ramaswamy H. Sarma