Discovery of novel BRD4-BD2 inhibitors via in silico approaches: QSAR techniques, molecular docking, and molecular dynamics simulations

Bromodomain-containing protein 4(BRD4) plays an important role in the occurrence and development of various malignant tumors, which has attracted the attention of scientific research institutions and pharmaceutical companies. The structural modification of most currently available BRD4 inhibitors is relatively simple, but the drug effectiveness is limited. Research has found that the inhibition of BD1 may promote the differentiation of oligodendrocyte progenitor cell; however, the inhibition of BD2 will not cause this outcome. Therefore, newly potential drugs which target BRD4-BD2 need further research. Herein, we initially built QSAR models out of 49 compounds using HQSAR, CoMFA, CoMSIA, and Topomer CoMFA technology. All of the models have shown suitable reliabilities (q² = 0.778, 0.533, 0.640, 0.702, respectively) and predictive abilities (r²_pred = 0.716, 0.6289, 0.6153, 0.7968, respectively) for BRD4-BD2 inhibitors. On the basis of QSAR results and the search of the R-group in the topomer search module, we designed 20 new compounds with high activity that showed appropriate docking score and suitable ADMET. Docking studies and MD simulation were carried out to reveal the amino acid residues (Asn351, Cys347, Tyr350, Pro293, and Asp299) at the active site of BRD4-BD2. Free energy calculations and free energy landscapes verified the stable binding results and indicated stable conformations of the complexes. These theoretical studies provide guidance and theoretical basis for designing and developing novel BRD4-BD2 inhibitors.

3D/2D-QSAR; ADMET; BRD4 inhibitors; Free energy landscape; MM/PBSA; Molecular docking; Molecular dynamics simulations.