Bromine-containing domain protein 4 (BRD4) plays a crucial role in regulating transcription and genome stability. Selective inhibitors of BRD4-BD1 can specifically target specific bromine domains to affect cell proliferation, apoptosis, and differentiation. In this work, 43 selective benzoazepinone BRD4-BD1 inhibitors were studied using molecular simulations and three-dimensional quantitative conformation relationships (3D-QSAR). A reliable 3D-QSAR model was established based on COMFA (Q2 = 0.532, R2 = 0.981) and COMSIA (S + E + H (Q2 = 0.536, R2 = 0.979) two different analysis methods. Through 3D-QSAR model prediction and quantum chemical analysis, 15 small molecules with stronger inhibitory activity than the template compounds were constructed, and 5 new compounds with higher predictive activity and binding affinity were screened by molecular docking and ADMET methods. According to the molecular dynamics simulation, the key residues that can interact with BRD4-BD1 protein and molecular docking results are consistent, including ASN140, MET132, GLN85, MET105, ASN135 and TYR97. From the MD trajectory, we calculated and analyzed RMSD, RMSF, free binding energy, FECM, DCCM and PCA, the loop region formed by amino acids VAL45∼PRO62 showed α-helix, β-folding and clustering towards the active center with the extension of simulation time. Further optimization of the structure of active candidate compounds A6, A11, A14, and A15 will provide the necessary theoretical basis for the synthesis and activity evaluation of novel BRD4-BD1 inhibitors.Communicated by Ramaswamy H. Sarma.
Keywords: 3D-QSAR; bromodomain-containing protein 4 (BRD4); molecular docking; molecular dynamics simulation.