Binding Thermodynamics and Dissociation Kinetics Analysis Uncover the Key Structural Motifs of Phenoxyphenol Derivatives as the Direct InhA Inhibitors and the Hotspot Residues of InhA

Given the current epidemic of multidrug-resistant tuberculosis, there is an urgent need to develop new drugs to combat drug-resistant tuberculosis. Direct inhibitors of the InhA target do not require activation and thus can overcome drug resistance caused by mutations in drug-activating enzymes. In this work, the binding thermodynamic and kinetic information of InhA to its direct inhibitors, phenoxyphenol derivatives, were explored through multiple computer-aided drug design (CADD) strategies. The results show that the van der Waals interactions were the main driving force for protein–ligand binding, among which hydrophobic residues such as Tyr158, Phe149, Met199 and Ile202 have high energy contribution. The AHRR pharmacophore model generated by multiple ligands demonstrated that phenoxyphenol derivatives inhibitors can form pi–pi stacking and hydrophobic interactions with InhA target. In addition, the order of residence time predicted by random acceleration molecular dynamics was consistent with the experimental values. The intermediate states of these inhibitors could form hydrogen bonds and van der Waals interactions with surrounding residues during dissociation. Overall, the binding and dissociation mechanisms at the atomic level obtained in this work can provide important theoretical guidance for the development of InhA direct inhibitors with higher activity and proper residence time.