Diving into the Water: Inducible Binding Conformations for BRD4, TAF1(2), BRD9, and CECR2 Bromodomains

Development and comparison of single FLT3-inhibitors to dual FLT3/TAF1-inhibitors as an anti-leukemic approach

Acute myeloid leukemia (AML) is characterized by several recurrent mutations that affect disease biology and phenotype, response to therapy and risk of subsequent relapse. Though tyrosine kinase inhibitors have gained regulatory approval for the treatment of AML, it is unclear whether single drugs targeting a specific genomic alteration will be sufficient to eradicate disease. Fortuitously, kinase/bromodomain inhibitors allow targeting of downstream transcriptional effectors of oncogenic pathways, allowing impediment of drug resistance at the transcriptional level. Successful development of combinatorial therapeutic strategies to inhibit both upstream oncogenic pathways and their downstream effectors could thus impede the onset of resistant disease. By using a combination of high-throughput cell-based screening assays and structure-based design, we have developed a novel anti-proliferative 3i-compound scaffold with a diverse range of single and dual FLT3/TAF1(2) activity against AML. Our novel approach to target both FLT3 kinase and TAF1(2) bromodomain efficiently maintained potency against haematological cancers. However, reference compounds and in vitro cell viability and cytotoxicity assays in cancer cell lines demonstrated superior effects of high affinity tyrosine kinase inhibition compared to inhibition of the TAF1 bromodomain. Our results highlight the feasibility of dual tyrosine kinase-bromodomain targeting to overcome disease mechanisms while also revealing the increased efficacy of FLT3-targeted compounds in AML.

Identification of Novel Bromodomain-Containing Protein 4 (BRD4) Binders through 3D Pharmacophore-Based Repositioning Screening Campaign

A 3D structure-based pharmacophore model built for bromodomain-containing protein 4 (BRD4) is reported here, specifically developed for investigating and identifying the key structural features of the (+)-JQ1 known inhibitor within the BRD4 binding site. Using this pharmacophore model, 273 synthesized and purchased compounds previously considered for other targets but yielding poor results were screened in a drug repositioning campaign. Subsequently, only six compounds showed potential as BRD4 binders and were subjected to further biophysical and biochemical assays. Compounds 2, 5, and 6 showed high affinity for BRD4, with IC50 values of 0.60 ± 0.25 µM, 3.46 ± 1.22 µM, and 4.66 ± 0.52 µM, respectively. Additionally, these compounds were tested against two other bromodomains, BRD3 and BRD9, and two of them showed high selectivity for BRD4. The reported 3D structure-based pharmacophore model proves to be a straightforward and useful tool for selecting novel BRD4 ligands.

Development of Potent Dual BET/HDAC Inhibitors via Pharmacophore Merging and Structure-Guided Optimization

Bromodomain and extra-terminal domain (BET) proteins and histone deacetylases (HDACs) are prime targets in cancer therapy. Recent research has particularly focused on the development of dual BET/HDAC inhibitors for hard-to-treat tumors, such as pancreatic cancer. Here, we developed a new series of potent dual BET/HDAC inhibitors by choosing starting scaffolds that enabled us to optimally merge the two functionalities into a single compound. Systematic structure-guided modification of both warheads then led to optimized binders that were superior in potency to both parent compounds, with the best molecules of this series binding to both BRD4 bromodomains as well as HDAC1/2 with EC50 values in the 100 nM range in cellular NanoBRET target engagement assays. For one of our lead molecules, we could also show the selective inhibition of HDAC1/2 over all other zinc-dependent HDACs. Importantly, this on-target activity translated into promising efficacy in pancreatic cancer and NUT midline carcinoma cells. Our lead molecules effectively blocked histone H3 deacetylation in pancreatic cancer cells and upregulated the tumor suppressor HEXIM1 and proapoptotic p57, both markers of BET inhibition. In addition, they have the potential to downregulate the oncogenic drivers of NUT midline carcinoma, as demonstrated for MYC and TP63 mRNA levels. Overall, this study expands the portfolio of available dual BET/class I HDAC inhibitors for future translational studies in different cancer models.

Exploring the chemical space of functionalized [1,2,4]triazolo[4,3-a]quinoxaline-based compounds targeting the bromodomain of BRD9

Here we report a detailed structure-activity relationship (SAR) study related to [1,2,4]triazolo[4,3-a]quinoxaline-based compounds targeting the reader module of bromodomain containing-protein 9 (BRD9). 3D structure-based pharmacophore models, previously introduced by us, were here employed to evaluate a second generation of compounds, exploring different substitution patterns on the heterocyclic core. Starting from the promising data obtained from our previously identified [1,2,4]triazolo[4,3-a]quinoxaline-based compounds 1-4, the combination of in silico studies, chemical synthesis, biophysical and in vitro assays led to the identification of a new set of derivatives, selected for thoroughly exploring the chemical space of the bromodomain binding site. In more details, the investigation of different linkers at C-4 position highlighted the amine spacer as mandatory for the binding with the protein counterpart and the crucial role of the alkyl substituents at C-1 for increasing the selectivity toward BRD9. Additionally, the importance of a hydrogen bond donor group, critical to anchor the ZA region and required for the interaction with Ile53 residue, was inferred from the analysis of our collected results. Herein we also propose an optimization and an update of our previously reported "pharm-druglike2" 3D structure-based pharmacophore model, introducing it as "pharm-druglike2.1". Compounds 24-26, 32, 34 and 36 were identified as new valuable BRD9 binders featuring IC50 values in the low micromolar range. Among them, 24 and 36 displayed an excellent selectivity towards BRD9 and a good antiproliferative effect on a panel of leukemia models, especially toward CCRF-CEM cell line, with no cytotoxicity on healthy cells. Notably, the interaction of 24 and 36 with the bromodomain and PHD finger-containing protein 1 (BRPF1) also emerged, disclosing them as new and unexplored dual inhibitors for these two proteins highly involved in leukemia. These findings highlight the potential for the identification of new attractive dual epidrugs as well as a promising starting point for the development of chemical degraders endowed with anticancer activities.

GNE-235: A Lead Compound Selective for the Second Bromodomain of PBRM1

Bromodomains are acetyl-lysine binding modules that are found in different classes of chromatin-interacting proteins. Among these are large chromatin remodeling complexes such as BAF and PBAF (variants of human SWI/SNF). Previous work has identified chemical probes targeting a subset of the bromodomains present in the BAF and PBAF complexes. Selective inhibitors of the individual bromodomains have proven challenging to discover, as the domains are highly similar. Here, elaboration of an aminopyridazine scaffold used previously to develop probes for the bromodomains of SMARCA2, SMARCA4, and the fifth bromodomain of PBRM1 yielded compounds with both potency and unusual selectivity for the second bromodomain of PBRM1. One of these, GNE-235, and its enantiomer control GNE-234 are suggested for initial cellular investigations of the function of the second bromodomain of PBRM1.

Reevaluation of bromodomain ligands targeting BAZ2A

BAZ2A promotes migration and invasion in prostate cancer. Two chemical probes, the specific BAZ2-ICR, and the BAZ2/BRD9 cross-reactive GSK2801, interfere with the recognition of acetylated lysines in histones by the bromodomains of BAZ2A and of its BAZ2B paralog. The two chemical probes were tested in prostate cancer cell lines with opposite androgen susceptibility. BAZ2-ICR and GSK2801 showed different cellular efficacies in accordance with their unequal selectivity profiles. Concurrent inhibition of BAZ2 and BRD9 did not reproduce the effects observed with GSK2801, indicating possible off-targets for this chemical probe. On the other hand, the single BAZ2 inhibition by BAZ2-ICR did not phenocopy genetic ablation, demonstrating that bromodomain interference is not sufficient to strongly affect BAZ2A functionality and suggesting a PROTAC-based chemical ablation as an alternative optimization strategy and a possible therapeutic approach. In this context, we also present the crystallographic structures of BAZ2A in complex with the above chemical probes. Binding poses of TP-238 and GSK4027, chemical probes for the bromodomain subfamily I, and two ligands of the CBP/EP300 bromodomains identify additional headgroups for the development of BAZ2A ligands.

TAF1 bromodomain inhibition as a candidate epigenetic driver of congenital heart disease

Heart formation requires transcriptional regulators that underlie congenital anomalies and the fetal gene program activated during heart failure. Attributing the effects of congenital heart disease (CHD) missense variants to disruption of specific protein domains allows for a mechanistic understanding of CHDs and improved diagnostics. A combined chemical and genetic approach was employed to identify novel CHD drivers, consisting of chemical screening during pluripotent stem cell (PSC) differentiation, gene expression analyses of native tissues and primary cell culture models, and the in vitro study of damaging missense variants from CHD patients. An epigenetic inhibitor of the TATA-Box Binding Protein Associated Factor 1 (TAF1) bromodomain was uncovered in an unbiased chemical screen for activators of atrial and ventricular fetal myosins in differentiating PSCs, leading to the development of a high affinity inhibitor (5.1 nM) of the TAF1 bromodomain, a component of the TFIID complex. TAF1 bromodomain inhibitors were tested for their effects on stem cell viability and cardiomyocyte differentiation, implicating a role for TAF1 in cardiogenesis. Damaging TAF1 missense variants from CHD patients were studied by mutational analysis of the TAF1 bromodomain, demonstrating a repressive role of TAF1 that can be abrogated by the introduction of damaging bromodomain variants or chemical TAF1 bromodomain inhibition. These results indicate that targeting the TAF1/TFIID complex with chemical compounds modulates cardiac transcription and identify an epigenetically-driven CHD mechanism due to damaging variants within the TAF1 bromodomain.

Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities

Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.

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).

E3 Ligases Meet Their Match: Fragment-Based Approaches to Discover New E3 Ligands and to Unravel E3 Biology

Ubiquitination is a key post-translational modification of proteins, affecting the regulation of multiple cellular processes. Cells are equipped with over 600 ubiquitin orchestrators, called E3 ubiquitin ligases, responsible for directing the covalent attachment of ubiquitin to substrate proteins. Due to their regulatory role in cells, significant efforts have been made to discover ligands for E3 ligases. The recent emergence of the proteolysis targeting chimera (PROTAC) and molecular glue degrader (MGD) modalities has further increased interest in E3 ligases as drug targets. This perspective focuses on how fragment based lead discovery (FBLD) methods have been used to discover new ligands for this important target class. In some cases these efforts have led to clinical candidates; in others, they have provided tools for deepening our understanding of E3 ligase biology. Recently, FBLD-derived ligands have inspired the design of PROTACs that are able to artificially modulate protein levels in cells.

A novel pan-selective bromodomain inhibitor for epigenetic drug design

For a long time, the development of bromodomain (BD) inhibitors (BDi) was almost exclusively related to the BET family. More recently, BDi for BDs outside the BET family have also been developed. Here we present a novel pan-BDi with micromolar affinities to various BDs, and nanomolar affinities to representatives of BD families I, II (Bromodomain and Extra-Terminal Domain (BET) family), III, and IV. The inhibitor shows a broad activity profile with nanomolar growth inhibition (GI50) values on various cancer cell lines. Subsequently, we were able to control the selectivity of the inhibitor by simple modifications and turned it into a highly selective BRD9 inhibitor.

Identification of 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based small molecules as selective BRD9 binders

Targeting bromodomain-containing protein 9 (BRD9) represents a promising strategy for the development of new agents endowed with anticancer properties. With this aim, a set of 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based compounds was investigated following a combined approach that relied on in silico studies, chemical synthesis, biophysical and biological evaluation of the most promising items. The protocol was initially based on molecular docking experiments, accounting a library of 1896 potentially synthesizable items tested in silico against the bromodomain of BRD9. A first set of 21 compounds (1-21) was selected and the binding on BDR9 was assessed through AlphaScreen assays. The obtained results disclosed compounds 17 and 20 able to bind BRD9 in the submicromolar range (IC₅₀ = 0.35 ± 0.18 μM and IC₅₀ = 0.14 ± 0.03 μM, respectively) showing a promising selectivity profile when tested against further nine bromodomains. Taking advantage of 3D structure-based pharmacophore models, additional 10 derivatives were selected in silico for the synthetic step and binding assessment, highlighting seven compounds (22, 23, 25, 26, 28, 29, 31) able to selectively bind BRD9 among different bromodomains. The ability of the identified BRD9 binders to cross artificial membranes in vitro was also assessed, revealing a very good passive permeability profile. Preliminary studies were carried out on a panel of healthy and cancer human cell lines to explore the biological behavior of the selected compounds, disclosing a moderate activity and significant selectivity profile towards leukaemia cells. These results highlighted the applicability of the reported multidisciplinary approach for accelerating the selection of promising items and for driving the chemical synthesis of novel selective BRD9 binders. Moreover, the low molecular weight of the reported 2,4,5-trisubstituted-2,4-dihydro-3H-1,2,4-triazol-3-one-based BRD9 binders suggests the possibility for further exploring the chemical space in order to obtain new analogues with improved potency.

Absolute Binding Free Energy Calculations for Buried Water Molecules

Water often plays a key role in mediating protein-ligand interactions. Understanding contributions from active-site water molecules to binding thermodynamics of a ligand is important in predicting binding free energies for ligand optimization. In this work, we tested a non-equilibrium switching method for absolute binding free energy calculations on water molecules in binding sites of 13 systems. We discuss the lessons we learned about identified issues that affected our calculations and ways to address them. This work fits with our larger focus on how to do accurate ligand binding free energy calculations when water rearrangements are very slow, such as rearrangements due to ligand modification (as in relative free energy calculations) or ligand binding (as in absolute free energy calculations). The method studied in this work can potentially be used to account for limited water sampling via providing endpoint corrections to free energy calculations using our calculated binding free energy of water.

Identification of Histone Peptide Binding Specificity and Small-Molecule Ligands for the TRIM33α and TRIM33β Bromodomains

TRIM33 is a member of the tripartite motif (TRIM) family of proteins, some of which possess E3 ligase activity and are involved in the ubiquitin-dependent degradation of proteins. Four of the TRIM family proteins, TRIM24 (TIF1α), TRIM28 (TIF1β), TRIM33 (TIF1γ) and TRIM66, contain C-terminal plant homeodomain (PHD) and bromodomain (BRD) modules, which bind to methylated lysine (KMe_n) and acetylated lysine (KAc), respectively. Here we investigate the differences between the two isoforms of TRIM33, TRIM33α and TRIM33β, using structural and biophysical approaches. We show that the N1039 residue, which is equivalent to N140 in BRD4(1) and which is conserved in most BRDs, has a different orientation in each isoform. In TRIM33β, this residue coordinates KAc, but this is not the case in TRIM33α. Despite these differences, both isoforms show similar affinities for H3_1-27K18Ac, and bind preferentially to H3_1-27K9Me₃K18Ac. We used this information to develop an AlphaScreen assay, with which we have identified four new ligands for the TRIM33 PHD-BRD cassette. These findings provide fundamental new information regarding which histone marks are recognized by both isoforms of TRIM33 and suggest starting points for the development of chemical probes to investigate the cellular function of TRIM33.

Molecular Insights into the Recognition of Acetylated Histone Modifications by the BRPF2 Bromodomain

HBO1 [HAT bound to the origin recognition complex (ORC)], a member of the MYST family of histone acetyltransferases (HATs), was initially identified as a binding partner of ORC that acetylates free histone H3, H4, and nucleosomal H3. It functions as a quaternary complex with the BRPF (BRPF1/2/3) scaffolding protein and two accessory proteins, ING4/5 and Eaf6. Interaction of BRPF2 with HBO1 has been shown to be important for regulating H3K14 acetylation during embryonic development. However, how BRPF2 directs the HBO1 HAT complex to chromatin to regulate its HAT activity toward nucleosomal substrates remains unclear. Our findings reveal novel interacting partners of the BRPF2 bromodomain that recognizes different acetyllysine residues on the N-terminus of histone H4, H3, and H2A and preferentially binds to H4K5ac, H4K8ac, and H4K5acK12ac modifications. In addition, mutational analysis of the BRPF2 bromodomain coupled with isothermal titration calorimetry binding and pull-down assays on the histone substrates identified critical residues responsible for acetyllysine binding. Moreover, the BRPF2 bromodomain could enrich H4K5ac mark-bearing mononucleosomes compared to other acetylated H4 marks. Consistent with this, ChIP-seq analysis revealed that BRPF2 strongly co-localizes with HBO1 at histone H4K5ac and H4K8ac marks near the transcription start sites in the genome. Our study provides novel insights into how the histone binding function of the BRPF2 bromodomain directs the recruitment of the HBO1 HAT complex to chromatin to regulate gene expression.

GNE-064: A Potent, Selective, and Orally Bioavailable Chemical Probe for the Bromodomains of SMARCA2 and SMARCA4 and the Fifth Bromodomain of PBRM1

Bromodomains are acetyllysine recognition domains present in a variety of human proteins. Bromodomains also bind small molecules that compete with acetyllysine, and therefore bromodomains have been targets for drug discovery efforts. Highly potent and selective ligands with good cellular permeability have been proposed as chemical probes for use in exploring the functions of many of the bromodomain proteins. We report here the discovery of a class of such inhibitors targeting the family VIII bromodomains of SMARCA2 (BRM) and SMARCA4 (BRG1), and PBRM1 (polybromo-1) bromodomain 5. We propose one example from this series, GNE-064, as a chemical probe for the bromodomains SMARCA2, SMARCA4, and PBRM1(5) with the potential for in vivo use.

Alchemical Free Energy Methods Applied to Complexes of the First Bromodomain of BRD4

Accurate and rapid predictions of the binding affinity of a compound to a target are one of the ultimate goals of computer aided drug design. Alchemical approaches to free energy estimations follow the path from an initial state of the system to the final state through alchemical changes of the energy function during a molecular dynamics simulation. Herein, we explore the accuracy and efficiency of two such techniques: relative free energy perturbation (FEP) and multisite lambda dynamics (MSλD). These are applied to a series of inhibitors for the bromodomain-containing protein 4 (BRD4). We demonstrate a procedure for obtaining accurate relative binding free energies using MSλD when dealing with a change in the net charge of the ligand. This resulted in an impressive comparison with experiment, with an average difference of 0.4 ± 0.4 kcal mol-1. In a benchmarking study for the relative FEP calculations, we found that using 20 lambda windows with 0.5 ns of equilibration and 1 ns of data collection for each window gave the optimal compromise between accuracy and speed. Overall, relative FEP and MSλD predicted binding free energies with comparable accuracy, an average of 0.6 kcal mol-1 for each method. However, MSλD makes predictions for a larger molecular space over a much shorter time scale than relative FEP, with MSλD requiring a factor of 18 times less simulation time for the entire molecule space.

Enhancing Sampling of Water Rehydration on Ligand Binding: A Comparison of Techniques

Water often plays a key role in protein structure, molecular recognition, and mediating protein-ligand interactions. Thus, free energy calculations must adequately sample water motions, which often proves challenging in typical MD simulation time scales. Thus, the accuracy of methods relying on MD simulations ends up limited by slow water sampling. Particularly, as a ligand is removed or modified, bulk water may not have time to fill or rearrange in the binding site. In this work, we focus on several molecular dynamics (MD) simulation-based methods attempting to help rehydrate buried water sites: BLUES, using nonequilibrium candidate Monte Carlo (NCMC); grand, using grand canonical Monte Carlo (GCMC); and normal MD. We assess the accuracy and efficiency of these methods in rehydrating target water sites. We selected a range of systems with varying numbers of waters in the binding site, as well as those where water occupancy is coupled to the identity or binding mode of the ligand. We analyzed the rehydration of buried water sites in binding pockets using both clustering of trajectories and direct analysis of electron density maps. Our results suggest both BLUES and grand enhance water sampling relative to normal MD and grand is more robust than BLUES, but also that water sampling remains a major challenge for all of the methods tested. The lessons we learned for these methods and systems are discussed.

Discovery of Dual TAF1-ATR Inhibitors and Ligand-Induced Structural Changes of the TAF1 Tandem Bromodomain

Bromodomains regulate chromatin remodeling and gene transcription through recognition of acetylated lysines on histones and other proteins. Bromodomain-containing protein TAF1, a subunit of general transcription factor TFIID, initiates preinitiation complex formation and cellular transcription. TAF1 serves as a cofactor for certain oncogenic transcription factors and is implicated in regulating the p53 tumor suppressor. Therefore, TAF1 is a potential target to develop small molecule therapeutics for diseases arising from dysregulated transcription, such as cancer. Here, we report the ATR kinase inhibitor AZD6738 (Ceralasertib) and analogues thereof as bona fide inhibitors of TAF1. Crystallographic and small-angle X-ray scattering studies established that newly identified and previously reported inhibitors stabilize distinct structural states of the TAF1 tandem bromodomain through "open-closed" transitions and dimerization. Combined with functional studies on p53 signaling in cancer cell lines, the data provide new insights into the feasibility and challenges of TAF1 inhibitors as chemical probes and therapeutics.

A Structure-based Design Approach for Generating High Affinity BRD4 D1-Selective Chemical Probes

Chemical probes for epigenetic proteins are essential tools for dissecting the molecular mechanisms for gene regulation and therapeutic development. The bromodomain and extra-terminal (BET) proteins are master transcriptional regulators. Despite promising therapeutic targets, selective small molecule inhibitors for a single bromodomain remain an unmet goal due to their high sequence similarity. Here, we address this challenge via a structure-activity relationship study using 1,4,5-trisubstituted imidazoles against the BRD4 N-terminal bromodomain (D1). Leading compounds 26 and 30 have 15 and 18 nM affinity against BRD4 D1 and over 500-fold selectivity against BRD2 D1 and BRD4 D2 via ITC. Broader BET selectivity was confirmed by fluorescence anisotropy, thermal shift, and CETSA. Despite BRD4 engagement, BRD4 D1 inhibition was unable to reduce c-Myc expression at low concentration in multiple myeloma cells. Conversely, for inflammation, IL-8 and chemokine downregulation were observed. These results provide new design rules for selective inhibitors of an individual BET bromodomain.

Introducing structure-based three-dimensional pharmacophore models for accelerating the discovery of selective BRD9 binders

A well-structured in silico workflow is here reported for disclosing structure-based pharmacophore models against bromodomain-containing protein 9 (BRD9), accelerating virtual screening campaigns and facilitating the identification of novel binders. Specifically, starting from 23 known ligands co-crystallized with BRD9, three-dimensional pharmacophore models, namely placed in a reference protein structure, were developed. Specifically, we here introduce a fragment-related pharmacophore model, useful for the identification of new promising small chemical probes targeting the protein region responsible of the acetyllysine recognition, and two further pharmacophore models useful for the selection of compounds featuring drug-like properties. A pharmacophore-driven virtual screening campaign was then performed to facilitate the selection of new selective BRD9 ligands, starting from a large library of commercially available molecules. The identification of a promising BRD9 binder (7) prompted us to re-iterate this computational workflow on a second focused in-house built library of synthesizable compounds and, eventually, three further novel BRD9 binders were disclosed (8-10). Moreover, all these compounds were tested among a panel comprising other nine bromodomains, showing a high selectivity for BRD9. Preclinical bioscreens for potential anticancer activity highlighted compound 7 as that showing the most promising biological effects, proving the reliability of this in silico pipeline and confirming the applicability of the here introduced structure-based three-dimensional (3D) pharmacophore models as straightforward tools for the selection of new BRD9 ligands.

Insights into the Molecular Mechanisms of Histone Code Recognition by the BRPF3 Bromodomain

Bromodomains are evolutionarily conserved reader modules that recognize acetylated lysine residues on the histone tails to facilitate gene transcription. The bromodomain and PHD finger containing protein 3 (BRPF3) is a scaffolding protein that forms a tetrameric complex with HBO1 histone acetyltransferase (HAT) and two other subunits, which is known to regulate the HAT activity and substrate specificity. However, its molecular mechanism, histone ligands, and biological functions remain unknown. Herein, we identify mono- (H4K5ac) and di- (H4K5acK12ac) acetylated histone peptides as novel interacting partners of the BRPF3 bromodomain. Consistent with this, pull-down assays on purified histones from human cells confirm the interaction of BRPF3 bromodomain with acetylated histone H4. Further, MD simulation studies highlight the binding mode of acetyllysine (Kac) and the stability of bromodomain-histone peptide complexes. Collectively, our findings provide a key insight into how histone targets of the BRPF3 bromodomain direct the recruitment of HBO1 complex to chromatin for downstream transcriptional regulation.

From Hit Seeking to Magic Bullets: The Successful Union of Epigenetic and Fragment Based Drug Discovery (EPIDD + FBDD)

We review progress in the application of fragment-based drug discovery (FBDD) to epigenetic drug discovery (EPIDD) targeted at epigenetic writer and eraser enzymes as well as reader domains over the last 15 years. The greatest successes to date are in prospecting for bromodomain binding ligands. From a diverse array of fragment hits, multiple potent and selective compounds ensued, including the oncology clinical candidates mivebresib, ABBV-744, pelabresib, and PLX51107.

Ligand binding: evaluating the contribution of the water molecules network using the Fragment Molecular Orbital method

Water molecules play a crucial role in protein-ligand binding, and many tools exist that aim to predict the position and relative energies of these important, but challenging participants of biomolecular recognition. The available tools are, in general, capable of predicting the location of water molecules. However, predicting the effects of their displacement is still very challenging. In this work, a linear-scaling quantum mechanics-based approach was used to assess water network energetics and the changes in network stability upon ligand structural modifications. This approach offers a valuable way to improve understanding of SAR data and help guide compound design.

Optimization of Naphthyridones into Selective TATA-Binding Protein Associated Factor 1 (TAF1) Bromodomain Inhibitors

Bromodomain containing proteins and the acetyl-lysine binding bromodomains contained therein are increasingly attractive targets for the development of novel epigenetic therapeutics. To help validate this target class and unravel the complex associated biology, there has been a concerted effort to develop selective small molecule bromodomain inhibitors. Herein we describe the structure-based efforts and multiple challenges encountered in optimizing a naphthyridone template into selective TAF1(2) bromodomain inhibitors which, while unsuitable as chemical probes themselves, show promise for the future development of small molecules to interrogate TAF1(2) biology. Key to this work was the introduction and modulation of the basicity of a pendant amine which had a substantial impact on not only bromodomain selectivity but also cellular target engagement.

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 the first over the second bromodomain of BRD4: theoretical insights from free energy calculations and multiple short molecular dynamics simulations

Bromodomain-containing protein 4 (BRD4) plays an important role in mediating gene transcription involved in cancers and non-cancer diseases such as acute heart failure and inflammatory diseases. In this work, multiple short molecular dynamics (MSMD) simulations are integrated with a molecular mechanics generalized Born surface area (MM-GBSA) approach to decipher binding selectivity of three inhibitors 8NS, 82Y, and 837 toward two domains BD1 and BD2 of BRD4. The results demonstrate that the enthalpy effects play critical roles in selectivity identification of inhibitors toward BD1 and BD2, determining that 8NS has better selectivity toward BD2 than BD1, while 82Y and 837 more favorably bind to BD1 than BD2. A residue-based free-energy decomposition method was used to calculate an inhibitor-residue interaction spectrum and unveil contributions of separate residues to binding selectivity. The results identify six common residues, containing (P82, P375), (V87, V380), (L92, L385), (L94, L387), (N140, N433), and (I146, V439) individually belonging to (BD1, BD2) of BRD4, and yield a considerable binding difference of inhibitors to BD1 and BD2, suggesting that these residues play key roles in binding selectivity of inhibitors toward BD1 and BD2 of BRD4. Therefore, these results provide useful dynamics information and a structure affinity relationship for the development of highly selective inhibitors targeting BD1 and BD2 of BRD4.

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.

Enhancing Water Sampling in Free Energy Calculations with Grand Canonical Monte Carlo

The prediction of protein-ligand binding affinities using free energy perturbation (FEP) is becoming increasingly routine in structure-based drug discovery. Most FEP packages use molecular dynamics (MD) to sample the configurations of proteins and ligands, as MD is well-suited to capturing coupled motion. However, MD can be prohibitively inefficient at sampling water molecules that are buried within binding sites, which has severely limited the domain of applicability of FEP and its prospective usage in drug discovery. In this paper, we present an advancement of FEP that augments MD with grand canonical Monte Carlo (GCMC), an enhanced sampling method, to overcome the problem of sampling water. We accomplished this without degrading computational performance. On both old and newly assembled data sets of protein-ligand complexes, we show that the use of GCMC in FEP is essential for accurate and robust predictions for ligand perturbations that disrupt buried water.

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.

Application of Atypical Acetyl-lysine Methyl Mimetics in the Development of Selective Inhibitors of the Bromodomain-Containing Protein 7 (BRD7)/Bromodomain-Containing Protein 9 (BRD9) Bromodomains

Non-BET bromodomain-containing proteins have become attractive targets for the development of novel therapeutics targeting epigenetic pathways. To help facilitate the target validation of this class of proteins, structurally diverse small-molecule ligands and methodologies to produce selective inhibitors in a predictable fashion are in high demand. Herein, we report the development and application of atypical acetyl-lysine (KAc) methyl mimetics to take advantage of the differential stability of conserved water molecules in the bromodomain binding site. Discovery of the n-butyl group as an atypical KAc methyl mimetic allowed generation of 31 (GSK6776) as a soluble, permeable, and selective BRD7/9 inhibitor from a pyridazinone template. The n-butyl group was then used to enhance the bromodomain selectivity of an existing BRD9 inhibitor and to transform pan-bromodomain inhibitors into BRD7/9 selective compounds. Finally, a solvent-exposed vector was defined from the pyridazinone template to enable bifunctional molecule synthesis, and affinity enrichment chemoproteomic experiments were used to confirm several of the endogenous protein partners of BRD7 and BRD9, which form part of the chromatin remodeling PBAF and BAF complexes, respectively.

The BRPF1 bromodomain is a molecular reader of di-acetyllysine

Bromodomain-containing proteins are often part of chromatin-modifying complexes, and their activity can lead to altered expression of genes that drive cancer, inflammation and neurological disorders in humans. Bromodomain-PHD finger protein 1 (BRPF1) is part of the MOZ (monocytic leukemic zinc-finger protein) HAT (histone acetyltransferase) complex, which is associated with chromosomal translocations known to contribute to the development of acute myeloid leukemia (AML). BRPF1 contains a unique combination of chromatin reader domains including two plant homeodomain (PHD) fingers separated by a zinc knuckle (PZP domain), a bromodomain, and a proline-tryptophan-tryptophan-proline (PWWP) domain. BRPF1 is known to recruit the MOZ HAT complex to chromatin by recognizing acetylated lysine residues on the N-terminal histone tail region through its bromodomain. However, histone proteins can contain several acetylation modifications on their N-terminus, and it is unknown how additional marks influence bromodomain recruitment to chromatin. Here, we identify the BRPF1 bromodomain as a selective reader of di-acetyllysine modifications on histone H4. We used ITC assays to characterize the binding of di-acetylated histone ligands to the BRPF1 bromodomain and found that the domain binds preferentially to histone peptides H4K5acK8ac and H4K5acK12ac. Analytical ultracentrifugation (AUC) experiments revealed that the monomeric state of the BRPF1 bromodomain coordinates di-acetylated histone ligands. NMR chemical shift perturbation studies, along with binding and mutational analyses, revealed non-canonical regions of the bromodomain-binding pocket that are important for histone tail recognition. Together, our findings provide critical information on how the combinatorial action of post-translational modifications can modulate BRPF1 bromodomain binding and specificity.

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.

BET proteins: Investigating BRDT as a potential target for male contraception

While many contraception options are available for women, birth control methods for men are limited to condoms and vasectomy. Past research into male contraceptives has focused on hormonal options but the associated side effects have thus far precluded this method from reaching the market. Non-hormonal male contraceptives and vas occlusion have also been explored, but to date no method has progressed past clinical testing. Recent interest in epigenetic research has unveiled a new potential non-hormonal male contraceptive target: the testis-specific bromodomain BRDT. Potent inhibitors for bromodomain-containing proteins are described in the literature, but a BRDT-specific compound has yet to be designed, prepared and tested. The high similarity between bromodomain proteins of the BET family makes development of selective and specific inhibitors both difficult and necessary. Selective inhibition of BRDT by a small molecule is an exciting new target in the search for a new non-hormonal male contraceptive.

In silico design and molecular basis for the selectivity of Olinone toward the first over the second bromodomain of BRD4

Bromodomains (BrDs), a conserved structural module in chromatin-associated proteins, are well known for recognizing ε-N-acetyl lysine residues on histones. One of the most relevant BrDs is BRD4, a tandem BrD containing protein (BrD1 and BrD2) that plays a critical role in numerous diseases including cancer. Growing evidence shows that the two BrDs of BRD4 have different biological functions; hence selective ligands that can be used to study their functions are of great interest. Here, as a follow-up of our previous work, we first provide a detailed characterization study of the in silico rational design of Olinone as part of a series of five tetrahydropyrido indole-based compounds as BRD4 BrD1 inhibitors. Additionally, we investigated the molecular basis for Olinone's selective recognition by BrD1 over BrD2. Molecular dynamics simulations, free energy calculations, and conformational analyses of the apo-BRD4-BrD1|2 and BRD4-BrD1|2/Olinone complexes showed that Olinone's selectivity is facilitated by five key residues: Leu92 in BrD1|385 in BrD2 of ZA loop, Asn140|433, Asp144|His437 and Asp145|Glu438 of BC loop, and Ile146|Val49 of helix C. Furthermore, the difference in hydrogen bonds number and in mobility of the ZA and BC loops of the acetyl-lysine binding site between BRD4 BrD1/Olinone and BrD2/Olinone complexes also contribute to the difference in Olinone's binding affinity and selectivity toward BrD1 over BrD2. Altogether, our computer-aided molecular design techniques can effectively guide the development of small-molecule BRD4 BrD1 inhibitors, explain their selectivity origin, and further open doors to the design of new therapeutically improved derivatives.

Dual Inhibition of TAF1 and BET Bromodomains from the BI-2536 Kinase Inhibitor Scaffold

Recent reports have highlighted the dual bromodomains of TAF1 (TAF1(1,2)) as synergistic with BET inhibition in cellular cancer models, engendering interest in TAF/BET polypharmacology. Here, we examine structure activity relationships within the BI-2536 PLK1 kinase inhibitor scaffold, previously reported to bind BRD4. We examine binding by this ligand to TAF1(2) and apply structure guided design strategies to discriminate binding to both the PLK1 kinase and BRD4(1) bromodomain while retaining activity on TAF1(2). Through this effort we discover potent dual inhibitors of TAF1(2)/BRD4(1), as well as biased derivatives showing marked TAF1 selectivity. We resolve X-ray crystallographic data sets to examine the mechanisms of the observed TAF1 selectivity and to provide a resource for further development of this scaffold.

A Qualified Success: Discovery of a New Series of ATAD2 Bromodomain Inhibitors with a Novel Binding Mode Using High-Throughput Screening and Hit Qualification

The bromodomain of ATAD2 has proved to be one of the least-tractable proteins within this target class. Here, we describe the discovery of a new class of inhibitors by high-throughput screening and show how the difficulties encountered in establishing a screening triage capable of finding progressible hits were overcome by data-driven optimization. Despite the prevalence of nonspecific hits and an exceptionally low progressible hit rate (0.001%), our optimized hit qualification strategy employing orthogonal biophysical methods enabled us to identify a single active series. The compounds have a novel ATAD2 binding mode with noncanonical features including the displacement of all conserved water molecules within the active site and a halogen-bonding interaction. In addition to reporting this new series and preliminary structure-activity relationship, we demonstrate the value of diversity screening to complement the knowledge-based approach used in our previous ATAD2 work. We also exemplify tactics that can increase the chance of success when seeking new chemical starting points for novel and less-tractable targets.

Bromodomains: a new target class for drug development

Less than a decade ago, it was shown that bromodomains, acetyl lysine 'reader' modules found in proteins with varied functions, were highly tractable small-molecule targets. This is an unusual property for protein-protein or protein-peptide interaction domains, and it prompted a wave of chemical probe discovery to understand the biological potential of new agents that targeted bromodomains. The original examples, inhibitors of the bromodomain and extra-terminal (BET) class of bromodomains, showed enticing anti-inflammatory and anticancer activities, and several compounds have since advanced to human clinical trials. Here, we review the current state of BET inhibitor biology in relation to clinical development, and we discuss the next wave of bromodomain inhibitors with clinical potential in oncology and non-oncology indications. The lessons learned from BET inhibitor programmes should affect efforts to develop drugs that target non-BET bromodomains and other epigenetic readers.

Advancements in the Development of non-BET Bromodomain Chemical Probes

The bromodomain and extra terminal (BET) family of bromodomain-containing proteins (BCPs) have been the subject of extensive research over the past decade, resulting in a plethora of high-quality chemical probes for their tandem bromodomains. In turn, these chemical probes have helped reveal the profound biological role of the BET bromodomains and their role in disease, ultimately leading to a number of molecules in active clinical development. However, the BET subfamily represents just 8/61 of the known human bromodomains, and attention has now expanded to the biological role of the remaining 53 non-BET bromodomains. Rapid growth of this research area has been accompanied by a greater understanding of the requirements for an effective bromodomain chemical probe and has led to a number of new non-BET bromodomain chemical probes being developed. Advances since December 2015 are discussed, highlighting the strengths/caveats of each molecule, and the value they add toward validating the non-BET bromodomains as tractable therapeutic targets.

Importance of a crystalline water network in docking-based virtual screening: a case study of BRD4

As a member of the bromodomain and extra terminal domain (BET) protein family, bromodomain-containing protein 4 (BRD4) is an epigenetic reader and can recognize acetylated lysine residues in histones.

Solvents to Fragments to Drugs: MD Applications in Drug Design

Simulations of molecular dynamics (MD) are playing an increasingly important role in structure-based drug discovery (SBDD). Here we review the use of MD for proteins in aqueous solvation, organic/aqueous mixed solvents (MDmix) and with small ligands, to the classic SBDD problems: Binding mode and binding free energy predictions. The simulation of proteins in their condensed state reveals solvent structures and preferential interaction sites (hot spots) on the protein surface. The information provided by water and its cosolvents can be used very effectively to understand protein ligand recognition and to improve the predictive capability of well-established methods such as molecular docking. The application of MD simulations to the study of the association of proteins with drug-like compounds is currently only possible for specific cases, as it remains computationally very expensive and labor intensive. MDmix simulations on the other hand, can be used systematically to address some of the common tasks in SBDD. With the advent of new tools and faster computers we expect to see an increase in the application of mixed solvent MD simulations to a plethora of protein targets to identify new drug candidates.

Structure investigation, enrichment analysis and structure-based repurposing of FDA-approved drugs as inhibitors of BET-BRD4

We report herein detailed structural insights into the ligand recognition modes guiding bromodomain selectivity, enrichment analysis and docking-based database screening for the identification of the FDA-approved drugs that have potential to be the human BRD4 inhibitors. Analysis of multiple X-ray structures prevailed that the lysine-recognition sites are highly conserved, and apparently, the dynamic ZA loop guides the specific ligand-recognition. The protein-ligand interaction profiling revealed that both BRD2 and BRD4 shared hydrophobic interaction of bound ligands with PRO-98/PRO-82, PHE-99/PHE-83, LEU-108/LEU-92 and direct H-bonding with ASN-156/ASN-140 (BRD2/BRD4), while on the other hand the water-mediated H-bonding of bound ligands with PRO-82, GLN-85, PRO-86, VAL-87, ASP-88, LEU-92, TYR-97 and MET-132, and aromatic π-π stacking with TRP-81 prevailed as unique interaction in BRD4, and were not observed in BRD2. Subsequently, through ROC curve analysis, the best enrichment was found with PDB-ID 4QZS of BRD4 structures. Finally, through docking-based database screening study, we found that several drugs have better binding affinity than the control candidate lead (+)-JQ1 (Binding affinity = -7.9 kcal/mol), a well-known BRD4 inhibitor. Among the top-ranked drugs, azelastine, a selective histamine H1 receptor antagonist, showed the best binding affinity of -9.3 kcal/mol and showed interactions with several key residues of the acetyl lysine binding pocket. Azelastine may serve as a promising template for further medicinal chemistry. These insights may serve as basis for structure-based drug design, drug repurposing and the discovery of novel BRD4 inhibitors. Communicated by Ramaswamy H. Sarma.