Drug Discovery in Low Data Regimes: Leveraging a Computational Pipeline for the Discovery of Novel SARS-CoV-2 Nsp14-MTase Inhibitors.

The COVID-19 pandemic, caused by the SARS-CoV-2 virus, has led to significant global morbidity and mortality. A crucial viral protein, the non-structural protein 14 (nsp14), catalyzes the methylation of viral RNA and plays a critical role in viral genome replication and transcription. Due to the low mutation rate in the nsp region among various SARS-CoV-2 variants, nsp14 has emerged as a promising therapeutic target. However, discovering potential inhibitors remains a challenge. In this work, we introduce a computational pipeline for the rapid and efficient identification of potential nsp14 inhibitors by leveraging virtual screening and the NCI open compound collection, which contains 250,000 freely available molecules for researchers worldwide. The introduced pipeline provides a cost-effective and efficient approach for early-stage drug discovery by allowing researchers to evaluate promising molecules without incurring synthesis expenses. Our pipeline successfully identified seven promising candidates after experimentally validating only 40 compounds. Notably, we discovered NSC620333, a compound that exhibits a strong binding affinity to nsp14 with a dissociation constant of 427 ± 84 nM. In addition, we gained new insights into the structure and function of this protein through molecular dynamics simulations. We identified new conformational states of the protein and determined that residues Phe367, Tyr368, and Gln354 within the binding pocket serve as stabilizing residues for novel ligand interactions. We also found that metal coordination complexes are crucial for the overall function of the binding pocket. Lastly, we present the solved crystal structure of the nsp14-MTase complexed with SS148 (PDB:8BWU), a potent inhibitor of methyltransferase activity at the nanomolar level (IC50 value of 70 ± 6 nM). Our computational pipeline accurately predicted the binding pose of SS148, demonstrating its effectiveness and potential in accelerating drug discovery efforts against SARS-CoV-2 and other emerging viruses.

The co-crystal structure of Cbl-b and a small-molecule inhibitor reveals the mechanism of Cbl-b inhibition.

Cbl-b is a RING-type E3 ubiquitin ligase that is expressed in several immune cell lineages, where it negatively regulates the activity of immune cells. Cbl-b has specifically been identified as an attractive target for cancer immunotherapy due to its role in promoting an immunosuppressive tumor environment. A Cbl-b inhibitor, Nx-1607, is currently in phase I clinical trials for advanced solid tumor malignancies. Using a suite of biophysical and cellular assays, we confirm potent binding of C7683 (an analogue of Nx-1607) to the full-length Cbl-b and its N-terminal fragment containing the TKBD-LHR-RING domains. To further elucidate its mechanism of inhibition, we determined the co-crystal structure of Cbl-b with C7683, revealing the compound's interaction with both the TKBD and LHR, but not the RING domain. Here, we provide structural insights into a novel mechanism of Cbl-b inhibition by a small-molecule inhibitor that locks the protein in an inactive conformation by acting as an intramolecular glue.

Structure-Based Discovery of Inhibitors of the SARS-CoV-2 Nsp14 N7-Methyltransferase.

An under-explored target for SARS-CoV-2 is the S-adenosyl methionine (SAM)-dependent methyltransferase Nsp14, which methylates the N7-guanosine of viral RNA at the 5'-end, allowing the virus to evade host immune response. We sought new Nsp14 inhibitors with three large library docking strategies. First, up to 1.1 billion lead-like molecules were docked against the enzyme's SAM site, leading to three inhibitors with IC50 values from 6 to 50 µM. Second, docking a library of 16 million fragments revealed 9 new inhibitors with IC50 values from 12 to 341 µM. Third, docking a library of 25 million electrophiles to covalently modify Cys387 revealed 7 inhibitors with IC50 values from 3.5 to 39 µM. Overall, 32 inhibitors encompassing 11 chemotypes had IC50 values < 50 µM and 5 inhibitors in 4 chemotypes had IC50 values < 10 µM. These molecules are among the first non-SAM-like inhibitors of Nsp14, providing starting points for future optimization.

Crystal structure of SARS-CoV-2 nsp10-nsp16 in complex with small molecule inhibitors, SS148 and WZ16.

SARS-CoV-2 nsp10-nsp16 complex is a 2'-O-methyltransferase (MTase) involved in viral RNA capping, enabling the virus to evade the immune system in humans. It has been considered a valuable target in the discovery of antiviral therapeutics, as the RNA cap formation is crucial for viral propagation. Through cross-screening of the inhibitors that we previously reported for SARS-CoV-2 nsp14 MTase activity against nsp10-nsp16 complex, we identified two compounds (SS148 and WZ16) that also inhibited nsp16 MTase activity. To further enable the chemical optimization of these two compounds towards more potent and selective dual nsp14/nsp16 MTase inhibitors, we determined the crystal structure of nsp10-nsp16 in complex with each of SS148 and WZ16. As expected, the structures revealed the binding of both compounds to S-adenosyl-L-methionine (SAM) binding pocket of nsp16. However, our structural data along with the biochemical mechanism of action determination revealed an RNA-dependent SAM-competitive pattern of inhibition for WZ16, clearly suggesting that binding of the RNA first may help the binding of some SAM competitive inhibitors. Both compounds also showed some degree of selectivity against human protein MTases, an indication of great potential for chemical optimization towards more potent and selective inhibitors of coronavirus MTases.

A chemical probe targeting the PWWP domain alters NSD2 nucleolar localization.

Nuclear receptor-binding SET domain-containing 2 (NSD2) is the primary enzyme responsible for the dimethylation of lysine 36 of histone 3 (H3K36), a mark associated with active gene transcription and intergenic DNA methylation. In addition to a methyltransferase domain, NSD2 harbors two proline-tryptophan-tryptophan-proline (PWWP) domains and five plant homeodomains (PHDs) believed to serve as chromatin reading modules. Here, we report a chemical probe targeting the N-terminal PWWP (PWWP1) domain of NSD2. UNC6934 occupies the canonical H3K36me2-binding pocket of PWWP1, antagonizes PWWP1 interaction with nucleosomal H3K36me2 and selectively engages endogenous NSD2 in cells. UNC6934 induces accumulation of endogenous NSD2 in the nucleolus, phenocopying the localization defects of NSD2 protein isoforms lacking PWWP1 that result from translocations prevalent in multiple myeloma (MM). Mutations of other NSD2 chromatin reader domains also increase NSD2 nucleolar localization and enhance the effect of UNC6934. This chemical probe and the accompanying negative control UNC7145 will be useful tools in defining NSD2 biology.

Probing the SAM Binding Site of SARS-CoV-2 Nsp14 In Vitro Using SAM Competitive Inhibitors Guides Developing Selective Bisubstrate Inhibitors.

The COVID-19 pandemic has clearly brought the healthcare systems worldwide to a breaking point, along with devastating socioeconomic consequences. The SARS-CoV-2 virus, which causes the disease, uses RNA capping to evade the human immune system. Nonstructural protein (nsp) 14 is one of the 16 nsps in SARS-CoV-2 and catalyzes the methylation of the viral RNA at N7-guanosine in the cap formation process. To discover small-molecule inhibitors of nsp14 methyltransferase (MTase) activity, we developed and employed a radiometric MTase assay to screen a library of 161 in-house synthesized S-adenosylmethionine (SAM) competitive MTase inhibitors and SAM analogs. Among six identified screening hits, SS148 inhibited nsp14 MTase activity with an IC50 value of 70 ± 6 nM and was selective against 20 human protein lysine MTases, indicating significant differences in SAM binding sites. Interestingly, DS0464 with an IC50 value of 1.1 ± 0.2 µM showed a bisubstrate competitive inhibitor mechanism of action. DS0464 was also selective against 28 out of 33 RNA, DNA, and protein MTases. The structure-activity relationship provided by these compounds should guide the optimization of selective bisubstrate nsp14 inhibitors and may provide a path toward a novel class of antivirals against COVID-19, and possibly other coronaviruses.

Development of LM98, a Small-Molecule TEAD Inhibitor Derived from Flufenamic Acid.

The YAP-TEAD transcriptional complex is responsible for the expression of genes that regulate cancer cell growth and proliferation. Dysregulation of the Hippo pathway due to overexpression of TEAD has been reported in a wide range of cancers. Inhibition of TEAD represses the expression of associated genes, demonstrating the value of this transcription factor for the development of novel anti-cancer therapies. We report herein the design, synthesis and biological evaluation of LM98, a flufenamic acid analogue. LM98 shows strong affinity to TEAD, inhibits its autopalmitoylation and reduces the YAP-TEAD transcriptional activity. Binding of LM98 to TEAD was supported by 19 F-NMR studies while co-crystallization experiments confirmed that LM98 is anchored within the palmitic acid pocket of TEAD. LM98 reduces the expression of CTGF and Cyr61, inhibits MDA-MB-231 breast cancer cell migration and arrests cell cycling in the S phase during cell division.

The Structure-Based Design of SARS-CoV-2 nsp14 Methyltransferase Ligands Yields Nanomolar Inhibitors.

In this study, we have focused on the structure-based design of the inhibitors of one of the two SARS-CoV-2 methyltransferases (MTases), nsp14. This MTase catalyzes the transfer of the methyl group from S-adenosyl-l-methionine (SAM) to cap the guanosine triphosphate moiety of the newly synthesized viral RNA, yielding the methylated capped RNA and S-adenosyl-l-homocysteine (SAH). As the crystal structure of SARS-CoV-2 nsp14 is unknown, we have taken advantage of its high homology to SARS-CoV nsp14 and prepared its homology model, which has allowed us to identify novel SAH derivatives modified at the adenine nucleobase as inhibitors of this important viral target. We have synthesized and tested the designed compounds in vitro and shown that these derivatives exert unprecedented inhibitory activity against this crucial enzyme. The docking studies nicely explain the contribution of an aromatic part attached by a linker to the position 7 of the 7-deaza analogues of SAH.

Discovery of the SMYD3 Inhibitor BAY-6035 Using Thermal Shift Assay (TSA)-Based High-Throughput Screening.

SMYD3 (SET and MYND domain-containing protein 3) is a protein lysine methyltransferase that was initially described as an H3K4 methyltransferase involved in transcriptional regulation. SMYD3 has been reported to methylate and regulate several nonhistone proteins relevant to cancer, including mitogen-activated protein kinase kinase kinase 2 (MAP3K2), vascular endothelial growth factor receptor 1 (VEGFR1), and the human epidermal growth factor receptor 2 (HER2). In addition, overexpression of SMYD3 has been linked to poor prognosis in certain cancers, suggesting SMYD3 as a potential oncogene and attractive cancer drug target. Here we report the discovery of a novel SMYD3 inhibitor. We performed a thermal shift assay (TSA)-based high-throughput screening (HTS) with 410,000 compounds and identified a novel benzodiazepine-based SMYD3 inhibitor series. Crystal structures revealed that this series binds to the substrate binding site and occupies the hydrophobic lysine binding pocket via an unprecedented hydrogen bonding pattern. Biochemical assays showed substrate competitive behavior. Following optimization and extensive biophysical validation with surface plasmon resonance (SPR) analysis and isothermal titration calorimetry (ITC), we identified BAY-6035, which shows nanomolar potency and selectivity against kinases and other PKMTs. Furthermore, BAY-6035 specifically inhibits methylation of MAP3K2 by SMYD3 in a cellular mechanistic assay with an IC50 <100 nM. Moreover, we describe a congeneric negative control to BAY-6035. In summary, BAY-6035 is a novel selective and potent SMYD3 inhibitor probe that will foster the exploration of the biological role of SMYD3 in diseased and nondiseased tissues.

Probing the SAM Binding Site of SARS-CoV-2 nsp14 in vitro Using SAM Competitive Inhibitors Guides Developing Selective bi-substrate Inhibitors.

The COVID-19 pandemic has clearly brought the healthcare systems world-wide to a breaking point along with devastating socioeconomic consequences. The SARS-CoV-2 virus which causes the disease uses RNA capping to evade the human immune system. Non-structural protein (nsp) 14 is one of the 16 nsps in SARS-CoV-2 and catalyzes the methylation of the viral RNA at N7-guanosine in the cap formation process. To discover small molecule inhibitors of nsp14 methyltransferase (MT) activity, we developed and employed a radiometric MT assay to screen a library of 161 in house synthesized S-adenosylmethionine (SAM) competitive methyltransferase inhibitors and SAM analogs. Among seven identified screening hits, SS148 inhibited nsp14 MT activity with an IC 50 value of 70 ± 6 nM and was selective against 20 human protein lysine methyltransferases indicating significant differences in SAM binding sites. Interestingly, DS0464 with IC 50 value of 1.1 ± 0.2 µM showed a bi-substrate competitive inhibitor mechanism of action. Modeling the binding of this compound to nsp14 suggests that the terminal phenyl group extends into the RNA binding site. DS0464 was also selective against 28 out of 33 RNA, DNA, and protein methyltransferases. The structure-activity relationship provided by these compounds should guide the optimization of selective bi-substrate nsp14 inhibitors and may provide a path towards a novel class of antivirals against COVID-19, and possibly other coronaviruses.

A First-in-Class, Highly Selective and Cell-Active Allosteric Inhibitor of Protein Arginine Methyltransferase 6.

Protein arginine methyltransferase 6 (PRMT6) catalyzes monomethylation and asymmetric dimethylation of arginine residues in various proteins, plays important roles in biological processes, and is associated with multiple cancers. To date, a highly selective PRMT6 inhibitor has not been reported. Here we report the discovery and characterization of a first-in-class, highly selective allosteric inhibitor of PRMT6, (R)-2 (SGC6870). (R)-2 is a potent PRMT6 inhibitor (IC50 = 77 ± 6 nM) with outstanding selectivity for PRMT6 over a broad panel of other methyltransferases and nonepigenetic targets. Notably, the crystal structure of the PRMT6-(R)-2 complex and kinetic studies revealed (R)-2 binds a unique, induced allosteric pocket. Additionally, (R)-2 engages PRMT6 and potently inhibits its methyltransferase activity in cells. Moreover, (R)-2's enantiomer, (S)-2 (SGC6870N), is inactive against PRMT6 and can be utilized as a negative control. Collectively, (R)-2 is a well-characterized PRMT6 chemical probe and a valuable tool for further investigating PRMT6 functions in health and disease.

Discovery of a First-in-Class Protein Arginine Methyltransferase 6 (PRMT6) Covalent Inhibitor.

Protein arginine methyltransferase 6 (PRMT6) plays important roles in several biological processes associated with multiple cancers. Well-characterized potent, selective, and cell-active PRMT6 inhibitors are invaluable tools for testing biological and therapeutic hypotheses. Although there are several known reversible PRMT6 inhibitors, covalent PRMT6 inhibitors have not been reported. Based on a cocrystal structure of PRMT6-MS023 (a type I PRMT inhibitor), we discovered the first potent and cell-active irreversible PRMT6 inhibitor, 4 (MS117). The covalent binding mode of compound 4 to PRMT6 was confirmed by mass spectrometry and kinetic studies and by a cocrystal structure. Compound 4 did not covalently modify other closely related PRMTs, potently inhibited PRMT6 in cells, and was selective for PRMT6 over other methyltransferases. We also developed two structurally similar control compounds, 5 (MS167) and 7 (MS168). We provide these valuable chemical tools to the scientific community for further studying PRMT6 physiological and pathophysiological functions.

A Chemical Probe for Tudor Domain Protein Spindlin1 to Investigate Chromatin Function.

Modifications of histone tails, including lysine/arginine methylation, provide the basis of a "chromatin or histone code". Proteins that contain "reader" domains can bind to these modifications and form specific effector complexes, which ultimately mediate chromatin function. The spindlin1 (SPIN1) protein contains three Tudor methyllysine/arginine reader domains and was identified as a putative oncogene and transcriptional coactivator. Here we report a SPIN1 chemical probe inhibitor with low nanomolar in vitro activity, exquisite selectivity on a panel of methyl reader and writer proteins, and with submicromolar cellular activity. X-ray crystallography showed that this Tudor domain chemical probe simultaneously engages Tudor domains 1 and 2 via a bidentate binding mode. Small molecule inhibition and siRNA knockdown of SPIN1, as well as chemoproteomic studies, identified genes which are transcriptionally regulated by SPIN1 in squamous cell carcinoma and suggest that SPIN1 may have a role in cancer related inflammation and/or cancer metastasis.

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.

Kinetic characterization of human histone H3 lysine 36 methyltransferases, ASH1L and SETD2.

BACKGROUND: Dysregulation of methylation of lysine 36 on histone H3 (H3K36) have been implicated in a variety of diseases including cancers. ASH1L and SETD2 are two enzymes among others that catalyze H3K36 methylation. H3K4 methylation has also been reported for ASH1L. METHODS: Radioactivity-based enzyme assays, Western and immunoblotting using specific antibodies and molecular modeling were used to characterize substrate specificity of ASH1L and SETD2. RESULTS: Here we report on the assay development and kinetic characterization of ASH1L and SETD2 and their substrate specificities in vitro. Both enzymes were active with recombinant nucleosome as substrate. However, SETD2 but not ASH1L methylated histone peptides as well indicating that the interaction of the basic post-SET extension with substrate may not be critical for SETD2 activity. Both enzymes were not active with nucleosome containing a H3K36A mutation indicating their specificity for H3K36. Analyzing the methylation state of the products of ASH1L and SETD2 reactions also confirmed that both enzymes mono- and dimethylate H3K36 and are inactive with H3K4 as substrate, and that only SETD2 is able to trimethylate H3K36 in vitro. CONCLUSIONS: We determined the kinetic parameters for ASH1L and SETD2 activity enabling screening for inhibitors that can be used to further investigate the roles of these two proteins in health and disease. Both ASH1L and SETD2 are H3K36 specific methyltransferases but only SETD2 can trimethylate this mark. The basic post-SET extension is critical for ASH1L but not SETD2 activity. GENERAL SIGNIFICANCE: We provide full kinetic characterization of ASH1L and SETD2 activity.

Trimethylation of histone H3 lysine 36 by human methyltransferase PRDM9 protein.

PRDM9 (PR domain-containing protein 9) is a meiosis-specific protein that trimethylates H3K4 and controls the activation of recombination hot spots. It is an essential enzyme in the progression of early meiotic prophase. Disruption of the PRDM9 gene results in sterility in mice. In human, several PRDM9 SNPs have been implicated in sterility as well. Here we report on kinetic studies of H3K4 methylation by PRDM9 in vitro indicating that PRDM9 is a highly active histone methyltransferase catalyzing mono-, di-, and trimethylation of the H3K4 mark. Screening for other potential histone marks, we identified H3K36 as a second histone residue that could also be mono-, di-, and trimethylated by PRDM9 as efficiently as H3K4. Overexpression of PRDM9 in HEK293 cells also resulted in a significant increase in trimethylated H3K36 and H3K4 further confirming our in vitro observations. Our findings indicate that PRDM9 may play critical roles through H3K36 trimethylation in cells.

Small-molecule inhibition of MLL activity by disruption of its interaction with WDR5.

WDR5 (WD40 repeat protein 5) is an essential component of the human trithorax-like family of SET1 [Su(var)3-9 enhancer-of-zeste trithorax 1] methyltransferase complexes that carry out trimethylation of histone 3 Lys4 (H3K4me3), play key roles in development and are abnormally expressed in many cancers. In the present study, we show that the interaction between WDR5 and peptides from the catalytic domain of MLL (mixed-lineage leukaemia protein) (KMT2) can be antagonized with a small molecule. Structural and biophysical analysis show that this antagonist binds in the WDR5 peptide-binding pocket with a Kd of 450 nM and inhibits the catalytic activity of the MLL core complex in vitro. The degree of inhibition was enhanced at lower protein concentrations consistent with a role for WDR5 in directly stabilizing the MLL multiprotein complex. Our data demonstrate inhibition of an important protein-protein interaction and form the basis for further development of inhibitors of WDR5-dependent enzymes implicated in MLL-rearranged leukaemias or other cancers.

Protein lysine methyltransferase G9a inhibitors: design, synthesis, and structure activity relationships of 2,4-diamino-7-aminoalkoxy-quinazolines.

Protein lysine methyltransferase G9a, which catalyzes methylation of lysine 9 of histone H3 (H3K9) and lysine 373 (K373) of p53, is overexpressed in human cancers. Genetic knockdown of G9a inhibits cancer cell growth, and the dimethylation of p53 K373 results in the inactivation of p53. Initial SAR exploration of the 2,4-diamino-6,7-dimethoxyquinazoline template represented by 3a (BIX01294), a selective small molecule inhibitor of G9a and GLP, led to the discovery of 10 (UNC0224) as a potent G9a inhibitor with excellent selectivity. A high resolution X-ray crystal structure of the G9a-10 complex, the first cocrystal structure of G9a with a small molecule inhibitor, was obtained. On the basis of the structural insights revealed by this cocrystal structure, optimization of the 7-dimethylaminopropoxy side chain of 10 resulted in the discovery of 29 (UNC0321) (Morrison K(i) = 63 pM), which is the first G9a inhibitor with picomolar potency and the most potent G9a inhibitor to date.