With tuberculosis still being one of leading causes of death in the world and the emergence of drug-resistant strains of Mycobacterium tuberculosis (Mtb), researchers have been seeking to find further therapeutic strategies or more specific molecular targets. PknB is one of the 11 Ser/Thr protein kinases of Mtb and is responsible for phosphorylation-mediated signaling, mainly involved in cell wall synthesis, cell division and metabolism. With the amount of structural information available and the great interest in protein kinases, PknB has become an attractive target for drug development. This work describes the optimization and application of an in silico computational protocol to find new PknB inhibitors. This multi-level computational approach combines protein-ligand docking, structure-based virtual screening, molecular dynamics simulations and free energy calculations. The optimized protocol was applied to screen a large dataset containing 129,650 molecules, obtained from the ZINC/FDA-Approved database, Mu.Ta.Lig Virtual Chemotheca and Chimiothèque Nationale. It was observed that the most promising compounds selected occupy the adenine-binding pocket in PknB, and the main interacting residues are Leu17, Val26, Tyr94 and Met155. Only one of the compounds was able to move the active site residues into an open conformation. It was also observed that the P-loop and magnesium position loops change according to the characteristics of the ligand. This protocol led to the identification of six compounds for further experimental testing while also providing additional structural information for the design of more specific and more effective derivatives.
Keywords: MM/GBSA; Mycobacterium tuberculosis; PknB; molecular docking; molecular dynamics simulations; serine/threonine protein kinases; virtual screening.