Many thermodynamic quantities can be extracted from computer simulations that generate an ensemble of microstates according to the principles of statistical mechanics. Among these quantities is the ...free energy of binding of a small molecule to a macromolecule, such as a protein. Here, we present an introductory overview of a protocol that allows for the estimation of ligand binding free energies via molecular dynamics simulations. While we focus on the binding of organic molecules to proteins, the approach is in principle transferable to any pair of molecules.
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% O2), the inhibition of the CREBBP/EP300 bromodomain results in the enhanced stabilization of HIF-1α.
Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking ...methods to more computationally expensive molecular‐dynamics (MD)‐based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end‐point method (MM/GBSA), and two MD‐based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD‐based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD‐based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.
Multiple physics‐based methods exist for ranking compounds using their predicted affinity. We probed the ranking ability of these methods against sixteen pharmaceutically relevant targets with a total of 423 known binders. Our results show that, despite system dependency, MD‐based free energy methods should be used for prioritization whenever feasible.
Display omitted
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).
Medicinal chemistry optimization and structural biology of a novel chemical series of covalent PPARG inverse agonists from screening hit compound 7a to BAY-5516, BAY-5094, and BAY-9683. Display ...omitted
PPAR gamma (PPARG) is a ligand activated transcription factor that regulates genes involved in inflammation, bone biology, lipid homeostasis, as well as a master regulator of adipogenesis and a potential lineage driver of luminal bladder cancer. While PPARG agonists lead to transcriptional activation of canonical target genes, inverse agonists have the opposite effect through inducing a transcriptionally repressive complex leading to repression of canonical target gene expression. While many agonists have been described and tested clinically, inverse agonists offer an underexplored avenue to modulate PPARG biology in vivo. Current inverse agonists lack favorable in vivo properties; herein we describe the discovery and characterization of a series of orally bioavailable 4-chloro-6-fluoroisophthalamides as covalent PPARG inverse-agonists, BAY-5516, BAY-5094, and BAY-9683. Structural studies of this series revealed distinct pre- and post-covalent binding positions, which led to the hypothesis that interactions in the pre-covalent conformation are primarily responsible for driving affinity, while interactions in the post-covalent conformation are more responsible for cellular functional effects by enhancing PPARG interactions with its corepressors. The need to simultaneously optimize for two distinct states may partially explain the steep SAR observed. Exquisite selectivity was achieved over related nuclear receptors in the subfamily due in part to a covalent warhead with low reactivity through an SNAr mechanism in addition to the specificity gained through covalent binding to a reactive cysteine uniquely positioned within the PPARG LBD. BAY-5516, BAY-5094, and BAY-9683 lead to pharmacodynamic regulation of PPARG target gene expression in vivo comparable to known inverse agonist SR10221 and represent new tools for future in vivo studies to explore their potential utility for treatment of disorders of hyperactivated PPARG including luminal bladder cancer and other disorders.
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.
TRIM33 is a member of the tripartite motif (TRIM) family of proteins, some of which possess E3 ligase activity and are involved in the ubiquitin-dependent degradation of proteins. Four of the TRIM ...family proteins, TRIM24 (TIF1α), TRIM28 (TIF1β), TRIM33 (TIF1γ) and TRIM66, contain C-terminal plant homeodomain (PHD) and bromodomain (BRD) modules, which bind to methylated lysine (KMe n ) and acetylated lysine (KAc), respectively. Here we investigate the differences between the two isoforms of TRIM33, TRIM33α and TRIM33β, using structural and biophysical approaches. We show that the N1039 residue, which is equivalent to N140 in BRD4(1) and which is conserved in most BRDs, has a different orientation in each isoform. In TRIM33β, this residue coordinates KAc, but this is not the case in TRIM33α. Despite these differences, both isoforms show similar affinities for H31–27K18Ac, and bind preferentially to H31–27K9Me3K18Ac. We used this information to develop an AlphaScreen assay, with which we have identified four new ligands for the TRIM33 PHD-BRD cassette. These findings provide fundamental new information regarding which histone marks are recognized by both isoforms of TRIM33 and suggest starting points for the development of chemical probes to investigate the cellular function of TRIM33.
The Trypanosoma cruzi (T. cruzi) parasite is the cause of Chagas disease, a neglected disease endemic in South America. The life cycle of the T. cruzi parasite is complex and includes transitions ...between distinct life stages. This change in phenotype (without a change in genotype) could be controlled by epigenetic regulation, and might involve the bromodomain-containing factors 1–5 (TcBDF1–5). However, little is known about the function of the TcBDF1–5. Here we describe a fragment-based approach to identify ligands for T. cruzi bromodomain-containing factor 3 (TcBDF3). We expressed a soluble construct of TcBDF3 in E. coli, and used this to develop a range of biophysical assays for this protein. Fragment screening identified 12 compounds that bind to the TcBDF3 bromodomain. On the basis of this screen, we developed functional ligands containing a fluorescence or 19F reporter group, and a photo-crosslinking probe for TcBDF3. These tool compounds will be invaluable in future studies on the function of TcBDF3 and will provide insight into the biology of T. cruzi.
PDF
