Structure-based discovery of selective BRPF1 bromodomain inhibitors

Inhibitors targeting the PWWP domain-containing proteins

PWWP domain-containing proteins play a pivotal role in chromatin-mediated biological processes, and their aberrant regulation is linked to various human diseases. Recent years have witnessed remarkable strides in unraveling the structural and functional features of PWWP domain-containing proteins, propelling significant advances in targeting the PWWP domain-containing proteins for drug discovery purposes. Several drugs have already been approved, while others are currently in clinical trials. This review offers a comprehensive overview of the latest developments on PWWP domain-containing proteins, including their structural characteristics and biological significance. It also provides detailed insights into the drug discovery process targeting these proteins, including screening, design, and structural optimization.

Accelerating BRPF1b hit identification with BioPhysical and Active Learning Screening (BioPALS)

We report the development of BioPhysical and Active Learning Screening (BioPALS); a rapid and versatile hit identification protocol combining AI-powered virtual screening with a GCI-driven biophysical confirmation workflow. Its application to the BRPF1b bromodomain afforded a range of novel micromolar binders with favorable ADMET properties. In addition to the excellent in silico/in vitro confirmation rate demonstrated with BRPF1b, binding kinetics were determined, and binding topologies predicted for all hits. BioPALS is a lean, data-rich, and standardized approach to hit identification applicable to a wide range of biological targets.

Integrated virtual screening and MD simulation approaches toward discovering potential inhibitors for targeting BRPF1 bromodomain in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most aggressive and life-threatening cancers. Although multiple treatment options are available, the prognosis of HCC patients is poor due to metastasis and drug resistance. Hence, discovering novel targets is essential for better therapeutic development for HCC. In this study, we used the cancer genome atlas (TCGA) dataset to analyze the expression of bromodomain-containing proteins in HCC, as bromodomains are emerging attractive therapeutic targets. Our analysis identified BRPF1 as the most highly upregulated gene in HCC among the 43 bromodomain-containing genes. Upregulation of BRPF1 was significantly associated with poorer patient survival. Therefore, targeting BRPF1 may be an approach for HCC treatment. Previously, several potential inhibitors of BRPF1 bromodomain have been discovered. However, due to the limited clinical success of the current inhibitors, we aim to search for new inhibitors with high affinity and specificity for the BRPF1 bromodomain. In this study, we utilized high-throughput virtual screening methods to screen synthetic and natural compound databases against the BRPF1 bromodomain. In addition, we used machine learning-based QSAR modeling to predict the IC50 values of the selected BRPF1 bromodomain inhibitors. Extensive MD simulations were used to calculate the binding free energies of BRPF1 bromodomain and inhibitor complexes. Using this approach, we identified four lead scaffolds with a similar or better binding affinity towards the BRPF1 bromodomain than the previously reported inhibitors. Overall, this study discovered some promising compounds that have the potential to act as potent BRPF1 bromodomain inhibitors.

Mosaicism in BRPF1-Related Neurodevelopmental Disorder: Report of Two Sisters and Literature Review

Bromodomain and PHD finger containing 1 (BRPF1)-related neurodevelopmental disorder is characterized by intellectual disability, developmental delay, hypotonia, dysmorphic facial features, ptosis, and blepharophimosis. Both de novo and inherited pathogenic variants have been previously reported in association with this disorder. We report two affected female siblings with a novel variant in BRPF1 c.2420_2433del (p.Q807Lfs∗27) identified through whole-exome sequencing. Their history of mild intellectual disability, speech delay, attention deficient hyperactivity disorder (ADHD), and ptosis align with the features previously reported in the literature. The absence of the BRPF1 variant in parental buccal samples provides evidence of a de novo frameshift pathogenic variant, most likely as a result of parental gonadal mosaicism, which has not been previously reported. The frameshift pathogenic variant reported here lends further support to haploinsufficiency as the underlying mechanism of disease. We review the literature, compare the clinical features seen in our patients with others reported, and explore the possibility of genotype-phenotype correlation based on the location of pathogenic variants in BRPF1. Our study helps to summarize available knowledge and report the first case of a de novo frameshift pathogenic variant in BRPF1 in two siblings with this neurodevelopmental disorder.

The role of loop dynamics in the prediction of ligand-protein binding enthalpy

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R² = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol^-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R² = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol^-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

Targeting bromodomain-containing proteins: research advances of drug discovery

Bromodomain (BD) is an evolutionarily conserved protein module found in 46 different BD-containing proteins (BCPs). BD acts as a specific reader for acetylated lysine residues (KAc) and serves an essential role in transcriptional regulation, chromatin remodeling, DNA damage repair, and cell proliferation. On the other hand, BCPs have been shown to be involved in the pathogenesis of a variety of diseases, including cancers, inflammation, cardiovascular diseases, and viral infections. Over the past decade, researchers have brought new therapeutic strategies to relevant diseases by inhibiting the activity or downregulating the expression of BCPs to interfere with the transcription of pathogenic genes. An increasing number of potent inhibitors and degraders of BCPs have been developed, some of which are already in clinical trials. In this paper, we provide a comprehensive review of recent advances in the study of drugs that inhibit or down-regulate BCPs, focusing on the development history, molecular structure, biological activity, interaction with BCPs and therapeutic potentials of these drugs. In addition, we discuss current challenges, issues to be addressed and future research directions for the development of BCPs inhibitors. Lessons learned from the successful or unsuccessful development experiences of these inhibitors or degraders will facilitate the further development of efficient, selective and less toxic inhibitors of BCPs and eventually achieve drug application in the clinic.

BRPF1-KAT6A/KAT6B Complex: Molecular Structure, Biological Function and Human Disease

The bromodomain and PHD finger-containing protein1 (BRPF1) is a member of family IV of the bromodomain-containing proteins that participate in the post-translational modification of histones. It functions in the form of a tetrameric complex with a monocytic leukemia zinc finger protein (MOZ or KAT6A), MOZ-related factor (MORF or KAT6B) or HAT bound to ORC1 (HBO1 or KAT7) and two small non-catalytic proteins, the inhibitor of growth 5 (ING5) or the paralog ING4 and MYST/Esa1-associated factor 6 (MEAF6). Mounting studies have demonstrated that all the four core subunits play crucial roles in different biological processes across diverse species, such as embryonic development, forebrain development, skeletal patterning and hematopoiesis. BRPF1, KAT6A and KAT6B mutations were identified as the cause of neurodevelopmental disorders, leukemia, medulloblastoma and other types of cancer, with germline mutations associated with neurodevelopmental disorders displaying intellectual disability, and somatic variants associated with leukemia, medulloblastoma and other cancers. In this paper, we depict the molecular structures and biological functions of the BRPF1-KAT6A/KAT6B complex, summarize the variants of the complex related to neurodevelopmental disorders and cancers and discuss future research directions and therapeutic potentials.

Computational studies with flavonoids and terpenoids as BRPF1 inhibitors: in silico biological activity prediction, molecular docking, molecular dynamics simulations, MM/PBSA calculations

The BRPF1 protein is encoded by the BRPF1 gene. In addition, the BRPF1 gene is known to be upregulated in leukaemia. Recent studies have shown that it is also overexpressed in hepatocellular carcinoma (HCC) as well. Therefore, BRPF1 is a significant target for anti-cancer drug development studies, especially on HCC. 40 terpenoids and flavonoids were chosen because of their anticancer properties given in the literature. In this study, the biological activity of molecules was also investigated with in silico structure-activity relationship analysis. In addition, interactions between a series of terpenoids and flavonoids and the BRPF1 protein were investigated by molecular docking and molecular dynamics simulations. The energy change caused by the interactions of BRPF1 with different compounds was also evaluated by MM/PBSA calculations. It has been revealed that compound 5 (-9.2 kcal/mol), a kind of secoclerodane type diterpenoid, has a higher affinity both compared to other flavonoids and terpenoids, and 9F9 (-7.9 kcal/mol), a selective BRPF1 inhibitor. The study presented in this article demonstrates that compound 5, as a natural product, could form a chemical scaffold for the development of selective BRPF1 bromodomain inhibitors.

Discovery, optimization and evaluation of 1-(indolin-1-yl)ethan-1-ones as novel selective TRIM24/BRPF1 bromodomain inhibitors

TRIM24 (tripartite motif-containing protein 24) and BRPF1 (bromodomain and PHD finger containing protein 1) are epigenetics "readers" and potential therapeutic targets for cancer and other diseases. Here we describe the structure-guided design of 1-(indolin-1-yl)ethan-1-ones as novel TRIM24/BRPF1 bromodomain inhibitors. The representative compound 20l (Y08624) is a new TRIM24/BRPF1 dual inhibitor, with IC₅₀ values of 0.98 and 1.16 μM, respectively. Cellular activity of 20l was validated by viability assay in prostate cancer (PC) cell lines. In PC xenograft models, 20l suppressed tumor growth (50 mg/kg/day, TGI = 53%) without exhibiting noticeable toxicity. Compound 20l represents a versatile starting point for the development of more potent TRIM24/BRPF1 inhibitors.

Insights into the post-translational modification and its emerging role in shaping the tumor microenvironment

More and more in-depth studies have revealed that the occurrence and development of tumors depend on gene mutation and tumor heterogeneity. The most important manifestation of tumor heterogeneity is the dynamic change of tumor microenvironment (TME) heterogeneity. This depends not only on the tumor cells themselves in the microenvironment where the infiltrating immune cells and matrix together forming an antitumor and/or pro-tumor network. TME has resulted in novel therapeutic interventions as a place beyond tumor beds. The malignant cancer cells, tumor infiltrate immune cells, angiogenic vascular cells, lymphatic endothelial cells, cancer-associated fibroblastic cells, and the released factors including intracellular metabolites, hormonal signals and inflammatory mediators all contribute actively to cancer progression. Protein post-translational modification (PTM) is often regarded as a degradative mechanism in protein destruction or turnover to maintain physiological homeostasis. Advances in quantitative transcriptomics, proteomics, and nuclease-based gene editing are now paving the global ways for exploring PTMs. In this review, we focus on recent developments in the PTM area and speculate on their importance as a critical functional readout for the regulation of TME. A wealth of information has been emerging to prove useful in the search for conventional therapies and the development of global therapeutic strategies.

The Discovery of Potential MDM2 Inhibitors: A Combination of Pharmacophore Modeling, Virtual Screening, Molecular Docking Studies, and in vitro/in vivo Biological Evaluation

Small-molecule inhibitors of MDM2 that block the MDM2-p53 protein-protein interaction have been considered as potential therapeutic agents for the treatment of cancer. Here, we identify five highly potent inhibitors of MDM2 (termed as WY 1-5) that display significant inhibitory effects on MDM2-p53 interaction by using a combined strategy of pharmacophore modeling, virtual screening, and molecular docking studies. Among them, WY-5 is the most active MDM2 inhibitor with an IC₅₀ value of 14.1±2.8 nM. Moreover, WY-5 significantly up-regulate the protein level of p53 in SK-Hep-1 cells harboring wild-type p53. In vitro anticancer study reveals that WY-5 markedly inhibits the survival of SK-Hep-1 cells. In vivo anticancer study suggests that WY-5 significantly inhibits the growth of SK-Hep-1 cells-derived xenograft in nude mice, with no observable toxicity. Our results demonstrate that WY-5 may be a promising candidate for the treatment of cancer harboring wild-type p53.

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.

Functional Roles of Bromodomain Proteins in Cancer

Histone acetylation is generally associated with an open chromatin configuration that facilitates many cellular processes including gene transcription, DNA repair, and DNA replication. Aberrant levels of histone lysine acetylation are associated with the development of cancer. Bromodomains represent a family of structurally well-characterized effector domains that recognize acetylated lysines in chromatin. As part of their fundamental reader activity, bromodomain-containing proteins play versatile roles in epigenetic regulation, and additional functional modules are often present in the same protein, or through the assembly of larger enzymatic complexes. Dysregulated gene expression, chromosomal translocations, and/or mutations in bromodomain-containing proteins have been correlated with poor patient outcomes in cancer. Thus, bromodomains have emerged as a highly tractable class of epigenetic targets due to their well-defined structural domains, and the increasing ease of designing or screening for molecules that modulate the reading process. Recent developments in pharmacological agents that target specific bromodomains has helped to understand the diverse mechanisms that bromodomains play with their interaction partners in a variety of chromatin processes, and provide the promise of applying bromodomain inhibitors into the clinical field of cancer treatment. In this review, we explore the expression and protein interactome profiles of bromodomain-containing proteins and discuss them in terms of functional groups. Furthermore, we highlight our current understanding of the roles of bromodomain-containing proteins in cancer, as well as emerging strategies to specifically target bromodomains, including combination therapies using bromodomain inhibitors alongside traditional therapeutic approaches designed to re-program tumorigenesis and metastasis.

Research progress of dual inhibitors targeting crosstalk between histone epigenetic modulators for cancer therapy

Abnormal epigenetics is a critical hallmark of human cancers. Anticancer drug discovery directed at histone epigenetic modulators has gained impressive advances with six drugs available for cancer therapy and numerous other candidates undergoing clinical trials. However, limited therapeutic profile, drug resistance, narrow safety margin, and dose-limiting toxicities pose intractable challenges for their clinical utility. Because histone epigenetic modulators undergo intricate crosstalk and act cooperatively to shape an aberrant epigenetic profile, co-targeting histone epigenetic modulators with a different mechanism of action has rapidly emerged as an attractive strategy to overcome the limitations faced by the single-target epigenetic inhibitors. In this review, we summarize in detail the crosstalk of histone epigenetic modulators in regulating gene transcription and the progress of dual epigenetic inhibitors targeting this crosstalk.

Exploring the Role of Asp1116 in Selective Drug Targeting of CREBcAMP- Responsive Element-binding Protein Implicated in Prostate Cancer

BACKGROUND

The selective targeting of CREB-cAMP-responsive element-binding protein (CBP) has recently evolved as a vital therapeutic approach for curtailing its aberrant upregulation associated with the development of prostate cancer. Inhibition of CBP has been discovered to be an important therapeutic option in androgen receptor signalling pathway mediated prostate cancer. Y08197, a novel selective inhibitor of CBP, has shown promising therapeutic outcome in prostate carcinogenesis over non-selective analogues such as CPI-637.

METHODS/RESULTS

Herein, we used molecular dynamics simulation to gain insights into the mechanistic and selective targeting of Y08197 at the bromodomain active site. Molecular Mechanics/ Poisson-Boltzmann Surface Area (MM/PBSA) analysis revealed a similar inhibitory effect between Y08197 and CPI-637. Furthermore, in exploring the selective affinity of Y08197 towards CBP in combination with Bromodomain and PHD finger-containing protein 1(BRPF1), our findings highlighted Asp1116 as the 'culprit' residue responsible for this selective targeting. Upon binding, Asp1116 assumed a conformation that altered the architecture of the bromodomain active site, thereby orienting the helices around the active site in a more compacted position. In addition to some specific structural perturbations mediated by Asp1116 on the dynamics of CBP, our study revealed that the strong hydrogen bond interaction (N-H...O) elicited in CBP-Y08197 sequestered Y08197 tightly into the CBP bromodomain active site.

CONCLUSION

Conclusively, the inhibition and selective pattern of Y08197 can be replicated in future structure-based CBP inhibitors and other bromodomain implicated in carcinogenesis.

Design, synthesis, and biological evaluation of dual targeting inhibitors of histone deacetylase 6/8 and bromodomain BRPF1

Histone modifying proteins, specifically histone deacetylases (HDACs) and bromodomains, have emerged as novel promising targets for anticancer therapy. In the current work, based on available crystal structures and docking studies, we designed dual inhibitors of both HDAC6/8 and the bromodomain and PHD finger containing protein 1 (BRPF1). Biochemical and biophysical tests showed that compounds 23a,b and 37 are nanomolar inhibitors of both target proteins. Detailed structure-activity relationships were deduced for the synthesized inhibitors which were supported by extensive docking and molecular dynamics studies. Cellular testing in acute myeloid leukemia (AML) cells showed only a weak effect, most probably because of the poor permeability of the inhibitors. We also aimed to analyse the target engagement and the cellular activity of the novel inhibitors by determining the protein acetylation levels in cells by western blotting (tubulin vs histone acetylation), and by assessing their effects on various cancer cell lines.

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.

Identification and characterization of fragment binding sites for allosteric ligand design using the site identification by ligand competitive saturation hotspots approach (SILCS-Hotspots)

BACKGROUND

Fragment-based ligand design is used for the development of novel ligands that target macromolecules, most notably proteins. Central to its success is the identification of fragment binding sites that are spatially adjacent such that fragments occupying those sites may be linked to create drug-like ligands. Current experimental and computational approaches that address this problem typically identify only a limited number of sites as well as use a limited number of fragment types.

METHODS

The site-identification by ligand competitive saturation (SILCS) approach is extended to the identification of fragment bindings sites, with the method termed SILCS-Hotspots. The approach involves precomputation of the SILCS FragMaps following which the identification of Hotspots, performed by identifying of all possible fragment binding sites on the full 3D structure of the protein followed by spatial clustering.

RESULTS

The SILCS-Hotspots approach identifies a large number of sites on the target protein, including many sites not accessible in experimental structures due to low binding affinities and binding sites on the protein interior. The identified sites are shown to recapitulate the location of known drug-like molecules in both allosteric and orthosteric binding sites on seven proteins including the androgen receptor, the CDK2 and Erk5 kinases, PTP1B phosphatase and three GPCRs; the β2-adrenergic, GPR40 fatty-acid binding and M2-muscarinic receptors. Analysis indicates the importance of considering all possible fragment binding sites, and not just those accessible to experimental methods, when identifying novel binding sites and performing ligand design versus just considering the most favorable sites. The approach is shown to identify a larger number of known binding sites of drug-like molecules versus the commonly used FTMap and Fpocket methods.

GENERAL SIGNIFICANCE

The present results indicate the potential utility of the SILCS-Hotspots approach for fragment-based rational design of ligands, including allosteric modulators.

Comprehensive evaluation of the MM-GBSA method on bromodomain-inhibitor sets

MM-PB/GBSA methods represent a higher-level scoring theory than docking. This study reports an extensive testing of different MM-GBSA scoring schemes on two bromodomain (BRD) datasets. The first set is composed of 24 BRPF1 complexes, and the second one is a nonredundant set constructed from the PDBbind and composed of 28 diverse BRD complexes. A variety of MM-GBSA schemes were analyzed to evaluate the performance of four protocols with different numbers of minimization and MD steps, 10 different force fields and three different water models. Results showed that neither additional MD steps nor unfixing the receptor atoms improved scoring or ranking power. On the contrary, our results underscore the advantage of fixing receptor atoms or limiting the number of MD steps not only for a reduction in the computational costs but also for boosting the prediction accuracy. Among Amber force fields tested, ff14SB and its derivatives rather than ff94 or polarized force fields provided the most accurate scoring and ranking results. The TIP3P water model yielded the highest scoring and ranking power compared to the others. Posing power was further evaluated for the BRPF1 set. A slightly better posing power for the protocol which uses both minimization and MD steps with a fixed receptor than the one which uses only minimization with a fully flexible receptor-ligand system was observed. Overall, this study provides insights into the usage of the MM-GBSA methods for screening of BRD inhibitors, substantiating the benefits of shorter protocols and latest force fields and maintaining the crystal waters for accuracy.

Molecular Basis for the PZP Domain of BRPF1 Association with Chromatin

The assembly of human histone acetyltransferase MOZ/MORF complexes relies on the scaffolding bromodomain plant homeodomain (PHD) finger 1 (BRPF1) subunit. The PHD-zinc-knuckle-PHD module of BRPF1 (BRPF1_PZP) has been shown to associate with the histone H3 tail and DNA; however, the molecular mechanism underlying recognition of H3 and the relationship between the histone and DNA-binding activities remain unclear. In this study, we report the crystal structure of BRPF1_PZP bound to the H3 tail and characterize the role of the bipartite interaction in the engagement of BRPF1_PZP with the nucleosome core particle (NCP). We find that although both interactions of BRPF1_PZP with the H3 tail and DNA are required for tight binding to NCP and for acetyltransferase function of the BRPF1-MORF-ING5-MEAF6 complex, binding to extranucleosomal DNA dominates. Our findings suggest that functionally active BRPF1_PZP might be important in stabilization of the MOZ/MORF complexes at chromatin with accessible DNA.

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.

Disulfide bridge formation influences ligand recognition by the ATAD2 bromodomain

The ATPase family, AAA domain-containing protein 2 (ATAD2) has a C-terminal bromodomain, which functions as a chromatin reader domain recognizing acetylated lysine on the histone tails within the nucleosome. ATAD2 is overexpressed in many cancers and its expression is correlated with poor patient outcomes, making it an attractive therapeutic target and potential biomarker. We solved the crystal structure of the ATAD2 bromodomain and found that it contains a disulfide bridge near the base of the acetyllysine binding pocket (Cys1057-Cys1079). Site-directed mutagenesis revealed that removal of a free C-terminal cysteine (C1101) residue greatly improved the solubility of the ATAD2 bromodomain in vitro. Isothermal titration calorimetry experiments in combination with the Ellman's assay demonstrated that formation of an intramolecular disulfide bridge negatively impacts the ligand binding affinities and alters the thermodynamic parameters of the ATAD2 bromodomain interaction with a histone H4K5ac peptide as well as a small molecule bromodomain ligand. Molecular dynamics simulations indicate that the formation of the disulfide bridge in the ATAD2 bromodomain does not alter the structure of the folded state or flexibility of the acetyllysine binding pocket. However, consideration of this unique structural feature should be taken into account when examining ligand-binding affinity, or in the design of new bromodomain inhibitor compounds that interact with this acetyllysine reader module.

From bench to bedside, via desktop. Recent advances in the application of cutting-edge in silico tools in the research of drugs targeting bromodomain modules

The discipline of drug discovery has greatly benefited by computational tools and in silico algorithms aiming at rationalization of many related processes, from the stage of early hit identification to the preclinical phases of drug candidate validation. The various methodologies referred to as molecular modeling tools span a broad spectrum of applications, from straightforward approaches such as virtual screening of compound libraries to more advanced techniques involving the precise estimation of free energy upon binding of the candidate drug to its macromolecular target. To this end, we report an overview of specific studies where implementation of such sophisticated modeling algorithms has shown to be indispensable for addressing challenging systems and biological questions otherwise difficult to answer. We focus our attention on the emerging field of bromodomain inhibitors. Bromodomains are small modules involved in epigenetic signaling and currently comprise high-priority targets for developing both drug candidates and chemical probes for basic biomedical research. We attempt a critical presentation of selected cases utilizing cutting-edge in silico methodologies, with our main emphasis being on absolute or relative free energy simulations, on implementation of quantum-mechanics level calculations and on characterization of solvent thermodynamics. We discuss the advantages and strengths as well as the drawbacks and weaknesses of computational tools utilized in those works and we attempt to comment on specific issues related to their integration into the regular medicinal chemistry practice. Our conclusion is that while such methods indeed represent highly promising resources for further advancing the discipline, their application is not always trivial.

Binding Motifs in the CBP Bromodomain: An Analysis of 20 Crystal Structures of Complexes with Small Molecules

We analyze 20 crystal structures of complexes between the CBP bromodomain and small-molecule ligands that belong to eight different chemotypes identified by docking. The binding motif of the moiety that mimics the natural ligand (acetylated side chain of lysine) at the bottom of the binding pocket is conserved. In stark contrast, the rest of the ligands form different interactions with different side chains and backbone polar groups on the outer rim of the binding pocket. Hydrogen bonds are direct or water-bridged. van der Waals contacts are optimized by rotations of hydrophobic side chains and a slight inward displacement of the ZA loop. Rare types of interactions are observed for some of the ligands.

In silico fragment-based drug design with SEED

We report on two fragment-based drug design protocols, SEED2XR and ALTA, which start by high-throughput docking. SEED2XR is a two-stage protocol for fragment-based drug design. The first stage is in silico and consists of the automatic docking of 10³-10⁴ fragments using SEED, which requires about 1 s per fragment. SEED is a docking software developed specifically for fragment docking and binding energy evaluation by a force field with implicit solvent. In the second stage of SEED2XR, the 10-10² fragments with the most favorable predicted binding energies are validated by protein X-ray crystallography. The recent applications of SEED2XR to bromodomains demonstrate that the whole SEED2XR protocol can be carried out in about a week of working time, with hit rates ranging from 10% to 40%. Information on fragment-target interactions generated by the SEED2XR protocol or directly from SEED docking has been used for the discovery of hundreds of hits. ALTA is a computational protocol for screening which identifies candidate ligands that preserve the interactions between the optimal SEED fragments and the protein target. Medicinal chemistry optimization of ligands predicted by ALTA has resulted in pre-clinical candidates for protein kinases and bromodomains. The high-throughput, very low cost, sustainability, and high hit rate of the SEED-based protocols, unreachable by purely experimental techniques, make them perfectly suitable for both academic and industrial drug discovery research.