Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Identification of novel natural inhibitors of Opisthorchis felineus cytochrome P450 using structure-based screening and molecular dynamic simulation

Opisthorchis felineus is the etiological agent of opisthorchiasis in humans. O. felineus cytochrome P450 (OfCYP450) is an important enzyme in the parasite xenobiotic metabolism. To identify the potential anti-opisthorchid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,65,869) against the OfCYP450. The ligands were screened against OfCYP450 in four sequential docking modes that resulted in 361 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 10 compounds were found to fit well with re-docking studies. After refinement by docking and drug-likeness analyses, four potential inhibitors (ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789) were identified. These ligands with reference compounds (itraconazole and fluconazole) were further subjected to molecular dynamics simulation (MDS) and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the OfCYP450 and bound complexes. The binding energy analyses were also calculated. The results suggested that the compounds had a negative binding energy with ?259.41, ?110.09, ?188.25, ?163.30, ?202.10, and ?158.79 kJ mol^?1 for itraconazole, fluconazole, and compounds with IDs ZINC2358298, ZINC8790946, ZINC70707116, and ZINC85878789, respectively. These lead compounds displayed significant pharmacological and structural properties to be drug candidates. On the basis of MDS results and binding energy analyses, we concluded that ZINC8790946, ZINC70707116, and ZINC85878789 have excellent potential to inhibit OfCYP450.     

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     

Molecular modelling of Mycobacterium tuberculosis acetolactate synthase catalytic subunit and its molecular docking study with inhibitors

Mycobacterium tuberculosis is a leading cause of infectious disease in the world today. This outlook is aggravated by a growing number of M. tuberculosis infections in individuals who are immunocompromised as a result of HIV infections. Thus, new and more potent anti-TB agents are necessary. Therefore, acetolactate synthase (mtALS) was selected as a target enzyme to combat M. tuberculosis. In this work, the three-dimensional molecular model of the hypothetical structure for the ALS catalytic subunit of M. tuberculosis was elucidated by homology modelling. In addition, the orientations and binding affinities of sulfonylurea inhibitors with the new structure was investigated. Our findings could be helpful for the design of new, more potent mtAHAS 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.     

The product complex of M. tuberculosis malate synthase revisited

Enzymes of the glyoxylate shunt have been implicated as virulence factors in several pathogenic organisms, notably Mycobacterium tuberculosis and Candida albicans. Malate synthase has thus emerged as a promising target for design of anti-microbial agents. For this effort, it is essential to have reliable models for enzyme:substrate complexes. A 2.7 Angstroms resolution crystal structure for M. tuberculosis malate synthase in the ternary complex with magnesium, malate, and coenzyme A has been previously described. However, some unusual aspects of malate and Mg(++) binding prompted an independent determination of the structure at 2.3 Angstroms resolution, in the presence of saturating concentrations of malate. The electron density map of the complex reveals the position and conformation of coenzyme A to be unchanged from that found in the previous study. However, the coordination of Mg(++) and orientation of bound malate within the active site are different. The revised position of bound malate is consistent with a reaction mechanism that does not require reorientation of the electrophilic substrate during the catalytic cycle, while the revised Mg(++) coordination is octahedral, as expected. The results should be useful in the design of malate synthase inhibitors.     

Structures of citrate synthase and malate dehydrogenase of Mycobacterium tuberculosis

The tricarboxylic acid (TCA) cycle is a central metabolic pathway of all aerobic organisms and is responsible for the synthesis of many important precursors and molecules. TCA cycle plays a key role in the metabolism of Mycobacterium tuberculosis and is involved in the adaptation process of the bacteria to the host immune response. We present here the first crystal structures of M. tuberculosis malate dehydrogenase and citrate synthase, two consecutive enzymes of the TCA, at 2.6 Å and 1.5 Å resolution, respectively. General analogies and local differences with the previously reported homologous protein structures are described. Proteins 2015; 83:389–394. © 2014 Wiley Periodicals, Inc.     

Allosteric inhibitors of Mycobacterium tuberculosis tryptophan synthase

Global dispersion of multidrug resistant bacteria is very common and evolution of antibiotic‐resistance is occurring at an alarming rate, presenting a formidable challenge for humanity. The development of new therapeuthics with novel molecular targets is urgently needed. Current drugs primarily affect protein, nucleic acid, and cell wall synthesis. Metabolic pathways, including those involved in amino acid biosynthesis, have recently sparked interest in the drug discovery community as potential reservoirs of such novel targets. Tryptophan biosynthesis, utilized by bacteria but absent in humans, represents one of the currently studied processes with a therapeutic focus. It has been shown that tryptophan synthase (TrpAB) is required for survival of Mycobacterium tuberculosis in macrophages and for evading host defense, and therefore is a promising drug target. Here we present crystal structures of TrpAB with two allosteric inhibitors of M. tuberculosis tryptophan synthase that belong to sulfolane and indole‐5‐sulfonamide chemical scaffolds. We compare our results with previously reported structural and biochemical studies of another, azetidine‐containing M. tuberculosis tryptophan synthase inhibitor. This work shows how structurally distinct ligands can occupy the same allosteric site and make specific interactions. It also highlights the potential benefit of targeting more variable allosteric sites of important metabolic enzymes.     

Identification of peptidomimetic compounds as potential inhibitors against MurA enzyme of Mycobacterium tuberculosis

Abstract

Increasing prevalence of resistance to anti-tubercular drugs has become the foremost challenge now. According to WHO, over half a million of multidrug resistance cases (rifampicin, isoniazid, etc.) were reported in 2017, mostly emerging from countries such as China, India, and Russia. Therefore, developing new drugs or repurposing existing ones is need of the hour. The Mycobacterium cell wall biogenesis pathway offers many attractive targets for drug discovery against Tuberculosis (TB). MurA, a transferase enzyme that catalyzes the initial step of peptidoglycan (PG) biosynthesis, is one among them. A peptidoglycan layer resides over the plasma membrane and is an integral component of the bacterial cell wall. Therefore, disruption of their formation through inhibition of MurA enzyme should lead to deficiency in Mycobacterium cell synthesis. Based on this strategy, we have designed this study where two libraries of peptidomimetic compounds (Asinex & ChemDiv) were first screened against our modeled MurA structure and then validated through molecular dynamic simulations. From our virtual screening, top four compounds (ChemDiv: D675-0102, D675-0217; Asinex: BDE25373574, BDE 26717803) were selected based on their docking scores, binding energies, and interactions with catalytic site residues, for further evaluation. Results revealed stable ligand-MurA interactions throughout 50?ns of MD simulation and also druggability acceptable pharmacokinetic profile for all four compounds. Thus, based on our findings, these compounds could be considered as potential inhibitors of Mycobacterium MurA enzyme and hence be further tested for in vitro experimental validation as TB therapeutic drug candidate.

Communicated by Ramaswamy H. Sarma     

In vitro anti-TB properties,in silico target validation,molecular docking and dynamics studies of substituted 1,2,4-oxadiazole analogues against Mycobacterium tuberculosis

The alarming increase in multi- and extensively drug-resistant (MDR and XDR) strains of Mycobacterium tuberculosis (MTB) has triggered the scientific community to search for novel, effective, and safer therapeutics. To this end, a series of 3,5-disubstituted-1,2,4-oxadiazole derivatives (3a–3i) were tested against H37Rv, MDR and XDR strains of MTB. Of which, compound 3a with para-trifluorophenyl substituted oxadiazole showed excellent activity against the susceptible H37Rv and MDR-MTB strain with a MIC values of 8 and 16 µg/ml, respectively.To understand the mechanism of action of these compounds (3a–3i) and identify their putative drug target, molecular docking and dynamics studies were employed against a panel of 20 mycobacterial enzymes reported to be essential for mycobacterial growth and survival. These computational studies revealed polyketide synthase (Pks13) enzyme as the putative target. Moreover, in silico ADMET predictions showed satisfactory properties for these compounds, collectively, making them, particularly compound 3a, promising leads worthy of further optimisation.     

以蛋白激酶B为靶点的新型抗结核药物高通量筛选模型的建立和应用

【目的】建立以结核分枝杆菌蛋白激酶B为靶点的高通量筛选模型,并运用此模型进行化合物的筛选。【方法】克隆和表达结核分枝杆菌蛋白激酶B,并以其为靶酶建立并优化PknB抑制剂高通量筛选模型,利用该模型对化合物样品进行筛选,并对筛选到的阳性化合物进行抗菌和抑酶活性评价。【结果】利用该模型筛选了化合物样品18 000个,得到具有抑酶活性的阳性化合物8个,其中3个化合物具有较好的对结核分枝杆菌、海分枝杆菌、耻垢分枝杆菌的抑菌活性。【结论】建立的以PknB为靶点的抗结核药物高通量筛选模型具有灵敏度高、稳定性强等优点,可成功用于化合物的高效筛选。筛选得到3个在抑酶水平和抗菌方面均具有良好活性的阳性化合物样品,值得进一步研究。     

3D-QSAR studies and shape based virtual screening for identification of novel hits to inhibit MbtA in Mycobacterium tuberculosis

Mycobacterium tuberculosis, the pathogen responsible for tuberculosis, uses various strategies to survive in a variety of host lesions. The re-emergence of multi-drug-resistant strains of M. tuberculosis underlines the necessity to discover new molecules. Inhibitors of aryl acid adenylating enzyme, MbtA, involved in siderophore biosynthesis in M. tuberculosis, are being explored as potential anti tubercular agents. In this study, we have used 3D-QSAR models and shape based virtual screening to identify novel MbtA inhibitors. 3D-QSAR studies were carried out on nucleoside bisubstrate derivatives. Both Comparative Molecular Field Analysis (r²?=?.944 and r²_pred?=?.938) and Comparative Molecular Similarity Indices Analysis (r²?=?.892 and r²_pred?=?.842) models, developed using Gasteiger charges with all fields, predicted efficiently. A total of 13 hits were identified as novel prospective inhibitors for MbtA by utilizing an insilico workflow. Out of 13 hits, five top ranked hits were used for further molecular dynamics studies to gain more insights about the stability of the complexes.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

以苯丙氨酰-tRNA合成酶为靶点的新型抗结核药物高通量筛选模型的建立和应用

【目的】建立结核分枝杆菌PheRS抑制剂高通量模型,并运用此模型筛选化合物和发酵液样品。【方法】克隆和表达结核分枝杆菌PheRS蛋白并优化其酶活测定方法,在此基础上建立结核分枝杆菌PheRS抑制剂高通量筛选模型,并通过耻垢分枝杆菌作为检定菌对筛选到的样品进行抗菌活性测定及细胞毒性评价。【结果】运用此模型筛选了化合物样品11 600个,发酵液样品5 200个,筛选得到阳性化合物9个,阳性发酵液37个。而后通过耻垢分枝杆菌作为检定菌的抗菌活性测定及细胞毒性评价后,得到了6个发酵液阳性样品。【结论】建立的PheRS抑制剂模型可成功用于化合物和微生物发酵液的高效筛选,得到的6个发酵液阳性样品在酶水平和抗分枝杆菌方面均具有良好活性且毒性较低,值得进一步研究。     

In silico discovery of potential drug molecules to improve the treatment of isoniazid-resistant Mycobacterium tuberculosis

The emergence of multidrug-resistant Mycobacterium tuberculosis (M.tb) has become one of the major hurdles in the treatment of tuberculosis (TB). Drug-resistant M.tb has evolved with various strategies to avoid killing by the anti-tubercular drugs. Thus, there is a rising need to develop effective anti-TB drugs to improve the treatment of these strains. Traditional drug design approach has earned little success due to time and the cost involved in the process of development of anti-infective drugs. Numerous reports have demonstrated that several mutations in the drug target sites cause emergence of drug-resistant M.tb strains. In this study, we performed computational mutational analysis of M.tb inhA, fabD, and ahpC genes, which are the primary targets for first-line isoniazid (INH) drug. In silico virtual drug screening was performed to identify the potent drugs from a ChEMBL compound library to improve the treatment of INH-resistant M.tb. Further, these compounds were analyzed for their binding efficiency against active drug binding cavity of M.tb wild-type and mutant InhA, FabD and AhpC proteins. The drug efficacy of predicted lead compounds was verified by molecular docking using M.tb wild-type and mutant InhA, FabD and AhpC protein template models. Different in silico and pharmacophore analysis predicted three potent lead compounds with better drug-like properties against both M.tb wild-type and mutant InhA, FabD, and AhpC proteins as compared to INH drug, and thus may be considered as effective drugs for the treatment of INH-resistant M.tb strains. We hypothesize that this work may accelerate drug discovery process for the treatment of drug-resistant TB.

Communicated by Ramaswamy H. Sarma     

A comparative modeling and molecular docking study on Mycobacterium tuberculosis targets involved in peptidoglycan biosynthesis

An alarming rise of multidrug-resistant Mycobacterium tuberculosis strains and the continuous high global morbidity of tuberculosis have reinvigorated the need to identify novel targets to combat the disease. The enzymes that catalyze the biosynthesis of peptidoglycan in M. tuberculosis are essential and noteworthy therapeutic targets. In this study, the biochemical function and homology modeling of MurI, MurG, MraY, DapE, DapA, Alr, and Ddl enzymes of the CDC1551 M. tuberculosis strain involved in the biosynthesis of peptidoglycan cell wall are reported. Generation of the 3D structures was achieved with Modeller 9.13. To assess the structural quality of the obtained homology modeled targets, the models were validated using PROCHECK, PDBsum, QMEAN, and ERRAT scores. Molecular dynamics simulations were performed to calculate root mean square deviation (RMSD) and radius of gyration (Rg) of MurI and MurG target proteins and their corresponding templates. For further model validation, RMSD and Rg for selected targets/templates were investigated to compare the close proximity of their dynamic behavior in terms of protein stability and average distances. To identify the potential binding mode required for molecular docking, binding site information of all modeled targets was obtained using two prediction algorithms. A docking study was performed for MurI to determine the potential mode of interaction between the inhibitor and the active site residues. This study presents the first accounts of the 3D structural information for the selected M. tuberculosis targets involved in peptidoglycan biosynthesis.     

Discovery of novel inhibitors of Mycobacterium tuberculosis MurG: homology modelling,structure based pharmacophore,molecular docking,and molecular dynamics simulations

MurG (Rv2153c) is a key player in the biosynthesis of the peptidoglycan layer in Mycobacterium tuberculosis (Mtb). This work is an attempt to highlight the structural and functional relationship of Mtb MurG, the three-dimensional (3D) structure of protein was constructed by homology modelling using Discovery Studio 3.5 software. The quality and consistency of generated model was assessed by PROCHECK, ProSA and ERRAT. Later, the model was optimized by molecular dynamics (MD) simulations and the optimized model complex with substrate Uridine-diphosphate-N-acetylglucosamine (UD1) facilitated us to employ structure-based virtual screening approach to obtain new hits from Asinex database using energy-optimized pharmacophore modelling (e-pharmacophore). The pharmacophore model was validated using enrichment calculations, and finally, validated model was employed for high-throughput virtual screening and molecular docking to identify novel Mtb MurG inhibitors. This study led to the identification of 10 potential compounds with good fitness, docking score, which make important interactions with the protein active site. The 25 ns MD simulations of three potential lead compounds with protein confirmed that the structure was stable and make several non-bonding interactions with amino acids, such as Leu290, Met310 and Asn167. Hence, we concluded that the identified compounds may act as new leads for the design of Mtb MurG inhibitors.