Insight into selective mechanism of class of I-BRD9 inhibitors toward BRD9 based on molecular dynamics simulations

Discovery of a Highly Potent and Selective BRD9 PROTAC Degrader Based on E3 Binder Investigation for the Treatment of Hematological Tumors

BRD9 is a pivotal epigenetic factor involved in cancers and inflammatory diseases. Still, the limited selectivity and poor phenotypic activity of targeted agents make it an atypically undruggable target. PROTAC offers an alternative strategy for overcoming the issue. In this study, we explored diverse E3 ligase ligands for the contribution of BRD9 PROTAC degradation. Through molecular docking, binding affinity analysis, and structure-activity relationship study, we identified a highly potent PROTAC E5, with excellent BRD9 degradation (DC50 = 16 pM) and antiproliferation in MV4-11 cells (IC50 = 0.27 nM) and OCI-LY10 cells (IC50 = 1.04 nM). E5 can selectively degrade BRD9 and induce cell cycle arrest and apoptosis. Moreover, the therapeutic efficacy of E5 was confirmed in xenograft tumor models, accompanied by further RNA-seq analysis. Therefore, these results may pave the way and provide the reference for the discovery and investigation of highly effective PROTAC degraders.

Binding Mechanism of Inhibitors to BRD4 and BRD9 Decoded by Multiple Independent Molecular Dynamics Simulations and Deep Learning

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.

Pan-cancer analyses of bromodomain containing 9 as a novel therapeutic target reveals its diagnostic, prognostic potential and biological mechanism in human tumours

BACKGROUND

Mutations in one or more genes responsible for encoding subunits within the SWItch/Sucrose Non-Fermentable (SWI/SNF) chromatin-remodelling complexes are found in approximately 25% of cancer patients. Bromodomain containing 9 (BRD9) is a more recently identified protein coding gene, which can encode SWI/SNF chromatin-remodelling complexes subunits. Although initial evaluations of the potential of BRD9-based targeted therapy have been explored in the clinical application of a small number of cancer types, more detailed study of the diagnostic and prognostic potential, as well as the detailed biological mechanism of BRD9 remains unreported.

METHODS

We used various bioinformatics tools to generate a comprehensive, pan-cancer analyses of BRD9 expression in multiple disease types described in The Cancer Genome Atlas (TCGA). Experimental validation was conducted in tissue microarrays and cell lines derived from lung and colon cancers.

RESULTS

Our study revealed that BRD9 exhibited elevated expression in a wide range of tumours. Analysis of survival data and DNA methylation for BRD9 indicated distinct conclusions for multiple tumours. mRNA splicing and molecular binding were involved in the functional mechanism of BRD9. BRD9 may affect cancer progression through different phosphorylation sites or N6 -methyladenosine site modifications. BRD9 could potentially serve as a novel biomarker for diagnosing different cancer types, especially could accurately forecast the prognosis of melanoma patients receiving anti-programmed cell death 1 immunotherapy. BRD9 has the potential to serve as a therapeutic target, when pairing with etoposide in patients with melanoma. The BRD9/SMARCD1 axis exhibited promising discriminative performance in forecasting the prognosis of patients afflicted with liver hepatocellular carcinoma (LIHC) and mesothelioma. Additionally, this axis appears to potentially influence the immune response in LIHC by regulating the programmed death-ligand 1 immune checkpoint. For experimental validation, high expression levels of BRD9 were observed in tumour tissue samples from both lung and colon cancer patients. Knocking down BRD9 led to the inhibition of lung and colon cancer development, likely via the Wnt/β-catenin signalling pathway.

CONCLUSIONS

These pan-cancer study revealed the diagnostic and prognostic potential, along with the biological mechanism of BRD9 as a novel therapeutic target in human tumours.

Conformational States of the GDP- and GTP-Bound HRAS Affected by A59E and K117R: An Exploration from Gaussian Accelerated Molecular Dynamics

The HRAS protein is considered a critical target for drug development in cancers. It is vital for effective drug development to understand the effects of mutations on the binding of GTP and GDP to HRAS. We conducted Gaussian accelerated molecular dynamics (GaMD) simulations and free energy landscape (FEL) calculations to investigate the impacts of two mutations (A59E and K117R) on GTP and GDP binding and the conformational states of the switch domain. Our findings demonstrate that these mutations not only modify the flexibility of the switch domains, but also affect the correlated motions of these domains. Furthermore, the mutations significantly disrupt the dynamic behavior of the switch domains, leading to a conformational change in HRAS. Additionally, these mutations significantly impact the switch domain's interactions, including their hydrogen bonding with ligands and electrostatic interactions with magnesium ions. Since the switch domains are crucial for the binding of HRAS to effectors, any alterations in their interactions or conformational states will undoubtedly disrupt the activity of HRAS. This research provides valuable information for the design of drugs targeting HRAS.

Targeting epigenetic deregulations for the management of esophageal carcinoma: recent advances and emerging approaches

Ranking from seventh in incidence to sixth in mortality, esophageal carcinoma is considered a severe malignancy of food pipe. Later-stage diagnosis, drug resistance, and a high mortality rate contribute to its lethality. Esophageal squamous cell carcinoma and esophageal adenocarcinoma are the two main histological subtypes of esophageal carcinoma, with squamous cell carcinoma alone accounting for more than eighty percent of its cases. While genetic anomalies are well known in esophageal cancer, accountability of epigenetic deregulations is also being explored for the recent two decades. DNA methylation, histone modifications, and functional non-coding RNAs are the crucial epigenetic players involved in the modulation of different malignancies, including esophageal carcinoma. Targeting these epigenetic aberrations will provide new insights into the development of biomarker tools for risk stratification, early diagnosis, and effective therapeutic intervention. This review discusses different epigenetic alterations, emphasizing the most significant developments in esophageal cancer epigenetics and their potential implication for the detection, prognosis, and treatment of esophageal carcinoma. Further, the preclinical and clinical status of various epigenetic drugs has also been reviewed.

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

Binding mechanism of anacardic acid, carnosol and garcinol with PCAF: A comprehensive study using molecular docking and molecular dynamics simulations and binding free energy analysis

The p300/CBP associated factor bromodomain (PCAF Brd) is emerged as one of the promising target proteins for different types of cancers. PCAF is one among the histone acetyltransferase enzymes which involved in the regulation of transcriptase process by modifying the chromatin structure. Anacardic acid, carnosol, garcinol are the experimentally reported inhibitors of PCAF Brd; however, their detailed binding mechanism these inhibitors are not yet known. The intermolecular interaction, binding energy, and the stability of these inhibitors with the active site of PCAF Brd are playing the key role in the binding of these inhibitors with PCAF. The in silico study incorporates the molecular docking and dynamics simulations; these molecular level simulations allow to understand the binding mechanism. In the present study, the induced fit molecular docking and molecular dynamics of anacardic acid, carnosol and garcinol molecules against the PCAF Brd have been performed. The docking score values of these molecules are -5.112 (anacardic acid), -5.141 (carnosol), -5.199 (garcinol) and -3.641 (L45) kcal/mol, respectively. Further, the molecular dynamics simulation was carried out for these docked complexes to understand their conformational their stability and binding energy from the roots means square deviation (RMSD) and root means square of fluctuation (RMSF), and molecular mechanics with the generalized born and surface area solvation (MM/GBSA) binding free energy calculations. The intermolecular interactions and binding free energy values confirm that garcinol forms key interactions and has high binding affinity towards PCAF Brd on compare with the other two inhibitors. Therefore, garcinol may be considered as a potential inhibitor of PCAF Brd.

Impacts of Mutations in the P-Loop on Conformational Alterations of KRAS Investigated with Gaussian Accelerated Molecular Dynamics Simulations

G12 mutations heavily affect conformational transformation and activity of KRAS. In this study, Gaussian accelerated molecular dynamics (GaMD) simulations were performed on the GDP-bound wild-type (WT), G12A, G12D, and G12R KRAS to probe mutation-mediated impacts on conformational alterations of KRAS. The results indicate that three G12 mutations obviously affect the structural flexibility and internal dynamics of the switch domains. The analyses of the free energy landscapes (FELs) suggest that three G12 mutations induce more conformational states of KRAS and lead to more disordered switch domains. The principal component analysis shows that three G12 mutations change concerted motions and dynamics behavior of the switch domains. The switch domains mostly overlap with the binding region of KRAS to its effectors. Thus, the high disorder states and concerted motion changes of the switch domains induced by G12 mutations affect the activity of KRAS. The analysis of interaction network of GDP with KRAS signifies that the instability in the interactions of GDP and magnesium ion with the switch domain SW1 drives the high disordered state of the switch domains. This work is expected to provide theoretical aids for understanding the function of KRAS.

Deciphering Selectivity Mechanism of BRD9 and TAF1(2) toward Inhibitors Based on Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

BRD9 and TAF1(2) have been regarded as significant targets of drug design for clinically treating acute myeloid leukemia, malignancies, and inflammatory diseases. In this study, multiple short molecular dynamics simulations combined with the molecular mechanics generalized Born surface area method were employed to investigate the binding selectivity of three ligands, 67B, 67C, and 69G, to BRD9/TAF1(2) with IC50 values of 230/59 nM, 1400/46 nM, and 160/410 nM, respectively. The computed binding free energies from the MM-GBSA method displayed good correlations with that provided by the experimental data. The results indicate that the enthalpic contributions played a critical factor in the selectivity recognition of inhibitors toward BRD9 and TAF1(2), indicating that 67B and 67C could more favorably bind to TAF1(2) than BRD9, while 69G had better selectivity toward BRD9 over TAF1(2). In addition, the residue-based free energy decomposition approach was adopted to calculate the inhibitor–residue interaction spectrum, and the results determined the gatekeeper (Y106 in BRD9 and Y1589 in TAF1(2)) and lipophilic shelf (G43, F44, and F45 in BRD9 and W1526, P1527, and F1528 in TAF1(2)), which could be identified as hotspots for designing efficient selective inhibitors toward BRD9 and TAF1(2). This work is also expected to provide significant theoretical guidance and insightful molecular mechanisms for the rational designs of efficient selective inhibitors targeting BRD9 and TAF1(2).

BRD9 Inhibition by Natural Polyphenols Targets DNA Damage/Repair and Apoptosis in Human Colon Cancer Cells

Epigenetic mechanisms play an important role in the etiology of colorectal cancer (CRC) and other malignancies due, in part, to deregulated bromodomain (BRD) functions. Inhibitors of the bromodomain and extraterminal (BET) family have entered into clinical trials as anticancer agents, and interest has grown in other acetyl 'reader' proteins as therapeutic targets, including non-BET member bromodomain-containing protein 9 (BRD9). We report here that overexpression of BRD9 is associated with poor prognosis in CRC patients, and that siRNA-mediated knockdown of BRD9 decreased cell viability and activated apoptosis in human colon cancer cells, coincident with increased DNA damage. Seeking natural compounds as BRD9 antagonists, molecular docking in silico identified several polyphenols such as Epigallocatechin-3-gallate (EGCG), Equol, Quercetin, and Aspalathin, with favorable binding energies, supported by BROMOscan^® (DiscoverX) and isothermal titration calorimetry experiments. Polyphenols mimicked BRD9 knockdown and iBRD9 treatment in reducing colon cancer cell viability, inhibiting colony formation, and enhancing DNA damage and apoptosis. Normal colonic epithelial cells were unaffected, signifying cancer-specific effects. These findings suggest that natural polyphenols recognize and target BRD9 for inhibition, and might serve as useful lead compounds for bromodomain therapeutics in the clinical setting.

Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations

Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.

Dynamics of the Histone Acetyltransferase Lysine-Rich Loop in the Catalytic Core of the CREB-Binding Protein

The tight control of transcriptional coactivators is a fundamental aspect of gene expression in cells. The regulation of the CREB-binding protein (CBP) and p300 coactivators, two paralog multidomain proteins, involves an autoinhibitory loop (AIL) of the histone acetyltransferase (HAT) domain. There is experimental evidence for the AIL engaging with the HAT binding site, thus interrupting the acetylation of histone tails or other proteins. Both CBP and p300 contain a domain of about 110 residues (called the bromodomain) that recognizes histone tails with one or more acetylated lysine side chains. Here, we investigate by molecular dynamics simulations whether the AIL of CBP (residues 1556-1618) acetylated at the side chain of Lys1595 can bind to the bromodomain. The structural instability and fast unbinding kinetics of the AIL from the bromodomain pocket suggest that the AIL is not a ligand of the bromodomain on the same protein chain. This is further supported by the absence of strong and persistent contacts at the binding interface. Furthermore, the simulations of unbinding show an initial fast detachment of the acetylated lysine and a slower phase necessary for complete AIL dissociation. We provide further evidence for the instability of the AIL intramolecular binding by comparison with a natural ligand, the histone peptide H3K56ac, which shows higher stability in the pocket.

Theoretical exploration of the binding selectivity of inhibitors to BRD7 and BRD9 with multiple short molecular dynamics simulations

Bromodomain-containing proteins 7 and 9 (BRD7 and BRD9) have been considered as potential targets of clinical drug design toward treatment of human cancers and other diseases. Multiple short molecular dynamics simulations and binding free energy predictions were carried out to decipher the binding selectivity of three inhibitors 4L2, 5U6, and 6KT toward BRD7 and BRD9. The results show that 4L2 has more favorable binding ability to BRD7 over BRD9 compared to 5U6 and 6KT, while 5U6 and 6KT possess more favorable associations with BRD9 than BRD7. Furthermore, estimations of residue-based free energy decompositions further identify that four common residue pairs, including (F155, F44), (V160, V49), (Y168, Y57) and (Y217, Y106) in (BRD7, BRD9) generate obvious binding differences with 4L2, 5U6, and 6KT, which mostly drives the binding selectivity of 4L2, 5U6, and 6KT to BRD7 and BRD9. Dynamic information arising from trajectory analysis also suggests that inhibitor bindings affect structural flexibility and motion modes, which is responsible for the partial selectivity of 4L2, 5U6, and 6KT toward BRD7 and BRD9. As per our expectation, this study theoretically provides useful hints for design of dual inhibitors with high selectivity on BRD7 and BRD9.

Bromodomains (BRDs) are structurally conserved epigenetic reader modules observed in numerous chromatin- and transcription-associated proteins that have a capability to identify acetylated lysine residues.

The binding mechanism of NHWD-870 to bromodomain-containing protein 4 based on molecular dynamics simulations and free energy calculation

Bromodomain and extra-terminal (BET) proteins (BRD2, BRD3, BRD4, and BRDT) are epigenetic readers with tandem bromodomains.

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

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.

Insights into the Ligand Binding to Bromodomain-Containing Protein 9 (BRD9): A Guide to the Selection of Potential Binders by Computational Methods

The estimation of the binding of a set of molecules against BRD9 protein was carried out through an in silico molecular dynamics-driven exhaustive analysis to guide the identification of potential novel ligands. Starting from eight crystal structures of this protein co-complexed with known binders and one apo form, we conducted an exhaustive molecular docking/molecular dynamics (MD) investigation. To balance accuracy and an affordable calculation time, the systems were simulated for 100 ns in explicit solvent. Moreover, one complex was simulated for 1 µs to assess the influence of simulation time on the results. A set of MD-derived parameters was computed and compared with molecular docking-derived and experimental data. MM-GBSA and the per-residue interaction energy emerged as the main indicators for the good interaction between the specific binder and the protein counterpart. To assess the performance of the proposed analysis workflow, we tested six molecules featuring different binding affinities for BRD9, obtaining promising outcomes. Further insights were reported to highlight the influence of the starting structure on the molecular dynamics simulations evolution. The data confirmed that a ranking of BRD9 binders using key parameters arising from molecular dynamics is advisable to discard poor ligands before moving on with the synthesis and the biological tests.

Selectively targeting individual bromodomain: Drug discovery and molecular mechanisms

Bromodomain-containing proteins include bromodomain and extra-terminal (BET) and non-BET families. Due to the conserved bromodomain (BD) module between BD-containing proteins, and especially BETs with each member having two BDs (BD1 and BD2), the high degree of structural similarity makes BD-selective inhibitors much difficult to be designed. However, increasing evidences emphasized that individual BDs had distinct functions and different cellular phenotypes after pharmacological inhibition, and selectively targeting one of the BDs could result in a different efficacy and tolerability profile. This review is to summarize the pioneering progress of BD-selective inhibitors targeting BET and non-BET proteins, focusing on their structural features, biological activity, therapeutic application and experimental/theoretical mechanisms. The present proteolysis targeting chimeras (PROTAC) degraders targeting BDs, and clinical status of BD-selective inhibitors were also analyzed, providing a new insight into future direction of bromodomain-selective drug discovery.

Structural variation of protein-ligand complexes of the first bromodomain of BRD4

The bromodomain-containing protein 4 (BRD4), a member of the bromodomain and extra-terminal domain (BET) family, plays a key role in several diseases, especially cancers. With increased interest in BRD4 as a therapeutic target, many X-ray crystal structures of the protein in complex with small molecule inhibitors are publicly available over the recent decade. In this study, we use this structural information to investigate the conformations of the first bromodomain (BD1) of BRD4. Structural alignment of 297 BRD4-BD1 complexes shows a high level of similarity between the structures of BRD4-BD1, regardless of the bound ligand. We employ WONKA, a tool for detailed analyses of protein binding sites, to compare the active site of over 100 of these crystal structures. The positions of key binding site residues show a high level of conformational similarity, with the exception of Trp81. A focused analysis on the highly conserved water network in the binding site of BRD4-BD1 is performed to identify the positions of these water molecules across the crystal structures. The importance of the water network is illustrated using molecular docking and absolute free energy perturbation simulations. 82% of the ligand poses were better predicted when including water molecules as part of the receptor. Our analysis provides guidance for the design of new BRD4-BD1 inhibitors and the selection of the best structure of BRD4-BD1 to use in structure-based drug design, an important approach for faster and more cost-efficient lead discovery.

Binding Selectivity of Inhibitors toward Bromodomains BAZ2A and BAZ2B Uncovered by Multiple Short Molecular Dynamics Simulations and MM-GBSA Calculations

Two Bromodomain-Containing proteins BAZ2A and BAZ2B are responsible for remodeling chromatin and regulating noncoding RNAs. As for our current studies, integration of multiple short molecular dynamics simulations (MSMDSs) with molecular mechanics generalized Born surface area (MM-GBSA) method is adopted for insights into binding selectivity of three small molecules D8Q, D9T and UO1 to BAZ2A against BAZ2B. The calculations of MM-GBSA unveil that selectivity of inhibitors toward BAZ2A and BAZ2B highly depends on the enthalpy changes and the details uncover that D8Q has better selectivity toward BAZ2A than BAZ2B, D9T more favorably bind to BAZ2B than BAZ2A, and UO1 does not show obvious selectivity toward these two proteins. The analysis of interaction network between residues and inhibitors indicates that seven residues are mainly responsible for the selectivity of D8Q, six residues for D9T and four residues provide significant contributions to associations of UO1 with two proteins. Moreover the analysis of interaction network not only reveals warm spots of inhibitor bindings to BAZ2A and BAZ2B but also unveils that common residue pairs, including (W1816, W1887), (P1817, P1888), (F1818, F1889), (V1822, V1893), (N1823, N1894),(L1826, L1897), (V1827, V1898), (F1872, F1943), (N1873, N1944) and (V1879, I1950) belonging to (BAZ2A, BAZ2B), induce mainly binding differences of inhibitors to BAZ2A and BAZ2B. Hence, insights from our current studies offer useful dynamics information relating with conformational alterations and structure-affinity relationship at atomistic levels for novel therapeutic strategies toward BAZ2A and BAZ2B.

Targeting BRD9 for Cancer Treatment: A New Strategy

Bromodomain-containing protein 9 (BRD9) is a newly identified subunit of the non-canonical barrier-to-autointegration factor (ncBAF) complex and a member of the bromodomain family IV. Studies have confirmed that BRD9 plays an oncogenic role in multiple cancer types, by regulating tumor cell growth. The tumor biological functions of BRD9 are mainly due to epigenetic modification mediated by its bromodomain. The bromodomain recruits the ncBAF complex to the promoter to regulate gene transcription. This review summarizes the potential mechanisms of action of BRD9 in carcinogenesis and the emerging strategies for targeting BRD9 for cancer therapeutics. Although the therapeutic potential of BRD9 has been exploited to some extent, research on the detailed biological mechanisms of BRD9 is still in its infancy. Therefore, targeting BRD9 to study its biological roles will be an attractive tool for cancer diagnosis and treatment, but it remains a great challenge.

BRD4/8/9 are prognostic biomarkers and associated with immune infiltrates in hepatocellular carcinoma

Bromodomain (BRD)-containing proteins are a class of epigenetic readers with unique recognition for N-acetyl-lysine in histones and functions of gene transcription and chromatin modification, known to be critical in various cancers. However, little is known about the roles of distinct BRD-containing protein genes in hepatocellular carcinoma (HCC). Most recently, we investigated the transcriptional and survival data of BRD1, BRD2, BRD3, BRD4, BRD7, BRD8, BRD9 in HCC patients through ONCOMINE, UALCAN, Human Protein Atlas, GEPIA, cBioPortal, STRING, TIMER databases. BRD1/2/3/4/7/8/9 were over-expressed in HCC and were significantly associated with clinical cancer stages and pathological tumor grades. High mRNA expressions of BRD4/8/9 were promising candidate biomarkers in HCC patients. The rate of sequence alternations in BRD1/2/3/4/7/8/9 was relatively high (52%) in HCC patients, and the genetic alternations were correlated with shorter overall survival and disease-free survival in HCC patients. Additionally, the mRNA expression levels of individual BRD genes were significantly positively associated with the immune infiltrating levels of B cells, CD8⁺ T cells, CD4⁺ T cells, macrophages, neutrophils, and dendritic cells. And the associations between BRD1/2/3/4/7/8/9 and diverse immune marker sets showed a significance. Overall, these results indicated that BRD4/8/9 could be potential prognostic markers and druggable epigenetic targets in HCC patients.

Hitting a Moving Target: Simulation and Crystallography Study of ATAD2 Bromodomain Blockers

Small molecule ligand binding to the ATAD2 bromodomain is investigated here through the synergistic combination of molecular dynamics and protein crystallography. A previously unexplored conformation of the binding pocket upon rearrangement of the gatekeeper residue Ile1074 has been found. Further, our investigations reveal how minor structural differences in the ligands result in binding with different plasticity of the ZA loop for this difficult-to-drug bromodomain.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

Revealing binding selectivity of inhibitors toward bromodomain-containing proteins 2 and 4 using multiple short molecular dynamics simulations and free energy analyses

Emerging evidences indicate bromodomain-containing proteins 2 and 4 (BRD2 and BRD4) play critical roles in cancers, inflammations, cardiovascular diseases and other pathologies. Multiple short molecular dynamics (MSMD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were applied to investigate the binding selectivity of three inhibitors 87D, 88M and 89G towards BRD2 over BRD4. The root-mean-square fluctuation (RMSF) analysis indicates that the structural flexibility of BRD4 is stronger than that of BRD2. Moreover the calculated distances between the C_α atoms in the centres of the ZA_loop and BC_loop of BRD4 are also bigger than that of BRD2. The rank of binding free energies calculated using MM-GBSA method agrees well with that determined by experimental data. The results show that 87D can bind more favourably to BRD2 than BRD4, while 88M has better selectivity on BRD4 over BRD2. Residue-based free-energy decomposition method was utilized to estimate the inhibitor-residue interaction spectrum and the results not only identify the hot interaction spots of inhibitors with BRD2 and BRD4, but also demonstrate that several common residues, including (W370, W374), (P371, P375), (V376, V380) and (L381, L385) belonging to (BRD2, BRD4), generate significant binding difference of inhibitors to BRD2 and BRD4.

Molecular mechanism of inhibitor bindings to bromodomain-containing protein 9 explored based on molecular dynamics simulations and calculations of binding free energies

Recently, bromodomain-containing protein 9 (BRD9) has been a prospective therapeutic target for anticancer drug design. Molecular dynamics (MD) simulations combined with molecular mechanics generalized Born surface area (MM-GBSA) method were adopted to explore binding modes of three inhibitors (5SW, 5U2, and 5U6) to BRD9 and identify the hot spot of the inhibitor-BRD9 binding. The results indicate that the inhibitor 5SW has the strongest binding ability to BRD9 among the current three inhibitors. Furthermore, the rank of the binding free energies predicted by MM-GBSA approach agrees with that determined by the experimental values. In addition, inhibitor-residue interactions were computed by using residue-based free-energy decomposition method and the results suggest that residue His42 produces the CH-H interactions, residues Asn100, Ile53 and Val49 produce the CH-[Formula: see text] interactions with three inhibitors and Tyr106, Phe45 and Phe44 generate the π-π interactions with inhibitors. Notably, the residue Asn140 forms hydrogen bonding interactions with three inhibitors. This research is expected to provide useful molecular basis and dynamics information at atomic levels for the design of potent inhibitors inhibiting the activity of BRD9.

Epigenetic Targets in Synovial Sarcoma: A Mini-Review

Synovial Sarcomas (SS) are a type of Soft Tissue Sarcoma (STS) and represent 8-10% of all STS cases. Although SS can arise at any age, it typically affects younger individuals aged 15-35 and is therefore part of both pediatric and adult clinical practices. SS occurs primarily in the limbs, often near joints, but can present anywhere. It is characterized by the recurrent pathognomonic chromosomal translocation t(X;18)(p11.2;q11.2) that most frequently fuses SSX1 or SSX2 genes with SS18. This leads to the expression of the SS18-SSX fusion protein, which causes disturbances in several interacting multiprotein complexes such as the SWItch/Sucrose Non-Fermentable (SWI/SNF) complex, also known as the BAF complex and the Polycomb Repressive Complex 1 and 2 (PRC1 and PRC2). Furthermore, this promotes widespread epigenetic rewiring, leading to aberrant gene expression that drives the pathogenesis of SS. Good prognoses are characterized predominantly by small tumor size and young patient age. Whereas, high tumor grade and an increased genomic complexity of the tumor constitute poor prognostic factors. The current therapeutic strategy relies on chemotherapy and radiotherapy, the latter of which can lead to chronic side effects for pediatric patients. We will focus on the known roles of SWI/SNF, PRC1, and PRC2 as the main effectors of the SS18-SSX-mediated genome modifications and we present existing biological rationale for potential therapeutic targets and treatment strategies.

Molecular Mechanism of Binding Selectivity of Inhibitors toward BACE1 and BACE2 Revealed by Multiple Short Molecular Dynamics Simulations and Free-Energy Predictions

The β-amyloid cleaving enzymes 1 and 2 (BACE1 and BACE2) have been regarded as the prospective targets for clinically treating Alzheimer's disease (AD) in the last two decades. Thus, insight into the binding differences of inhibitors to BACE1 and BACE2 is of significance for designing highly selective inhibitors toward the two proteins. In this work, multiple short molecular dynamics (MSMD) simulations are coupled with the molecular mechanics generalized Born surface area (MM-GBSA) method to probe the binding selectivity of three inhibitors DBO, CS9, and SC7 on BACE1 over BACE2. The results show that the entropy effect plays a key role in selectivity identification of inhibitors toward BACE1 and BACE2, which determines that DBO has better selectivity toward BACE2 over BACE1, while CS9 and CS7 can more favorably bind to BACE1 than BACE2. The hierarchical clustering analysis based on energetic contributions of residues suggests that BACE1 and BACE2 share the common hot interaction spots. The residue-based free-energy decomposition method was applied to compute the inhibitor-residue interaction spectrum, and the results recognize four common binding subpockets corresponding to the different groups of inhibitors, which can be used as efficient targets for designing highly selective inhibitors toward BACE1 and BACE2. Therefore, these results provide a useful molecular basis and dynamics information for development of highly selective inhibitors targeting BACE1 and BACE2.