Abstract
Bromodomain 4 and 9 (BRD4 and BRD9) have been regarded as important targets of drug designs in regard to the treatment of multiple diseases. In our current study, molecular dynamics (MD) simulations, deep learning (DL) and binding free energy calculations are integrated to probe the binding modes of three inhibitors (H1B, JQ1 and TVU) to BRD4 and BRD9. The MD trajectory-based DL successfully identify significant functional function domains, such as BC-loop and ZA-loop. The information from the post-processing analysis of MD simulations indicates that inhibitor binding highly influences the structural flexibility and dynamic behavior of BRD4 and BRD9. The results of the MM-GBSA calculations not only suggest that the binding ability of H1B, JQ1 and TVU to BRD9 are stronger than to BRD4, but they also verify that van der Walls interactions are the primary forces responsible for inhibitor binding. The hot spots of BRD4 and BRD9 revealed by residue-based free energy estimation provide target sites of drug design in regard to BRD4 and BRD9. This work is anticipated to provide useful theoretical aids for the development of selective inhibitors over BRD family members.