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.

BET proteins: Biological functions and therapeutic interventions

Bromodomain and extra-terminal (BET) family member proteins (BRD2, BRD3, BRD4 and BRDT) play a pivotal role in interpreting the epigenetic information of histone Kac modification, thus controlling gene expression, remodeling chromatin structures and avoid replicative stress-induced DNA damages. Abnormal activation of BET proteins is tightly correlated to various human diseases, including cancer. Therefore, BET bromodomain inhibitors (BBIs) were considered as promising therapeutics to treat BET-related diseases, raising >70 clinical trials in the past decades. Despite preliminary effects achieved, drug resistance and adverse events represent two major challenges for current BBIs development. In this review, we will introduce the biological functions of BET proteins in both physiological and pathological conditions; and summarize the progress in current BBI drug development. Moreover, we will also discuss the major challenges in the front of BET inhibitor development and provide rational strategies to overcome these obstacles.

Bromodomain Interactions with Acetylated Histone 4 Peptides in the BRD4 Tandem Domain: Effects on Domain Dynamics and Internal Flexibility

The bromodomain and extra-terminal (BET) protein BRD4 regulates gene expression via recruitment of transcriptional regulatory complexes to acetylated chromatin. Like other BET proteins, BRD4 contains two bromodomains, BD1 and BD2, that can interact cooperatively with target proteins and designed ligands, with important implications for drug discovery. Here, we used nuclear magnetic resonance (NMR) spectroscopy to study the dynamics and interactions of the isolated bromodomains, as well as the tandem construct including both domains and the intervening linker, and investigated the effects of binding a tetra-acetylated peptide corresponding to the tail of histone 4. The peptide affinity is lower for both domains in the tandem construct than for the isolated domains. Using ¹⁵N spin relaxation, we determined the global rotational correlation times and residue-specific order parameters for BD1 and BD2. Isolated BD1 is monomeric in the apo state but apparently dimerizes upon binding the tetra-acetylated peptide. Isolated BD2 partially dimerizes in both the apo and peptide-bound states. The backbone order parameters reveal marked differences between BD1 and BD2, primarily in the acetyl-lysine binding site where the ZA loop is more flexible in BD2. Peptide binding reduces the order parameters of the ZA loop in BD1 and the ZA and BC loops in BD2. The AB loop, located distally from the binding site, shows variable dynamics that reflect the different dimerization propensities of the domains. These results provide a basis for understanding target recognition by BRD4.

Bromodomain and Extra-Terminal Protein Inhibitors: Biologic Insights and Therapeutic Potential in Pediatric Brain Tumors

Pediatric brain tumors have surpassed leukemia as the leading cause of cancer-related death in children. Several landmark studies from the last two decades have shown that many pediatric brain tumors are driven by epigenetic dysregulation within specific developmental contexts. One of the major determinants of epigenetic control is the histone code, which is orchestrated by a number of enzymes categorized as writers, erasers, and readers. Bromodomain and extra-terminal (BET) proteins are reader proteins that bind to acetylated lysines in histone tails and play a crucial role in regulating gene transcription. BET inhibitors have shown efficacy in a wide range of cancers, and a number have progressed to clinical phase testing. Here, we review the evidence for BET inhibitors in pediatric brain tumor experimental models, as well as their translational potential.

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.

Selectively targeting individual bromodomain: Drug discovery and molecular mechanisms

Bromodomain-containing proteins include bromodomain and extra-terminal (BET) and non-BET families. Due to the conserved bromodomain (BD) module between BD-containing proteins, and especially BETs with each member having two BDs (BD1 and BD2), the high degree of structural similarity makes BD-selective inhibitors much difficult to be designed. However, increasing evidences emphasized that individual BDs had distinct functions and different cellular phenotypes after pharmacological inhibition, and selectively targeting one of the BDs could result in a different efficacy and tolerability profile. This review is to summarize the pioneering progress of BD-selective inhibitors targeting BET and non-BET proteins, focusing on their structural features, biological activity, therapeutic application and experimental/theoretical mechanisms. The present proteolysis targeting chimeras (PROTAC) degraders targeting BDs, and clinical status of BD-selective inhibitors were also analyzed, providing a new insight into future direction of bromodomain-selective drug discovery.

Dasatinib-SIK2 Binding Elucidated by Homology Modeling, Molecular Docking, and Dynamics Simulations

Salt-inducible kinases (SIKs) are calcium/calmodulin-dependent protein kinase (CAMK)-like (CAMKL) family members implicated in insulin signal transduction, metabolic regulation, inflammatory response, and other processes. Here, we focused on SIK2, which is a target of the Food and Drug Administration (FDA)-approved pan inhibitor N-(2-chloro-6-methylphenyl)-2-(6-(4-(2-hydroxyethyl)piperazin-1-yl)-2-methylpyrimidin-4-ylamino)thiazole-5-carboxamide (dasatinib), and constructed four representative SIK2 structures by homology modeling. We investigated the interactions between dasatinib and SIK2 via molecular docking, molecular dynamics simulation, and binding free energy calculation and found that dasatinib showed strong binding affinity for SIK2. Binding free energy calculations suggested that the modification of various dasatinib regions may provide useful information for drug design and to guide the discovery of novel dasatinib-based SIK2 inhibitors.

Efficacy of pH-Sensitive Nanomedicines in Tumors with Different c-MYC Expression Depends on the Intratumoral Activation Profile

Effective inhibition of the protein derived from cellular myelocytomatosis oncogene (c-Myc) is one of the most sought-after goals in cancer therapy. While several c-Myc inhibitors have demonstrated therapeutic potential, inhibiting c-Myc has proven challenging, since c-Myc is essential for normal tissues and tumors may present heterogeneous c-Myc levels demanding contrasting therapeutic strategies. Herein, we developed tumor-targeted nanomedicines capable of treating both tumors with high and low c-Myc levels by adjusting their ability to spatiotemporally control drug action. These nanomedicines loaded homologues of the bromodomain and extraterminal (BET) motif inhibitor JQ1 as epigenetic c-Myc inhibitors through pH-cleavable bonds engineered for fast or slow drug release at intratumoral pH. In tumors with high c-Myc expression, the fast-releasing (FR) nanomedicines suppressed tumor growth more effectively than the slow-releasing (SR) ones, whereas, in the low c-Myc tumors, the efficacy of the nanomedicines was the opposite. By studying the tumor distribution and intratumoral activation of the nanomedicines, we found that, despite SR nanomedicines achieved higher accumulation than the FR counterparts in both c-Myc high and low tumors, the antitumor activity profiles corresponded with the availability of activated drugs inside the tumors. These results indicate the potential of engineered nanomedicines for c-Myc inhibition and spur the idea of precision pH-sensitive nanomedicine based on cancer biomarker levels.

Selective inhibition mechanism of nitroxoline to the BET family: Insight from molecular simulations

Although the proteins in bromodomain and extra-terminal domain (BET) family are promising therapy drug targets for numerous human diseases, the binding effectiveness is interfered by the competition from non-BET protein BRD9. In this study, molecular docking, molecular dynamics simulations, binding free energy calculations and per-residue energy decomposition methods were employed to clarify the selective inhibition mechanism of nitroxoline. The results showed that the different cavity volume of effective embedding inhibitor and the changes in conserved residues were associated with the significant higher selectivity of inhibitor nitroxoline for BET family than non-BET protein (BRD9). In addition, the non-polar interactions occurred in Phe83, Val87 at ZA loop, and the polar interaction appeared in Met132, Asn135 at BC loop. Therefore, when designing a new inhibitor, it could better improve the inhibitor activity by introducing the heteroatom conjugated pyridine-like moiety and the strong electron-withdrawing nitro-like moiety. Overall, this study not only clarified the molecular mechanism of the selective inhibition of nitroxoline, but also provided insight into designing more effective BET inhibitors in next step.

BET bromodomains as novel epigenetic targets for brain health and disease

Epigenetic pharmacotherapy for CNS-related diseases is a burgeoning area of research. In particular, members of the bromodomain and extra-terminal domain (BET) family of proteins have emerged as intriguing therapeutic targets due to their putative involvement in an array of brain diseases. With their ability to bind to acetylated histones and act as a scaffold for chromatin modifying complexes, BET proteins were originally thought of as passive epigenetic 'reader' proteins. However, new research depicts a more complex reality where BET proteins act as key nodes in lineage-specific and signal-dependent transcriptional mechanisms to influence disease-relevant functions. Amid a recent wave of drug development efforts from basic scientists and pharmaceutical companies, BET inhibitors are currently being studied in several CNS-related disease models, but safety and tolerability remain a concern. Here we review the progress in understanding the neurobiological mechanisms of BET proteins and the therapeutic potential of targeting BET proteins for brain health and disease.

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.

Computational study on the selective inhibition mechanism of MS402 to the first and second bromodomains of BRD4

As a member of the bromodomain and extraterminal domain (BET) family, BRD4 is considered as a potential target for cancer treatment. However, because of the highly conservation of its two homologous bromodomains (BD1/BD2), selective inhibition of each bromodomain remains a challenge. MS402 is a domain-selective inhibitor of BRD4-BD1 over BRD4-BD2 reported recently. Understanding the selectivity mechanism would be very useful for the further design of more potent BD1-selectivity inhibitors. Molecular dynamics simulation, adaptive biasing force and multiple-walker adaptive biasing force were performed to study the inhibition and domain-selective mechanism of MS402 toward BRD4-BD1 over BRD4-BD2 here. Results demonstrate BRD4-BD1 binds to MS402 with lower binding free energy than BRD4-BD2. Residues Gln85, Pro86, Asn140, and Ile146 are crucial for MS402's selectively binding to BRD4-BD1. MS402 needs to overcome more energy barrier to dissociate from BD1 than from BD2 pocket. These findings will be helpful for rational structural modification of existing inhibitors to increase their BD1-selectivity.

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.

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.