Pharmacophore modeling, 3D-QSAR and docking study of 2-phenylpyrimidine analogues as selective PDE4B inhibitors

Pharmacophore modeling, molecular docking, and molecular dynamics (MD) simulation studies have been performed, to explore the putative binding modes of 2-phenylpyrimidine series as PDE4B selective inhibitors. A five point pharmacophore model was developed using 87 molecules having pIC50 ranging from 8.52 to 5.07. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a high correlation coefficient (R(2)=0.918), cross validation coefficient (Q(2)=0.852), and F value (175) at 4 component PLS factor. The external validation indicated that our QSAR model possessed high predictive power (R(2)=0.70). The generated model was further validated by enrichment studies using the decoy test. To evaluate the effectiveness of docking protocol in flexible docking, we have selected crystallographic bound compound to validate our docking procedure as evident from root mean square deviation. A 10ns molecular dynamics simulation confirmed the docking results of both stability of the 1XMU-ligand complex and the presumed active conformation. Further, similar orientation was observed between the superposition of the conformations of 85 after MD simulation and best XP-docking pose; MD simulation and 3D-QSAR pose; best XP-docking and 3D-QSAR poses. Outcomes of the present study provide insight in designing novel molecules with better PDE4B selective inhibitory activity.