Binding Mode of Acetylated Histones to Bromodomains: Variations on a Common Motif

Deciphering the structure, function, and mechanism of lysine acetyltransferase cGNAT2 in cyanobacteria

Lysine acetylation is a conserved regulatory posttranslational protein modification that is performed by lysine acetyltransferases (KATs). By catalyzing the transfer of acetyl groups to substrate proteins, KATs play critical regulatory roles in all domains of life; however, no KATs have yet been identified in cyanobacteria. Here, we tested all predicted KATs in the cyanobacterium Synechococcus sp. PCC 7002 (Syn7002) and demonstrated that A1596, which we named cyanobacterial Gcn5-related N-acetyltransferase (cGNAT2), can catalyze lysine acetylation in vivo and in vitro. Eight amino acid residues were identified as the key residues in the putative active site of cGNAT2, as indicated by structural simulation and site-directed mutagenesis. The loss of cGNAT2 altered both growth and photosynthetic electron transport in Syn7002. In addition, quantitative analysis of the lysine acetylome identified 548 endogenous substrates of cGNAT2 in Syn7002. We further demonstrated that cGNAT2 can acetylate NAD(P)H dehydrogenase J (NdhJ) in vivo and in vitro, with the inability to acetylate K89 residues, thus decreasing NdhJ activity and affecting both growth and electron transport in Syn7002. In summary, this study identified a KAT in cyanobacteria and revealed that cGNAT2 regulates growth and photosynthesis in Syn7002 through an acetylation-mediated mechanism.

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.

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.

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

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

Emerging oncogene ATAD2: Signaling cascades and therapeutic initiatives

ATAD2 is a promising oncoprotein with tumor-promoting functions in many cancers. It is a valid cancer drug-target and a potential cancer-biomarker for multiple malignancies. As a cancer/testis antigen (CTA), ATAD2 could also be a probable candidate for immunotherapy. It is a unique CTA that belongs to both AAA+ ATPase and bromodomain family proteins. Since 2007, several research groups have been reported on the pleiotropic oncogenic functions of ATAD2 in diverse signaling pathways, including Rb/E2F-cMyc pathway, steroid hormone signaling pathway, p53 and p38-MAPK-mediated apoptotic pathway, AKT pathway, hedgehog signaling pathway, HIF1α signaling pathway, and Epithelial to Mesenchymal Transition (EMT) pathway in various cancers. In all these pathways, ATAD2 participates in chromatin dynamics, DNA replication, and gene transcription, demonstrating its role as an epigenetic reader and transcription factor or coactivator to promote tumorigenesis. However, despite the progress, an overall mechanism of ATAD2-mediated oncogenesis in diverse origin is elusive. In this review, we summarize the accumulated evidence to envision the overall ATAD2 signaling networks during carcinogenesis and highlight the area where missing links await further research. Besides, the structure-function aspect of ATAD2 is also discussed. Since the efforts have already been initiated to explore targeted drug molecules and RNA-based therapeutic alternatives against ATAD2, their potency and prospects have been elucidated. Together, we believe this is a well-rounded review on ATAD2, facilitating a new drift in ATAD2 research, essential for its clinical implication as a biomarker and/or cancer drug-target.

An ABSINTH-Based Protocol for Predicting Binding Affinities between Proteins and Small Molecules

The core task in computational drug discovery is to accurately predict binding free energies in receptor-ligand systems for large libraries of putative binders. Here, the ABSINTH implicit solvent model and force field are extended to describe small, organic molecules and their interactions with proteins. We show that an automatic pipeline based on partitioning arbitrary molecules into substructures corresponding to model compounds with known free energies of solvation can be combined with the CHARMM general force field into a method that is successful at the two important challenges a scoring function faces in virtual screening work flows: it ranks known binders with correlation values rivaling that of comparable state-of-the-art methods and it enriches true binders in a set of decoys. Our protocol introduces innovative modifications to common virtual screening workflows, notably the use of explicit ions as competitors and the integration over multiple protein and ligand species differing in their protonation states. We demonstrate the value of modifications to both the protocol and ABSINTH itself. We conclude by discussing the limitations of high-throughput implicit methods such as the one proposed here.

Substrate Affinity and Specificity of the ScSth1p Bromodomain Are Fine-Tuned for Versatile Histone Recognition

Bromodomains recognize a wide range of acetylated lysines in histones and other nuclear proteins. Substrate specificity is critical for their biological function and arises from unique acetyl-lysine binding sites formed by variable loop regions. Here, we analyzed substrate affinity and specificity of the yeast ScSth1p bromodomain, an essential component of the "Remodels the Structure of Chromatin" complex, and found that the wild-type bromodomain preferentially recognizes H3K14ac and H4K20ac peptides. Mutagenesis studies-guided by our crystal structure determined at 2.7-Å resolution-revealed loop residues Ser1276 and Trp1338 as key determinants for such interactions. Strikingly, point mutations of each of these residues substantially increased peptide binding affinity and selectivity, respectively. Our data demonstrate that the ScSth1p bromodomain is not optimized for binding to an individual acetylation mark, but fine-tuned for interactions with several such modifications, consistent with the versatile and multivalent nature of histone recognition by reader modules such as bromodomains.

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.

HDAC stimulates gene expression through BRD4 availability in response to IFN and in interferonopathies

In contrast to the common role of histone deacetylases (HDACs) for gene repression, HDAC activity provides a required positive function for IFN-stimulated gene (ISG) expression. Here, we show that HDAC1/2 as components of the Sin3A complex are required for ISG transcriptional elongation but not for recruitment of RNA polymerase or transcriptional initiation. Transcriptional arrest by HDAC inhibition coincides with failure to recruit the epigenetic reader Brd4 and elongation factor P-TEFb due to sequestration of Brd4 on hyperacetylated chromatin. Brd4 availability is regulated by an equilibrium cycle between opposed acetyltransferase and deacetylase activities that maintains a steady-state pool of free Brd4 available for recruitment to inducible promoters. An ISG expression signature is a hallmark of interferonopathies and other autoimmune diseases. Combined inhibition of HDAC1/2 and Brd4 resolved the aberrant ISG expression detected in cells derived from patients with two inherited interferonopathies, ISG15 and USP18 deficiencies, defining a novel therapeutic approach to ISG-associated autoimmune diseases.

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.

Focused conformational sampling in proteins

A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.

Leveraging Epigenetics to Enhance the Cellular Response to Chemotherapies and Improve Tumor Immunogenicity

Cancer chemotherapeutic drugs have greatly advanced our ability to successfully treat a variety of human malignancies. The different forms of stress produced by these agents in cancer cells result in both cell autonomous and cell nonautonomous effects. Desirable cell autonomous effects include reduced proliferative potential, cellular senescence, and cell death. More recently recognized cell nonautonomous effects, usually in the form of stimulating an antitumor immune response, have significant roles in therapeutic efficiency for a select number of chemotherapies. Unfortunately, the success of these therapeutics is not universal as not all tumors respond to treatment, and those that do respond will frequently relapse into therapy-resistant disease. Numerous strategies have been developed to sensitize tumors toward chemotherapies as a means to either improve initial responses, or serve as a secondary treatment strategy for therapy-resistant disease. Recently, targeting epigenetic regulators has emerged as a viable method of sensitizing tumors to the effects of chemotherapies, many of which are cytotoxic. In this review, we summarize these strategies and propose a path for future progress.

Mechanistic insights into peptide and ligand binding of the ATAD2-bromodomain via atomistic simulations disclosing a role of induced fit and conformational selection

Atomistic simulations of the ATAD2-bromodomain disclose a role of induced fit and conformational selection upon ligand and peptide binding.

Discovery of Inhibitors of Four Bromodomains by Fragment-Anchored Ligand Docking

The high-throughput docking protocol called ALTA-VS (anchor-based library tailoring approach for virtual screening) was developed in 2005 for the efficient in silico screening of large libraries of compounds by preselection of only those molecules that have optimal fragments (anchors) for the protein target. Here we present an updated version of ALTA-VS with a broader range of potential applications. The evaluation of binding energy makes use of a classical force field with implicit solvent in the continuum dielectric approximation. In about 2 days per protein target on a 96-core compute cluster (equipped with Xeon E3-1280 quad core processors at 2.5 GHz), the screening of a library of nearly 77 000 diverse molecules with the updated ALTA-VS protocol has resulted in the identification of 19, 3, 3, and 2 μM inhibitors of the human bromodomains ATAD2, BAZ2B, BRD4(1), and CREBBP, respectively. The success ratio (i.e., number of actives in a competition binding assay in vitro divided by the number of compounds tested) ranges from 8% to 13% in dose-response measurements. The poses predicted by fragment-based docking for the three ligands of the BAZ2B bromodomain were confirmed by protein X-ray crystallography.

The ATAD2 bromodomain binds different acetylation marks on the histone H4 in similar fuzzy complexes

Bromodomains are protein modules adopting conserved helix bundle folds. Some bromodomain-containing proteins, such as ATPase family AAA domain-containing protein 2 (ATAD2), isoform A, have attracted much interest because they are overexpressed in many types of cancer. Bromodomains bind to acetylated lysine residues on histone tails and thereby facilitate the reading of the histone code. Epigenetic regulators in general have been implicated as indicators, mediators, or causes of a large number of diseases and disorders. To interfere with or modulate these processes, it is therefore of fundamental interest to understand the molecular mechanisms by which epigenetic regulation occurs. Here, we present results from molecular dynamics simulations of a doubly acetylated histone H4 peptide bound to the bromodomain of ATAD2 (hereafter referred to as ATAD2A). These simulations revealed how the flexibility of ATAD2A's major loop, the so-called ZA loop, creates an adaptable interface that preserves the disorder of both peptide and loop in the bound state. We further demonstrate that the binding involves an almost identical average pattern of interactions irrespective of which acetyl mark is inserted into the pocket. In conjunction with a likely mechanism of electrostatically driven recruitment, our simulation results highlight how the bromodomain is built toward promiscuous binding with low specificity. In conclusion, the simulations indicate that disorder and electrostatic steering function jointly to recruit ATAD2A to the histone core and that these fuzzy interactions may promote cooperativity between nearby epigenetic marks.

tBRD-1 and tBRD-2 regulate expression of genes necessary for spermatid differentiation

Male germ cell differentiation proceeds to a large extent in the absence of active gene transcription. In Drosophila, hundreds of genes whose proteins are required during post-meiotic spermatid differentiation (spermiogenesis) are transcribed in primary spermatocytes. Transcription of these genes depends on the sequential action of the testis meiotic arrest complex (tMAC), Mediator complex, and testis-specific TFIID (tTFIID) complex. How the action of these protein complexes is coordinated and which other factors are involved in the regulation of transcription in spermatocytes is not well understood. Here, we show that the bromodomain proteins tBRD-1 and tBRD-2 regulate gene expression in primary spermatocytes and share a subset of target genes. The function of tBRD-1 was essential for the sub-cellular localization of endogenous tBRD-2 but dispensable for its protein stability. Our comparison of different microarray data sets showed that in primary spermatocytes, the expression of a defined number of genes depends on the function of the bromodomain proteins tBRD-1 and tBRD-2, the tMAC component Aly, the Mediator component Med22, and the tTAF Sa.

Derivatives of 3-Amino-2-methylpyridine as BAZ2B Bromodomain Ligands: In Silico Discovery and in Crystallo Validation

The 3-amino-2-methylpyridine derivative 1 was identified as ligand of the BAZ2B bromodomain by automatic docking of nearly 500 compounds, selected on the basis of previous fragment hits. Hit expansion by two in silico approaches, pharmacophore search followed by docking, and substructure search resulted in five additional ligands. The predicted binding mode of the six 3-amino-2-methylpyridine derivatives was validated by protein crystallography. A small displacement of residues 1894-1899 of the ZA loop is observed for two of the six ligands. In all structures, the pyridine head is involved in a water-mediated hydrogen bond with the side chain of the conserved Tyr1901 while the 3-amino linker acts as hydrogen bond donor for the backbone carbonyl of Pro1888. Heterogeneous orientations are observed for the tail groups (i.e., the 3-amino substituents). The sulfonyl group in the tail of compounds 1 and 2 is involved in a hydrogen bond with the backbone amide of Asn1894.

Exploring the role of water in molecular recognition: predicting protein ligandability using a combinatorial search of surface hydration sites

The interaction between any two biological molecules must compete with their interaction with water molecules. This makes water the most important molecule in medicine, as it controls the interactions of every therapeutic with its target. A small molecule binding to a protein is able to recognize a unique binding site on a protein by displacing bound water molecules from specific hydration sites. Quantifying the interactions of these water molecules allows us to estimate the potential of the protein to bind a small molecule. This is referred to as ligandability. In the study, we describe a method to predict ligandability by performing a search of all possible combinations of hydration sites on protein surfaces. We predict ligandability as the summed binding free energy for each of the constituent hydration sites, computed using inhomogeneous fluid solvation theory. We compared the predicted ligandability with the maximum observed binding affinity for 20 proteins in the human bromodomain family. Based on this comparison, it was determined that effective inhibitors have been developed for the majority of bromodomains, in the range from 10 to 100 nM. However, we predict that more potent inhibitors can be developed for the bromodomains BPTF and BRD7 with relative ease, but that further efforts to develop inhibitors for ATAD2 will be extremely challenging. We have also made predictions for the 14 bromodomains with no reported small molecule K d values by isothermal titration calorimetry. The calculations predict that PBRM1(1) will be a challenging target, while others such as TAF1L(2), PBRM1(4) and TAF1(2), should be highly ligandable. As an outcome of this work, we assembled a database of experimental maximal K d that can serve as a community resource assisting medicinal chemistry efforts focused on BRDs. Effective prediction of ligandability would be a very useful tool in the drug discovery process.

Selectivity on-target of bromodomain chemical probes by structure-guided medicinal chemistry and chemical biology

Targeting epigenetic proteins is a rapidly growing area for medicinal chemistry and drug discovery. Recent years have seen an explosion of interest in developing small molecules binding to bromodomains, the readers of acetyl-lysine modifications. A plethora of co-crystal structures has motivated focused fragment-based design and optimization programs within both industry and academia. These efforts have yielded several compounds entering the clinic, and many more are increasingly being used as chemical probes to interrogate bromodomain biology. High selectivity of chemical probes is necessary to ensure biological activity is due to an on-target effect. Here, we review the state-of-the-art of bromodomain-targeting compounds, focusing on the structural basis for their on-target selectivity or lack thereof. We also highlight chemical biology approaches to enhance on-target selectivity.

Fragment-based in silico screening of bromodomain ligands

We review the results of fragment-based high-throughput docking to the N-terminal bromodomain of BRD4 and the CREBBP bromodomain. In both docking campaigns the ALTA (anchor-based library tailoring) procedure was used to reduce the size of the initial library by selecting for flexible docking only the molecules that contain a fragment with favorable predicted binding energy. Ranking by a force field-based energy with solvation has resulted in small-molecule hits with low-micromolar affinity and favorable ligand efficiency. Importantly, the binding modes predicted by docking have been validated by X-ray crystallography. One of the hits for the CREBBP bromodomain has been optimized by medicinal chemistry into a series of potent and selective ligands.

The bromodomain inhibitor N-methyl pyrrolidone reduced fat accumulation in an ovariectomized rat model

Background

Weight gain is one of the consequences of estrogen deficiency and constitutes a major health problem. The present study highlights the effects of N-methyl pyrrolidone (NMP) on adipogenesis in osteoporosis induced by estrogen deficiency in an ovariectomized rat model.

Results

Ovariectomy resulted in body weight gain, increased femoral marrow adipocytes, and hypertrophic adipocytes in white adipose tissue, distorted serum leptin, and TNF-α and PPARγ levels. Treatment with NMP normalized these parameters similar to the control group. In vitro, NMP inhibited the differentiation of 3T3-L1 pre-adipocytes and hMSCs, indicating its anti-adipogenic effect. Moreover, PPARγ was significantly reduced with NMP treatment in in vivo and in vitro experiments. NMP inhibited BRD2 and BRD4 binding in an AlphaScreen assay, with an IC50 of 3 and 4 mM, respectively. The effect of NMP was consistent with its role as a bromodomain inhibitor.

Conclusions

Our data indicates that NMP inhibits the adipogenic effect of estrogen deficiency at the level of PPARγ expression by BRD4 inhibition.

Electronic supplementary material

The online version of this article (doi:10.1186/s13148-016-0209-2) contains supplementary material, which is available to authorized users.