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

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.

Identification of novel natural product inhibitors of BRD4 using high throughput virtual screening and MD simulation

Bromodomains are evolutionarily conserved structural motifs that recognize acetylated lysine residues on histone tails. They play a crucial role in shaping chromatin architecture and regulating gene expression in various biological processes. Mutations in bromodomains containing proteins lead to multiple human diseases, which makes them attractive target for therapeutic intervention. Extensive studies have been done on BRD4 as a target for several cancers, such as Acute Myeloid Leukemia (AML) and Burkitt Lymphoma. Several potential inhibitors have been identified against the BRD4 bromodomain. However, most of these inhibitors have drawbacks such as non-specificity and toxicity, decreasing their appeal and necessitating the search for novel non-toxic inhibitors. This study aims to address this need by virtually screening natural compounds from the NPASS database against the Kac binding site of BRD4-BD1 using high throughput molecular docking followed by similarity clustering, pharmacokinetic screening, MD simulation and MM-PBSA binding free energy calculations. Using this approach, we identified five natural product inhibitors having a similar or better binding affinity to the BRD4 bromodomain compared to JQ1 (previously reported inhibitor of BRD4). Further systematic analysis of these inhibitors resulted in the top three hits: NPC268484 (Palodesangren-B), NPC295021 (Candidine) and NPC313112 (Buxifoliadine-D). Collectively, our in silico results identified some promising natural products that have the potential to act as potent BRD4-BD1 inhibitors and can be considered for further validation through future in vitro and in vivo studies.Communicated by Ramaswamy H. Sarma.

A comprehensive atlas of epigenetic regulators reveals tissue-specific epigenetic regulation patterns

Epigenetic machinery contributes to gene regulation in eukaryotic species. However, the machinery including more than 600 epigenetic regulator (ER) genes responsible for reading, writing, and erasing histone modifications and DNA modifications remains largely uncharacterized across species. We compile a comprehensive list of ERs based on an evolutionary analysis across 23 species, which is the most comprehensive ER list in various species until recently. We further perform comparative transcriptomic analyses across different tissues in humans, mice, as well as other amniote species. We observe a consistent tissue-of-origin expression specificity pattern of duplicated ER genes across species and suggest links between expression specificity and ER gene evolution as well as ER function. Additional analyses further suggest that ER duplication can generate tissue-specific ER genes with the same epigenetic substrates, which may be closely related to their regulatory specificity in tissue development. Our work can serve as a foundation to better comprehend the tissue-specific expression patterns of ER genes from an evolutionary perspective and also the functional implications of ERs in tissue-specific epigenetic regulation.

SIN-3 acts in distinct complexes to regulate the germline transcriptional program in Caenorhabditis elegans

The transcriptional co-regulator SIN3 influences gene expression through multiple interactions that include histone deacetylases. Haploinsufficiency and mutations in SIN3 are the underlying cause of Witteveen-Kolk syndrome and related intellectual disability and autism syndromes, emphasizing its key role in development. However, little is known about the diversity of its interactions and functions in developmental processes. Here, we show that loss of SIN-3, the single SIN3 homolog in Caenorhabditis elegans, results in maternal-effect sterility associated with de-regulation of the germline transcriptome, including de-silencing of X-linked genes. We identify at least two distinct SIN3 complexes containing specific histone deacetylases and show that they differentially contribute to fertility. Single-cell, single-molecule fluorescence in situ hybridization reveals that in sin-3 mutants the X chromosome becomes re-expressed prematurely and in a stochastic manner in individual germ cells, suggesting a role for SIN-3 in its silencing. Furthermore, we identify histone residues whose acetylation increases in the absence of SIN-3. Together, this work provides a powerful framework for the in vivo study of SIN3 and associated proteins.

Uncovering the Domain-Specific Interactome of the TAF1 Tandem Reader Using Site-Specific Azide-Acetyllysine Photochemistry

Combinatorial readout of histone post-translational modifications by tandem reader modules mediates crosstalk among different histone modifications. To identify the domain-specific interactome of the tandem reader, we engineered the dual bromodomain of TATA-binding protein-associated factor-1 (TAF1) to carry a photoactivatable unnatural amino acid, 4-azido-l-phenylalanine (AzF), via amber suppressor mutagenesis. Using computational approaches, we modeled the targeted residues of TAF1 with AzF to predict the cross-linking distance between the reactive arylazide and its interacting partner. We developed three photoactivatable TAF1 tandem-bromodomain analogues, viz., Y1403AzF in bromodomain 1 (BD1), W1526AzF in bromodomain 2 (BD2), and Y1403AzF/W1526AzF in both BD1 and BD2. Circular dichroism and a thermal shift assay were used to confirm the structural integrity of the engineered readers. Using the TAF1 tandem-bromodomain analogues, we characterized their histone ligand binding properties by isothermal titration calorimetry and photo-cross-linking experiments. We found that the dual bromodomain of TAF1 independently binds and cross-links to different acetylated histone ligands. We further used the engineered BD1 and BD2 analogues of the TAF1 tandem readers to identify their domain-specific interacting partners at the cellular level. Both BD1 and BD2 independently cross-link to a unique interactome, and importantly, the dual cross-linker carrying TAF1 analogue could capture both BD1- and BD2-specific interactomes. Our work concludes that BD1 and BD2 of the TAF1 tandem reader independently recognize their interacting partners to regulate downstream cellular functions.

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.

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.

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

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