Discovery of CREBBP Bromodomain Inhibitors by High-Throughput Docking and Hit Optimization Guided by Molecular Dynamics

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.

Enriching screening libraries with bioactive fragment space

By deconvoluting 238,073 bioactive molecules in the ChEMBL library into extended Murcko ring systems, we identified a set of 2245 ring systems present in at least 10 molecules. These ring systems belong to 2221 clusters by ECFP4 fingerprints with a minimum intracluster similarity of 0.8. Their overlap with ring systems in commercial libraries was further quantified. Our findings suggest that success of a small fragment library is driven by the convergence of effective coverage of bioactive ring systems (e.g., 10% coverage by 1000 fragments vs. 40% by 2million HTS compounds), high enrichment of bioactive ring systems, and low molecular complexity enhancing the probability of a match with the protein targets. Reconciling with the previous studies, bioactive ring systems are underrepresented in screening libraries. As such, we propose a library of virtual fragments with key functionalities via fragmentation of bioactive molecules. Its utility is exemplified by a prospective application on protein kinase CK2, resulting in the discovery of a series of novel inhibitors with the most potent compound having an IC50 of 0.5μM and a ligand efficiency of 0.41kcal/mol per heavy atom.

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

The biological role played by non-BET bromodomains remains poorly understood, and it is therefore imperative to identify potent and highly selective inhibitors to effectively explore the biology of individual bromodomain proteins. A ligand-efficient nonselective bromodomain inhibitor was identified from a 6-methyl pyrrolopyridone fragment. Small hydrophobic substituents replacing the N-methyl group were designed directing toward the conserved bromodomain water pocket, and two distinct binding conformations were then observed. The substituents either directly displaced and rearranged the conserved solvent network, as in BRD4(1) and TAF1(2), or induced a narrow hydrophobic channel adjacent to the lipophilic shelf, as in BRD9 and CECR2. The preference of distinct substituents for individual bromodomains provided selectivity handles useful for future lead optimization efforts for selective BRD9, CECR2, and TAF1(2) inhibitors.

High-Throughput Fragment Docking into the BAZ2B Bromodomain: Efficient in Silico Screening for X-Ray Crystallography

Bromodomains are protein modules that bind to acetylated lysine side chains in histones and other proteins. The bromodomain adjacent to zinc finger domain protein 2B (BAZ2B) has been reported to be poorly druggable. Here, we screened an in-house library of 350 fragments by automatic docking to the BAZ2B bromodomain. The top 12 fragments according to the predicted binding energy were selected for experiments of soaking into apo crystals of BAZ2B which yielded the structure of the complex for four of them, which is a hit rate of 33%. Additional crystal structures were solved for BAZ2B and two scaffolds identified by analogy. For three topologically similar fragments, the crystal structures reveal binding modes with different penetration, i.e., with zero, one, and two water molecules, respectively, located between the fragment and the side chain of a conserved tyrosine (Tyr1901) in the bottom of the acetyl lysine pocket of BAZ2B. Furthermore, a remarkable stereoselectivity of the acetyl lysine pocket emerges from the crystal structures of the bromodomains of BAZ2B and SMARCA4 in complex with the chiral diol MPD (2-methyl-2,4-pentanediol).

Docking Screens for Novel Ligands Conferring New Biology

It is now plausible to dock libraries of 10 million molecules against targets over several days or weeks. When the molecules screened are commercially available, they may be rapidly tested to find new leads. Although docking retains important liabilities (it cannot calculate affinities accurately nor even reliably rank order high-scoring molecules), it can often can distinguish likely from unlikely ligands, often with hit rates above 10%. Here we summarize the improvements in libraries, target quality, and methods that have supported these advances, and the open access resources that make docking accessible. Recent docking screens for new ligands are sketched, as are the binding, crystallographic, and in vivo assays that support them. Like any technique, controls are crucial, and key experimental ones are reviewed. With such controls, docking campaigns can find ligands with new chemotypes, often revealing the new biology that may be docking's greatest impact over the next few years.

4-Acyl Pyrrole Derivatives Yield Novel Vectors for Designing Inhibitors of the Acetyl-Lysine Recognition Site of BRD4(1)

Several human diseases, including cancer, show altered signaling pathways resulting from changes in the activity levels of epigenetic modulators. In the past few years, small-molecule inhibitors against specific modulators, including the bromodomain and extra-terminal (BET) bromodomain family of acetylation readers, have shown early promise in the treatment of the genetically defined midline carcinoma and hematopoietic malignancies. We have recently developed a novel potent inhibitor of BET proteins, 1 (XD14[ Angew. Chem., Int. Ed. 2013, 52, 14055]), which exerts a strong inhibitory potential on the proliferation of specific leukemia cell lines. In the study presented here, we designed analogues of 1 to study the potential of substitutions on the 4-acyl pyrrole backbone to occupy additional sites within the substrate recognition site of BRD4(1). The compounds were profiled using ITC, DSF, and X-ray crystallography. We could introduce several substitutions that address previously untargeted areas of the substrate recognition site. This work may substantially contribute to the development of therapeutics with increased target specificity against BRD4-related malignancies.

Discovery of Benzo[cd]indol-2(1H)-ones as Potent and Specific BET Bromodomain Inhibitors: Structure-Based Virtual Screening, Optimization, and Biological Evaluation

The discovery of inhibitors of bromodomain and extra terminal domain (BET) has achieved great progress, and at least seven inhibitors have progressed into clinical trials for the treatment of cancer or inflammatory diseases. Here, we describe the identification, optimization, and evaluation of benzo[cd]indol-2(1H)-one containing compounds as a new class of BET bromodomain inhibitors, starting from structure-based virtual screening (SBVS). Through structure-based optimization, potent compounds were obtained with significantly improved activity. The two most potent compounds bind to the BRD4 bromodomain, with Kd values of 124 and 137 nM. Selected compounds exhibited high selectivity over other non-BET subfamily members. Notably, compound 85 demonstrated a reasonable antiproliferation effect on MV4;11 leukemia cells and exhibited a good pharmacokinetic profile with high oral bioavailability (75.8%) and moderate half-life (T1/2 = 3.95 h). The resulting lead molecule 85 represents a new, potent, and selective class of BET bromodomain inhibitors for the development of therapeutics to treat cancer and inflammatory diseases.

Cation–π interactions in CREBBP bromodomain inhibition: an electrostatic model for small-molecule binding affinity and selectivity

A new model is presented to explain and predict binding affinity of aromatic and heteroaromatic ligands for the CREBBP bromodomain based on cation–π interaction strength.

Disrupting Acetyl-Lysine Recognition: Progress in the Development of Bromodomain Inhibitors

Bromodomains, small protein modules that recognize acetylated lysine on histones, play a significant role in the epigenome, where they function as "readers" that ultimately determine the functional outcome of the post-translational modification. Because the initial discovery of selective BET inhibitors have helped define the role of that protein family in oncology and inflammation, BET bromodomains have continued to garner the most attention of any other bromodomain. More recently, non-BET bromodomain inhibitors that are potent and selective have been disclosed for ATAD2, CBP, BRD7/9, BRPF, BRPF/TRIM24, CECR2, SMARCA4, and BAZ2A/B. Such novel inhibitors can be used to probe the physiological function of these non-BET bromodomains and further understanding of their role in certain disease states. Here, we provide an update to the progress in identifying selective bromodomain inhibitors and their use as biological tools, as well as our perspective on the field.

Fragment-Based Design of Selective Nanomolar Ligands of the CREBBP Bromodomain

Novel ligands of the CREBBP bromodomain were identified by fragment-based docking. The in silico discovered hits have been optimized by chemical synthesis into selective nanomolar compounds, thereby preserving the ligand efficiency. The selectivity for the CREBBP bromodomain over other human bromodomain subfamilies has achieved by a benzoate moiety which was predicted by docking to be involved in favorable electrostatic interactions with the Arg1173 side chain, a prediction that could be verified a posteriori by the high-resolution crystal structure of the CREBBP bromodomain in complex with ligand 6 and also by MD simulations (see Xu, M.; Unzue, A.; Dong, J.; Spiliotopoulos, D.; Nevado, C.; Caflisch, A. Discovery of CREBBP bromodomain inhibitors by high-throughput docking and hit optimization guided by molecular dynamics. J. Med. Chem. 2015, DOI: 10.1021/acs.jmedchem.5b00171).