Research and Evaluation of the Allosteric Protein-Specific Force Field Based on a Pre-Training Deep Learning Model.

Allosteric modulators are important regulation elements that bind the allosteric site beyond the active site, leading to the changes in dynamic and/or thermodynamic properties of the protein. Allosteric modulators have been a considerable interest as potential drugs with high selectivity and safety. However, current experimental methods have limitations to identify allosteric sites. Therefore, molecular dynamics simulation based on empirical force field becomes an important complement of experimental methods. Moreover, the precision and efficiency of current force fields need improvement. Deep learning and reweighting methods were used to train allosteric protein-specific precise force field (named APSF). Multiple allosteric proteins were used to evaluate the performance of APSF. The results indicate that APSF can capture different types of allosteric pockets and sample multiple energy-minimum reference conformations of allosteric proteins. At the same time, the efficiency of conformation sampling for APSF is higher than that for ff14SB. These findings confirm that the newly developed force field APSF can be effectively used to identify the allosteric pocket that can be further used to screen potential allosteric drugs based on these pockets.

Targeting a cryptic allosteric site of SIRT6 with small-molecule inhibitors that inhibit the migration of pancreatic cancer cells.

SIRT6 belongs to the conserved NAD+-dependent deacetylase superfamily and mediates multiple biological and pathological processes. Targeting SIRT6 by allosteric modulators represents a novel direction for therapeutics, which can overcome the selectivity problem caused by the structural similarity of orthosteric sites among deacetylases. Here, developing a reversed allosteric strategy AlloReverse, we identified a cryptic allosteric site, Pocket Z, which was only induced by the bi-directional allosteric signal triggered upon orthosteric binding of NAD+. Based on Pocket Z, we discovered an SIRT6 allosteric inhibitor named JYQ-42. JYQ-42 selectively targets SIRT6 among other histone deacetylases and effectively inhibits SIRT6 deacetylation, with an IC50 of 2.33 µmol/L. JYQ-42 significantly suppresses SIRT6-mediated cancer cell migration and pro-inflammatory cytokine production. JYQ-42, to our knowledge, is the most potent and selective allosteric SIRT6 inhibitor. This study provides a novel strategy for allosteric drug design and will help in the challenging development of therapeutic agents that can selectively bind SIRT6.

Discovery of the Clinical Candidate MAK683: An EED-Directed, Allosteric, and Selective PRC2 Inhibitor for the Treatment of Advanced Malignancies.

Polycomb Repressive Complex 2 (PRC2) plays an important role in transcriptional regulation during animal development and in cell differentiation, and alteration of PRC2 activity has been associated with cancer. On a molecular level, PRC2 catalyzes methylation of histone H3 lysine 27 (H3K27), resulting in mono-, di-, or trimethylated forms of H3K27, of which the trimethylated form H3K27me3 leads to transcriptional repression of polycomb target genes. Previously, we have shown that binding of the low-molecular-weight compound EED226 to the H3K27me3 binding pocket of the regulatory subunit EED can effectively inhibit PRC2 activity in cells and reduce tumor growth in mouse xenograft models. Here, we report the stepwise optimization of the tool compound EED226 toward the potent and selective EED inhibitor MAK683 (compound 22) and its subsequent preclinical characterization. Based on a balanced PK/PD profile, efficacy, and mitigated risk of forming reactive metabolites, MAK683 has been selected for clinical development.

Computational elucidation of allosteric communication in proteins for allosteric drug design.

Allosteric modulators target topologically distal allosteric sites in order to modulate orthosteric sites, providing enhanced specificity and physiochemical properties. Harnessing allostery for drug discovery is an emerging paradigm in modern pharmaceutics. Allosteric regulation substantially depends on the propagation of allosteric signaling. Delineating allosteric signaling pathways is therefore one of the leading prerequisites for allosteric drug discovery. Allosteric signal transduction is subtle and dynamic, posing challenges for characterization through traditional experimental techniques, but computational strategies promise to provide a solution to this problem. Here, we comprehensively review bioinformatic methods for elucidating allosteric communication, along with their successful applications in allosteric drug design. Current challenges and future perspectives are also discussed. We aim to provide guidance for the future application and optimization of these computational strategies, thereby promoting rational allosteric drug discovery.

Harnessing Reversed Allosteric Communication: A Novel Strategy for Allosteric Drug Discovery.

Allostery is a fundamental and extensive mechanism of intramolecular signal transmission. Allosteric drugs possess several unique pharmacological advantages over traditional orthosteric drugs, including greater selectivity, better physicochemical properties, and lower off-target toxicity. However, owing to the complexity of allosteric regulation, experimental approaches for the development of allosteric modulators are traditionally serendipitous. Recently, the reversed allosteric communication theory has been proposed, providing a feasible tool for the unbiased detection of allosteric sites. Herein, we review the latest research on the reversed allosteric communication effect using the examples of sirtuin 6, epidermal growth factor receptor, 3-phosphoinositide-dependent protein kinase 1, and Related to A and C kinases (RAC) serine/threonine protein kinase B and recapitulate the methodologies of reversed allosteric communication strategy. The novel reversed allosteric communication strategy greatly expands the horizon of allosteric site identification and allosteric mechanism exploration and is expected to accelerate an end-to-end framework for drug discovery.

Allosteric binding sites at the receptor-lipid bilayer interface: novel targets for GPCR drug discovery.

As a superfamily of membrane receptors, G-protein-coupled receptors (GPCRs) have significant roles in human physiological processes, including cell proliferation, metabolism, and neuromodulation. GPCRs are vital targets of therapeutic drugs, and their allosteric regulation represents a novel direction for drug discovery. Given the numerous breakthroughs in structural biology, diverse allosteric sites on GPCRs have been identified within the extracellular and intracellular loops, and the seven core transmembrane helices. However, a unique type of allosteric site has also been discovered at the interface of the receptor-lipid bilayer, similar to the ß2-adrenergic receptor. Here, we review recent identifications of these allosteric sites and the detailed modulator-target interactions within the interface for each modulator to highlight the role of lipids in GPCR allosteric drug discovery.

Discovery of Small-Molecule Antagonists of the H3K9me3 Binding to UHRF1 Tandem Tudor Domain.

Ubiquitin-like with PHD and RING finger domains 1 (UHRF1) is a multidomain protein that plays a critical role in maintaining DNA methylation patterns through concurrent recognition of hemimethylated DNA and histone marks by various domains, and recruitment of DNA methyltransferase 1 (DNMT1). UHRF1 is overexpressed in various cancers, including breast cancer. The tandem tudor domain (TTD) of UHRF1 specifically and tightly binds to histone H3 di- or trimethylated at lysine 9 (H3K9me2 or H3K9me3, respectively), and this binding is essential for UHRF1 function. We developed an H3K9me3 peptide displacement assay, which was used to screen a library of 44,000 compounds for small molecules that disrupt the UHRF1-H3K9me3 interaction. This screen resulted in the identification of NV01, which bound to UHRF1-TTD with a Kd value of 5 µM. The structure of UHRF1-TTD in complex with NV01 confirmed binding to the H3K9me3-binding pocket. Limited structure-based optimization of NV01 led to the discovery of NV03 (Kd of 2.4 µM). These well-characterized small-molecule antagonists of the UHRF1-H3K9me2/3 interaction could be valuable starting chemical matter for developing more potent and cell-active probes toward further characterizing UHRF1 function, with possible applications as anticancer therapeutics.

Discovery of Potent and Selective Allosteric Inhibitors of Protein Arginine Methyltransferase 3 (PRMT3).

PRMT3 catalyzes the asymmetric dimethylation of arginine residues of various proteins. It is crucial for maturation of ribosomes and has been implicated in several diseases. We recently disclosed a highly potent, selective, and cell-active allosteric inhibitor of PRMT3, compound 4. Here, we report comprehensive structure-activity relationship studies that target the allosteric binding site of PRMT3. We conducted design, synthesis, and evaluation of novel compounds in biochemical, selectivity, and cellular assays that culminated in the discovery of 4 and other highly potent (IC50 values: ∼10-36 nM), selective, and cell-active allosteric inhibitors of PRMT3 (compounds 29, 30, 36, and 37). In addition, we generated compounds that are very close analogs of these potent inhibitors but displayed drastically reduced potency as negative controls (compounds 49-51). These inhibitors and negative controls are valuable chemical tools for the biomedical community to further investigate biological functions and disease associations of PRMT3.

An allosteric PRC2 inhibitor targeting the H3K27me3 binding pocket of EED.

Polycomb repressive complex 2 (PRC2) consists of three core subunits, EZH2, EED and SUZ12, and plays pivotal roles in transcriptional regulation. The catalytic subunit EZH2 methylates histone H3 lysine 27 (H3K27), and its activity is further enhanced by the binding of EED to trimethylated H3K27 (H3K27me3). Small-molecule inhibitors that compete with the cofactor S-adenosylmethionine (SAM) have been reported. Here we report the discovery of EED226, a potent and selective PRC2 inhibitor that directly binds to the H3K27me3 binding pocket of EED. EED226 induces a conformational change upon binding EED, leading to loss of PRC2 activity. EED226 shows similar activity to SAM-competitive inhibitors in blocking H3K27 methylation of PRC2 target genes and inducing regression of human lymphoma xenograft tumors. Interestingly, EED226 also effectively inhibits PRC2 containing a mutant EZH2 protein resistant to SAM-competitive inhibitors. Together, we show that EED226 inhibits PRC2 activity via an allosteric mechanism and offers an opportunity for treatment of PRC2-dependent cancers.

Discovery and Molecular Basis of a Diverse Set of Polycomb Repressive Complex 2 Inhibitors Recognition by EED.

Polycomb repressive complex 2 (PRC2), a histone H3 lysine 27 methyltransferase, plays a key role in gene regulation and is a known epigenetics drug target for cancer therapy. The WD40 domain-containing protein EED is the regulatory subunit of PRC2. It binds to the tri-methylated lysine 27 of the histone H3 (H3K27me3), and through which stimulates the activity of PRC2 allosterically. Recently, we disclosed a novel PRC2 inhibitor EED226 which binds to the K27me3-pocket on EED and showed strong antitumor activity in xenograft mice model. Here, we further report the identification and validation of four other EED binders along with EED162, the parental compound of EED226. The crystal structures for all these five compounds in complex with EED revealed a common deep pocket induced by the binding of this diverse set of compounds. This pocket was created after significant conformational rearrangement of the aromatic cage residues (Y365, Y148 and F97) in the H3K27me3 binding pocket of EED, the width of which was delineated by the side chains of these rearranged residues. In addition, all five compounds interact with the Arg367 at the bottom of the pocket. Each compound also displays unique features in its interaction with EED, suggesting the dynamics of the H3K27me3 pocket in accommodating the binding of different compounds. Our results provide structural insights for rational design of novel EED binder for the inhibition of PRC2 complex activity.

Discovery of First-in-Class, Potent, and Orally Bioavailable Embryonic Ectoderm Development (EED) Inhibitor with Robust Anticancer Efficacy.

Overexpression and somatic heterozygous mutations of EZH2, the catalytic subunit of polycomb repressive complex 2 (PRC2), are associated with several tumor types. EZH2 inhibitor, EPZ-6438 (tazemetostat), demonstrated clinical efficacy in patients with acceptable safety profile as monotherapy. EED, another subunit of PRC2 complex, is essential for its histone methyltransferase activity through direct binding to trimethylated lysine 27 on histone 3 (H3K27Me3). Herein we disclose the discovery of a first-in-class potent, selective, and orally bioavailable EED inhibitor compound 43 (EED226). Guided by X-ray crystallography, compound 43 was discovered by fragmentation and regrowth of compound 7, a PRC2 HTS hit that directly binds EED. The ensuing scaffold hopping followed by multiparameter optimization led to the discovery of 43. Compound 43 induces robust and sustained tumor regression in EZH2MUT preclinical DLBCL model. For the first time we demonstrate that specific and direct inhibition of EED can be effective as an anticancer strategy.

Histone methyltransferase SETDB1 regulates liver cancer cell growth through methylation of p53.

SETDB1 is a histone H3K9 methyltransferase that has a critical role in early development. It is located within a melanoma susceptibility locus and facilitates melanoma formation. However, the mechanism by which SETDB1 regulates tumorigenesis remains unknown. Here we report the molecular interplay between SETDB1 and the well-known hotspot gain-of-function (GOF) TP53 R249S mutation. We show that in hepatocellular carcinoma (HCC) SETDB1 is overexpressed with moderate copy number gain, and GOF TP53 mutations including R249S associate with this overexpression. Inactivation of SETDB1 in HCC cell lines bearing the R249S mutation suppresses cell growth. The TP53 mutation status renders cancer cells dependent on SETDB1. Moreover, SETDB1 forms a complex with p53 and catalyses p53K370 di-methylation. SETDB1 attenuation reduces the p53K370me2 level, which subsequently leads to increased recognition and degradation of p53 by MDM2. Together, we provide both genetic and biochemical evidence for a mechanism by which SETDB1 regulates cancer cell growth via methylation of p53.

Conformational Preferences of π-π Stacking Between Ligand and Protein, Analysis Derived from Crystal Structure Data Geometric Preference of π-π Interaction.

π-π Interaction is a direct attractive non-covalent interaction between aromatic moieties, playing an important role in DNA stabilization, drug intercalation, etc. Aromatic rings interact through several different conformations including face-to-face, T-shaped, and offset stacked conformation. Previous quantum calculations indicated that T-shaped and offset stacked conformations are preferred for their smaller electron repulsions. However, substitution group on aromatic ring could have a great impact on π-π interaction by changing electron repulsion force between two rings. To investigate π-π interaction between ligand and aromatic side chain of protein, Brookhaven Protein Data Bank was analyzed. We extracted isolated dimer pairs with the aim of excluding multiple π-π stacking effects and found that T-shaped conformation is prevalent among aromatic interaction between phenyl ring of ligand and protein, which corresponds with the phenomenon of Phe-Phe interactions in small peptide. Specifically, for the non-substitution model, both Phe-Phe and Phenyl-Phe exhibit a favored T-shaped conformation whose dihedral angle is around 50°-70° and centroid distance is between 5.0 and 5.6 Å. However, it could be changed by substituent effect. The hydroxyl group could contact in the case of Tyr-Tyr pairs, while they point away from phenyl plane in Phe-Tyr pairs.

A potent, selective and cell-active allosteric inhibitor of protein arginine methyltransferase†3 (PRMT3).

PRMT3 catalyzes the asymmetric dimethylation of arginine residues of various proteins. It is essential for maturation of ribosomes, may have a role in lipogenesis, and is implicated in several diseases. A potent, selective, and cell-active PRMT3 inhibitor would be a valuable tool for further investigating PRMT3 biology. Here we report the discovery of the first PRMT3 chemical probe, SGC707, by structure-based optimization of the allosteric PRMT3 inhibitors we reported previously, and thorough characterization of this probe in biochemical, biophysical, and cellular assays. SGC707 is a potent PRMT3 inhibitor (IC50 =31±2 nM, KD =53±2 nM) with outstanding selectivity (selective against 31 other methyltransferases and more than 250 non-epigenetic targets). The mechanism of action studies and crystal structure of the PRMT3-SGC707 complex confirm the allosteric inhibition mode. Importantly, SGC707 engages PRMT3 and potently inhibits its methyltransferase activity in cells. It is also bioavailable and suitable for animal studies. This well-characterized chemical probe is an excellent tool to further study the role of PRMT3 in health and disease.

Global chromatin profiling reveals NSD2 mutations in pediatric acute lymphoblastic leukemia.

Epigenetic dysregulation is an emerging hallmark of cancers. We developed a high-information-content mass spectrometry approach to profile global histone modifications in human cancers. When applied to 115 lines from the Cancer Cell Line Encyclopedia, this approach identified distinct molecular chromatin signatures. One signature was characterized by increased histone 3 lysine 36 (H3K36) dimethylation, exhibited by several lines harboring translocations in NSD2, which encodes a methyltransferase. A previously unknown NSD2 p.Glu1099Lys (p.E1099K) variant was identified in nontranslocated acute lymphoblastic leukemia (ALL) cell lines sharing this signature. Ectopic expression of the variant induced a chromatin signature characteristic of NSD2 hyperactivation and promoted transformation. NSD2 knockdown selectively inhibited the proliferation of NSD2-mutant lines and impaired the in vivo growth of an NSD2-mutant ALL xenograft. Sequencing analysis of >1,000 pediatric cancer genomes identified the NSD2 p.E1099K alteration in 14% of t(12;21) ETV6-RUNX1-containing ALLs. These findings identify NSD2 as a potential therapeutic target for pediatric ALL and provide a general framework for the functional annotation of cancer epigenomes.

Selective inhibition of Ezh2 by a small molecule inhibitor blocks tumor cells proliferation.

Ezh2 (Enhancer of zeste homolog 2) protein is the enzymatic component of the Polycomb repressive complex 2 (PRC2), which represses gene expression by methylating lysine 27 of histone H3 (H3K27) and regulates cell proliferation and differentiation during embryonic development. Recently, hot-spot mutations of Ezh2 were identified in diffused large B-cell lymphomas and follicular lymphomas. To investigate if tumor growth is dependent on the enzymatic activity of Ezh2, we developed a potent and selective small molecule inhibitor, EI1, which inhibits the enzymatic activity of Ezh2 through direct binding to the enzyme and competing with the methyl group donor S-Adenosyl methionine. EI1-treated cells exhibit genome-wide loss of H3K27 methylation and activation of PRC2 target genes. Furthermore, inhibition of Ezh2 by EI1 in diffused large B-cell lymphomas cells carrying the Y641 mutations results in decreased proliferation, cell cycle arrest, and apoptosis. These results provide strong validation of Ezh2 as a potential therapeutic target for the treatment of cancer.

Molecular dynamics simulations of CYP2E1.

CYP2E1, as a member of the cytochrome P450s (CYPs) super-family, is in charge of six percent drug metabolism involving a diversity of drugs distinct in structures and chemical properties, such as alcohols, monocyclic compounds (e.g., acetaminophen, benzene, p-nitrophenol), bicyclic heterocycles (e.g., coumarin, caffeine) and even fatty acids. The aromatic molecules form a vital species catalyzed by CYP2E1. To investigate the mechanism of metabolizing a diversity of aromatic molecules, five representative aromatic substrates were selected: (1) benzene, the non-polar simple ring; (2) aniline, the monocyclic substrate with smallest substitution on the phenyl ring; (3) acetaminophen, a large monocyclic substrate with highly active reactivity; (4) chlorzoxazone, and (5) theophylline, the bicyclic substrates with low or high catalytic activities. They were docked into X-ray structure of CYP2E1, after which all-atom molecular dynamics simulations of 5 ns were performed on each model. It was found that the active site interact with the aromatic substrates mainly through π-π stacking, supplied by five hydrophobic phenylalanines in the active site. Our simulations also illustrated the specific movement of different kinds of aromatic substrates in the pocket. Small monocyclic substrates show highly frequent self-rotation and limited translation movement. Substrates with single catalytic position are less movable in the pocket than substrates with multiple products. All these findings are quite useful for understanding the catalytic mechanism of CYP2E1, stimulating novel strategies for conducting further mutagenesis studies for specific drug design.

Identification of a selective thieno[2,3-c]pyridine inhibitor of COT kinase and TNF-alpha production.

COT (Tpl2 in mice) is a serine/threonine MAP3 kinase that regulates production of TNF-alpha and other pro-inflammatory cytokines such as IL-1beta via the ERK/MAP kinase pathway. As TNF-alpha and IL-1beta are clinically validated targets for therapeutic intervention in rheumatoid arthritis (RA), blocking COT provides a potential avenue for amelioration of disease. Herein we describe identification of a cellular active selective small molecule inhibitor of COT kinase.

Discovery of thieno[2,3-c]pyridines as potent COT inhibitors.

Evaluation of hit chemotypes from high throughput screening identified a novel series of 2,4-disubstituted thieno[2,3-c]pyridines as COT kinase inhibitors. Structural modifications exploring SAR at the 2- and 4-positions resulting in inhibitors with improved enzyme potency and cellular activity are disclosed.