Theoretically exploring selective-binding mechanisms of BRD4 through integrative computational approaches.

The origin of cancer is related to the dysregulation of multiple signal pathways and of physiological processes. Bromodomain-containing protein 4 (BRD4) has become an attractive target for the development of anticancer and anti-inflammatory agents since it can epigenetically regulate the transcription of growth-promoting genes. The synthesized BRD4 inhibitors with new chemical structures can reduce the drug resistance, but their binding modes and the inhibitory mechanism remain unclear. Here, we initially constructed robust QSAR models based on 68 reported tetrahydropteridin analogues using topomer CoMFA and HQSAR. On the basis of QSAR results, we designed 16 novel tetrahydropteridin analogues with modified structures and carried out docking studies. Instead of significant hydrogen bondings with amino acid residue Asn140 as reported in previous research, the molecular docking modelling suggested a novel docking pose that involves the amino acid residues (Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asp144, and Ile146) at the active site of BRD4. The MD simulations, free energy calculations, and residual energy contributions all indicate that hydrophobic interactions are decisive factors affecting bindings between inhibitors and BRD4. The current study provides new insights that can aid the discovery of BRD4 inhibitors with enhanced anti-cancer ability.