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

Structure-Guided Design and Synthesis of Pyridinone-Based Selective Bromodomain and Extra-Terminal Domain (BET)-First Bromodomain (BD1) Inhibitors

The bromodomain and extra-terminal domain (BET) proteins are epigenetic readers, regulating transcription via two highly homologous tandem bromodomains, BD1 and BD2. Clinical development of nonselective pan-BD BET inhibitors has been challenging, partly due to dose-limiting side effects such as thrombocytopenia. This has prompted the push for domain-selective BET inhibitors to achieve a more favorable therapeutic window. We report a structure-guided drug design campaign that led to the development of a potent BD1-selective BET inhibitor, 33 (XL-126), with a Kd of 8.9 nM and 185-fold BD1/BD2 selectivity. The high selectivity was first assayed by SPR, validated by a secondary time-resolved fluorescence energy transfer assay, and further corroborated by BROMOscan (∼57-373 fold selectivity). The cocrystal of 33 with BRD4 BD1 and BD2 demonstrates the source of selectivity: repulsion with His437 and lost binding with the leucine clamp. Notably, the BD1 selectivity of BET inhibitor 33 leads to both the preservation of platelets and potent anti-inflammatory efficacy.

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.

Mutate and Conjugate: A Method to Enable Rapid In-Cell Target Validation

Target validation remains a challenge in drug discovery, which leads to a high attrition rate in the drug discovery process, particularly in Phase II clinical trials. Consequently, new approaches to enhance target validation are valuable tools to improve the drug discovery process. Here, we report the combination of site-directed mutagenesis and electrophilic fragments to enable the rapid identification of small molecules that selectively inhibit the mutant protein. Using the bromodomain-containing protein BRD4 as an example, we employed a structure-based approach to identify the L94C mutation in the first bromodomain of BRD4 [BRD4(1)] as having a minimal effect on BRD4(1) function. We then screened a focused, KAc mimic-containing fragment set and a diverse fragment library against the mutant and wild-type proteins and identified a series of fragments that showed high selectivity for the mutant protein. These compounds were elaborated to include an alkyne click tag to enable the attachment of a fluorescent dye. These clickable compounds were then assessed in HEK293T cells, transiently expressing BRD4(1)WT or BRD4(1)L94C, to determine their selectivity for BRD4(1)L94C over other possible cellular targets. One compound was identified that shows very high selectivity for BRD4(1)L94C over all other proteins. This work provides a proof-of-concept that the combination of site-directed mutagenesis and electrophilic fragments, in a mutate and conjugate approach, can enable rapid identification of small molecule inhibitors for an appropriately mutated protein of interest. This technology can be used to assess the cellular phenotype of inhibiting the protein of interest, and the electrophilic ligand provides a starting point for noncovalent ligand development.

Bioactivation and reactivity research advances - 2022 year in review‡

With the 50th year mark since the launch of Drug Metabolism and Disposition journal, the field of drug metabolism and bioactivation has advanced exponentially in the past decades (Guengerich 2023).This has, in a major part, been due to the continued advances across the whole spectrum of applied technologies in hardware, software, machine learning (ML), and artificial intelligence (AI). LC-MS platforms continue to evolve to support key applications in the field, and automation is also improving the accuracy, precision, and throughput of these supporting assays. In addition, sample generation and processing is being aided by increased diversity and quality of reagents and bio-matrices so that what is being analyzed is more relevant and translatable. The application of in silico platforms (applied software, ML, and AI) is also making great strides, and in tandem with the more traditional approaches mentioned previously, is significantly advancing our understanding of bioactivation pathways and how these play a role in toxicity. All of this continues to allow the area of bioactivation to evolve in parallel with associated fields to help bring novel or improved medicines to patients with urgent or unmet needs.Shuai Wang and Cyrus Khojasteh, on behalf of the authors.

Selectivity Mechanism of Pyrrolopyridone Analogues Targeting Bromodomain 2 of Bromodomain-Containing Protein 4 from Molecular Dynamics Simulations

Bromodomain and extra-terminal domain (BET) proteins play an important role in epigenetic regulation and are linked to several diseases; therefore, they are interesting targets. BET has two bromodomains: bromodomain 1 (BD1) and BD2. Selective targeting of BD1 or BD2 may produce different activities and greater effects than pan-BD inhibitors. However, the selective mechanism of the specific core must be studied at the atomic level. This study determined the effectiveness of pyrrolopyridone analogues to selectively inhibit BD2 using a pan-BD inhibitor (ABBV-075) and a selective-BD2 inhibitor (ABBV-744). Molecular dynamics simulations and calculations of binding free energies were used to systematically study the selectivity of BD2 inhibition by the pyrrolopyridone analogues. Overall, the pyrrolopyridone analogue inhibitors targeting BD2 interacted mainly with the following amino acid pairs between bromodomain-containing protein 4 (BRD4)-BD1 and BRD4-BD2 complexes: I146/V439, N140/N433, D144/H437, P82/P375, V87/V380, D88/D381, and Y139/Y432. The pyrrolopyridone analogues targeting BRD4-BD2 were divided into five regions based on selectivity mechanism. These results suggest that the R3 and R5 regions of pyrrolopyridone analogues can be modified to improve the selectivity between BRD4-BD1 and BRD4-BD2. The selectivity of BD2 inhibition by pyrrolopyridone analogues can be used to design novel BD2 inhibitors based on a pyrrolopyridone core.

The role of loop dynamics in the prediction of ligand-protein binding enthalpy

The enthalpic and entropic components of ligand-protein binding free energy reflect the interactions and dynamics between ligand and protein. Despite decades of study, our understanding and hence our ability to predict these individual components remains poor. In recent years, there has been substantial effort and success in the prediction of relative and absolute binding free energies, but the prediction of the enthalpic (and entropic) contributions in biomolecular systems remains challenging. Indeed, it is not even clear what kind of performance in terms of accuracy could currently be obtained for such systems. It is, however, relatively straight-forward to compute the enthalpy of binding. We thus evaluated the performance of absolute enthalpy of binding calculations using molecular dynamics simulation for ten inhibitors against a member of the bromodomain family, BRD4-1, against isothermal titration calorimetry data. Initial calculations, with the AMBER force-field showed good agreement with experiment (R² = 0.60) and surprisingly good accuracy with an average of root-mean-square error (RMSE) = 2.49 kcal mol^-1. Of the ten predictions, three were obvious outliers that were all over-predicted compared to experiment. Analysis of various simulation factors, including parameterization, buffer concentration and conformational dynamics, revealed that the behaviour of a loop (the ZA loop on the periphery of the binding site) strongly dictates the enthalpic prediction. Consistent with previous observations, the loop exists in two distinct conformational states and by considering one or the other or both states, the prediction for the three outliers can be improved dramatically to the point where the R² = 0.95 and the accuracy in terms of RMSE improves to 0.90 kcal mol^-1. However, performance across force-fields is not consistent: if OPLS and CHARMM are used, different outliers are observed and the correlation with the ZA loop behaviour is not recapitulated, likely reflecting parameterization as a confounding problem. The results provide a benchmark standard for future study and comparison.

From Hit Seeking to Magic Bullets: The Successful Union of Epigenetic and Fragment Based Drug Discovery (EPIDD + FBDD)

We review progress in the application of fragment-based drug discovery (FBDD) to epigenetic drug discovery (EPIDD) targeted at epigenetic writer and eraser enzymes as well as reader domains over the last 15 years. The greatest successes to date are in prospecting for bromodomain binding ligands. From a diverse array of fragment hits, multiple potent and selective compounds ensued, including the oncology clinical candidates mivebresib, ABBV-744, pelabresib, and PLX51107.

Discovery of a Novel Bromodomain and Extra Terminal Domain (BET) Protein Inhibitor, I-BET282E, Suitable for Clinical Progression

The functions of the bromodomain and extra terminal (BET) family of proteins have been implicated in a wide range of diseases, particularly in the oncology and immuno-inflammatory areas, and several inhibitors are under investigation in the clinic. To mitigate the risk of attrition of these compounds due to structurally related toxicity findings, additional molecules from distinct chemical series were required. Here we describe the structure- and property-based optimization of the in vivo tool molecule I-BET151 toward I-BET282E, a molecule with properties suitable for progression into clinical studies.

The dioxasilepanyl group as a versatile organometallic unit: studies on stability, reactivity, and utility

Organic synthesis is performed based on precise choices of functional groups and reactions employed. In a multistep synthesis, an ideal functional group should be compatible with various reaction conditions and unaltered until it is subjected to a selective conversion. The current study was set out to search for a silicon functionality that meets these criteria. Here we have established a new silicon-based synthetic methodology centred on a bulky 7-membered dialkoxysilyl group (2,4,4,7,7-pentamethyl-1,3,2-dioxasilepan-2-yl) that uniquely has both stability and on-demand reactivity. The exceptional stability of this functional group was corroborated by both experimental and computational studies which demonstrated that key factors for its stability were a 7-membered structure and steric hindrance. In turn, the dioxasilepanyl group was found to become reactive and to be easily transformed in the presence of appropriate activators. Combined with the development of easy and robust methods to introduce the dioxasilepanyl group onto aryl rings, these findings have allowed a shorter and more efficient synthesis of a bioactive molecule, thus demonstrating the potential utility of the easily accessible dioxasilepanyl group in organic synthesis.

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.

Rigorous Computational Study Reveals What Docking Overlooks: Double Trouble from Membrane Association in Protein Kinase C Modulators

Increasing protein kinase C (PKC) activity is of potential therapeutic value. Its activation involves an interaction between the C1 domain and diacylglycerol (DAG) at intracellular membrane surfaces; DAG mimetics hold promise as new drugs. We previously developed the isophthalate derivative HMI-1a3, an effective but highly lipophilic (clogP = 6.46) DAG mimetic. Although a less lipophilic pyrimidine analog, PYR-1gP (clogP = 3.30), gave positive results in computational docking, it unexpectedly presented greatly diminished binding to PKC in vitro. Through more rigorous computational molecular modeling, we reveal that, unlike HMI-1a3, PYR-1gP forms an intramolecular hydrogen bond, which both obstructs binding and reorients PYR-1gP in the membrane in a fashion that prevents it from correctly accessing the PKC C1 domain. Our results highlight the great value of molecular dynamics simulations as a key component for the drug design process of ligands targeting weakly membrane-associated proteins, where simulation in the relevant membrane environment is crucial for obtaining biologically applicable results.

Design, synthesis and biological evaluation of imidazolopyridone derivatives as novel BRD4 inhibitors

Bromodomain containing protein 4 (BRD4) has been demonstrated to play critical roles in cellular proliferation and cell cycle progression. In this study, using the BRD4 inhibitor Fragment 9 as a lead compound, a series of imidazolopyridone derivatives were designed and tested for their inhibitory activity against BRD4 protein in vitro. Among them, HB100-A7 showed excellent BRD4(1) inhibitory activities with an IC₅₀ value of 0.035 μM in amplified luminescent proximity homogeneous assay (Alphascreen). The result of MTT assay showed that HB100-A7 could suppress the proliferation of pancreatic cancer cells. In addition, flow cytometry further illustrated that HB100-A7 treatment resulted in G0/G1 phase arrest and promoted apoptosis of BxPc3 cells. Furthermore, the in vivo study found that HB100-A7 displayed significant tumor growth inhibition in a pancreatic mouse tumor model (Panc-02). Moreover, IHC staining suggested that HB100-A7 induce cell apoptosis in pancreatic cancer tumor tissue. Together, this study revealed, for the first time, HB100-A7 is a promising lead compound for further development as a new generation of small molecule inhibitors targeting the BRD4 protein.

Structure-Activity Relationship in Pyrazolo[4,3-c]pyridines, First Inhibitors of PEX14-PEX5 Protein-Protein Interaction with Trypanocidal Activity

Trypanosoma protists are pathogens leading to a spectrum of devastating infectious diseases. The range of available chemotherapeutics against Trypanosoma is limited, and the existing therapies are partially ineffective and cause serious adverse effects. Formation of the PEX14-PEX5 complex is essential for protein import into the parasites' glycosomes. This transport is critical for parasite metabolism and failure leads to mislocalization of glycosomal enzymes, with fatal consequences for the parasite. Hence, inhibiting the PEX14-PEX5 protein-protein interaction (PPI) is an attractive way to affect multiple metabolic pathways. Herein, we have used structure-guided computational screening and optimization to develop the first line of compounds that inhibit PEX14-PEX5 PPI. The optimization was driven by several X-ray structures, NMR binding data, and molecular dynamics simulations. Importantly, the developed compounds show significant cellular activity against Trypanosoma, including the human pathogen Trypanosoma brucei gambiense and Trypanosoma cruzi parasites.

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.

BET bromodomain inhibitors: fragment-based in silico design using multi-target QSAR models

Epigenetics has become a focus of interest in drug discovery. In this sense, bromodomain-containing proteins have emerged as potential epigenetic targets in cancer research and other therapeutic areas. Several computational approaches have been applied to the prediction of bromodomain inhibitors. Nevertheless, such approaches have several drawbacks such as the fact that they predict activity against only one bromodomain-containing protein, using structurally related compounds. Also, there are no reports focused on meaningfully analyzing the physicochemical/structural features that are necessary for the design of a bromodomain inhibitor. This work describes the development of two different multi-target models based on quantitative structure-activity relationships (mt-QSAR) for the prediction and in silico design of multi-target bromodomain inhibitors against the proteins BRD2, BRD3, and BRD4. The first model relied on linear discriminant analysis (LDA) while the second focused on artificial neural networks. Both models exhibited accuracies higher than 85% in the dataset. Several molecular fragments were extracted, and their contributions to the inhibitory activity against the three BET proteins were calculated by the LDA model. Six molecules were designed by assembling the fragments with positive contributions, and they were predicted as multi-target BET bromodomain inhibitors by the two mt-QSAR models. Molecular docking calculations converged with the predictions performed by the mt-QSAR models, suggesting that the designed molecules can exhibit potent activity against the three BET proteins. These molecules complied with the Lipinski's rule of five.