Discovery of new anti-tuberculosis drugs with novel mode of action is urgently needed. The tryptophan synthase is a genetically validated enzyme that catalyzes last step of tryptophan biosynthetic pathway required for growth and survival of Mycobacterium tuberculosis. Here, a ligand-based pharmacophore model was built using molecular operating environment (MOE) software (version 2010.12) and validation of generated pharmacophoric features was done using active, inactive and decoy set of molecules. The generated pharmacophore model was used for screening of 7,523,972 drug-like molecules of ZINC database. The best matches (RMSD < 1) retrieved as a result of screening were subjected to molecular docking studies into active pocket of α-subunit of tryptophan synthase from M. tuberculosis. The five hits were selected and validated through anti-tuberculosis activity analysis. Finally, a new inhibitor ZINC09150898 has been identified with best binding score -32.07 kcal/mol, showing 100% growth inhibition of M. tuberculosis (H37Rv strain) at 50 µg/mL. This identified inhibitor-protein complex was further subjected to MD simulations studies (50 ns) involving root mean square deviation, root mean square fluctuation, secondary structure analysis and pocket interaction analysis to explore its binding mode stability inside active pocket. The binding free energies of inhibitor-protein complex through MM-PBSA analysis suggested that van der Waals interactions play a vital role for retention of identified inhibitor inside the protein pocket. All these analyses confirmed retention of ligand inside pocket and no unfolding in protein structure was observed over explored time scale.Communicated by Ramaswamy H. Sarma.
Keywords: MD simulations; MMPBSA; Pharmacophore; RMSD; anti-tuberculosis; molecular docking.