Next Article in Journal
Surface Interaction of Ionic Liquids: Stabilization of Polyethylene Terephthalate-Degrading Enzymes in Solution
Next Article in Special Issue
Accurate Physical Property Predictions via Deep Learning
Previous Article in Journal
Screening of Plant Pollen Sources, Polyphenolic Compounds, Fatty Acids and Antioxidant/Antimicrobial Activity from Bee Pollen
Previous Article in Special Issue
High-Throughput Screening Campaign Identified a Potential Small Molecule RXFP3/4 Agonist
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 27
Issue 1
10.3390/molecules27010118
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors

by
Siao Chen
,
Yi He
,
Yajiao Geng
,
Zhi Wang
,
Lu Han
* and
Weiwei Han
*
Key Laboratory for Molecular Enzymology and Engineering of Ministry of Education, National Engineering Laboratory for AIDS Vaccine, School of Life Science, Jilin University, Changchun 130012, China
*
Authors to whom correspondence should be addressed.
Molecules 2022, 27(1), 118; https://doi.org/10.3390/molecules27010118
Submission received: 11 November 2021 / Revised: 22 December 2021 / Accepted: 22 December 2021 / Published: 26 December 2021
(This article belongs to the Special Issue Advances in Computational Chemistry for Drug Design, Discovery and Screening)
Browse Figures
Versions Notes

Abstract

:
Bromodomain and extra-terminal domain (BET) subfamily is the most studied subfamily of bromodomain-containing proteins (BCPs) family which can modulate acetylation signal transduction and produce diverse physiological functions. Thus, the BET family can be treated as an alternative strategy for targeting androgen-receptor (AR)-driven cancers. In order to explore the effect of inhibitors binding to BRD4 (the most studied member of BET family), four 150 ns molecular dynamic simulations were performed (free BRD4, Cpd4-BRD4, Cpd9-BRD4 and Cpd19-BRD4). Docking studies showed that Cpd9 and Cpd19 were located at the active pocket, as well as Cpd4. Molecular dynamics (MD) simulations indicated that only Cpd19 binding to BRD4 can induce residue Trp81-Ala89 partly become α-helix during MD simulations. MM-GBSA calculations suggested that Cpd19 had the best binding effect with BRD4 followed by Cpd4 and Cpd9. Computational alanine scanning results indicated that mutations in Phe83 made the greatest effects in Cpd9-BRD4 and Cpd19-BRD4 complexes, showing that Phe83 may play crucial roles in Cpd9 and Cpd19 binding to BRD4. Our results can provide some useful clues for further BCPs family search.

1. Introduction

Bromodomain-containing proteins (BCPs) can decipher the acetylation code of proteins via binding to their acetylated lysine (KAc) residues, which further modulates acetylation signal transduction and produces diverse physiological functions [1,2,3,4,5,6,7,8,9,10]. BCPs are divided into eight subfamilies containing the bromodomain and extra-terminal domain (BET) subfamily [11,12]. The BET subfamily [12] consists of BRD2 [13], BRD3 [14], BRD4 [12,15] and BRD testis-specific protein BRDT [16,17]. The BET protein family is highly similar in structure and function and is mainly composed of BD1, BD2 and ET domains and CHAT and CTM functional domains. BD1 and BD2 are responsible for the recognition and binding of histone acetylation sites at the H3 and H4 ends. For BRD4 protein, BD1 and BD2 also have partial kinase properties. BRD4-BD1 protein consists of A, B, C, Z four α-helix, ZA and BC two loops, in addition to ZA channel, WPF pocket and the cavity between four α-helix. The active pocket of BRD4 was presented in Figure 1. Recognition of histones with BRD4 is mediated through these cavities. Androgen regulated BRDs increased chromatin accessibility by enhancing AR recruitment to chromatin, which may drive cancer progression. BET inhibitors could be an alternative strategy for targeting AR-driven cancers due to the interaction between the BET family proteins and AR. In particular, the effect of BRD inhibitors is focused on cell growth and resistance in cancer cells and on inflammatory signaling in immune cells [18]. The use of such inhibitors has helped to understand the cell/tissue-specific functions that these proteins perform.
Many BRD inhibitors, such as (+)-JQ-1 [19], OTX-015 [20], I-BET762 [21], CPI0610 [22] and ABBV-07 [23] have been exploited. For example, Hu et al. pointed out some novel BRD4-selective inhibitors with an adihydroquinoxalin-2(1H)-one scaffold from the PLK1-BRD4 inhibitor [24,25]. Watts et al. designed a series of ALK and BRD4 dual inhibitors from BI-2536 [26]. Liu et al. described a series of compounds with different selectivity between PLK1 and BRD4 using BI-2536 as the starting compound [27] and so on.
In 2019, Hu et al. performed further lead optimization to generate a novel potent series of bromodomain and extra-terminal (BET) inhibitors with a (1,2,4-triazol-5-yl)-3,4-dihydroquinoxalin-2(1H)-one structure, and they found that compound 19 is a potential BET protein drug candidate for the treatment of cancer [12]. However, what effects do the potent and selective inhibitors binding to BRD4 have on the active site? MD simulations were a widely used powerful tool, simulating the physical movements of atoms and molecules by Newton’s equations of motion and qualitatively and quantitatively analyze conformational changes related to its function. Fernando D. Prieto-Martínez et al. performed the MD simulations to explore the conformational changes for Flavonoids binding to BRD4. The results showed that amentoflavone makes numerous contacts in the ZA channel [28]. In order to explore the effect of inhibitors binding to BRD4, four 150 ns molecular dynamic simulations were performed (free BRD4, compound 4 (Cpd4)-BRD4, compound 9 (Cpd9)-BRD4, compound 19 (Cpd19)-BRD4). Our theoretical results may be useful for designing excellent inhibitors for BRD4.

2. Methods

2.1. Protein Prepare

BRD4 was obtained from the Protein Data Bank (https://www.rcsb.org, accessed on 20 Febraury 2019) and used as starting structures (PDB code: 6JI3) [12]. The 3D structure of BRD4-Cpd4 was download from PDB (PDB code: 6JI4) [12]. The 3D structure of Cpd9 and Cpd19 were obtained by Gaussian View [29] and then optimized with Gaussian 09 software [29] at B3LYP 6-31G* set [30]. Then, Cpd9 and Cpd19 were docked to BRD4 with AutoDock 4.2.6 software [31]. The Lamarckian genetic algorithm was used to predict the stereo conformation and select the lowest energy system for further MD.

2.2. Molecular Dynamics Simulations

The complexes of the three inhibitors and BRD4 were performed by Amber 16 software with 150 ns molecular dynamics simulations (The sufficiency to really reach the proper equilibrium of the MD simulation needs to be explored according to the characteristics of the system. We have found these papers for the BET family with MD simulations for 100–200 ns [32,33,34,35,36,37,38]. In this study, the 150 ns MD simulations were performed). The system setting and the 2D chemical structure of the potent inhibitors were shown in Table 1. The force field of BRD4 was ff99SB force field [39] and the inhibitors were generated with GAFF parameters and RESP partial charges through the Antechamber program.
The protein-ligand complexes were placed in a cube box with a boundary of 12 Å filled with periodic boundary conditions and a TIP3P water model [40]. Cl− was added to maintain the electrical neutrality of the system [41]. The energy of the system was minimized by the steepest descent method and conjugate gradient algorithm, and then the simulated system was heated (from 0–310 K), density balanced (50 ps), and slowly released (500 ps at 310 K), respectively.
At last, four 150 ns MD unconstrained simulations were conducted by Amber 16 packages using the particle-mesh Ewald (PME) method in each system [42]. The SHAKE algorithm is used to constrain the hydrogen bonds [43]. Langevin dynamics method (1 atm constant pressure and 310 K temperature) was treated as Langevin thermostat during MD simulations [44]. The Berendsen pressurizer was used for constant pressure simulation and the cut-off radius was set to 8 Å. All of the coordinates were saved every 2 ps. We have performed three times MD simulations for four systems (see Figures S1–S3).
In order to analyze the conformational changes of BRD4-inhibitors during simulations, its backbone was superimposed and then clustered with a distance cut-off of 2.5 Å. Root-mean-square deviation (RMSD), root-mean-square fluctuation (RMSF), hydrogen bonding, dihedral angle distribution, and secondary structure form action were analyzed by using the Cpptraj module in Amber 16. Moreover, the Bio3D software was employed to perform principal component fluctuation analysis and dynamic residue cross-correlation analysis. The popular VMD software was used to visualize the MD trajectories.

2.3. Principal Components Analysis (PCA) and Cross-Correlation Matrix

PCA is a popular statistical multivariate method that analyzes trajectories data from MD simulations [45]. Dynamics cross-correlation matrix (DCCM) among the atoms of protein I and j (cij) was selected as a matrix from MD trajectories using the Bio3d package by R program [46].

2.4. MM/GBSA Calculations

Free energy calculations play an important role in computational biology, such as drug design and protein structure prediction, which require the participation of free energy calculations of ligands binding to protein [47,48]. There are many ways to calculate the free energy of ligands binding, such as free energy perturbation, thermodynamic integration and so on. These methods are strict in theory, they require very high computation, and become very expensive and time-consuming with the increase of system size and have poor convergence for complex systems. However, MM/GBSA is an analytical technique widely used in MD simulations calculated by the MMPBSA.py program in the Amber 16 package using Generalized Born (GB) implicit solvent models [49,50]. The sum and the free energy are divided into the molecular mechanics and solvation energies, the energy of the acceptor, the Ligand and the complex in the solution were calculated. Then the difference energy was calculated, and the combined free energy data and energy decomposition data were obtained. A total of 1000 frames were extracted from the 150 ns in each MD simulation for the calculation of binding free energies.

3. Result and Discussion

3.1. Docking Study

The 3D structure of the complex Cpd4-BRD4 was downloaded from Protein Data Bank (PDB code 6JI3) [12]. Figure S4 indicated that the active residues around Cpd4 binding to BRD4. It can be seen that Asn140, Tyr97, Ile146, Phe83, Leu94, Pro82, Leu92 and Tyr139. In particular, Asn140 made a hydrogen bond with BRD4, which suggested that Asn140 may play important role in Cpd4 binding to BRD4.
We redocked Cpd4 to BRD4 with AutoDock 4.2.6 software to compare with the docking pose with Cpd4 to BRD4 [30]. The best docking score was shown in Figure S5. The RMSD plot between the docking score and the reference Cpd4 was 0.88Å, showing that AutoDock 4.2.6 software was perfect and can be used for further docking study. Figure 2a,b showed that Cpd9 and Cpd19 were located at the active pocket. As seen in Figure 2a, Ile146, Asn140, Tyr97, Pro82, Phe83, Val87, Leu92 and Trp81 were the key residues for Cpd9 binding to BRD4. As seen from Figure 2b, Thr131, Asn135, Cys136, Met132, Met105, Asp136 Asn140, Tyr97, Ile146, Phe83, Val87, Pro82, Gln85, Trp81 and Lys91 were the key residues for Cpd19 binding to BRD4. In particular, Cys136, Met105 and Lys91 made hydrogen bonds with Cpd19. We also estimated the free energies of binding for BRD4-Cpd4, BRD4-Cpd9 and BRD4-Cpd19, which were −5.37, −5.80 and −7.16 kcal/mol, respectively. To sum up, the three complexes were all located at the active pocket and all were stable; thence the three complexes can be used for further study.

3.2. The Stability of the Four Systems

The RMSD plots during the 150 ns MD simulations were shown in Figure 3a,b. As seen in the RMSD plots, the three systems have reached equilibrium after 120 ns. Average RMSD values of free-BRD4, BRD4-Cpd4, BRD4-Cpd9, BRD4-Cpd19 were RMSD 3.51, 2.76, 2.41 and 2.11 Å, respectively.
As shown in Figure 3c,d, the Rg values of Free-BRD4, BRD4-Cpd4 and BRD4-Cpd9 were stabilized at 14.88 Å, 14.71 Å and 14.87 Å, respectively, while the BRD4-Cpd19 system was stabilized at 14.87 Å. Figure 3e,f showed the hydrophilicity of the four systems throughout the 150 ns MD simulations. In Figure 3e,f, it is obvious that the solvent accessible surface area (SASA) plot of free-BRD4 was similar to those of the other three systems in the 150 ns. All the results indicated that the four systems were stable and can be used for further study.

3.3. The Conformational Changes for Three Inhibitors Binding

After 150 ns MD simulations, the RMSF of the backbone atoms of the proteins in the four systems was calculated and the results were shown in Figure 4. It can be seen that residues 81 to 89 only in the Cpd19-BRD4 complex exhibited distinct atom-positional fluctuation amplitudes during 150 ns MD simulations. The results indicated that Cpd19 bound to BRD4 can make it less flexible and maintain an ordered structure of BRD4, thereby allowing the tighter fitting of Cpd19 into the enzyme active site.
To investigate the conformational changes, the difference in the secondary structure (DSSP) plot (Figure 5) was obtained. The DSSP of residue 81 to 83 in the BRD4-Cpd19 system significantly differed from the other three systems. As shown in Figure 5b, residue 81 to 83 in the BRD4-Cpd19 system partly formed α-helix. However, these residues of the other three systems showed a turn of BRD4. It was reported that the WPF shelf (W81, P82, F83), as well as the N140, are essential for acetylated lysine binding [43]. Our results also indicated that WPF shelves were important for inhibitors binding.
We calculated the probability to form α-helix of the four systems. The probabilities of residues Trp81, ProP82 and Phe83 in four systems were listed in Table 2. The probabilities in BRD4-Cpd19 (33.25 %) were considerably higher than those in the three other systems, suggesting that residue 81 to 83 in the BRD4-Cpd19 system partly formed α-helix during MD simulations. Therefore, the binding of Cpd19 was beneficial to maintain the order of the structure of BRD4 and it is useful to Cpd19 slide into the active pocket. Cpd19 is an excellent inhibitor binding to BRD4, which is consistent with the experimental data (IC50: 5.3 ± 0.4 nM) [12].

3.4. Cross-Correlation Analysis

The DCCM for all pairs of Cα atoms was shown in Figure 6a–d. Seen from Figure 6d, it can be concluded that the BRD4-Cpd19 system has the lightest color on the cross-correlation matrix map, indicating that the BRD4-Cpd19 system has experienced the weakest flexibility and the most stability during the MD simulation, which showed that BRD4 was an excellent inhibitor binding to BRD4 with the experimental data (IC50: 5.3 ± 0.4 nM) [12].
Figure 7a–d showed the B factor value during MD simulations. It was found that the protease catalytic active sites were usually located in the low B value region, Figure 7a is consistent with the conclusion. B factor reflects protein flexibility. The higher the B factor in the protein structure is, the better the mobility of the protein. It can be seen in Figure 7d, the active region of the BRD4-Cpd19 system experienced the weakest flexibility and the most stability during the MD simulation.
In this study, the similar structures of the trajectories of the four systems were divided into different groups using the RMSD-based clustering method (Figure 8a–d). Through cluster analysis, the most representative structure was discovered and was chosen to compare the binding affinities of different ligands.
For BRD4-Cpd4 complex (Figure 9a), Pro82, Phe83, Gln85, Val87, Pro86, Asp106, Met107, Met132, Gln135, Cys136 and Tyr139 had interactions with Cpd4. For BRD4-Cpd9 (Figure 9b), it can be seen that Phe83, Pro82, Lys91, Leu92, Tyr139, Cys136, Ile146 and Asn140 had interactions with Cpd9. From Figure 9c, it can be seen that Pro82, Phe83, Gln85, Val87, Leu92, Leu94, Met105, Asp106, Cys136, Tyr139, Ile146 and Asp145 had interactions with Cpd19. Among three systems, BRD4 is bound most closely to Cpd19 (12 residues), and weakly bound to Cpd4 (11 residues) and Cpd9 (eight residues).

3.5. The Free Energy Binding of Three Systems

The free energy between BRD4 and the three inhibitors was calculated in AmberTools17 package [41]. The results were listed in Table 3. In Table 3, the energy in the BRD4-Cpd19 system was lower than that in the two complexes. The results were consistent with experimental data [12].
In this analysis, the residues surrounding the inhibitors (Cpd4, Cpd9 and Cpd19) were mutated to alanine to explore their influence on the stability of protein structures. The results were presented in Figure 10a–c. Among them (Figure 10b,c), mutations in Phe83 indicated the greatest effects on the protein-inhibitor complex, indicating that Phe83 may play crucial roles in binding to Cpd9 and Cpd19.

4. Conclusions

The addition of different inhibitors had different local effects on the protein. Molecular dynamics simulations showed that the binding of three inhibitors influenced the secondary structure changes in BRD4. Only Cpd19 binding to BRD4 induced residues Trp81-Ala89 to partly become a helix but the effects of Cpd4 and Cpd9 were not significant. MM-GBSA calculations suggested that, among three inhibitors, Cpd19 had the lowest binding energy with BRD4 followed by Cpd4, Cpd9.
The above results showed that the Cpd19 combined with BRD4 is the best and consisted with the highly efficient experiment results. Cpd4 and Cpd9, as mild inhibitors, had little effect on the motions of protein and lower binding energy to protein. Residue Trp81-Ala89 in the computer-simulated inhibitor-protein complex may provide some useful clues for further BCPs family search.

Supplementary Materials

The following are available online. Figure S1: (a) Root-mean-square deviation (RMSD) for the backbone atoms, (b) Radius-of-gyration (Rg) plot, (c) Solvent-accessible-surface-area (SASA) plot of the second MD simulation; Figure S2: (a) Root-mean-square deviation (RMSD) for the backbone atoms, (b) Radius-of-gyration (Rg) plot, (c) Solvent-accessible-surface-area (SASA) plot of the third time MD simulation; Figure S3: (a) average Root-mean-square deviation (RMSD) for the backbone atoms, (b) average Radius-of-gyration (Rg) plot of the third time MD simulation, (c) average Solvent-accessible-surface-area (SASA) plot; Figure S4: Binging pocket of The active residues around Cpd 4 binding to BRD4; Figure S5: Cpd 4 redocked to the pocket of BRD4 with AutoDock 4.2.6 software to compared with the docking pose of 6JI3. (blue is Cpd4 of 6JI3 and green is redock Cpd4).

Author Contributions

S.C. wrote this paper. Y.H. and Y.G. prepared the tables and figures. Z.W. made the Supplementary Materials. L.H. and W.H. revised this paper. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (31870201), Overseas Cooperation Project of Jilin Province (20200801069GH), Science and technology development program of Jilin Province (20200301029RQ).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank D.A. Case for providing free access to the AMBER software.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Urbanucci, A.; Mills, I.G. Bromodomain-containing proteins in prostate cancer. Mol. Cell. Endocrinol. 2018, 462, 31–40. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Fujisawa, T.; Filippakopoulos, P. Functions of bromodomain-containing proteins and their roles in homeostasis and cancer. Nat. Rev. Mol. Cell Biol. 2017, 18, 246–262. [Google Scholar] [CrossRef] [PubMed]
  3. Zhou, Q.; Huang, J.; Zhang, C.; Zhao, F.; Kim, W.; Tu, X.; Zhang, Y.; Nowsheen, S.; Zhu, Q.; Deng, M.; et al. The bromodomain containing protein BRD-9 orchestrates RAD51-RAD54 complex formation and regulates homologous recombination-mediated repair. Nat. Commun. 2020, 11, 1–13. [Google Scholar] [CrossRef]
  4. Wu, S.-Y.; Lee, C.-F.; Lai, H.-T.; Yu, C.-T.; Lee, J.-E.; Zuo, H.; Tsai, S.Y.; Tsai, M.-J.; Ge, K.; Wan, Y.; et al. Opposing functions of BRD4 isoforms in breast cancer. Mol. Cell 2020, 78, 1114–1132. [Google Scholar] [CrossRef] [PubMed]
  5. Wen, N.; Guo, B.; Zheng, H.; Xu, L.; Liang, H.; Wang, Q.; Wang, D.; Chen, X.; Zhang, S.; Li, Y.; et al. Bromodomain inhibitor jq1 induces cell cycle arrest and apoptosis of glioma stem cells through the VEGF/PI3K/AKT signaling pathway. Int. J. Oncol. 2019, 55, 879–895. [Google Scholar] [CrossRef]
  6. Wang, D.; Kon, N.; Lasso, G.; Jiang, L.; Leng, W.; Zhu, W.-G.; Qin, J.; Hong, B.; Gu, W. Acetylation-regulated interaction between p53 and SET reveals a widespread regulatory mode. Nature 2016, 538, 118–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Zanconato, F.; Battilana, G.; Forcato, M.; Filippi, L.; Azzolin, L.; Manfrin, A.; Quaranta, E.; Di Biagio, D.; Sigismondo, G.; Guzzardo, V.; et al. Transcriptional addiction in cancer cells is mediated by YAP/TAZ through BRD4. Nat. Med. 2018, 24, 1599–1610. [Google Scholar] [CrossRef] [PubMed]
  8. Jain, A.K.; Michelle, C.B. Bromodomain histone readers and cancer. J. Mol. Biol. 2017, 429, 2003–2010. [Google Scholar] [CrossRef]
  9. Lu, L.; Chen, Z.; Lin, X.; Tian, L.; Su, Q.; An, P.; Li, W.; Wu, Y.; Du, J.; Shan, H.; et al. Inhibition of BRD4 suppresses the malignancy of breast cancer cells via regulation of Snail. Cell Death Differ. 2020, 27, 255–268. [Google Scholar] [CrossRef] [Green Version]
  10. Zheng, J.; Zhang, Y.; Cai, S.; Dong, L.; Hu, X.; Chen, M.-B.; Zhu, Y.-H. MicroRNA-4651 targets bromodomain-containing protein 4 to inhibit non-small cell lung cancer cell progression. Cancer Lett. 2020, 476, 129–139. [Google Scholar] [CrossRef]
  11. Sanchez, R.; Meslamani, J.; Zhou, M.-M. The bromodomain: From epigenome reader to druggable target. Biochim. et Biophys. Acta 2014, 1839, 676–685. [Google Scholar] [CrossRef] [Green Version]
  12. Hu, J.; Tian, C.-Q.; Damaneh, M.S.; Li, Y.; Cao, D.; Lv, K.; Yu, T.; Meng, T.; Chen, D.; Wang, X.; et al. Structure-based discovery and development of a series of potent and selective bromodomain and extra-terminal protein inhibitors. J. Med. Chem. 2019, 62, 8642–8663. [Google Scholar] [CrossRef] [PubMed]
  13. Tian, X.; Cai, J.; Ma, S.; Fang, Y.; Huang, H.; Lin, T.; Rao, H.; Li, M.; Xia, Z.; Kang, T.; et al. BRD2 induces drug resistance through activation of the RasGRP1/Ras/ERK signaling pathway in adult T-cell lymphoblastic lymphoma. Cancer Commun. 2020, 40, 245–259. [Google Scholar] [CrossRef] [PubMed]
  14. Daneshvar, K.; Ardehali, M.B.; Klein, I.A.; Hsieh, F.-K.; Kratkiewicz, A.J.; Mahpour, A.; Cancelliere, S.O.L.; Zhou, C.; Cook, B.M.; Li, W.; et al. lncRNA DIGIT and BRD3 protein form phase-separated condensates to regulate endoderm differentiation. Nat. Cell Biol. 2020, 22, 1211–1222. [Google Scholar] [CrossRef]
  15. Donati, B.; Lorenzini, E.; Ciarrocchi, A. BRD4 and cancer: Going beyond transcriptional regulation. Mol. Cancer 2018, 17, 164. [Google Scholar] [CrossRef] [PubMed]
  16. Chen, L.; Cai, S.; Wang, J.-M.; Huai, Y.-Y.; Lu, P.-H.; Chu, Q. BRDT promotes ovarian cancer cell growth. Cell Death Dis. 2020, 11, 1021. [Google Scholar] [CrossRef]
  17. Filippakopoulos, P.; Stefan, K. Targeting bromodomains: Epigenetic readers of lysine acetylation. Nat. Rev. Drug Discov. 2014, 13, 337–356. [Google Scholar] [CrossRef] [PubMed]
  18. Müller, A.-C.; Eckert, F.; Heinrich, V.; Bamberg, M.; Brucker, S.; Hehr, T. Resurgery and chest wall re-irradiation for recurrent breast cancer: A second curative approach. BMC Cancer 2011, 11, 197. [Google Scholar] [CrossRef] [Green Version]
  19. Filippakopoulos, P.; Qi, J.; Picaud, S.; Shen, Y.; Smith, W.B.; Fedorov, O.; Morse, E.M.; Keates, T.; Hickman, T.T.; Felletar, I.; et al. Selective inhibition of BET bromodomains. Nature 2010, 468, 1067–1073. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  20. Noel, J.K.; Iwata, K.; Ooike, S.; Sugahara, K.; Nakamura, H.; Daibata, M. Abstract C244: Development of the BET bromodomain inhibitor OTX015. Mol. Cancer Ther. 2013, 12, C244. [Google Scholar] [CrossRef]
  21. Nicodeme, E.; Jeffrey, K.; Schaefer, U.; Beinke, S.; Dewell, S.; Chung, C.-W.; Chandwani, R.; Marazzi, I.; Wilson, P.; Coste, H.; et al. Suppression of inflammation by a synthetic histone mimic. Nature 2010, 468, 1119–1123. [Google Scholar] [CrossRef]
  22. Albrecht, B.K.; Gehling, V.S.; Hewitt, M.C.; Vaswani, R.G.; Côté, A.; Leblanc, Y.; Nasveschuk, C.G.; Bellon, S.; Bergeron, L.; Campbell, R.; et al. Identification of a benzoisoxazoloazepine inhibitor (CPI-0610) of the bromodomain and extra-terminal (BET) family as a candidate for human clinical trials. J. Med. Chem. 2016, 59, 1330–1339. [Google Scholar] [CrossRef] [Green Version]
  23. Bui, M.H.; Lin, X.; Albert, D.H.; Li, L.; Lam, L.T.; Faivre, E.J.; Warder, S.E.; Huang, X.; Wilcox, D.; Donawho, C.K.; et al. Preclinical characterization of BET family bromodomain inhibitor ABBV-075 suggests combination therapeutic strategies. Cancer Res. 2017, 77, 2976–2989. [Google Scholar] [CrossRef] [Green Version]
  24. Hu, J.; Wang, Y.; Li, Y.; Cao, D.; Xu, L.; Song, S.; Damaneh, M.S.; Li, J.; Chen, Y.; Wang, X.; et al. Structure-based optimization of a series of selective BET inhibitors containing aniline or indoline groups. Eur. J. Med. Chem. 2018, 150, 156–175. [Google Scholar] [CrossRef] [PubMed]
  25. Hu, J.; Wang, Y.; Li, Y.; Xu, L.; Cao, D.; Song, S.; Damaneh, M.S.; Wang, X.; Meng, T.; Chen, Y.-L.; et al. Discovery of a series of dihydroquinoxalin-2(1H)-ones as selective BET inhibitors from a dual PLK1-BRD4 inhibitor. Eur. J. Med. Chem. 2017, 137, 176–195. [Google Scholar] [CrossRef]
  26. Watts, E.; Heidenreich, D.; Tucker, E.; Raab, M.; Strebhardt, K.; Chesler, L.; Knapp, S.; Bellenie, B.; Hoelder, S. Designing dual inhibitors of anaplastic lymphoma kinase (ALK) and bromodomain-4 (BRD4) by tuning kinase selectivity. J. Med. Chem. 2019, 62, 2618–2637. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Liu, S.; Yosief, H.O.; Dai, L.; Huang, H.; Dhawan, G.; Zhang, X.; Muthengi, A.M.; Roberts, J.; Buckley, D.L.; Perry, J.A.; et al. Structure-guided design and development of potent and selective dual bromodomain 4 (BRD4)/polo-like kinase 1 (PLK1) inhibitors. J. Med. Chem. 2018, 61, 7785–7795. [Google Scholar] [CrossRef]
  28. Prieto-Martínez, F.D.; Medina-Franco, J.L. Flavonoids as Putative Epi-Modulators: Insight into Their Binding Mode with BRD4 Bromodomains Using Molecular Docking and Dynamics. Biomolecules 2018, 8, 61. [Google Scholar] [CrossRef] [Green Version]
  29. Kadir, M.A.; Abdul Razak, F.I.; Haris, N.S.H. Experimental and DFT data of p-chlorocalix[4]arene as drugs receptor. Data Brief 2020, 32, 106263. [Google Scholar] [CrossRef]
  30. Tirado-Rives, J.; Jorgensen, W.L. Performance of B3LYP density functional methods for a large set of organic molecules. J. Chem. Theory Comput. 2008, 4, 297–306. [Google Scholar] [CrossRef] [PubMed]
  31. Tahara, T.; Watanabe, A.; Yutani, M.; Yamano, Y.; Sagara, M.; Nagai, S.; Saito, K.; Yamashita, M.; Ihara, M.; Iida, A. STAT3 inhibitory activity of naphthoquinones isolated from Tabebuia avellanedae. Bioorg. Med. Chem. 2020, 28, 115347. [Google Scholar] [CrossRef] [PubMed]
  32. Su, J.; Liu, X.; Zhang, S.; Yan, F.; Zhang, Q.; Chen, J. Insight into selective mechanism of class of I-BRD9 inhibitors toward BRD9 based on molecular dynamics simulations. Chem. Biol. Drug Des. 2019, 93, 163–176. [Google Scholar] [CrossRef] [PubMed]
  33. Su, J.; Liu, X.; Zhang, S.; Yan, F.; Zhang, Q.; Chen, J. A computational insight into binding modes of inhibitors XD29, XD35, and XD28 to bromodomain-containing protein 4 based on molecular dynamics simulations. J. Biomol. Struct. Dyn. 2018, 36, 1212–1224. [Google Scholar] [CrossRef]
  34. Wu, S.L.; Wang, L.F.; Sun, H.B.; Wang, W.; Yu, Y.X. Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis. SAR QSAR Environ. Res. 2020, 31, 547–570. [Google Scholar] [CrossRef] [PubMed]
  35. Tumdam, R.; Kumar, A.; Subbarao, N.; Balaji, B.S. In silico study directed towards identification of novel high-affinity inhibitors targeting an oncogenic protein: BRD4-BD1. SAR QSAR Environ. Res. 2018, 29, 975–996. [Google Scholar] [CrossRef]
  36. Zhu, M.; Song, L.T.; Liu, R.R.; Zhai, H.L.; Meng, Y.J.; Ren, C.L. Selective inhibition mechanism of nitroxoline to the BET family: Insight from molecular simulations. Life Sci. 2021, 270, 119141. [Google Scholar] [CrossRef]
  37. Zhong, H.; Wang, Z.; Wang, X.; Liu, H.; Li, D.; Liu, H.; Yao, X.; Hou, T. Importance of a crystalline water network in docking-based virtual screening: A case study of BRD4. Phys. Chem. Chem. Phys. 2019, 21, 25276–25289. [Google Scholar] [CrossRef]
  38. Mathur, S.; Deshmukh, P.; Tripathi, S.; Marimuthu, P.; Padmanabhan, B. Insights into the crystal structure of BRD2-BD2—Phenanthridinone complex and theoretical studies on phenanthridinone analogs. J. Biomol. Struct. Dyn. 2018, 36, 2342–2360. [Google Scholar] [CrossRef]
  39. Weber, O.C.; Uversky, V.N. How accurate are your simulations? Effects of confined aqueous volume and AMBER FF99SB and CHARMM22/CMAP force field parameters on structural ensembles of intrinsically disordered proteins: Amyloid-β42 in water. Intrinsically Disord. Proteins 2017, 5, e1377813. [Google Scholar] [CrossRef] [Green Version]
  40. Janowski, P.A.; Liu, C.; Deckman, J.; Case, D.A. Molecular dynamics simulation of triclinic lysozyme in a crystal lattice. Protein Sci. A Publ. Protein Soc. 2016, 25, 87–102. [Google Scholar] [CrossRef] [Green Version]
  41. Nguyen, H.T.; Pabit, S.A.; Pollack, L.; Case, D.A. Extracting water and ion distributions from solution x-ray scattering experiments. J. Chem. Phys. 2016, 144, 214105. [Google Scholar] [CrossRef] [Green Version]
  42. Linse, B.; Linse, P. Tuning the smooth particle mesh Ewald sum: Application on ionic solutions and dipolar fluids. J. Chem. Phys. 2014, 141, 184114. [Google Scholar] [CrossRef] [PubMed]
  43. Miyamoto, S.; Kollman, P.A. Absolute and relative binding free energy calculations of the interaction of biotin and its analogs with streptavidin using molecular dynamics/free energy perturbation approaches. Proteins 1993, 16, 226–245. [Google Scholar] [CrossRef] [PubMed]
  44. Hijazi, M.; Wilkins, D.M.; Ceriotti, M. Fast-forward Langevin dynamics with momentum flips. J. Chem. Phys. 2018, 148, 184109. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. David, C.C.; Jacobs, D.J. Principal component analysis: A method for determining the essential dynamics of proteins. Methods Mol. Biol. 2014, 1084, 193–226. [Google Scholar] [CrossRef] [Green Version]
  46. Kumar, D.T.; Shaikh, N.; Kumar, S.U.; Doss, C.G.P.; Zayed, H. Structure-Based Virtual Screening to Identify Novel Potential Compound as an Alternative to Remdesivir to Overcome the RdRp Protein Mutations in SARS-CoV-2. Front. Mol. Biosci. 2021, 8, 645216. [Google Scholar] [CrossRef]
  47. Liu, Y.; Yu, Z.; Zhu, J.; Wang, S.; Xu, D.; Han, W. Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study. Front. Chem. 2018, 6, 223. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Liu, Y.; Zhu, J.; Guo, X.; Huang, T.; Han, J.; Gao, J.; Xu, D.; Han, W. How oncogenic mutations activate human MAP kinase 1 (MEK1): A molecular dynamics simulation study. J. Biomol. Struct. Dyn. 2020, 38, 3942–3958. [Google Scholar] [CrossRef]
  49. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
  50. Rahnamoun, A.; Kaymak, M.C.; Manathunga, M.; Götz, A.W.; van Duin, A.C.T.; Merz, K.M., Jr.; Aktulga, H.M. ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations. J. Chem. Theory Comput. 2020, 16, 7645–7654. [Google Scholar] [CrossRef]
Molecules 27 00118 g001 550
Figure 1. The structure of BRD4 protein and the binding pocket of BRD4.
Figure 1. The structure of BRD4 protein and the binding pocket of BRD4.
Molecules 27 00118 g001
Molecules 27 00118 g002 550
Figure 2. The binding pocket of (a) the active residues around Cpd9 binding to BRD4. (b) the active residues around Cpd19 binding to BRD4.
Figure 2. The binding pocket of (a) the active residues around Cpd9 binding to BRD4. (b) the active residues around Cpd19 binding to BRD4.
Molecules 27 00118 g002
Molecules 27 00118 g003 550
Figure 3. (a) Root-mean-square deviation (RMSD) for the backbone atoms, (b) the average RMSD of the four systems, (c) Radius-of-gyration (Rg) plot and (d) the average Rg plot, (e) Solvent-accessible-surface-area (SASA) plot (f) SASA average plot. The free-BRD4 is represented in black, BRD4-Cpd4 is represented in red, BRD4-Cpd9 is represented in blue and BRD4-Cpd19 is represented in green.
Figure 3. (a) Root-mean-square deviation (RMSD) for the backbone atoms, (b) the average RMSD of the four systems, (c) Radius-of-gyration (Rg) plot and (d) the average Rg plot, (e) Solvent-accessible-surface-area (SASA) plot (f) SASA average plot. The free-BRD4 is represented in black, BRD4-Cpd4 is represented in red, BRD4-Cpd9 is represented in blue and BRD4-Cpd19 is represented in green.
Molecules 27 00118 g003
Molecules 27 00118 g004 550
Figure 4. RMSF for the backbone atoms of the four systems. The free-BRD4 is represented in black, BRD4-Cpd4 is represented in red, BRD4-Cpd9 is represented in blue and BRD4-Cpd19 is represented in green.
Figure 4. RMSF for the backbone atoms of the four systems. The free-BRD4 is represented in black, BRD4-Cpd4 is represented in red, BRD4-Cpd9 is represented in blue and BRD4-Cpd19 is represented in green.
Molecules 27 00118 g004
Molecules 27 00118 g005 550
Figure 5. (a) Differences in the secondary structures of residues Trp81 to Ala89 of the four systems (b) Structure of Trp81 to Ala89.
Figure 5. (a) Differences in the secondary structures of residues Trp81 to Ala89 of the four systems (b) Structure of Trp81 to Ala89.
Molecules 27 00118 g005
Molecules 27 00118 g006 550
Figure 6. Cross-correlation matrix of the fluctuations of each of the x, y and z coordinates of the Cα atoms from their average during 150 ns MD for (a) the Free-BRD4 with the prosthetic group, (b)BRD4-Cpd4, (c)BRD4-Cpd9, and (d) BRD4-Cpd19.
Figure 6. Cross-correlation matrix of the fluctuations of each of the x, y and z coordinates of the Cα atoms from their average during 150 ns MD for (a) the Free-BRD4 with the prosthetic group, (b)BRD4-Cpd4, (c)BRD4-Cpd9, and (d) BRD4-Cpd19.
Molecules 27 00118 g006
Molecules 27 00118 g007 550
Figure 7. Motions of the four systems (a) Free-BRD4, (b) BRD4-Cpd4, (c) BRD4-Cpd9, and (d) BRD4-Cpd19. A thick red line represents a high B value and a thin blue line represents a low B value.
Figure 7. Motions of the four systems (a) Free-BRD4, (b) BRD4-Cpd4, (c) BRD4-Cpd9, and (d) BRD4-Cpd19. A thick red line represents a high B value and a thin blue line represents a low B value.
Molecules 27 00118 g007
Molecules 27 00118 g008 550
Figure 8. The cluster analysis for (a) Free-BRD4, (b) BRD4-Cpd4, (c) BRD4-Cpd9 and (d) BRD4-Cpd19.
Figure 8. The cluster analysis for (a) Free-BRD4, (b) BRD4-Cpd4, (c) BRD4-Cpd9 and (d) BRD4-Cpd19.
Molecules 27 00118 g008
Molecules 27 00118 g009 550
Figure 9. Subnetwork analysis of the protein−ligand interaction. (a) Subnetwork betweenBRD4 and Cpd4; (b) subnetwork between BRD4 and Cpd9; (c) subnetwork between BRD4 and Cpd19.
Figure 9. Subnetwork analysis of the protein−ligand interaction. (a) Subnetwork betweenBRD4 and Cpd4; (b) subnetwork between BRD4 and Cpd9; (c) subnetwork between BRD4 and Cpd19.
Molecules 27 00118 g009
Molecules 27 00118 g010 550
Figure 10. Computational alanine scanning of the active site residues in three complexes. (a) BRD4-Cpd4; (b) BRD4-Cpd9; (c) BRD4-Cpd19.
Figure 10. Computational alanine scanning of the active site residues in three complexes. (a) BRD4-Cpd4; (b) BRD4-Cpd9; (c) BRD4-Cpd19.
Molecules 27 00118 g010
Table
Table 1. Four systems in 150 ns molecular dynamics simulations.
Table 1. Four systems in 150 ns molecular dynamics simulations.
SystemProteinLigand
Free-BRD4-BD1BRD4-BD1None
BRD4-BD1-Cpd4BRD4-BD1 Molecules 27 00118 i001
BRD4-BD1-Cpd9BRD4-BD1 Molecules 27 00118 i002
BRD4-BD1-Cpd19BRD4-BD1 Molecules 27 00118 i003
Table
Table 2. The probability of the α-helix in the four systems.
Table 2. The probability of the α-helix in the four systems.
Residue No.Free-BRD4BRD4-Cpd4BRD4-Cpd9BRD4-Cpd19
W8114.9414.0913.2033.25
P8214.9414.0913.2033.25
F8314.9414.0913.2033.25
Table
Table 3. MM-GBSA results (kcal/mol).
Table 3. MM-GBSA results (kcal/mol).
SystemBRD4-Cpd4BRD4-Cpd9BRD4-Cpd19
△Gvdw−32.94 ± 5.14−38.17 ± 4.10−46.04 ± 6.40
△Gele−12.47 ± 6.03−13.72 ± 5.40−9.58 ± 8.74
△Gpolar24.22 ± 4.8727.73 ± 4.9728.64 ± 8.62
△Gnonpolar−4.12 ± 0.48−4.69 ± 0.47−5.48 ± 0.76
△Ggas−45.41 ± 9.18−51.89 ± 7.81−55.62 ± 11.58
△Gsolv20.11 ± 4.6323.04 ± 4.7423.15 ± 8.40
△Gtotal−25.31 ± 5.48−28.85 ± 4.23−32.46 ± 5.60
IC50303.0 nM142.0 ± 3.0 nM5.3 ± 0.4 nM
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Chen, S.; He, Y.; Geng, Y.; Wang, Z.; Han, L.; Han, W. Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors. Molecules 2022, 27, 118. https://doi.org/10.3390/molecules27010118

AMA Style

Chen S, He Y, Geng Y, Wang Z, Han L, Han W. Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors. Molecules. 2022; 27(1):118. https://doi.org/10.3390/molecules27010118

Chicago/Turabian Style

Chen, Siao, Yi He, Yajiao Geng, Zhi Wang, Lu Han, and Weiwei Han. 2022. "Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors" Molecules 27, no. 1: 118. https://doi.org/10.3390/molecules27010118

APA Style

Chen, S., He, Y., Geng, Y., Wang, Z., Han, L., & Han, W. (2022). Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors. Molecules, 27(1), 118. https://doi.org/10.3390/molecules27010118

Article Metrics

Back to TopTop