Structure-Based Design and in Silico Screening of Virtual Combinatorial Library of Benzamides Inhibiting 2-trans Enoyl-Acyl Carrier Protein Reductase of

backgroundDuring the previous decade a new class of benzamide-based inhibitors of 2-trans enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (Mt) with unusual binding mode have emerged. Here we report in silico design and evaluation of novel benzamide InhA-Mt inhibitors with favorable predicted pharmacokinetic profiles.

methodsBy using in situ modifications of the crystal structure of N-benzyl-4-((heteroaryl)methyl) benzamide (BHMB)-InhA complex (PDB entry 4QXM), 3D models of InhA-BHMBx complexes were prepared for a training set of 19 BHMBs with experimentally determined inhibitory potencies (half-maximal inhibitory concentrations IC50exp). In the search for active conformation of the BHMB1-19, linear QSAR model was prepared, which correlated computed gas phase enthalpies of formation (∆∆HMM) of InhA-BHMBx complexes with the IC50exp. Further, taking into account the solvent effect and entropy changes upon ligand, binding resulted in a superior QSAR model correlating computed complexation Gibbs free energies (∆∆Gcom). The successive pharmacophore model (PH4) generated from the active conformations of BHMBs served as a virtual screening tool of novel analogs included in a virtual combinatorial library (VCL) of compounds containing benzamide scaffolds. The VCL filtered by Lipinski's rule-of-five was screened by the PH4 model to identify new BHMB analogs.

resultsGas phase QSAR model: -log

conclusionsCombined use of QSAR models that considered binding of the BHMBs to InhA, pharmacophore model, and ADME properties helped to recognize bound active conformation of the benzamide inhibitors, permitted in silico screening of VCL of compounds sharing benzamide scaffold and identification of new analogs with predicted high inhibitory potencies and favorable pharmacokinetic profiles.