The Pan-Canadian Chemical Library: A Mechanism to Open Academic Chemistry to High-Throughput Virtual Screening.

Computationally screening chemical libraries to discover molecules with desired properties is a common technique used in early-stage drug discovery. Recent progress in the field now enables the efficient exploration of billions of molecules within days or hours, but this exploration remains confined within the boundaries of the accessible chemistry space. While the number of commercially available compounds grows rapidly, it remains a limited subset of all druglike small molecules that could be synthesized. Here, we present a workflow where chemical reactions typically developed in academia and unconventional in drug discovery are exploited to dramatically expand the chemistry space accessible to virtual screening. We use this process to generate a first version of the Pan-Canadian Chemical Library, a collection of nearly 150 billion diverse compounds that does not overlap with other ultra-large libraries such as Enamine REAL or SAVI and could be a resource of choice for protein targets where other libraries have failed to deliver bioactive molecules.

Subpocket Similarity-Based Hit Identification for Challenging Targets: Application to the WDR Domain of LRRK2.

We herewith applied a priori a generic hit identification method (POEM) for difficult targets of known three-dimensional structure, relying on the simple knowledge of physicochemical and topological properties of a user-selected cavity. Searching for local similarity to a set of fragment-bound protein microenvironments of known structure, a point cloud registration algorithm is first applied to align known subpockets to the target cavity. The resulting alignment then permits us to directly pose the corresponding seed fragments in a target cavity space not typically amenable to classical docking approaches. Last, linking potentially connectable atoms by a deep generative linker enables full ligand enumeration. When applied to the WD40 repeat (WDR) central cavity of leucine-rich repeat kinase 2 (LRRK2), an unprecedented binding site, POEM was able to quickly propose 94 potential hits, five of which were subsequently confirmed to bind in vitro to LRRK2-WDR.

A chemical probe to modulate human GID4 Pro/N-degron interactions.

The C-terminal to LisH (CTLH) complex is a ubiquitin ligase complex that recognizes substrates with Pro/N-degrons via its substrate receptor Glucose-Induced Degradation 4 (GID4), but its function and substrates in humans remain unclear. Here, we report PFI-7, a potent, selective and cell-active chemical probe that antagonizes Pro/N-degron binding to human GID4. Use of PFI-7 in proximity-dependent biotinylation and quantitative proteomics enabled the identification of GID4 interactors and GID4-regulated proteins. GID4 interactors are enriched for nucleolar proteins, including the Pro/N-degron-containing RNA helicases DDX21 and DDX50. We also identified a distinct subset of proteins whose cellular levels are regulated by GID4 including HMGCS1, a Pro/N-degron-containing metabolic enzyme. These data reveal human GID4 Pro/N-degron targets regulated through a combination of degradative and nondegradative functions. Going forward, PFI-7 will be a valuable research tool for investigating CTLH complex biology and facilitating development of targeted protein degradation strategies that highjack CTLH E3 ligase activity.

Development of Peptide Displacement Assays to Screen for Antagonists of DDB1 Interactions.

The DNA damage binding protein 1 (DDB1) is an essential component of protein complexes involved in DNA damage repair and the ubiquitin-proteasome system (UPS) for protein degradation. As an adaptor protein specific to Cullin-RING E3 ligases, DDB1 binds different receptors that poise protein substrates for ubiquitination and subsequent degradation by the 26S proteasome. Examples of DDB1-binding protein receptors are Cereblon (CRBN) and the WD-repeat containing DDB1- and CUL4-associated factors (DCAFs). Cognate substrates of CRBN and DCAFs are involved in cancer-related cellular processes or are mimicked by viruses to reprogram E3 ligases for the ubiquitination of antiviral host factors. Thus, disrupting interactions of DDB1 with receptor proteins might be an effective strategy for anticancer and antiviral drug discovery. Here, we developed fluorescence polarization (FP)-based peptide displacement assays that utilize full-length DDB1 and fluorescein isothiocyanate (FITC)-labeled peptide probes derived from the specific binding motifs of DDB1 interactors. A general FP-based assay condition applicable to diverse peptide probes was determined and optimized. Mutagenesis and biophysical analyses were then employed to identify the most suitable peptide probe. The FITC-DCAF15 L49A peptide binds DDB1 with a dissociation constant of 68 nM and can be displaced competitively by unlabeled peptides at sub-µM to low nM concentrations. These peptide displacement assays can be used to screen small molecule libraries to identify novel modulators that could specifically antagonize DDB1 interactions toward development of antiviral and cancer therapeutics.

Emerging structure-based computational methods to screen the exploding accessible chemical space.

Structure-based virtual screening can be a valuable approach to computationally select hit candidates based on their predicted interaction with a protein of interest. The recent explosion in the size of chemical libraries increases the chances of hitting high-quality compounds during virtual screening exercises but also poses new challenges as the number of chemically accessible molecules grows faster than the computing power necessary to screen them. We review here two novel approaches rapidly gaining in popularity to address this problem: machine learning-accelerated and synthon-based library screening. We summarize the results from seminal proof-of-concept studies, highlight the latest developments, and discuss limitations and future directions.

Small Molecule Screen Identifies Non-catalytic USP3 Chemical Handle.

Zinc-finger ubiquitin-binding domains (ZnF-UBDs) are noncatalytic domains mostly found in deubiquitylases (DUBs) such as USP3. They represent an underexplored opportunity for the development of deubiquitylase-targeting chimeras (DUBTACs) to pharmacologically induce the deubiquitylation of target proteins. We previously showed that ZnF-UBDs are ligandable domains. Here, a focused small molecule library screen against a panel of 11 ZnF-UBDs led to the identification of compound 59, a ligand engaging the ZnF-UBD of USP3 with a KD of 14 µM. The compound binds the expected C-terminal ubiquitin binding pocket of USP3 as shown by hydrogen-deuterium exchange mass spectrometry experiments and does not inhibit the cleavage of K48-linked diubiquitin by USP3. As such, this molecule is a chemical starting point toward chemical tools that could be used to interrogate the function of the USP3 Znf-UBD and the consequences of recruiting USP3 to ubiquitylated proteins.

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.

Methods for computer-assisted PROTAC design.

Proximity-induced pharmacology is an emerging field in chemical biology and drug discovery where a small molecule induces non-natural interactions between two proteins, leading to specific phenotypic responses. Proteolysis targeting chimeras (PROTACs) are the most mature examples, where ligands for an E3 ligase and a target protein are linked to induce the ubiquitination and subsequent degradation of the target. The discovery of PROTACs typically relies on a trial-and-error approach where chemical handles and linker chemistry, length and attachment points are systematically varied in the hope that one of the combinations will produce an active molecule. Novel computational methods and tools are developed in an attempt to rationalize and accelerate this process and differ significantly from traditional structure-based drug design approaches. In this chapter, we review three different solutions for computer-assisted PROTAC design: MOE, ICM and PRosettaC. Specifically, we describe protocols to predict the structure of ternary complexes (E3 ligase-PROTAC-target protein) and to screen virtually libraries of PROTAC candidates. We also provide troubleshooting tips. Rational PROTAC design is still in its infancy. By opening this space to users and developers, we hope that this methods article will contribute to much needed advancement in the field.

Discovery and Characterization of a Chemical Probe Targeting the Zinc-Finger Ubiquitin-Binding Domain of HDAC6.

Histone deacetylase 6 (HDAC6) inhibition is an attractive strategy for treating numerous cancers, and HDAC6 catalytic inhibitors are currently in clinical trials. The HDAC6 zinc-finger ubiquitin-binding domain (UBD) binds free C-terminal diglycine motifs of unanchored ubiquitin polymer chains and protein aggregates, playing an important role in autophagy and aggresome assembly. However, targeting this domain with small molecule antagonists remains an underdeveloped avenue of HDAC6-focused drug discovery. We report SGC-UBD253 (25), a chemical probe potently targeting HDAC6-UBD in vitro with selectivity over nine other UBDs, except for weak USP16 binding. In cells, 25 is an effective antagonist of HDAC6-UBD at 1 µM, with marked proteome-wide selectivity. We identified SGC-UBD253N (32), a methylated derivative of 25 that is 300-fold less active, serving as a negative control. Together, 25 and 32 could enable further exploration of the biological function of the HDAC6-UBD and investigation of the therapeutic potential of targeting this domain.

Discovery of a Novel DCAF1 Ligand Using a Drug-Target Interaction Prediction Model: Generalizing Machine Learning to New Drug Targets.

DCAF1 functions as a substrate recruitment subunit for the RING-type CRL4DCAF1 and the HECT family EDVPDCAF1 E3 ubiquitin ligases. The WDR domain of DCAF1 serves as a binding platform for substrate proteins and is also targeted by HIV and SIV lentiviral adaptors to induce the ubiquitination and proteasomal degradation of antiviral host factors. It is therefore attractive both as a potential therapeutic target for the development of chemical inhibitors and as an E3 ligase that could be recruited by novel PROTACs for targeted protein degradation. In this study, we used a proteome-scale drug-target interaction prediction model, MatchMaker, combined with cheminformatics filtering and docking to identify ligands for the DCAF1 WDR domain. Biophysical screening and X-ray crystallographic studies of the predicted binders confirmed a selective ligand occupying the central cavity of the WDR domain. This study shows that artificial intelligence-enabled virtual screening methods can successfully be applied in the absence of previously known ligands.

Target 2035 - an update on private sector contributions.

Target 2035, an international federation of biomedical scientists from the public and private sectors, is leveraging 'open' principles to develop a pharmacological tool for every human protein. These tools are important reagents for scientists studying human health and disease and will facilitate the development of new medicines. It is therefore not surprising that pharmaceutical companies are joining Target 2035, contributing both knowledge and reagents to study novel proteins. Here, we present a brief progress update on Target 2035 and highlight some of industry's contributions.

RAS and PP2A activities converge on epigenetic gene regulation.

RAS-mediated human cell transformation requires inhibition of the tumor suppressor protein phosphatase 2A (PP2A). However, the phosphoprotein targets and cellular processes in which RAS and PP2A activities converge in human cancers have not been systematically analyzed. Here, we discover that phosphosites co-regulated by RAS and PP2A are enriched on proteins involved in epigenetic gene regulation. As examples, RAS and PP2A co-regulate the same phosphorylation sites on HDAC1/2, KDM1A, MTA1/2, RNF168, and TP53BP1. We validate RAS- and PP2A-elicited regulation of HDAC1/2 chromatin recruitment, of RNF168-TP53BP1 interaction, and of gene expression. Consistent with their known synergistic effects in cancer, RAS activation and PP2A inhibition resulted in epigenetic reporter derepression and activation of oncogenic transcription. Transcriptional derepression by PP2A inhibition was associated with an increase in euchromatin and a decrease in global DNA methylation. Collectively, the results indicate that epigenetic protein complexes constitute a significant point of convergence for RAS hyperactivity and PP2A inhibition in cancer. Furthermore, the work provides an important resource for future studies focusing on phosphoregulation of epigenetic gene regulation in cancer and in other RAS/PP2A-regulated cellular processes.

ProxyBind: A compendium of binding sites for proximity-induced pharmacology.

Proximity-induced pharmacology (ProxPharm) is a novel paradigm in drug discovery where a small molecule brings two proteins in close proximity to elicit a signal, generally from one protein onto another. The potential of ProxPharm compounds as a new therapeutic modality is firmly established by proteolysis targeting chimeras (PROTACs) that bring an E3 ubiquitin ligase in proximity to a target protein to induce ubiquitination and subsequent degradation of the target. The concept can be expanded to induce other post-translational modifications via the recruitment of different types of protein-modifying enzymes. To survey the human proteome for opportunities in proximity pharmacology, we systematically mapped non-catalytic drug binding pockets on the structure of protein-modifying enzymes available from the Protein Databank. In addition to binding sites exploited by previously reported ProxPharm compounds, we identified putative ligandable non-catalytic pockets in 236 kinases, 45 phosphatases, 37 deubiquitinases, 14 methyltransferases, 11 acetyltransferases, 13 glycosyltransferases, 4 deacetylases, 7 demethylases and 2 glycosidases, including cavities occupied by chemical matter that may serve as starting points for future ProxPharm compounds. This systematic survey confirms that proximity pharmacology is a versatile modality with largely unexplored and promising potential and reveals novel opportunities to pharmacologically rewire molecular circuitries. All data is available from the ProxyBind database at https://polymorph.sgc.utoronto.ca/proxybind/index.php.

Combinatorial Anticancer Drug Screen Identifies Off-Target Effects of Epigenetic Chemical Probes.

Anticancer drug response is determined by genetic and epigenetic mechanisms. To identify the epigenetic regulators of anticancer drug response, we conducted a chemical epigenetic screen using chemical probes that target different epigenetic modulators. In this screen, we tested 31 epigenetic probes in combination with 14 mechanistically diverse anticancer agents and identified 8 epigenetic probes that significantly potentiate the cytotoxicity of TAK-243, a first-in-class ubiquitin-activating enzyme (UBA1) inhibitor evaluated in several solid and hematologic malignancies. These probes are TP-472, GSK864, A-196, UNC1999, SGC-CBP30, and PFI-4 (and its related analogues GSK6853 and GSK5959), and they target BRD9/7, mutant IDH1, SUV420H1/2, EZH2/1, p300/CBP, and BRPF1B, respectively. In contrast to epigenetic probes, negative control compounds did not have a significant impact on TAK-243 cytotoxicity. Potentiation of TAK-243 cytotoxicity was associated with reduced ubiquitylation and induction of apoptosis. Mechanistically, these epigenetic probes exerted their potentiation by inhibiting the efflux transporter ATP-binding cassette subfamily G member 2 (ABCG2) without inducing significant changes in the ubiquitylation pathways or ABCG2 expression levels. As assessed by docking analysis, the identified probes could potentially interact with ABCG2. Based on these data, we have developed a cell-based assay that can quantitatively evaluate ABCG2 inhibition by drug candidates. In conclusion, our study identifies epigenetic probes that profoundly potentiate TAK-243 cytotoxicity through off-target ABCG2 inhibition. We also provide experimental evidence that several negative control compounds cannot exclude a subset of off-target effects of chemical probes. Finally, potentiation of TAK-243 cytotoxicity can serve as a quantitative measure of ABCG2-inhibitory activity.

ChemBioPort: an online portal to navigate the structure, function and chemical inhibition of the human proteome.

Chemical probes are important tools to investigate the function of proteins, evaluate their potential as therapeutic targets and provide chemical starting points for drug discovery. As a result, a growing federation of scientists aims to generate chemical probes for all human druggable proteins. A diverse array of data typically guides target selection and chemical probe discovery: information on protein function can help prioritize targets, domain architecture can provide insight on druggability, structural data enables molecular design and existing chemical ligands can serve as foundation or inspiration for chemical probe development. But these heterogenous data types are dispersed across a variety of public repositories that are difficult to cross-reference by non-experts. We developed ChemBioPort, an online resource that allows users to combine queries related to the ontology, domain architecture or name of human proteins to produce downloadable tables that integrate information on function, disease association, essentiality, tissue enrichment, domain architecture, structure and chemical ligands of proteins. Users can convert these tables into dendrograms reflecting sequence similarity, onto which they can graphically project all data types, linked via a mouse-click to their original repositories or published articles. This interface will support the growing community of chemical biologists, chemists, cell and structural biologists on their perilous journey from genes to medicines. Database URL: https://chembioport.thesgc.org.

CACHE (Critical Assessment of Computational Hit-finding Experiments): A public-private partnership benchmarking initiative to enable the development of computational methods for hit-finding.

One aspirational goal of computational chemistry is to predict potent and drug-like binders for any protein, such that only those that bind are synthesized. In this Roadmap, we describe the launch of Critical Assessment of Computational Hit-finding Experiments (CACHE), a public benchmarking project to compare and improve small molecule hit-finding algorithms through cycles of prediction and experimental testing. Participants will predict small molecule binders for new and biologically relevant protein targets representing different prediction scenarios. Predicted compounds will be tested rigorously in an experimental hub, and all predicted binders as well as all experimental screening data, including the chemical structures of experimentally tested compounds, will be made publicly available, and not subject to any intellectual property restrictions. The ability of a range of computational approaches to find novel binders will be evaluated, compared, and openly published. CACHE will launch 3 new benchmarking exercises every year. The outcomes will be better prediction methods, new small molecule binders for target proteins of importance for fundamental biology or drug discovery, and a major technological step towards achieving the goal of Target 2035, a global initiative to identify pharmacological probes for all human proteins.

Target 2035 - update on the quest for a probe for every protein.

Twenty years after the publication of the first draft of the human genome, our knowledge of the human proteome is still fragmented. The challenge of translating the wealth of new knowledge from genomics into new medicines is that proteins, and not genes, are the primary executers of biological function. Therefore, much of how biology works in health and disease must be understood through the lens of protein function. Accordingly, a subset of human proteins has been at the heart of research interests of scientists over the centuries, and we have accumulated varying degrees of knowledge about approximately 65% of the human proteome. Nevertheless, a large proportion of proteins in the human proteome (∼35%) remains uncharacterized, and less than 5% of the human proteome has been successfully targeted for drug discovery. This highlights the profound disconnect between our abilities to obtain genetic information and subsequent development of effective medicines. Target 2035 is an international federation of biomedical scientists from the public and private sectors, which aims to address this gap by developing and applying new technologies to create by year 2035 chemogenomic libraries, chemical probes, and/or biological probes for the entire human proteome.

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.