ABSTRACT
We have carried out a computational structure-based design of new potent pyrrolidine carboxamide (PCAMs) inhibitors of enoyl-acyl carrier protein reductase (InhA) of Mycobacterium tuberculosis (MTb). Three-dimensional (3D) models of InhA-PCAMx complexes were prepared by in situ modification of the crystal structure of InhA-PCAM1 (Protein Data Bank (PDB) entry code: 4U0J), the reference compound of a training set of 20 PCAMs with known experimental inhibitory potencies (IC50(exp)). First, we built a gas phase quantitative structure-activity relationships (QSAR) model, linearly correlating the computed enthalpy of the InhA-PCAM complex formation and the IC50(exp). Further, taking into account the solvent effect and loss of inhibitor entropy upon enzyme binding led to a QSAR model with a superior linear correlation between computed Gibbs free energies (ΔΔGcom) of InhA-PCAM complex formation and IC50(exp) (pIC50(exp) = -0.1552·ΔΔGcom + 5.0448, R² = 0.94), which was further validated with a 3D-QSAR pharmacophore model generation (PH4). Structural information from the models guided us in designing of a virtual combinatorial library (VL) of more than 17 million PCAMs. The VL was adsorption, distribution, metabolism and excretion (ADME) focused and reduced down to 1.6 million drug like orally bioavailable analogues and PH4 in silico screened to identify new potent PCAMs with predicted IC50(pre) reaching up to 5 nM. Combining molecular modeling and PH4 in silico screening of the VL resulted in the proposed novel potent antituberculotic agent candidates with favorable pharmacokinetic profiles.
