Volume 11 Supplement 1
Virtual Screening of potential drug-like inhibitors against Lysine/DAP pathway of Mycobacterium tuberculosis
© Garg et al; licensee BioMed Central Ltd. 2010
Published: 18 January 2010
An explosive global spreading of multidrug resistant Mycobacterium tuberculosis (Mtb) is a catastrophe, which demands an urgent need to design or develop novel/potent antitubercular agents. The Lysine/DAP biosynthetic pathway is a promising target due its specific role in cell wall and amino acid biosynthesis. Here, we report identification of potential antitubercular candidates targeting Mtb dihydrodipicolinate synthase (DHDPS) enzyme of the pathway using virtual screening protocols.
In the present study, we generated three sets of drug-like molecules in order to screen potential inhibitors against Mtb drug target DHDPS. The first set of compounds was a combinatorial library, which comprised analogues of pyruvate (substrate of DHDPS). The second set of compounds consisted of pyruvate-like molecules i.e. structurally similar to pyruvate, obtained using 3D flexible similarity search against NCI and PubChem database. The third set constituted 3847 anti-infective molecules obtained from PubChem. These compounds were subjected to Lipinski's rule of drug-like five filters. Finally, three sets of drug-like compounds i.e. 4088 pyruvate analogues, 2640 pyruvate-like molecules and 1750 anti-infective molecules were docked at the active site of Mtb DHDPS (PDB code: 1XXX used in the molecular docking calculations) to select inhibitors establishing favorable interactions.
The above-mentioned virtual screening procedures helped in the identification of several potent candidates that possess inhibitory activity against Mtb DHDPS. Therefore, these novel scaffolds/candidates which could have the potential to inhibit Mtb DHDPS enzyme would represent promising starting points as lead compounds and certainly aid the experimental designing of antituberculars in lesser time.
Causing massacre especially in Asia and Africa, Tuberculosis (TB) prevalence and mortality rates have probably been mounting globally for last several years . Further, association of TB with HIV patients and emergence of multiple drug-resistant Mycobacterium tuberculosis (Mtb) to isoniazid and rifampicin and extensive drug-resistant Mtb to any floroquinolone, amikacin and capreomycin is a growing alarm. Despondently, more than two million people happen to be victim of TB annually and globally [2–4]. World Health Organization (WHO) 2008 report has mentioned the statistics regarding the occurrence of 9.2 million new cases and 1.7 million deaths from TB in 2006, out of which 0.7 million cases and 0.2 million deaths were in HIV-positive patients . These numbers observed to be boosted compared with those reported by the WHO for the previous years. Therefore, discovery of novel unexploited drug target enzymes and their inhibitors besides generating analogues of existing drugs is a major challenge in the field of drug discovery and designing.
The amino acids play a major role in defining the cellular growth, cell wall and protein synthesis of bacterial system. Importantly, the absence of de novo synthesis of protein building blocks and requirement of amino acids as dietary components in mammals implies that specific inhibitors of amino acid biosynthetic pathways would display a novel class of antibacterial agents through inhibition of cell wall and protein synthesis with no mammalian toxicity. For past few years, Lysine/DAP biosynthetic pathway has been gaining high attention due to its foremost feature in the synthesis of D, L-diaminopimelic acid (meso-DAP) and lysine. Both components are essential for cross-linking peptidoglycan chains to provide strength and rigidity to the bacterial cell wall [6–8]. It has been observed that Mycobacterium cell walls are characterized by an unusual high DAP content. Moreover, gene-knockout experiments with Mycobacterium smegmatis has demonstrated the essentiality of DAP pathway for the bacteria, where the absence of DAP results in cell lysis and death . In view of its importance, the designing of potential inhibitors against any enzyme of this pathway may display a novel classes of antitubercular agents.
The present study mainly focused on dihydrodipicolinate synthase (DHDPS) enzyme of the pathway, catalyses the first committed step towards meso-DAP formation by condensation of substrate pyruvate with active site residue (LYS-171), which results in the formation of a Schiff-base [10, 11]. Next, tautomerisation and aldol type reaction with aspartate β-semialdehyde generates an enzyme-tethered acyclic intermediate that undergoes transimination to form heterocyclic [(4S)-4-hydroxy-2,3,4,5-tetrahydro-(2S)-dipicolinate] (HTPA). The release of HTPA from the active site with elimination of water molecules provides product dihydrodipicolinate (DHDP) . The three-dimensional crystal structures of DHDPS from Escherichia coli, Nicotiana sylvestris, Staphylococcus aureus, Mtb, Salmonella typhimurium, Bacillus anthracis, Clostridium botulinum, Corynebacterium glutamicum, Thermotoga maritime and Bacillus clausii are available at PDB database. Previously, various structural studies have reported the conservation of active site residues from different bacterial species [13–21].
Till date, designing of inhibitors against DHDPS (mainly from E. coli) is being carried out using experimental procedure; however, no potent inhibitors have been reported. However, analogues of pyruvate such as α-ketobutyrate, α-ketoglutarate, glyoxylate and fluoropyruvate have been shown to be competitive inhibitors of DHDPS with respect to pyruvate. Additionally, few inhibitors based on DHDP or HTPA structures showing weak to moderate inhibitory activity is also reported [22–24]. Recently, Mitsakos et al has demonstrated that several experimentally known inhibitors displayed a clear differentiation in inhibition of DHDPS enzymes from different bacterial species, hence, suggested that designing of inhibitors against DHDPS should be specific to bacterial species rather than a broad-spectrum inhibitor.
Keeping in view, the importance of DAP pathway in Mtb and low outcome of DHDPS inhibitors using experimental procedures, we have made an attempt to screen inhibitors of Mtb DHDPS using virtual screening procedures. The present work of screening of DHDPS inhibitors is reported for the first time, hence, would be a great help in aiding the experimental studies and rational development of novel drugs against Mtb.
Generation of combinatorial library of pyruvate analogues
Since the similarities in structures are indicative of similarities in bioactivities, therefore, structure based searching of databases/libraries has been gaining high demand nowadays [28, 29]. In the present study, VLifeMDS package was used for carrying out 3D flexible search (using pyruvate as a template) against most popular databases i.e. National cancer Institute (NCI) Open database with ~2,60,071 compounds, and a library of sub/superstructures of pyruvate from PubChem database which constituted 21,061 compounds. These databases were first subjected to Lipinski's rule of five constraints (mentioned earlier), which reduced their size to ~160,000 and 17,157 compounds respectively. The searching was based on 3D superimposition, where the query pyruvate and the compounds present in the databases were structurally superimposed for rmsd calculations, in order to check whether the atoms in the match mapping meets the spatial constrains (distance, angle, dihedral) in a query or not.
The databases generally store only a single low-energy conformation or a limited number of conformations for each compound, which may lead to the reduction in the hit rates. Therefore to make the search more effective, the 3D flexible search was carried out in the present study. Here, the conformers for each compound was generated using Metropolis Monte Carlo simulations, which explores the compound's conformational search space using random moves by altering torsion angle values. Here, the tolerance limit (which defines the rmsd cut-off) was set to 50%, such that the hits with rmsd value greater than 50% were discarded. Following the 3D search, 291 and 2349 compounds from NCI and PubChem libraries were retrieved. Finally, we got total 2640 molecules, which are structurally similar to pyruvate molecules. In this study, these molecules will be called pyruvate-like molecules.
Additionally, 3847 anti-infective compounds, consisting of 1743 antibacterials were retrieved from PubChem database. Out of 3847, only 1750 anti-infectives satisfied the Lipinski's rule of five constraints. These compounds were highly diverse from the pyruvate such that none of the anti-infectives showed 2D/3D similarity with pyruvate. Hence, the docking of these compounds would help to screen the diverse classes of antitubercular agents against Mtb DHDPS.
Ligand-receptor flexible docking
To find the binding affinities between target receptor and screened out compounds, an automated flexible docking of ligands at the flexible active site of receptor was carried out using AutoDock (v.4.0) software . The software facilitates the internal degree of freedom along with the values of translation and rotation for the side chains of selected active residues as well as for the ligands in search of its suitable bound conformations. Undoubtedly, introduction of flexibility makes the docking process computationally more expensive but more superior than rigid ligand-receptor docking. Before docking process, several separate pre-docking steps: ligand preparation, receptor preparation and grid map calculations were performed. The ligand and receptor preparation stage involved the addition of hydrogen atoms, computing charges, merging non-polar hydrogen atoms and defining AD4 atom types to ensure that atoms conformed to the AutoDock atom types. Next, information about rotatable torsion bonds that defines the bond flexibility was acquired. The ligands and receptor molecule preparation was followed by grid construction using AutoGrid module. During grid construction, atom types of the ligand, which acted as probes in the calculation of grid maps, were identified. The grid with default volume of 40 × 40 × 40 Å with a spacing of 0.375 Å centered on the receptor was prepared. For conformational search, the docking calculations using the genetic algorithm (GA) procedure with default parameters was performed. The GA computed the fitness of a docked candidate every time by measuring the minimum values of free energy binding (ΔG) based on different types of energy evaluations. In the present study, the python scripts were used for carrying out automated docking process.
Results and discussion
In the present study, different virtual screening approaches were used for selecting potential inhibitors against Mtb DHDPS. The first approach employed the generation of combinatorial library i.e. analogues of pyruvate. In the second approach pyruvate-like molecules generated using 3D flexible similarity search against available databases/libraries. Thirdly, to screen diverse classes of antitubercular agents, drug-like 1750 anti-infectives available at PubChem database were retrieved. Finally, these three sets of compounds i.e. generated pyruvate analogues, pyruvate-like molecules and anti-infectives were docked into the active site of receptor with the purpose of sorting potential inhibitors of Mtb DHDPS.
The three-dimensional crystal structure of Mtb DHDPS stored in the PDB file (code: 1XXX)  was obtained from protein databank. Mtb DHDPS is a homotetramer and each subunit with 300 amino acid residues comprises: N-terminal (β/α)8-barrel domain (residues 1-233) and a C-terminal domain (residues 234-300), which consists of three α-helices. The residues bounded the active site are THR54, THR55 TYR143, ARG148 and LYS171. Particularly, LYS171 responsible for substrate binding and catalysis are located at the centre of each monomer in the (b/α)8-barrel domain facing the central cavity of the tetramer. In E. coli DHDPS enzyme the equivalent residue is LYS161. Using PYMOL software, all the 1587 water molecules, eight DTT molecules, eight Mg2+ and eight Cl- ions were removed.
Docking of pyruvate analogues
Top 10 pyruvate analogues with better free binding energy values in comparison with substrate
ΔG values (kcal/mol)
Docking of pyruvate-like molecules
The virtual screening using 3D similarity based search could provide two main advantages- i) helps to narrow down the size of large databases/libraries to be screened, which eventually reduces the computational time required to dock each library compound and ii) 3D searching retrieves the compounds, out of which some are the same class of compounds to which query belongs, but some others may be entirely new classes of compounds, which may directly lead to the discovery of novel lead compounds.
Free binding energy values and structures for top ten NCI hits
ΔG values (kcal/mol)
Free binding energy values and structures for top five PubChem compounds
ΔG values (kcal/mol)
Docking of anti-infectives
Docking results for top five anti-infectives retrieved from PubChem database
ΔG values (kcal/mol)
To conclude, we have employed several virtual screening protocols such as generation of combinatorial library, 3D flexible search and molecular docking to identify potential inhibitors against Mtb DHDPS. Several potential drug-like inhibitors have been screened out showing strong binding affinity to Mtb DHDPS. Additionally, few anti-infectives with highly diverse topology from the pyruvate also displayed strong binding. Though experimental studies are indispensable to mark them as lead compound for the development of novel drugs against Mtb, however, screened out inhibitors would undoubtedly aid the experimental designing of antitubercular agents expeditiously.
List of abbreviations used
- The abbreviations used are:
DAP: Diaminopimelic Acid
- Mtb :
National cancer Institute
AG is thankful to CSIR for providing SRF. The authors are thankful to the CSIR and Department of Biotechnology, Government of India for financial assistance. The authors are also thankful to Mammon Rashid for correcting grammatical mistakes.
This article has been published as part of BMC Bioinformatics Volume 11 Supplement 1, 2010: Selected articles from the Eighth Asia-Pacific Bioinformatics Conference (APBC 2010). The full contents of the supplement are available online at http://www.biomedcentral.com/1471-2105/11?issue=S1.
- Bloom BR, Murray CJ: Tuberculosis: commentary on a reemergent killer. Science 1992, 257: 1055–64. 10.1126/science.257.5073.1055View ArticlePubMedGoogle Scholar
- LoBue P: Extensively drug-resistant tuberculosis. Curr Opin Infect Dis 2009, 22: 167–173. 10.1097/QCO.0b013e3283229fabView ArticlePubMedGoogle Scholar
- Jain A, Dixit P: Multidrug-resistant to extensively drug resistant tuberculosis: what is next? J Biosci 2008, 33: 605–616. 10.1007/s12038-008-0078-8View ArticlePubMedGoogle Scholar
- Rattan A, Kalia A, Ahmad N: Multidrug-resistant Mycobacterium tuberculosis : molecular perspectives. Emerg Infect Dis 1998, 4: 195–209. 10.3201/eid0402.980207PubMed CentralView ArticlePubMedGoogle Scholar
- WHO Report[http://www.who.int/tb/publications/global_report/2008/summary/en/index.html]
- Hutton CA, Perugini MA, Gerrard JA: Inhibition of lysine biosynthesis: an evolving antibiotic strategy. Mol BioSyst 2007, 3: 458–465. 10.1039/b705624aView ArticlePubMedGoogle Scholar
- Cox RJ, Sutherland A, Vederas JC: Bacterial diaminopimelate metabolism as a target for antibiotic design. Bioorg Med Chem 2000, 8: 843–871. 10.1016/S0968-0896(00)00044-4View ArticlePubMedGoogle Scholar
- Hutton CA, Southwood TJ, Turner JJ: Inhibitors of lysine biosynthesis as antibacterial agents. Mini Rev Med Chem 2003, 3: 115–127. 10.2174/1389557033405359View ArticlePubMedGoogle Scholar
- Pavelka MS Jr, Jacobs WR Jr: Biosynthesis of diaminopimelate, the precursor of lysine and a component of peptidoglycan, is an essential function of Mycobacterium smegmatis . J Bacteriol 1996, 178: 6496–507.PubMed CentralPubMedGoogle Scholar
- Shedlarski JG, Gilvarg C: The pyruvate-aspartic semialdehyde condensing enzyme of Escherichia coli . J Biol Chem 1970, 245: 1362–1373.PubMedGoogle Scholar
- Laber B, Gomis-Rüth FX, Romão MJ, Huber R: Escherichia coli dihydrodipicolinate synthase. Identification of the active site and crystallization. Biochem J 1992, 288: 691–695.PubMed CentralView ArticlePubMedGoogle Scholar
- Blickling S, Renner C, Laber B, Pohlenz H, Holak T, Huber R: Reaction mechanism of Escherichia coli dihydrodipicolinate synthase investigated by X-ray crystallography and NMR spectroscopy. Biochem J 1997, 36: 24–33. 10.1021/bi962272dView ArticleGoogle Scholar
- Karsten WE: Dihydrodipicolinate Synthase from Escherichia coli : pH dependent changes in the kinetic mechanism and kinetic mechanism of allosteric inhibition by L-Lysine. Biochem J 1997, 36: 1730–1739. 10.1021/bi962264xView ArticleGoogle Scholar
- Blagova E, Levdikov V, Milioti N, Fogg MJ, Kalliomaa AK, Brannigan JA, Wilson KS, Wilkinson AJ: Crystal structure of dihydrodipicolinate synthase (BA3935) from Bacillus anthracis at 1.94 A resolution. Proteins 2006, 62: 297. 10.1002/prot.20684View ArticlePubMedGoogle Scholar
- Burgess BR, Dobson RC, Dogovski C, Jameson GB, Parker MW, Perugini MA: Purification, crystallization and preliminary X-ray diffraction studies to near-atomic resolution of dihydrodipicolinate synthase from methicillin-resistant Staphylococcus aureus . Acta Crystallogr Sect F Struct Biol Cryst Commun 2008, 64: 659–61. 10.1107/S1744309108016746PubMed CentralView ArticlePubMedGoogle Scholar
- Pearce FG, Perugini MA, McKerchar HJ, Gerrard JA: Dihydrodipicolinate synthase from Thermotoga maritima . Biochem J 2006, 400: 359–66. 10.1042/BJ20060771PubMed CentralView ArticlePubMedGoogle Scholar
- Dobson RCJ, Valegard K, Gerrard JA: The crystal structure of three site-directed mutants of Escherichia coli dihydrodipicolinate synthase: further evidence for a catalytic triad. J Mol Biol 2004, 338: 329–339. 10.1016/j.jmb.2004.02.060View ArticlePubMedGoogle Scholar
- Blickling S, Beisel HG, Bozic D, Knablein J, Laber B, Huber R: Structure of dihydrodipicolinate synthase of Nicotiana sylvestris reveals novel quaternary structure. J Mol Biol 1997, 274: 608–621. 10.1006/jmbi.1997.1393View ArticlePubMedGoogle Scholar
- Mirwaldt C, Korndorfer I, Huber R: The crystal structure of dihydrodipicolinate synthase from Escherichia coli at 2.5 resolution. J Mol Biol 1995, 246: 227–239. 10.1006/jmbi.1994.0078View ArticlePubMedGoogle Scholar
- Kefala G, Evans GL, Griffin MD, Devenish SR, Pearce FG, Perugini MA, Gerrard JA, Weiss MS, Dobson RC: Crystal structure and kinetic study of dihydrodipicolinate synthase from Mycobacterium tuberculosis . Biochem J 2008, 411: 351–60. 10.1042/BJ20071360View ArticlePubMedGoogle Scholar
- Kefala G, Weiss MS: Cloning, expression, purification, crystallization and preliminary X-ray diffraction analysis of DapA (Rv2753c) from Mycobacterium tuberculosis . Acta Crystallograph Sect F Struct Biol Crystal Commun 2006, 62: 1116–1119. 10.1107/S1744309106039844View ArticleGoogle Scholar
- Turner JJ, Gerrard JA, Hutton CA: Heterocyclic inhibitors of dihydrodipicolinate synthase are not competitive. Bioorg Med Chem 2005, 13: 2133–40. 10.1016/j.bmc.2005.01.001View ArticlePubMedGoogle Scholar
- Turner JJ, Healy JP, Dobson RC, Gerrard JA, Hutton CA: Two new irreversible inhibitors of dihydrodipicolinate synthase: diethyl (E,E)-4-oxo-2,5-heptadienedioate and diethyl (E)-4-oxo-2-heptenedioate. Bioorg Med Chem Lett 2005, 15: 995–8. 10.1016/j.bmcl.2004.12.043View ArticlePubMedGoogle Scholar
- Turner JJ, Gerrard JA, Hutton CA: Conformationally constrained diketopimelic acid analogues as inhibitors of dihydrodipicolinate synthase. Bioorg Med Chem 2005, 13: 2133–40. 10.1016/j.bmc.2005.01.001View ArticlePubMedGoogle Scholar
- Mitsakos V, Dobson RCJ, Pearce FG, Devenish SR, Evans GL, Burgess BR, Perugini MA, Gerrard JA, Hutton CA: Inhibiting dihydrodipicolinate synthase across species: towards specificity for pathogens? Bioorg Med Chem Lett 2008, 18: 842–844. 10.1016/j.bmcl.2007.11.026View ArticlePubMedGoogle Scholar
- VLife Sciences technologies: MDS 2.0 Molecular Design Suite. Pvt Ltd, Pune, India; 2003.Google Scholar
- Lipinski CA: Drug-like properties and the causes of poor solubility and poor permeability. J Pharmacol Toxicol Methods 2000, 44: 235–249. 10.1016/S1056-8719(00)00107-6View ArticlePubMedGoogle Scholar
- Muegge I, Oloff S: Advances in virtual screening. Drug Discov Today 2006, 3: 405–411. 10.1016/j.ddtec.2006.12.002View ArticleGoogle Scholar
- Lyne PD: Structure-based virtual screening: an overview. Drug Discov Today 2002, 7: 1047–1055. 10.1016/S1359-6446(02)02483-2View ArticlePubMedGoogle Scholar
- Morris GM, Goodsell DS, Halliday RS, Huey R, Hart WE, Belew RK, Olson AJ: Automated docking using a Lamarckian Genetic algorithm and empirical binding free energy function. J Computational Chemistry 1998, 19: 1639–1662. 10.1002/(SICI)1096-987X(19981115)19:14<1639::AID-JCC10>3.0.CO;2-BView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.