In silico study directed towards identification of novel high-affinity inhibitors targeting an oncogenic protein:...

Bromodomain-containing protein 4 (BRD4) is a member of the bromodomain and extra-terminal domain (BET) family of proteins. It epigentically regulates the transcription of growth-promoting genes and has become an attractive target for the development of anticancer and anti-inflammatory agents. In the current study, we performed an in silico screening of a small-molecule chemical library against the acetyl-lysine binding site of the first bromodomain (BD1) in BRD4 protein. Potential inhibitors identified through virtual screening were further studied through molecular dynamics simulations, water entrapment analysis and Molecular Mechanics (MM)/Poisson-Boltzmann surface area (PBSA) binding free energy calculations. Many of the identified compounds exhibit better G-score (-11.64 kcal∙mol^-1 to -10.31 kcal∙mol^-1) and predicted binding affinity (-9.66 kcal∙mol^-1 to -6.63 kcal∙mol^-1) values towards BRD4-BD1 than that of the reference compound (+)-JQ1. Molecular dynamics simulation studies show that in free-form BRD4 the reported conserved water molecules are not retained at their specific positoins due to flexibiliy in the ZA-loop. In BRD4-ligand complexes the number and positions of conserved water molecules depends on the bound ligand. Identified potential inhibitors bind stably at the acetyl-lysine binding pocket of BRD4 and form direct and water-mediated hydrogen bonds with higher occupancy which may contribute to ligand specificity towards BRD4-BD1. Further, through MM/PBSA we calculated the binding free energies of selected compounds, which shows that they have comparable energies to that of (+)-JQ1, while NSC744713 shows better binding free energy.