Structural Insights into the Recognition of Mono- and Diacetylated Histones by the ATAD2B Bromodomain

Impact of Combinatorial Histone Modifications on Acetyllysine Recognition by the ATAD2 and ATAD2B Bromodomains

The ATPase family AAA+ domain containing 2 (ATAD2) protein and its paralog ATAD2B have a C-terminal bromodomain (BRD) that functions as a reader of acetylated lysine residues on histone proteins. Using a structure-function approach, we investigated the ability of the ATAD2/B BRDs to select acetylated lysine among multiple histone post-translational modifications. The ATAD2B BRD can bind acetylated histone ligands that also contain adjacent methylation or phosphorylation marks, while the presence of these modifications significantly weakened the acetyllysine binding activity of the ATAD2 BRD. Our structural studies provide mechanistic insights into how ATAD2/B BRD-binding pocket residues coordinate the acetyllysine group in the context of adjacent post-translational modifications. Furthermore, we investigated how sequence changes in amino acids of the histone ligands impact the recognition of an adjacent acetyllysine residue. Our study highlights how the interplay between multiple combinations of histone modifications influences the reader activity of the ATAD2/B BRDs, resulting in distinct binding modes.

A high-throughput molecular dynamics screening (HTMDS) approach to the design of novel cyclopeptide inhibitors of ATAD2B based on the non-canonical combinatorial library

Cyclic peptides (CPs) are a promising class of drugs because of their high biological activity and specificity. However, the design of CP remains challenging due to their conformational flexibility and difficulties in designing stable binding conformation. Herein, we present a high-throughput MD screening (HTMDS) process for the iterative design of stable CP binders with a combinatorial CP library composed of canonical and non-canonical amino acids. As a proof of concept, we apply our methods to design CP inhibitors for the bromodomain (BrD) of ATAD2B. 698,800 CP candidates with a total of 25,570 ns MD simulations were performed to study the protein-ligand binding interactions. The binding free energies (ΔGbind) estimated by MM/PBSA approach for eight lead CP designs were found to be low. CP-1st.43 was the best CP candidate with an estimated ΔGbind of -28.48 kcal/mol when compared to the standard inhibitor C-38 which has been experimentally validated and shown to exhibit ΔGbind of -17.11 kcal/mol. The major contribution of binding sites for BrD of ATAD2B involved the hydrogen-bonding anchor within the Aly-binding pocket, salt bridging, and hydrogen-bonding mediated stabilization of the ZA loop and BC loop, and the complementary Van der Waals attraction. Our methods demonstrate encouraging results by yielding conformationally stable and high-potential CP binders that should have potential applicability in future CP drug development.Communicated by Ramaswamy H. Sarma.

The Impact of Mutational Hotspots on Cancer Survival

BACKGROUND

Cofactors, biomarkers, and the mutational status of genes such as TP53, EGFR, IDH1/2, or PIK3CA have been used for patient stratification. However, many genes exhibit recurrent mutational positions known as hotspots, specifically linked to varying degrees of survival outcomes. Nevertheless, few hotspots have been analyzed (e.g., TP53 and EGFR). Thus, many other genes and hotspots remain unexplored.

METHODS

We systematically screened over 1400 hotspots across 33 TCGA cancer types. We compared the patients carrying a hotspot against (i) all cases, (ii) gene-mutated cases, (iii) other mutated hotspots, or (iv) specific hotspots. Due to the limited number of samples in hotspots and the inherent group imbalance, besides Cox models and the log-rank test, we employed VALORATE to estimate their association with survival precisely.

RESULTS

We screened 1469 hotspots in 6451 comparisons, where 314 were associated with survival. Many are discussed and linked to the current literature. Our findings demonstrate associations between known hotspots and survival while also revealing more potential hotspots. To enhance accessibility and promote further investigation, all the Kaplan-Meier curves, the log-rank tests, Cox statistics, and VALORATE-estimated null distributions are accessible on our website.

CONCLUSIONS

Our analysis revealed both known and putatively novel hotspots associated with survival, which can be used as biomarkers. Our web resource is a valuable tool for cancer research.

Nucleosome conformation dictates the histone code

Histone post-translational modifications (PTMs) play a critical role in chromatin regulation. It has been proposed that these PTMs form localized 'codes' that are read by specialized regions (reader domains) in chromatin-associated proteins (CAPs) to regulate downstream function. Substantial effort has been made to define [CAP: histone PTM] specificities, and thus decipher the histone code and guide epigenetic therapies. However, this has largely been done using the reductive approach of isolated reader domains and histone peptides, which cannot account for any higher-order factors. Here, we show that the [BPTF PHD finger and bromodomain: histone PTM] interaction is dependent on nucleosome context. The tandem reader selectively associates with nucleosomal H3K4me3 and H3K14ac or H3K18ac, a combinatorial engagement that despite being in cis is not predicted by peptides. This in vitro specificity of the BPTF tandem reader for PTM-defined nucleosomes is recapitulated in a cellular context. We propose that regulatable histone tail accessibility and its impact on the binding potential of reader domains necessitates we refine the 'histone code' concept and interrogate it at the nucleosome level.

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.

Genome-wide survey of D/E repeats in human proteins uncovers their instability and aids in identifying their role in the chromatin regulator ATAD2

D/E repeats are stretches of aspartic and/or glutamic acid residues found in over 150 human proteins. We examined genomic stability of D/E repeats and functional characteristics of D/E repeat-containing proteins vis-à-vis the proteins with poly-Q or poly-A repeats, which are known to undergo pathologic expansions. Mining of tumor sequencing data revealed that D/E repeat-coding regions are similar to those coding poly-Qs and poly-As in increased incidence of trinucleotide insertions/deletions but differ in types and incidence of substitutions. D/E repeat-containing proteins preferentially function in chromatin metabolism and are the more likely to be nuclear and interact with core histones, the longer their repeats are. One of the longest D/E repeats of unknown function is in ATAD2, a bromodomain family ATPase frequently overexpressed in tumors. We demonstrate that D/E repeat deletion in ATAD2 suppresses its binding to nascent and mature chromatin and to the constitutive pericentromeric heterochromatin, where ATAD2 represses satellite transcription.

•

Many human proteins contain runs of aspartic/glutamic acid residues (D/E repeats)

•

D/E repeats show increased incidence of in-frame insertions/deletions in tumors

•

Nuclear and histone-interacting proteins often have long D/E repeats

•

D/E repeat of the oncogene ATAD2 controls its binding to pericentric chromatin

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.

ATAD2 drives colorectal cancer progression by regulating TRIM25 expression via a positive feedback loop with E2F transcriptional factors

ATPase family AAA domain-containing protein 2 (ATAD2) is highly expressed in a variety of cancer types, and acts as a co-activator of androgen and estrogen receptors, as well as MYC and E2F transcription factors, to promote tumor cell proliferation. However, the regulation of ATAD2 and its related mechanisms are still elusive. Here, we show that ATAD2 protein was stabilized during DNA damage response in colorectal cancer (CRC) cells. TRIM25, an oncogenic ubiquitin E3 ligase, can interact with ATAD2 and stabilize ATAD2 upon genotoxic insult. We further demonstrated that ATAD2 played a tumor promoting role in CRC and acted as a transcriptional co-activator of E2Fs to promote the expression of TRIM25. Thus, our results revealed an unknown ATAD2-E2Fs-TRIM25 positive feedback loop that drove CRC progression.

The dCypher Approach to Interrogate Chromatin Reader Activity Against Posttranslational Modification-Defined Histone Peptides and Nucleosomes

Bulk chromatin encompasses complex sets of histone posttranslational modifications (PTMs) that recruit (or repel) the diverse reader domains of Chromatin-Associated Proteins (CAPs) to regulate genome processes (e.g., gene expression, DNA repair, mitotic transmission). The binding preference of reader domains for their PTMs mediates localization and functional output, and are often dysregulated in disease. As such, understanding chromatin interactions may lead to novel therapeutic strategies, However the immense chemical diversity of histone PTMs, combined with low-throughput, variable, and nonquantitative methods, has defied accurate CAP characterization. This chapter provides a detailed protocol for dCypher, a novel approach for the rapid, quantitative interrogation of CAPs (as mono- or multivalent Queries) against large panels (10s to 100s) of PTM-defined histone peptide and semisynthetic nucleosomes (the potential Targets). We describe key optimization steps and controls to generate robust binding data. Further, we compare the utility of histone peptide and nucleosome substrates in CAP studies, outlining important considerations in experimental design and data interpretation.

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.

Engineering an acetyllysine reader with a photocrosslinking amino acid for interactome profiling

The site-specific installation of light-activable crosslinker unnatural amino acids offers a powerful approach to trap transient protein-protein interactions both in vitro and in vivo. Herein, we engineer a bromodomain to introduce 4-benzoyl-L-phenylalanine (BzF) using amber suppressor mutagenesis without compromising its ability to recognize the acetylated histone proteins. We demonstrate the high crosslinking efficiency of the engineered reader towards the interacting partners and its suitability for profiling the transient bromodomain interactome.

Coordination of Di-Acetylated Histone Ligands by the ATAD2 Bromodomain

The ATPase Family, AAA domain-containing protein 2 (ATAD2) bromodomain (BRD) has a canonical bromodomain structure consisting of four α-helices. ATAD2 functions as a co-activator of the androgen and estrogen receptors as well as the MYC and E2F transcription factors. ATAD2 also functions during DNA replication, recognizing newly synthesized histones. In addition, ATAD2 is shown to be up-regulated in multiple forms of cancer including breast, lung, gastric, endometrial, renal, and prostate. Furthermore, up-regulation of ATAD2 is strongly correlated with poor prognosis in many types of cancer, making the ATAD2 bromodomain an innovative target for cancer therapeutics. In this study, we describe the recognition of histone acetyllysine modifications by the ATAD2 bromodomain. Residue-specific information on the complex formed between the histone tail and the ATAD2 bromodomain, obtained through nuclear magnetic resonance spectroscopy (NMR) and X-ray crystallography, illustrates key residues lining the binding pocket, which are involved in coordination of di-acetylated histone tails. Analytical ultracentrifugation, NMR relaxation data, and isothermal titration calorimetry further confirm the monomeric state of the functionally active ATAD2 bromodomain in complex with di-acetylated histone ligands. Overall, we describe histone tail recognition by ATAD2 BRD and illustrate that one acetyllysine group is primarily engaged by the conserved asparagine (N1064), the "RVF" shelf residues, and the flexible ZA loop. Coordination of a second acetyllysine group also occurs within the same binding pocket but is essentially governed by unique hydrophobic and electrostatic interactions making the di-acetyllysine histone coordination more specific than previously presumed.

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.