A series of potent CREBBP bromodomain ligands reveals an induced-fit pocket stabilized by a cation-π interaction

Bioisoteres for carboxylic acids: From ionized isosteres to novel unionized replacements

Carboxylic acids are key pharmacophoric elements in many molecules. They can be seen as a problem by some, due to perceived permeability challenges, potential for high plasma protein binding and the risk of forming reactive metabolites due to acyl-glucuronidation. By others they are viewed more favorably as they can decrease lipophilicity by adding an ionizable center which can be beneficial for solubility, and can add enthalpic interactions with the target protein. However, there are many instances where the replacement of a carboxylic acid with a bioisosteric group is required. This has led to the development of a number of ionizable groups which sufficiently mimic the carboxylic acid functionality whilst improving, for example, the metabolic profile of the molecule in question. An alternative strategy involves replacement of the carboxylate by neutral functional groups. This review initially details carefully selected examples whereby tetrazoles, acyl sulfonamides or isoxazolols have been beneficially utilized as carboxylic acid bioisosteres altering physicohemical properties, interactions with the target and metabolism and/or pharmacokinetics, before delving further into the binding mode of carboxylic acid derivatives with their target proteins. This analysis highlights new ways to consider the replacement of carboxylic acids by neutral bioisosteric groups which either rely on hydrogen bonds or cation-π interactions. It should serve as a useful guide for scientists working in drug discovery.

Synthesis of Dihydroquinoxalinones from Biomass-Derived Keto Acids and o-Phenylenediamines

A novel catalyst-free cascade amination/cyclization/reduction reaction was developed for the synthesis of various Dihydroquinoxalinones under mild conditions from accessible biomass-derived keto acids and 1,2-phenylenediamines with ammonia borane as a hydrogen donor. This single-step approach enables a simple and eco-friendly route toward the direct synthesis of 12 kinds of Dihydroquinoxalinones in moderate to excellent yields in the green solvent dimethyl carbonate. The results of deuterium-labeling experiments and density function calculations demonstrate that the reductive process proceeds along a double hydrogen transfer pathway. An acceptable yield of Dihydroquinoxalinone can be afforded in a gram-scale experiment, illustrating the practicality of the as-reported reaction system.

Discovery of a potent and selective CBP bromodomain inhibitor (Y08262) for treating acute myeloid leukemia

The bromodomain of CREB (cyclic-AMP response element binding protein) binding protein (CBP) is an epigenetic "reader" and plays a key role in transcriptional regulation. CBP bromodomain is considered to be a promising therapeutic target for acute myeloid leukemia (AML). Herein, we report the discovery of a series of 1-(indolizin-3-yl)ethan-1-one derivatives as potent, and selective CBP bromodomain inhibitors focused on improving cellular potency. One of the most promising compounds, 7e (Y08262), inhibits the CBP bromodomain at the nanomolar level (IC50 = 73.1 nM) with remarkable selectivity. In addition, the new inhibitor also displays potent inhibitory activities in AML cell lines. Collectively, this study provides a new lead compound for further validation of CBP bromodomain as a molecular target for AML drug development.

Small Molecule Ligands of the BET-like Bromodomain, SmBRD3, Affect Schistosoma mansoni Survival, Oviposition, and Development

Schistosomiasis is a disease affecting >200 million people worldwide, but its treatment relies on a single agent, praziquantel. To investigate new avenues for schistosomiasis control, we have conducted the first systematic analysis of bromodomain-containing proteins (BCPs) in a causative species, Schistosoma mansoni. Having identified 29 putative bromodomains (BRDs) in 22 S. mansoni proteins, we selected SmBRD3, a tandem BRD-containing BCP that shows high similarity to the human bromodomain and extra terminal domain (BET) family, for further studies. Screening 697 small molecules identified the human BET BRD inhibitor I-BET726 as a ligand for SmBRD3. An X-ray crystal structure of I-BET726 bound to the second BRD of SmBRD3 [SmBRD3(2)] enabled rational design of a quinoline-based ligand (15) with an ITC Kd = 364 ± 26.3 nM for SmBRD3(2). The ethyl ester pro-drug of compound 15 (compound 22) shows substantial effects on sexually immature larval schistosomula, sexually mature adult worms, and snail-infective miracidia in ex vivo assays.

Divergent Synthesis of Trifluoromethyl-Substituted 1,2-Dihydroquinoxalines and Diimines by Cascade Reactions of CF3-Imidoyl Sulfoxonium Ylides with Azo Compounds

A base-mediated cascade reaction of CF3-imidoyl sulfoxonium ylides and azo compounds has been achieved, allowing for facile access to trifluoromethyl-substituted 1,2-dihydroquinoxalines and diimines in moderate to excellent yields. Noteworthy is that the unusual N-N bond cleavage and rearrangement of azo compounds are involved in the transformations.

Recent Advances on Small-Molecule Bromodomain-Containing Histone Acetyltransferase Inhibitors

In recent years, substantial research has been conducted on molecular mechanisms and inhibitors targeting bromodomains (BRDs) and extra-terminal (BET) family proteins. On this basis, non-BET BRD is gradually becoming a research hot spot. BRDs are abundant in histone acetyltransferase (HAT)-associated activating transcription factors, and BRD-containing HATs have been linked to cancer, inflammation, and viral replication. Therefore, the development of BRD-containing HATs as chemical probes is useful for understanding the specific biological roles of BRDs in diseases and drug discovery. Several types of BRD-containing HATs, including CBP/P300, PCAF/GCN5, and TAF1, are discussed in this context in terms of their structures, functions, and small-molecule inhibitors. Additionally, progress in BRD inhibitors/chemical probes and proteolysis targeting chimeras in terms of drug design, biological activity, and disease application are summarized. These findings provide insights into the development of BRD inhibitors as potential drug candidates for various diseases.

Structural insights revealed by the cocrystal structure of CCS1477 in complex with CBP bromodomain

Inhibition of the bromodomain of the CREB (cyclic-AMP response element-binding protein) binding protein (CBP) is a particularly promising new therapeutic approach for cancer. Benzimidazole derivatives CCS1477 and its analogues (8 and 9) are highly potent and selective CBP bromodomain inhibitors, with K_d values of 26.4, 37.0, and 34.3 nM in ITC assay, respectively. Among these compounds, CCS1477 is undergoing phase Ib/IIa clinical trials for the treatment of various cancers. Thus, we determined the co-crystal structures of CCS1477 and its analogues in complex with CBP bromodomain and revealed the detailed binding modes. Furthermore, overlapping with other reported co-crystal structures allowed us to identify that interaction with Arg1173, LPF shelf, and ZA channel was critical for keeping strong biological activity and selectivity. Collectively, this study provided a structural basis for CBP bromodomain inhibitors design.

A formal [4 + 2] annulation of diamines and prop-2-ynyl sulfonium salts for the synthesis of tetrahydroquinoxalines

A formal [4 + 2] annulation of diamines and prop-2-ynyl sulfonium salts was developed.

Design, Synthesis, and Biological Evaluation of 1-(Indolizin-3-yl)ethan-1-ones as CBP Bromodomain Inhibitors for the Treatment of Prostate Cancer

CREB (cyclic-AMP responsive element binding protein) binding protein (CBP) is a potential target for prostate cancer treatment. Herein, we report the structural optimization of a series of 1-(indolizin-3-yl)ethan-1-one compounds as new selective CBP bromodomain inhibitors, aiming to improve cellular potency and metabolic stability. This process led to compound 9g (Y08284), which possesses good liver microsomal stability and pharmacokinetic properties (F = 25.9%). Furthermore, the compound is able to inhibit CBP bromodomain as well as the proliferation, colony formation, and migration of prostate cancer cells. Additionally, the new inhibitor shows promising antitumor efficacy in a 22Rv1 xenograft model (TGI = 88%). This study provides new lead compounds for further development of drugs for the treatment of prostate cancer.

α-Keto Acids as Triggers and Partners for the Synthesis of Quinazolinones, Quinoxalinones, Benzooxazinones, and Benzothiazoles in Water

A general and efficient method for the synthesis of quinazolinones, quinoxalinones, benzooxazinones, and benzothiazoles from the reactions of α-keto acids with 2-aminobenzamides, benzene-1,2-diamines, 2-aminophenols, and 2-aminobenzenethiols, respectively, is described. The reactions were conducted under catalyst-free conditions, using water as the sole solvent with no additive required, and successfully applied to the synthesis of sildenafil. More importantly, these reactions can be conducted on a mass scale, and the products can be easily purified through filtration and washing with ethanol (or crystallized).

Catalytic Enantioselective Access to Dihydroquinoxalinones via Formal α-Halo Acyl Halide Synthon in One Pot

An enantioselective one-pot catalytic strategy to dihydroquinoxalinones, featuring novel 1-phenylsulfonyl-1-cyano enantioenriched epoxides as masked α-halo acyl halide synthons, followed by a domino ring-opening cyclization (DROC), is documented. A popular quinine-derived urea served as the catalyst in two out of the three steps performed in the same solvent using commercially available aldehydes, (phenylsulfonyl)acetonitrile, cumyl hydroperoxide and 1,2-phenylendiamines. Medicinally relevant 3-aryl/alkyl-substituted heterocycles are isolated in generally good to high overall yield and high enantioselectivity (up to 99 % ee). A rare example of excellent reusability of an organocatalyst at higher scale, subjected to oxidative conditions, is demonstrated. Mechanistically, labile α-ketosulfone has been detected as the intermediate involved in the DROC process. Theoretical calculations on the key epoxidation step rationalize the observed stereocontrol, highlighting the important role played by the sulfone group.

Controlling Intramolecular Interactions in the Design of Selective, High-Affinity Ligands for the CREBBP Bromodomain

CREBBP (CBP/KAT3A) and its paralogue EP300 (KAT3B) are lysine acetyltransferases (KATs) that are essential for human development. They each comprise 10 domains through which they interact with >400 proteins, making them important transcriptional co-activators and key nodes in the human protein–protein interactome. The bromodomains of CREBBP and EP300 enable the binding of acetylated lysine residues from histones and a number of other important proteins, including p53, p73, E2F, and GATA1. Here, we report a work to develop a high-affinity, small-molecule ligand for the CREBBP and EP300 bromodomains [(−)-OXFBD05] that shows >100-fold selectivity over a representative member of the BET bromodomains, BRD4(1). Cellular studies using this ligand demonstrate that the inhibition of the CREBBP/EP300 bromodomain in HCT116 colon cancer cells results in lowered levels of c-Myc and a reduction in H3K18 and H3K27 acetylation. In hypoxia (<0.1% O₂), the inhibition of the CREBBP/EP300 bromodomain results in the enhanced stabilization of HIF-1α.

Bioactivation of Isoxazole-Containing Bromodomain and Extra-Terminal Domain (BET) Inhibitors

The 3,5-dimethylisoxazole motif has become a useful and popular acetyl-lysine mimic employed in isoxazole-containing bromodomain and extra-terminal (BET) inhibitors but may introduce the potential for bioactivations into toxic reactive metabolites. As a test, we coupled deep neural models for quinone formation, metabolite structures, and biomolecule reactivity to predict bioactivation pathways for 32 BET inhibitors and validate the bioactivation of select inhibitors experimentally. Based on model predictions, inhibitors were more likely to undergo bioactivation than reported non-bioactivated molecules containing isoxazoles. The model outputs varied with substituents indicating the ability to scale their impact on bioactivation. We selected OXFBD02, OXFBD04, and I-BET151 for more in-depth analysis. OXFBD’s bioactivations were evenly split between traditional quinones and novel extended quinone-methides involving the isoxazole yet strongly favored the latter quinones. Subsequent experimental studies confirmed the formation of both types of quinones for OXFBD molecules, yet traditional quinones were the dominant reactive metabolites. Modeled I-BET151 bioactivations led to extended quinone-methides, which were not verified experimentally. The differences in observed and predicted bioactivations reflected the need to improve overall bioactivation scaling. Nevertheless, our coupled modeling approach predicted BET inhibitor bioactivations including novel extended quinone methides, and we experimentally verified those pathways highlighting potential concerns for toxicity in the development of these new drug leads.

CBP/p300: Critical Co-Activators for Nuclear Steroid Hormone Receptors and Emerging Therapeutic Targets in Prostate and Breast Cancers

The CREB-binding protein (CBP) and p300 are two paralogous lysine acetyltransferases (KATs) that were discovered in the 1980s-1990s. Since their discovery, CBP/p300 have emerged as important regulatory proteins due to their ability to acetylate histone and non-histone proteins to modulate transcription. Work in the last 20 years has firmly established CBP/p300 as critical regulators for nuclear hormone signaling pathways, which drive tumor growth in several cancer types. Indeed, CBP/p300 are critical co-activators for the androgen receptor (AR) and estrogen receptor (ER) signaling in prostate and breast cancer, respectively. The AR and ER are stimulated by sex hormones and function as transcription factors to regulate genes involved in cell cycle progression, metabolism, and other cellular functions that contribute to oncogenesis. Recent structural studies of the AR/p300 and ER/p300 complexes have provided critical insights into the mechanism by which p300 interacts with and activates AR- and ER-mediated transcription. Breast and prostate cancer rank the first and forth respectively in cancer diagnoses worldwide and effective treatments are urgently needed. Recent efforts have identified specific and potent CBP/p300 inhibitors that target the acetyltransferase activity and the acetytllysine-binding bromodomain (BD) of CBP/p300. These compounds inhibit AR signaling and tumor growth in prostate cancer. CBP/p300 inhibitors may also be applicable for treating breast and other hormone-dependent cancers. Here we provide an in-depth account of the critical roles of CBP/p300 in regulating the AR and ER signaling pathways and discuss the potential of CBP/p300 inhibitors for treating prostate and breast cancer.

Stereo- and regiodefined DNA-encoded chemical libraries enable efficient tumour-targeting applications

The encoding of chemical compounds with amplifiable DNA tags facilitates the discovery of small-molecule ligands for proteins. To investigate the impact of stereo- and regiochemistry on ligand discovery, we synthesized a DNA-encoded library of 670,752 derivatives based on 2-azido-3-iodophenylpropionic acids. The library was selected against multiple proteins and yielded specific ligands. The selection fingerprints obtained for a set of protein targets of pharmaceutical relevance clearly showed the preferential enrichment of ortho-, meta- or para-regioisomers, which was experimentally verified by affinity measurements in the absence of DNA. The discovered ligands included novel selective enzyme inhibitors and binders to tumour-associated antigens, which enabled conditional chimeric antigen receptor T-cell activation and tumour targeting.

Development of Dimethylisoxazole-Attached Imidazo[1,2-a]pyridines as Potent and Selective CBP/P300 Inhibitors

The use of epigenetic bromodomain inhibitors as anticancer therapeutics has transitioned from targeting bromodomain extraterminal domain (BET) proteins into targeting non-BET bromodomains. The two most relevant non-BET bromodomain oncology targets are cyclic AMP response element-binding protein (CBP) and E1A binding protein P300 (EP300). To explore the growing CBP/EP300 interest, we developed a highly efficient two-step synthetic route for dimethylisoxazole-attached imidazo[1,2-a]pyridine scaffold-containing inhibitors. Our efficient two-step reactions enabled high-throughput synthesis of compounds designed by molecular modeling, which together with structure-activity relationship (SAR) studies facilitated an overarching understanding of selective targeting of CBP/EP300 over non-BET bromodomains. This led to the identification of a new potent and selective CBP/EP300 bromodomain inhibitor, UMB298 (compound 23, CBP IC₅₀ 72 nM and bromodomain 4, BRD4 IC₅₀ 5193 nM). The SAR we established is in good agreement with literature-reported CBP inhibitors, such as CBP30, and demonstrates the advantage of utilizing our two-step approach for inhibitor development of other bromodomains.

A Single-Stranded DNA-Encoded Chemical Library Based on a Stereoisomeric Scaffold Enables Ligand Discovery by Modular Assembly of Building Blocks

A versatile and Lipinski-compliant DNA-encoded library (DEL), comprising 366 600 glutamic acid derivatives coupled to oligonucleotides serving as amplifiable identification barcodes is designed, constructed, and characterized. The GB-DEL library, constructed in single-stranded DNA format, allows de novo identification of specific binders against several pharmaceutically relevant proteins. Moreover, hybridization of the single-stranded DEL with a set of known protein ligands of low to medium affinity coupled to a complementary DNA strand results in self-assembled selectable chemical structures, leading to the identification of affinity-matured compounds.

Current development of CBP/p300 inhibitors in the last decade

CBP/p300, functioning as histone acetyltransferases and transcriptional co-factors, represents an attractive target for various diseases, including malignant tumor. The development of small-molecule inhibitors targeting the bromodomain and HAT domains of CBP/p300 has aroused broad interests of medicinal chemist in expectation of providing new hope for anti-cancer treatment. In particular, the CBP/p300 bromodomain inhibitor CCS1477, identified by CellCentric, is currently undergone clinical evaluation for the treatment of haematological malignancies and prostate cancer. In this review, we depict the development of CBP/p300 inhibitors reported from 2010 to 2020 and particularly highlight their structure-activity relationships (SARs), binding modes, selectivity and pharmacological functions with the aim to facilitate rational design and development of CBP/p300 inhibitors.

Discovery of selective inhibitors for cyclic AMP response element-binding protein: a combined ligand and structure-based resources pipeline

Bromodomains are epigenetic readers of acetyl-lysine involved in chromatin remodeling and transcriptional regulations. Over the past few years, extensive research has been carried out to discover small-molecule inhibitors against bromodomains to treat various diseases. Cyclic AMP response element-binding protein (CREBBP) bromodomain has emerged as a hot target for cancer therapy. This study aims at discovering new inhibitors against CREBBP bromodomain using ligand-based molecular docking. A library of 2168 lead-like compounds were docked into the Kac binding site of CREBBP bromodomain. On the basis of the energy score and interaction analysis, six compounds were selected. In order to validate the stability of these six protein-ligand complexes 20 ns molecular dynamics simulations and principal component analyses were carried out. Based on the different analyses these six compounds may provide valuable information for developing CREBBP selective inhibitors.

Selective Fragments for the CREBBP Bromodomain Identified from an Encoded Self-assembly Chemical Library

DNA-encoded chemical libraries (DECLs) are collections of chemical moieties individually coupled to distinctive DNA barcodes. Compounds can be displayed either at the end of a single DNA strand (i. e., single-pharmacophore libraries) or at the extremities of two complementary DNA strands (i. e., dual-pharmacophore libraries). In this work, we describe the use of a dual-pharmacophore encoded self-assembly chemical (ESAC) library for the affinity maturation of a known 4,5-dihydrobenzodiazepinone ring (THBD) acetyl-lysine (KAc) mimic for the cyclic-AMP response element binding protein (CREB) binding protein (CREBBP or CBP) bromodomain. The new pair of fragments discovered from library selection showed a sub-micromolar affinity for the CREBBP bromodomain in fluorescence polarization and ELISA assays, and selectivity against BRD4(1).

Structure-Based Discovery and Development of a Series of Potent and Selective Bromodomain and Extra-Terminal Protein Inhibitors

BRD4 has recently emerged as a promising drug target. Therefore, identifying novel inhibitors with distinct properties could enrich their use in anticancer treatment. Guided by the cocrystal structure of hit compound 4 harboring a five-membered-ring linker motif, we quickly identified lead compound 7, which exhibited good antitumor effects in an MM.1S xenograft model by oral administration. Encouraged by its high potency and interesting scaffold, we performed further lead optimization to generate a novel potent series of bromodomain and extra-terminal (BET) inhibitors with a (1,2,4-triazol-5-yl)-3,4-dihydroquinoxalin-2(1H)-one structure. Among them, compound 19 was found to have the best balance of activity, stability, and antitumor efficacy. After confirming its low brain penetration, we conducted comprehensive preclinical studies, including a multiple-species pharmacokinetics profile, extensive cellular mechanism studies, hERG assay, and in vivo antitumor growth effect testing, and we found that compound 19 is a potential BET protein drug candidate for the treatment of cancer.

Underestimated Noncovalent Interactions in Protein Data Bank

Noncovalent interactions (NCIs) play essential roles in the structure and function of biomacromolecules. There are various NCIs, e.g., hydrogen bonds (HBs), cation-π and π-π interactions, and ionic bonds, among which HBs are the most widespread and well-studied. By utilizing the ratio of the observed HBs over pseudo HBs (1.0 Å longer than the HB distance criteria without angle constraints), we demonstrated that HBs in both protein-ligand and protein-protein interfaces are overlooked in structures deposited in PDB. After the QM/MM optimization of 12 protein-ligand complexes, we showed that the overlooked HBs could be recovered. With a systematic search in the PDB, we found that the HB number per residue (N_HB/R) in proteins decreases as structural resolution becomes lower, implying that HBs are overlooked even today, regardless of the type of refinement approach used. Similarly, cation-π, π-π, and ionic interactions were found to be significantly lost, manifesting the universal underestimation of various NCIs. Considering the vital role of NCIs, it is important to recover the NCIs to facilitate drug design, to explore protein-protein interaction, and to study protein structure and function.

Carbene-Catalyzed α-Carbon Amination of Chloroaldehydes for Enantioselective Access to Dihydroquinoxaline Derivatives

An NHC-catalyzed α-carbon amination of chloroaldehydes was developed. Cyclohexadiene-1,2-diimines are used as amination reagents and four-atom synthons. Our reaction affords optically enriched dihydroquinoxalines that are core structures in natural products and synthetic bioactive molecules.

Rewiring the Epigenetic Networks in MLL-Rearranged Leukemias: Epigenetic Dysregulation and Pharmacological Interventions

Leukemias driven by chromosomal translocation of the mixed-lineage leukemia gene (MLL or KMT2A) are highly prevalent in pediatric oncology. The poor survival rate and lack of an effective targeted therapy for patients with MLL-rearranged (MLL-r) leukemias emphasize an urgent need for improved knowledge and novel therapeutic approaches for these malignancies. The resulting chimeric products of MLL gene rearrangements, i.e., MLL-fusion proteins (MLL-FPs), are capable of transforming hematopoietic stem/progenitor cells (HSPCs) into leukemic blasts. The ability of MLL-FPs to reprogram HSPCs toward leukemia requires the involvement of multiple chromatin effectors, including the histone 3 lysine 79 methyltransferase DOT1L, the chromatin epigenetic reader BRD4, and the super elongation complex. These epigenetic regulators constitute a complicated network that dictates maintenance of the leukemia program, and therefore represent an important cluster of therapeutic opportunities. In this review, we will discuss the role of MLL and its fusion partners in normal HSPCs and hematopoiesis, including the links between chromatin effectors, epigenetic landscapes, and leukemia development, and summarize current approaches to therapeutic targeting of MLL-r leukemias.

qFit-ligand Reveals Widespread Conformational Heterogeneity of Drug-Like Molecules in X-Ray Electron Density Maps

Proteins and ligands sample a conformational ensemble that governs molecular recognition, activity, and dissociation. In structure-based drug design, access to this conformational ensemble is critical to understand the balance between entropy and enthalpy in lead optimization. However, ligand conformational heterogeneity is currently severely underreported in crystal structures in the Protein Data Bank, owing in part to a lack of automated and unbiased procedures to model an ensemble of protein-ligand states into X-ray data. Here, we designed a computational method, qFit-ligand, to automatically resolve conformationally averaged ligand heterogeneity in crystal structures, and applied it to a large set of protein receptor-ligand complexes. In an analysis of the cancer related BRD4 domain, we found that up to 29% of protein crystal structures bound with drug-like molecules present evidence of unmodeled, averaged, relatively isoenergetic conformations in ligand-receptor interactions. In many retrospective cases, these alternate conformations were adventitiously exploited to guide compound design, resulting in improved potency or selectivity. Combining qFit-ligand with high-throughput screening or multitemperature crystallography could therefore augment the structure-based drug design toolbox.

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.

BET bromodomain ligands: Probing the WPF shelf to improve BRD4 bromodomain affinity and metabolic stability

Ligands for the bromodomain and extra-terminal domain (BET) family of bromodomains have shown promise as useful therapeutic agents for treating a range of cancers and inflammation. Here we report that our previously developed 3,5-dimethylisoxazole-based BET bromodomain ligand (OXFBD02) inhibits interactions of BRD4(1) with the RelA subunit of NF-κB, in addition to histone H4. This ligand shows a promising profile in a screen of the NCI-60 panel but was rapidly metabolised (t½ = 39.8 min). Structure-guided optimisation of compound properties led to the development of the 3-pyridyl-derived OXFBD04. Molecular dynamics simulations assisted our understanding of the role played by an internal hydrogen bond in altering the affinity of this series of molecules for BRD4(1). OXFBD04 shows improved BRD4(1) affinity (IC50 = 166 nM), optimised physicochemical properties (LE = 0.43; LLE = 5.74; SFI = 5.96), and greater metabolic stability (t½ = 388 min).

Computer-Aided Drug Design in Epigenetics

Epigenetic dysfunction has been widely implicated in several diseases especially cancers thus highlights the therapeutic potential for chemical interventions in this field. With rapid development of computational methodologies and high-performance computational resources, computer-aided drug design has emerged as a promising strategy to speed up epigenetic drug discovery. Herein, we make a brief overview of major computational methods reported in the literature including druggability prediction, virtual screening, homology modeling, scaffold hopping, pharmacophore modeling, molecular dynamics simulations, quantum chemistry calculation, and 3D quantitative structure activity relationship that have been successfully applied in the design and discovery of epi-drugs and epi-probes. Finally, we discuss about major limitations of current virtual drug design strategies in epigenetics drug discovery and future directions in this field.

In silico quest of selective naphthyl-based CREBBP bromodomain inhibitor

The reader proteins like bromodomains have recently gained increased attentions in the area of epigenetic drug discovery, as they are the potent regulators in gene transcription process. Among the other bromodomains, cAMP response element-binding protein (CREB) binding protein or CREBBP bomodomain involved in various cancer progressions and therefore, efforts to develop specific inhibitors of CREBBP bomodomain are of clinical value. In this study, we tried to identify selective CREBBP bromodomain inhibitor, which was accomplished by using molecular docking, free energy calculation and molecular dynamics (MD) simulation studies, considering a series of naphthyl based compounds. The docking procedure was validated by comparing root mean square deviations (RMSDs) of crystallographic complex to docked complex. Favorable electrostatic interactions with the Arg1173 side chain were considered to attain selectivity for CREBBP bromodomain over other human bromodomain subfamilies. We found that naphthyl-based compounds have greater binding affinities towards the CREBBP bromodomain, and formed non-bonded interactions with various side chain residues that are important for bromodomain inhibition. From detailed investigation by induced fit docking, compound 31 was found to have favorable electrostatic interactions with the Arg1173 side chain by forming conventional hydrogen bonds. This result was further confirmed by analyzing hydrogen bond occupancy and bonding distance during the molecular dynamics simulation. We believe that these findings offer useful insight for the designing of target specific new bromodomain inhibitor and also promote further structure guided synthesis of analogues for identification of potent CREBBP bromodomain inhibitors as well as detailed in vitro and in vivo analyses.

Chemical Space Expansion of Bromodomain Ligands Guided by in Silico Virtual Couplings (AutoCouple)

Expanding the chemical space and simultaneously ensuring synthetic accessibility is of upmost importance, not only for the discovery of effective binders for novel protein classes but, more importantly, for the development of compounds against hard-to-drug proteins. Here, we present AutoCouple, a de novo approach to computational ligand design focused on the diversity-oriented generation of chemical entities via virtual couplings. In a benchmark application, chemically diverse compounds with low-nanomolar potency for the CBP bromodomain and high selectivity against the BRD4(1) bromodomain were achieved by the synthesis of about 50 derivatives of the original fragment. The binding mode was confirmed by X-ray crystallography, target engagement in cells was demonstrated, and antiproliferative activity was showcased in three cancer cell lines. These results reveal AutoCouple as a useful in silico coupling method to expand the chemical space in hit optimization campaigns resulting in potent, selective, and cell permeable bromodomain ligands.

Structure-based optimization of a series of selective BET inhibitors containing aniline or indoline groups

Recently, several kinase inhibitors were found to also inhibit bromodomains, providing a new strategy for the discovery of bromodomain inhibitors. Along this line, starting from PLK1-BRD4 dual inhibitor BI-2536, we discovered a new series of dihydroquinoxalin-2(1H)-one with aniline and indoline WPF binders as selective BRD4 inhibitors. They showed better BRD4-BD1 potency and negligible PLK1 kinase activity comparing with BI-2536. Additionally, dihydroquinoxalin-2(1H)-ones containing indoline group showed profound activities in molecular and cellular based assays. Throughout the study, compounds 9, 28 and 37 showed significant inhibitory activity for c-Myc or PD-L1 protein expression and mRNA transcription both at concentration of 0.2 and 1 μM. Compound 9 was found possessing the best balance of binding affinity, in vitro metabolic stability and in vivo pharmacokinetic properties. Therefore, it was selected for in vivo pharmacological study. By using MM.1S cell derived xenograft model, we confirmed compound 9 showed comparable in vivo tumor inhibition to phase II investigation drug I-BET762, which, together with the novel WPF binder, further indicated the utility of this series of BRD4 inhibitors.

Discovery and optimization of 1-(1H-indol-1-yl)ethanone derivatives as CBP/EP300 bromodomain inhibitors for the treatment of castration-resistant prostate cancer

The CREB (cAMP responsive element binding protein) binding protein (CBP) and its homolog EP300 have emerged as new therapeutic targets for the treatment of cancer and inflammatory diseases. Here we report the identification, optimization and evaluation of 1-(1H-indol-1-yl)ethanone derivatives as CBP/EP300 inhibitors starting from fragment-based virtual screening (FBVS). A cocrystal structure of the inhibitor (22e) in complex with CBP provides a solid structural basis for further optimization. The most potent compound 32h binds to the CBP bromodomain and has an IC₅₀ value of 0.037 μM in the AlphaScreen assay which was 2 times more potent than the reported CBP bromodomain inhibitor SGC-CBP30 in our hands. 32h also exhibit high selectivity for CBP/EP300 over other bromodomain-containing proteins. Notably, the ester derivative (29h) of compound 32h markedly inhibits cell growth in several prostate cancer cell lines including LNCaP, 22Rv1 and LNCaP derived C4-2B. Compound 29h suppresses the mRNA expression of full length AR (AR-FL), AR target genes and other oncogene in LNCaP cells. 29h also reduces the expression of PSA, the biomarker of prostate cancer. CBP/EP300 inhibitor 29h represents a promising lead compound for the development of new therapeutics for the treatment of castration-resistant prostate cancer.

Cation-π interactions in protein-ligand binding: theory and data-mining reveal different roles for lysine and arginine

The interactions of neutral aromatic ligands with cationic arginine, histidine and lysine amino acid residues have been studied with ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of X-ray structures.

We have studied the cation–π interactions of neutral aromatic ligands with the cationic amino acid residues arginine, histidine and lysine using ab initio calculations, symmetry adapted perturbation theory (SAPT), and a systematic meta-analysis of all available Protein Data Bank (PDB) X-ray structures. Quantum chemical potential energy surfaces (PES) for these interactions were obtained at the DLPNO-CCSD(T) level of theory and compared against the empirical distribution of 2012 unique protein–ligand cation–π interactions found in X-ray crystal structures. We created a workflow to extract these structures from the PDB, filtering by interaction type and residue pK_a. The gas phase cation–π interaction of lysine is the strongest by more than 10 kcal mol^–1, but the empirical distribution of 582 X-ray structures lies away from the minimum on the interaction PES. In contrast, 1381 structures involving arginine match the underlying calculated PES with good agreement. SAPT analysis revealed that underlying differences in the balance of electrostatic and dispersion contributions are responsible for this behavior in the context of the protein environment. The lysine–arene interaction, dominated by electrostatics, is greatly weakened by a surrounding dielectric medium and causes it to become essentially negligible in strength and without a well-defined equilibrium separation. The arginine–arene interaction involves a near equal mix of dispersion and electrostatic attraction, which is weakened to a much smaller degree by the surrounding medium. Our results account for the paucity of cation–π interactions involving lysine, even though this is a more common residue than arginine. Aromatic ligands are most likely to interact with cationic arginine residues as this interaction is stronger than for lysine in higher polarity surroundings.

Classical molecular dynamics and metadynamics simulations decipher the mechanism of CBP30 selectively inhibiting CBP/p300 bromodomains

The selective inhibition mechanism of CBP30 towards CBP/p300 over BRD4-BD1/BD2 bromodomains was revealed by conventional molecular dynamics and metadynamics simulations.

A Unique Approach to Design Potent and Selective Cyclic Adenosine Monophosphate Response Element Binding Protein, Binding Protein (CBP) Inhibitors

The epigenetic regulator CBP/P300 presents a novel therapeutic target for oncology. Previously, we disclosed the development of potent and selective CBP bromodomain inhibitors by first identifying pharmacophores that bind the KAc region and then building into the LPF shelf. Herein, we report the "hybridization" of a variety of KAc-binding fragments with a tetrahydroquinoline scaffold that makes optimal interactions with the LPF shelf, imparting enhanced potency and selectivity to the hybridized ligand. To demonstrate the utility of our hybridization approach, two analogues containing unique Asn binders and the optimized tetrahydroquinoline moiety were rapidly optimized to yield single-digit nanomolar inhibitors of CBP with exquisite selectivity over BRD4(1) and the broader bromodomain family.

GNE-781, A Highly Advanced Potent and Selective Bromodomain Inhibitor of Cyclic Adenosine Monophosphate Response Element Binding Protein, Binding Protein (CBP)

Inhibition of the bromodomain of the transcriptional regulator CBP/P300 is an especially interesting new therapeutic approach in oncology. We recently disclosed in vivo chemical tool 1 (GNE-272) for the bromodomain of CBP that was moderately potent and selective over BRD4(1). In pursuit of a more potent and selective CBP inhibitor, we used structure-based design. Constraining the aniline of 1 into a tetrahydroquinoline motif maintained potency and increased selectivity 2-fold. Structure-activity relationship studies coupled with further structure-based design targeting the LPF shelf, BC loop, and KAc regions allowed us to significantly increase potency and selectivity, resulting in the identification of non-CNS penetrant 19 (GNE-781, TR-FRET IC₅₀ = 0.94 nM, BRET IC₅₀ = 6.2 nM; BRD4(1) IC₅₀ = 5100 nΜ) that maintained good in vivo PK properties in multiple species. Compound 19 displays antitumor activity in an AML tumor model and was also shown to decrease Foxp3 transcript levels in a dose dependent manner.