Lysine-specific demethylase 1 (LSD1) is a histone-modifying enzyme, which has been proposed as a promising target for anticancer drug development. Extensive research on LSD1 inhibitors has been performed since its discovery. In order to get more information for lead identification and optimization, we carried out a molecular modeling study on a set of 43 thieno[3,2-b]pyrrole competitive inhibitors of LSD1 using three-dimensional quantitative structure-activity relationship (3D-QSAR), molecular docking and molecular dynamics (MD) simulations. Based on the co-crystallized conformer-based alignment (CCBA) method, 3D-QSAR model of thieno[3,2-b]pyrrole derivatives as LSD1 inhibitors was established. The significant statistics (q2 = 0.595, r2 = 0.959, r2pred = 0.846) of the 3D-QSAR indicated the good predictive power and statistical reliability of this model. Based on the corresponding contour maps six LSD1 inhibitors were designed and their activities were predicted by 3D-QSAR model. Meanwhile, molecular docking was performed to simulate the probable binding modes between ligands and LSD1 protein. The molecular interactions mainly contributions to the binding affinity for LSD1 inhibitions were further supplemented by 100 ns MD simulations and binding free energy calculation.
Keywords: 3D-QSAR; LSD1; MM-GBSA; epigenetic; thieno[3,2-b]pyrrole derivatives.