Binding Kinetics versus Affinities in BRD4 Inhibition

ClickRNA-PROTAC for Tumor-Selective Protein Degradation and Targeted Cancer Therapy

Proteolysis-targeting chimeras (PROTACs) show promise in tumor treatment. However, the E3 ligases VHL and CRBN, commonly used in PROTAC, are highly expressed in only a few tumors, thus limiting the application scope and efficacy of PROTAC drugs. Furthermore, the lack of tumor specificity in PROTAC drugs can result in toxic side effects. Therefore, there is an urgent need to develop tumor-selective PROTAC drugs that do not rely on endogenous E3 ligases. In this study, we introduce the ClickRNA-PROTAC system, which involves the expression of a fusion protein of the E3 ubiquitin ligase SIAH1 and SNAPTag through mRNA transfection and recruits the protein of interest (POI) using bio-orthogonal click chemistry. ClickRNA-PROTAC can effectively degrade various proteins such as BRD4, KRAS, and NFκB simply by replacing the warhead molecules. By employing a tumor-specific mRNA-responsive translation strategy, ClickRNA-PROTAC can selectively degrade POIs in tumor cells. Furthermore, ClickRNA-PROTAC demonstrated strong efficacy in targeted cancer therapy in a xenograft mouse model of adrenocortical carcinoma. In conclusion, this approach offers several advantages, including independence from endogenous E3 ubiquitin ligases, tumor specificity, and programmability, thereby paving the way for the development of PROTAC drugs.

Addressing Suboptimal Poses in Nonequilibrium Alchemical Calculations

Alchemical transformations can be used to quantitatively estimate absolute binding free energies at a reasonable computational cost. However, most of the approaches currently in use require knowledge of the correct (crystallographic) pose. In this paper, we present a combined Hamiltonian replica exchange nonequilibrium alchemical method that allows us to reliably calculate absolute binding free energies, even when starting from suboptimal initial binding poses. Performing a preliminary Hamiltonian replica exchange enhances the sampling of slow degrees of freedom of the ligand and the target, allowing the system to populate the correct binding pose when starting from an approximate docking pose. We apply the method on 6 ligands of the first bromodomain of the BRD4 bromodomain-containing protein. For each ligand, we start nonequilibrium alchemical transformations from both the crystallographic pose and the top-scoring docked pose that are often significantly different. We show that the method produces statistically equivalent binding free energies, making it a useful tool for computational drug discovery pipelines.

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.

Regulation of inflammatory genes in decidual cells: Involvement of the bromodomain and extra-terminal family proteins

The decidua undergoes proinflammatory activation in late pregnancy, promoting labor. Bromodomain and Extra-Terminal (BET) family proteins interact with acetylated histones and may control gene expression in inflammation. Here, we assessed whether BETs are involved in inflammatory gene regulation in human decidual cells. We have treated primary cultures of decidual stromal cells (DSCs) from term pregnancies with endotoxin (LPS) and measured the expression of a panel of pro-and anti-inflammatory genes. BET involvement was assessed using the selective BET inhibitors (+)-JQ1 and I-BET-762 or the negative control compound (-)-JQ1. Histone 3 and -4 acetylation and BETs binding at the target gene promoters were determined to assess whether these processes are involved in the actions of LPS, BETs, and BET inhibitors. LPS increased the expression of the proinflammatory (PTGS2, IL6, CXCL8/IL8, TNF) and the anti-inflammatory (IL10, IDO1) genes of the panel. The constitutively expressed inflammatory genes (PTGS1, PTGES) were unaffected. The BET inhibitors, but not the control compound, reduced the basal and LPS-induced expression of PTGS1, PTGS2, IL6, CXCL8/IL8, IL10, and IDO1. TNF expression was not changed by BET inhibition. The dominant BETs were Bromodomain-containing protein -2 (BRD2) and -4L (BRD4L) in DSCs. LPS increased histone 4 acetylation at the CXCL8/IL8 and TNF promoters and histone 3 and -4 acetylation at the IDO1 promoter, while (+)-JQ1 abrogated histone acetylation at several promoters. Overall, histone acetylation and promoter binding of BETs showed no consistent relationship with gene expression across the gene panel and the treatments. BET proteins, predominantly BRD2 and BRD4L, control critical pro- and anti-inflammatory genes in DSCs. TNF induction exemplifies a BET-independent pathway. Changing histone acetylation at the promoters is not a general obligatory requirement for inflammatory gene expression in response to LPS. BETs likely act at chromatin loci separate from the examined promoters. BET inhibitors may block decidual activation at labor.

Mechanistic Analysis of Chemically Diverse Bromodomain-4 Inhibitors Using Balanced QSAR Analysis and Supported by X-ray Resolved Crystal Structures

Bromodomain-4 (BRD-4) is a key enzyme in post-translational modifications, transcriptional activation, and many other cellular processes. Its inhibitors find their therapeutic usage in cancer, acute heart failure, and inflammation to name a few. In the present study, a dataset of 980 molecules with a significant diversity of structural scaffolds and composition was selected to develop a balanced QSAR model possessing high predictive capability and mechanistic interpretation. The model was built as per the OECD (Organisation for Economic Co-operation and Development) guidelines and fulfills the endorsed threshold values for different validation parameters (R²_tr = 0.76, Q²_LMO = 0.76, and R²_ex = 0.76). The present QSAR analysis identified that anti-BRD-4 activity is associated with structural characters such as the presence of saturated carbocyclic rings, the occurrence of carbon atoms near the center of mass of a molecule, and a specific combination of planer or aromatic nitrogen with ring carbon, donor, and acceptor atoms. The outcomes of the present analysis are also supported by X-ray-resolved crystal structures of compounds with BRD-4. Thus, the QSAR model effectively captured salient as well as unreported hidden pharmacophoric features. Therefore, the present study successfully identified valuable novel pharmacophoric features, which could be beneficial for the future optimization of lead/hit compounds for anti-BRD-4 activity.

BRD4: quantum mechanical protein–ligand binding free energies using the full-protein DFT-based QM-PBSA method

Fully quantum mechanical approaches to calculating protein–ligand free energies of binding have the potential to reduce empiricism and explicitly account for all physical interactions responsible for protein–ligand binding.

Recent progress on cheminformatics approaches to epigenetic drug discovery

The ability of epigenetic markers to affect genome function has enabled transformative changes in drug discovery, especially in cancer and other emerging therapeutic areas. Concordant with the introduction of the term 'epi-informatics', the size of the epigenetically relevant chemical space has grown substantially and so did the number of applications of cheminformatic methods to epigenetics. Recent progress in epi-informatics has improved our understanding of the structure-epigenetic activity relationships and boosted the development of models predicting novel epigenetic agents. Herein, we review the advances in computational approaches to drug discovery of small molecules with epigenetic modulation profiles, summarize the current chemogenomics data available for epigenetic targets, and provide a perspective on the greater utility of biomedical knowledge mining as a means to advance the epigenetic drug discovery.

Current advances on the development of BET inhibitors: insights from computational methods

Epigenetics was coined almost 70 years ago for the description of heritable phenotype without altering DNA sequences. Research on the field has uncovered significant roles of such mechanisms, that account for the biogenesis of several diseases. Further studies have led the way for drug development which targets epi-enzymes, mainly for cancer treatment. Of the numerous epi-targets involved with histone acetylation, bromodomains have captured the spotlight of drug discovery focused on novel therapies. However, due to high sequence identity, the development of potent and selective inhibitors poses a significant challenge. Herein, we discuss recent computational developments on BET inhibitors and other methods that may be applied for drug discovery in general. As a proof-of-concept, we discuss a virtual screening to identify novel BET inhibitors based on coumarin derivatives. From public data, we identified putative structure-activity relationships of coumarin scaffold and propose R-group modifications for BET selectivity. Results showed that the optimization and design of novel coumarins could be further explored.

BRD4 as a therapeutic target for nonfunctioning and growth hormone pituitary adenoma

BACKGROUND

Nonfunctioning pituitary adenoma (NFPA) and growth hormone pituitary adenoma (GHPA) are major subtypes of pituitary adenomas (PAs). The primary treatment is surgical resection. However, radical excision remains challenging, and few effective medical therapies are available. It is urgent to find novel targets for the treatment. Bromodomain-containing protein 4 (BRD4) is an epigenetic regulator that leads to aberrant transcriptional activation of oncogenes. Herein, we investigated the pathological role of BRD4 and evaluated the effectiveness of BRD4 inhibitors in the treatment of NFPA and GHPA.

METHODS

The expression of BRD4 was detected in NFPA, GHPA, and normal pituitary tissues. The efficacies of BRD4 inhibitors were evaluated in GH3 and MMQ cell lines, patient-derived tumor cells, and in vivo mouse xenograft models of PA. Standard western blots, real-time PCR, and flow cytometry experiments were performed to investigate the effect of BRD4 inhibitors on cell cycle progression, apoptosis, and the expression patterns of downstream genes.

RESULTS

Immunohistochemistry studies demonstrated the overexpression of BRD4 in NFPA and GHPA. In vitro and in vivo studies showed that treatment with the BRD4 inhibitor ZBC-260 significantly inhibited the proliferation of PA cells. Further mechanistic studies revealed that ZBC-260 could downregulate the expression of c-Myc, B-cell lymphoma 2 (Bcl2), and related genes, which are vital factors in pituitary tumorigenesis.

CONCLUSION

In this study, we determined the overexpression of BRD4 in NFPA and GHPA and assessed the effects of BRD4 inhibitors on PA cells in vitro and in vivo. Our findings suggest that BRD4 is a promising therapeutic target for NFPA and GHPA.

Probing molecular mechanism of inhibitor bindings to bromodomain-containing protein 4 based on molecular dynamics simulations and principal component analysis

It is well known that bromodomain-containing protein 4 (BRD4) has been thought as a promising target utilized for treating various human diseases, such as inflammatory disorders, malignant tumours, acute myelogenous leukaemia (AML), bone diseases, etc. For this study, molecular dynamics (MD) simulations, binding free energy calculations, and principal component analysis (PCA) were integrated together to uncover binding modes of inhibitors 8P9, 8PU, and 8PX to BRD4(1). The results obtained from binding free energy calculations show that van der Waals interactions act as the main regulator in bindings of inhibitors to BRD4(1). The information stemming from PCA reveals that inhibitor associations extremely affect conformational changes, internal dynamics, and movement patterns of BRD4(1). Residue-based free energy decomposition method was wielded to unveil contributions of independent residues to inhibitor bindings and the data signify that hydrogen bonding interactions and hydrophobic interactions are decisive factors affecting bindings of inhibitors to BRD4(1). Meanwhile, eight residues Trp81, Pro82, Val87, Leu92, Leu94, Cys136, Asn140, and Ile146 are recognized as the common hot interaction spots of three inhibitors with BRD4(1). The results from this work are expected to provide a meaningfully theoretical guidance for design and development of effective inhibitors inhibiting of the activity of BRD4.

In silico design and molecular basis for the selectivity of Olinone toward the first over the second bromodomain of BRD4

Bromodomains (BrDs), a conserved structural module in chromatin-associated proteins, are well known for recognizing ε-N-acetyl lysine residues on histones. One of the most relevant BrDs is BRD4, a tandem BrD containing protein (BrD1 and BrD2) that plays a critical role in numerous diseases including cancer. Growing evidence shows that the two BrDs of BRD4 have different biological functions; hence selective ligands that can be used to study their functions are of great interest. Here, as a follow-up of our previous work, we first provide a detailed characterization study of the in silico rational design of Olinone as part of a series of five tetrahydropyrido indole-based compounds as BRD4 BrD1 inhibitors. Additionally, we investigated the molecular basis for Olinone's selective recognition by BrD1 over BrD2. Molecular dynamics simulations, free energy calculations, and conformational analyses of the apo-BRD4-BrD1|2 and BRD4-BrD1|2/Olinone complexes showed that Olinone's selectivity is facilitated by five key residues: Leu92 in BrD1|385 in BrD2 of ZA loop, Asn140|433, Asp144|His437 and Asp145|Glu438 of BC loop, and Ile146|Val49 of helix C. Furthermore, the difference in hydrogen bonds number and in mobility of the ZA and BC loops of the acetyl-lysine binding site between BRD4 BrD1/Olinone and BrD2/Olinone complexes also contribute to the difference in Olinone's binding affinity and selectivity toward BrD1 over BrD2. Altogether, our computer-aided molecular design techniques can effectively guide the development of small-molecule BRD4 BrD1 inhibitors, explain their selectivity origin, and further open doors to the design of new therapeutically improved derivatives.

Discovery of alkoxy benzamide derivatives as novel BPTF bromodomain inhibitors via structure-based virtual screening

Bromodomain PHD finger transcription factor (BPTF), a bromodomain-containing protein, plays a crucial role in the regulation of downstream gene expression through the specific recognition of lysine acetylation on bulk histones. The dysfunction of BPTF is closely involved with the development and progression of many human diseases, especially cancer. Therefore, BPTF bromodomain has become a promising drug target for epigenetic cancer therapy. However, unlike BET family inhibitors, few BPTF bromodomain inhibitors have been reported. In this study, by integrating docking-based virtual screening with biochemical analysis, we identified a novel selective BPTF bromodomain inhibitor DCB29 with the IC50 value of 13.2 ± 1.6 μM by homogenous time-resolved fluorescence resonance energy transfer (HTRF) assays. The binding between DCB29 and BPTF was confirmed by NMR and SPR. Molecular docking disclosed that DCB29 occupied the pocket of acetylated H4 peptide substrate and provided detailed SAR explanations for its derivatives. Collectively, DCB29 presented great potential as a powerful tool for BPTF-related biological research and further medicinal chemistry optimization.

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.

Selective inhibition mechanism of RVX-208 to the second bromodomain of bromo and extraterminal proteins: insight from microsecond molecular dynamics simulations

RVX-208 is a recently reported inhibitor of bromo and extraterminal (BET) family proteins (including BRD2-4 and BRDT) with selectivity for the second bromodomain (BD2), currently in phase III clinical trials. Despite of its promising antitumor activity, due to the conserved folds of the first and second bromodomains (BD1 and BD2), the detailed selectivity mechanism of RVX-208 towards BD2 over BD1 is still unknown. To elucidate selective inhibition mechanism of RVX-208 to BD2, microsecond molecular dynamics simulations were performed in this study for BRD2-BD1, BRD2-BD2 and BRD4-BD1 with and without RVX-208, respectively. Binding free energy calculations show that there exists strongest interaction between RVX-208 and BRD2-BD2. Leu383 and Asn429 are two most important residues of BRD2-BD2 for binding to RVX-208. Structural network analysis reveals that RVX-208 can shorten the communication path of ZA and BC loops in BRD2-BD2 pocket, making pocket more suitable to accommodate RVX-208. Additionally, different behaviors of His433 (Asp160 in BRD2-BD1) and Val435 (Ile162 in BRD2-BD1) in BRD2-BD2 are key factors responsible for selective binding of RVX-208 to BRD2-BD2. The proposed selective inhibition mechanism of RVX-208 to BRD2-BD2 can be helpful for rational design of novel selective inhibitors of the second bromodomain of BET family proteins.

Bisphosphonate-Functionalized Hydroxyapatite Nanoparticles for the Delivery of the Bromodomain Inhibitor JQ1 in the Treatment of Osteosarcoma

Osteosarcoma (OS) is one of the most common neoplasia among children, and its survival statistics have been stagnating since the combinatorial anticancer therapy triad was first introduced. Here, we report on the assessment of the effect of hydroxyapatite (HAp) nanoparticles loaded with medronate, the simplest bisphosphonate, as a bone-targeting agent and JQ1, a small-molecule bromodomain inhibitor, as a chemotherapeutic in different 2D and 3D K7M2 OS in vitro models. Both additives decreased the crystallinity of HAp, but the effect was more intense for medronate because of its higher affinity for HAp. As the result of PO₄^3--NH⁺ binding, JQ1 shielded the surface phosphates of HAp and pushed its surface charge to more positive values, exhibiting the opposite effect from calcium-blocking medronate. In contrast to the faster and more exponential release of JQ1 from monetite, its release from HAp nanoparticles followed a zero-order kinetics, but 98% of the payload was released after 48 h. The apoptotic effect of HAp nanoparticles loaded with JQ1, with medronate and with both JQ1 and medronate, was selective in 2D culture: pronounced against the OS cells and nonexistent against the healthy fibroblasts. While OS cell invasion was significantly inhibited by all of the JQ1-containing HAp formulations, that is, with and without medronate, all of the combinations of the targeting compound, medronate, and the chemotherapeutic, JQ1, delivered using HAp, but not HAp alone, inhibited OS cell migration from the tumor spheroids. JQ1 delivered using HAp had an effect on tumor migration, invasion, and apoptosis even at extremely low, subnanomolar concentrations, at which no effect of JQ1 per se was observed, meaning that this form of delivery could help achieve a multifold increase of this drug's efficacy. More than 80% of OS cells internalized JQ1-loaded HAp nanoparticles after 24 h of coincubation, suggesting that this augmentation of the activity of JQ1 may be due to the intracellular delivery and sustained release of the drug enabled by HAp. In addition to the reduction of the OS cell viability, the reduction of the migration and invasion radii was observed in OS tumor spheroids challenged with even JQ1-free medronate-functionalized HAp nanoparticles, demonstrating a definite anticancer activity of medronate alone when combined with HAp. The effect of medronate-functionalized JQ1-loaded HAp nanoparticles was most noticeable against OS cells differentiated into an osteoblastic lineage, in which case they surpassed in effect pure JQ1 and medronate-free compositions. The activity of JQ1 was mediated through increased Ezrin expression and decreased RUNX2 expression and was MYC and FOSL1 independent, but these patterns of gene expression changed in cells challenged with the nanoparticulate form of delivery, having been accompanied by the upregulation of RUNX2 and downregulation of Ezrin in OS cells treated with medronate-functionalized JQ1-loaded HAp nanoparticles.

Machine-Learning-Assisted Approach for Discovering Novel Inhibitors Targeting Bromodomain-Containing Protein 4

Bromodomain-containing protein 4 (BRD4) is implicated in the pathogenesis of a number of different cancers, inflammatory diseases and heart failure. Much effort has been dedicated toward discovering novel scaffold BRD4 inhibitors (BRD4is) with different selectivity profiles and potential antiresistance properties. Structure-based drug design (SBDD) and virtual screening (VS) are the most frequently used approaches. Here, we demonstrate a novel, structure-based VS approach that uses machine-learning algorithms trained on the priori structure and activity knowledge to predict the likelihood that a compound is a BRD4i based on its binding pattern with BRD4. In addition to positive experimental data, such as X-ray structures of BRD4-ligand complexes and BRD4 inhibitory potencies, negative data such as false positives (FPs) identified from our earlier ligand screening results were incorporated into our knowledge base. We used the resulting data to train a machine-learning model named BRD4LGR to predict the BRD4i-likeness of a compound. BRD4LGR achieved a 20-30% higher AUC-ROC than that of Glide using the same test set. When conducting in vitro experiments against a library of previously untested, commercially available organic compounds, the second round of VS using BRD4LGR generated 15 new BRD4is. Moreover, inverting the machine-learning model provided easy access to structure-activity relationship (SAR) interpretation for hit-to-lead optimization.

Attach-Pull-Release Calculations of Ligand Binding and Conformational Changes on the First BRD4 Bromodomain

Bromodomains, protein domains involved in epigenetic regulation, are able to bind small molecules with high affinity. In the present study, we report free energy calculations for the binding of seven ligands to the first BRD4 bromodomain, using the attach-pull-release (APR) method to compute the reversible work of removing the ligands from the binding site and then allowing the protein to relax conformationally. We test three different water models, TIP3P, TIP4PEw, and SPC/E, as well as the GAFF and GAFF2 parameter sets for the ligands. Our simulations show that the apo crystal structure of BRD4 is only metastable, with a structural transition happening in the absence of the ligand typically after 20 ns of simulation. We compute the free energy change for this transition with a separate APR calculation on the free protein and include its contribution to the ligand binding free energies, which generally causes an underestimation of the affinities. By testing different water models and ligand parameters, we are also able to assess their influence in our results and determine which one produces the best agreement with the experimental data. Both free energies associated with the conformational change and ligand binding are affected by the choice of water model, with the two sets of ligand parameters affecting their binding free energies to a lesser degree. Across all six combinations of water model and ligand potential function, the Pearson correlation coefficients between calculated and experimental binding free energies range from 0.55 to 0.83, and the root-mean-square errors range from 1.4-3.2 kcal/mol. The current protocol also yields encouraging preliminary results when used to assess the relative stability of ligand poses generated by docking or other methods, as illustrated for two different ligands. Our method takes advantage of the high performance provided by graphics processing units and can readily be applied to other ligands as well as other protein systems.

Deciphering the mechanisms of selective inhibition for the tandem BD1/BD2 in the BET-bromodomain family

The bromodomain and extra terminal domain (BET) family of bromodomains (BRDs) are well-known drug targets for many human diseases.

Exploring Selective Inhibition of the First Bromodomain of the Human Bromodomain and Extra-terminal Domain (BET) Proteins

A midthroughput screening follow-up program targeting the first bromodomain of the human BRD4 protein, BRD4(BD1), identified an acetylated-mimic xanthine derivative inhibitor. This compound binds with an affinity in the low micromolar range yet exerts suitable unexpected selectivity in vitro against the other members of the bromodomain and extra-terminal domain (BET) family. A structure-based program pinpointed a role of the ZA loop, paving the way for the development of potent and selective BET-BRDi probes.

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.