RESUMEN
DprE1 is involved in the synthesis of Mycobacterium tuberculosis cell wall and is a potent drug target for Tuberculosis (TB) treatment. The structure and dynamics of the loops L-I and L-II flanking the inhibitor binding site was studied using molecular dynamics (MD) simulation and MMPBSA in Amber v18. Docking and three-dimensional quantitative structure-activity relationship (3D-QSAR) of 55 Morpholino-pyrimidine (MP) inhibitors was carried out using Autodock v1.2.0 and Forge v10. ADMET analysis was done using SwissADME and pkCSM. All MP inhibitors docked in the DprE1 binding pocket, making contacts with L-II residues. MD studies showed that L-I and L-II unfold in the absence of the inhibitor but fold stably structure with reduced protein motions in the presence of MP-38, the highest affinity inhibitor. This was confirmed by k-means clustering and secondary structure analysis. L-II residues, L317, F320 and R325 contributed most towards the MMPBSA binding free energy of MP-38. A robust field-based 3D-QSAR model showed values of r2train = 0.982, r2test = 0.702 and q2 = 0.516. The MP inhibitor field points were broadly divided into negative electrostatics near the A, B rings and hydrophobic electrostatics near the D, E rings. Addition of negative groups at methanone position and ring B as well as addition of hydrophobic and bulky groups at ring E will improve activity. Highly active compounds 47, 49 and 50 of MP series exhibited highly favourable drug-like properties. SAR and ADMET insights attained from this model will help in the development of active DprE1 inhibitors in future.Communicated by Ramaswamy H. Sarma.