The third Critical Assessment of Computational Hit-finding Experiments (CACHE) challenged computational teams to identify chemically novel ligands targeting the macrodomain 1 of SARS-CoV-2 Nsp3, a promising coronavirus drug target. Twenty-three groups deployed diverse design strategies to collectively select 1739 ligand candidates. While over 85% of the designed molecules were chemically novel, the best experimentally confirmed hits were structurally similar to previously published compounds. Confirming a trend observed in CACHE #1 and #2, two of the best-performing workflows used compounds selected by physics-based computational screening methods to train machine learning models able to rapidly screen large chemical libraries, while four others used exclusively physics-based approaches. Three pharmacophore searches and one fragment growing strategy were also part of the seven winning workflows. While active molecules discovered by CACHE #3 participants largely mimicked the adenine ring of the endogenous substrate, ADP-ribose, preserving the canonical chemotype commonly observed in previously reported Nsp3-Mac1 ligands, they still provide novel structure-activity relationship insights that may inform the development of future antivirals. Collectively, these results show that multiple molecular design strategies can efficiently converge on similar potent molecules.
Abstract Huntington’s disease is a fatal neurodegenerative disorder caused by expansion of a cytosine-adenosine-guanine repeat in exon 1 of the Huntingtin (HTT) gene, resulting in a polyglutamine-expanded HTT protein. Although the genetic cause of Huntington’s is well defined, the molecular functions of HTT and the mechanisms linking polyglutamine expansion to neurodegeneration remain incompletely understood. Huntingtin Associated Protein 40 kDa (HAP40) is a key HTT interaction partner that forms a stable complex with HTT and is increasingly recognized as an important player in the HTT structure-function paradigm. However, investigation of HAP40 biology has been limited by a lack of tools capable of directly targeting endogenous protein. Here, we report the discovery and characterization of a panel of nanomolar-affinity macrocyclic peptides targeting HAP40 identified using Random nonstandard Peptide Integrated discovery platform. We characterized macrocycle binding in vitro using surface plasmon resonance, fluorescence polarization, and hydrogen-deuterium exchange mass spectrometry, revealing selective, high-affinity engagement of distinct epitopes on HAP40. We further demonstrate that these macrocycles engage endogenous HAP40 in cellular lysates, enable selective isolation of HAP40-containing protein complexes using macrocycle precipitations and insights into interaction partners of distinct HTT and HAP40 proteoforms. Together, these macrocycles establish a new toolkit for investigating HTT-HAP40 biology and provide a framework for dissecting HAP40-specific functions relevant to Huntington’s disease pathogenesis.
Hydrogen-Deuterium eXchange Mass Spectrometry (HDX-MS) is a dynamics-sensitive structural method that has rapidly achieved widespread adoption in the biopharmaceuticals industry, owing to its ability to quickly identify binding sites and to provide a molecular mechanism of action for drug candidates. However, a central limitation of conventional HDX-MS is that it has substantially lower spatial resolution than other structural techniques, typically providing information averaged over "segments" of five amino acids or more. Here, we demonstrate a sensitive, broadly applicable method for single amino acid resolution, i.e., site-specific HDX-MS measurements. Using a set of five therapeutic candidates targeting the highly druggable cancer target WDR5, we explore the greatly enhanced analytical power that arises from site specificity, including binding mode characterization, affinity ranking, and the detection of features that are "silent" in conventional peptide-level HDX-MS experiments.
Hydrogen-Deuterium eXchange Mass Spectrometry (HDX-MS) is a dynamics-sensitive structural method that has rapidly achieved widespread adoption in the biopharmaceuticals industry, owing to its ability to quickly identify binding sites and to provide a molecular mechanism of action for drug candidates. However, a central limitation of conventional HDX-MS is that it has substantially lower spatial resolution than other structural techniques, typically providing information averaged over "segments" of five amino acids or more. Here, we demonstrate a sensitive, broadly applicable method for single amino acid resolution, i.e., site-specific HDX-MS measurements. Using a set of five therapeutic candidates targeting the highly druggable cancer target WDR5, we explore the greatly enhanced analytical power that arises from site specificity, including binding mode characterization, affinity ranking, and the detection of features that are "silent" in conventional peptide-level HDX-MS experiments.
Huntington's disease (HD) is a fatal neurodegenerative disorder caused by a Cytosine-Adenosine-Guanine (CAG) repeat expansion in the Huntingtin (HTT) gene, with no disease-modifying therapies currently available. The precise molecular function of the HTT protein is unclear, and the lack of selective chemical tools has limited functional studies. We have identified and characterized macrocyclic peptide binders targeting HTT. These binders exhibit low-nanomolar affinity in vitro and engage distinct HTT and HTT-HAP40 interfaces, as revealed by hydrogen-deuterium exchange mass spectrometry and cryoelectron microscopy. Chemoproteomics confirmed selective binding in cell extracts from wildtype but not HTT-null cell lines. HAP40 consistently and stoichiometrically copurified with HTT across cell lines, including with HTT variants containing different CAG repeat lengths, highlighting the broad presence of the HTT-HAP40 complex.
Significant efforts have been recently dedicated to the discovery of small molecule inhibitors against the Macrodomain 1 (Mac1) of nonstructural protein 3 (NSP3) as potential antivirals for SARS-CoV-2. Thus, Mac1 has also been selected as the target for the Critical Assessment of Hit-finding Experiments (CACHE) challenge #3. As contestants in that challenge, we developed a computational strategy that ranked on the top among all 23 participants in the competition and resulted in the discovery of a novel chemical series of non-charged Mac1 inhibitors. Those have been identified through the combination of machine learning-accelerated virtual screening of Enamine REAL Diversity Subset of approximately 25 million compounds and consequent hit expansion into the entire Enamine REAL Space library. In particular, the initially identified hit compound CACHE3-HI_1706_56 (KD = 20 μM) was explored by probing 17 close analogues from a library of 44 billion molecules from the Enamine REAL. All those analogues effectively displaced the Mac1-binding ADP-ribose peptide, and 12 were confirmed to engage with Mac1 by the Surface Plasmon Resonance experiments, revealing a new chemical series of compounds for hit-to-lead optimization. The structure of the CACHE3-HI_1706_56-Mac1 complex was further determined at high resolution with crystallography, confirming initial computational predictions. Our results illustrate the effectiveness of ML-accelerated docking to rapidly identify novel chemical series and provide a strong foundation for the development of SARS-CoV-2 NSP3 Mac1 inhibitors.
We report an enantioselective protein affinity selection mass spectrometry screening approach (E-ASMS) that enables the detection of weak binders, informs on selectivity, and generates orthogonal confirmation of binding. After method development with control proteins, we screen 31 human proteins against a designed library of 8,217 chiral compounds. We identify 16 binders to 12 targets, including many proteins predicted to be "challenging to ligand", and confirm their interactions through orthogonal biophysical assays. Seven binders to six targets display enantioselective binding, with KD values ranging from 3 to 20 µM. Binders for four targets (DDB1, WDR91, WDR55, and HAT1) are selected for in-depth characterization using X-ray crystallography. In all four cases, the mechanisms underlying enantioselectivity are readily explained. These results demonstrate that E-ASMS enables the identification and characterization of selective and weakly binding ligands for novel protein targets with unprecedented throughput and sensitivity.
A promising drug target, SETDB1, is a dual Kme reader and methyltransferase, which has been implicated in cancer and neurodegenerative disease progression. To help understand the role of the triple Tudor domain (3TD) of SETDB1, its Kme reader, we first identified a low micromolar small molecule ligand, UNC6535, which occupies simultaneously both the TD2 and TD3 reader binding sites. Further optimization led to the discovery of UNC10013, the first covalent 3TD ligand targeting Cys385 of SETDB1. UNC10013 is potent with a k inact /K I of 1.0 x 10 6 M -1 s -1 and demonstrated proteome-wide selectivity. In cells, negative allosteric modulation of SETDB1-mediated Akt methylation was observed after treatment with UNC10013. Therefore, UNC10013 is a potent, selective and cell-active covalent ligand for the 3TD of SETDB1, demonstrating negative allosteric modulator properties and making it a promising tool to study the biological role of SETDB1 in disease progression.
Fungal pathogens pose a serious threat to human health, with Candida and Aspergillus spp. representing some of the most significant opportunistic invaders. Aspergillus fumigatus causes aspergillosis, one of the most prevalent fungal diseases of humans. There is a limited number of drugs available to combat these infections, and antifungal drug resistance is on the rise. In this manuscript, we show 4-[4-(4-Piperidinyl) phenyl]-7-[4-(-(trifluoromethyl) phenyl]-2-naphthalenecarboxylic acid (PPTN), a highly specific antagonist of the human P2Y14 receptor, is a promising antifungal adjuvant against diverse fungal pathogens. PPTN interacts with caspofungin (CAS), ibrexafungerp, voriconazole (VOR), and amphotericin against A. fumigatus CAS- and VOR-resistant clinical isolates, and also CAS against Candida spp and Cryptococcus neoformans. The combination of PPTN and CAS increases cell death in A. fumigatus. In the model yeast Saccharomyces cerevisiae, heterozygous deletion of genes involved in chromatin remodeling results in PPTN hypersensitivity, and in A. fumigatus, PPTN can have increased fungicidal activity when combined with the histone deacetylase inhibitor trichostatin A and the DNA methyltransferase inhibitor 5-azacytidine. Finally, PPTN has reduced toxicity to human immortalized cell lineages and partially clears A. fumigatus conidia infection in A549 pulmonary epithelial cells. Our results indicate that PPTN is a novel adjuvant antifungal drug against fungal diseases caused by A. fumigatus and Candida spp. IMPORTANCE:Invasive fungal infections have a high mortality rate, causing more deaths annually than tuberculosis or malaria. Aspergillus fumigatus is the main etiological agent of aspergillosis, one of the most prevalent and deadly fungal diseases. There are few therapeutic options for treating this disease, and treatment commonly fails due to host complications or the emergence of antifungal resistance. Drug repurposing, where existing drugs are deployed for other clinical indications, has increasingly been used in the process of drug discovery. Here, we show that 4-[4-(4-Piperidinyl) phenyl]-7-[4-(-(trifluoromethyl) phenyl]-2-naphthalenecarboxylic acid (PPTN), a highly specific antagonist of the human P2Y14 receptor, when combined with caspofungin (CAS), ibrexafungerp, voriconazole (VOR), and amphotericin can increase the fungicidal activity against not only A. fumigatus CAS- and VOR-resistant clinical isolates but also CAS against Candida spp.
The leucine-rich repeat kinase 2 (LRRK2) is the most mutated gene in familial Parkinson's disease, and its mutations lead to pathogenic hallmarks of the disease. The LRRK2 WDR domain is an understudied drug target for Parkinson's disease, with no known inhibitors prior to the first phase of the Critical Assessment of Computational Hit-Finding Experiments (CACHE) Challenge. A unique advantage of the CACHE Challenge is that the predicted molecules are experimentally validated in-house. Here, we report the design and experimental confirmation of LRRK2 WDR inhibitor molecules. We used an active learning (AL) machine learning (ML) workflow based on optimized free-energy molecular dynamics (MD) simulations utilizing the thermodynamic integration (TI) framework to expand a chemical series around two of our previously confirmed hit molecules. We identified 8 experimentally verified novel inhibitors out of 35 experimentally tested (23% hit rate). These results demonstrate the efficacy of our free-energy-based active learning workflow to explore large chemical spaces quickly and efficiently while minimizing the number and length of expensive simulations. This workflow is widely applicable to screening any chemical space for small-molecule analogs with increased affinity, subject to the general constraints of RBFE calculations. The mean absolute error of the TI MD calculations was 2.69 kcal/mol, with respect to the measured KD of hit compounds.
Critical Assessment of Computational Hit-Finding Experiments (CACHE) Challenges emerged as real-life stress tests for computational hit-finding strategies. In CACHE Challenge #1, 23 participants contributed their original workflows to identify small-molecule ligands for the WD40 repeat (WDR) of LRRK2, a promising Parkinson's target. We applied the FRASE-based hit-finding robot (FRASE-bot), a platform for interaction-based screening allowing a drastic reduction of the explorable chemical space and a concurrent detection of putative ligand-binding sites. In two screening rounds, 84 compounds were procured for experimental testing and 8 were confirmed to bind LRRK2-WDR with dissociation constants (K d) ranging from 3 to 41 μM. To investigate the functional effect of WDR ligands, they were tested for their ability to modify the LRRK2 activity markers in HEK293T cells. Two compounds showed statistically significant increases in the kinase activity of WT LRRK2, and two compounds affected the conformation and kinase activity of major LRRK2 mutants.
Conventional bottom-up HDX-MS experiments are highly suitable for use in drug development; however, a major limitation of this approach is that it generally provides only peptide-level structural resolution. Site specific (i.e., single amino acid-resolved) HDX-MS measurements have been achieved using ECD/ETD, but the low efficiency of these fragmentation techniques, combined with poor ion transmission associated with 'detuning' the instrument to fully prevent deuterium scrambling, results in sensitivity losses that make ligand binding measurements impractical in a 'real-world' (e.g., drug development) context. Here we apply a recently developed method for zero scrambling, high efficiency ECD in the challenging context of ligand binding differential HDX experiments, demonstrating a wealth of additional information that can be acquired when HDX-MS analyses are conducted at the amino acid level.
Human DCAF1 is a multidomain protein that plays a critical role in protein homeostasis. Its WDR domain functions as a substrate recruitment module for RING-type CRL4 and HECT family EDVP E3 ubiquitin ligases, enabling the ubiquitination and proteasomal degradation of specific substrates. DCAF1's activity has been implicated in cell proliferation and is documented to promote tumorigenesis. Additionally, the DCAF1 WDR domain is hijacked by lentiviral accessory proteins to induce the degradation of host antiviral factors, such as SAMHD1 and UNG2. These diverse roles make DCAF1 an attractive target for therapeutic development in oncology and antiviral strategies. It is also a promising candidate for use in targeted protein degradation. We previously reported a novel ligand, OICR-8268, that targets the DCAF1 WDR domain. In this study, we present the development of OICR-41103, a potent, selective, and cell-active small molecule chemical probe for DCAF1, derived from OICR-8268. The co-crystal structure of the DCAF1-OICR-41103 complex reveals the ligand's binding mode within the WDR central pocket, demonstrating its potential for PROTAC design and development. Notably, OICR-41103 effectively displaces the lentiviral Vpr protein from DCAF1 in both biochemical and cellular settings, highlighting its potential for the development of HIV therapeutics.
Huntington’s disease (HD) is caused by a CAG repeat expansion mutation, giving rise to a polyglutamine expansion in the huntingtin (HTT). However, the explicit molecular functions of HTT and opportunities for direct pharmacological modulation remain incompletely understood. Here, we report the discovery of a small molecule ligand for the full-length HTT protein in complex with its partner, HAP40. Using affinity selection mass spectrometry (AS-MS), we identified a stereoselective binder, whose binding was characterized by surface plasmon resonance, hydrogen-deuterium exchange mass spectrometry, and cryo-electron microscopy at 2.3 Å resolution. The ligand binds HTT-HAP40 in vitro with single-digit micromolar affinity and one-to-one stoichiometry at a druggable interface previously predicted computationally. In silico studies predicted and experimental analyses confirmed the (R)-enantiomer as the eutomer and initial structure activity relationship was established experimentally. This work details a structurally-validated chemical scaffold and highlights a ligandable pocket which could enable development of chemical probes for probing HTT biology, as well as therapeutics such as degraders and imaging agents for HD. ### Competing Interest Statement The authors have declared no competing interest.
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.
Target 2035 is a global initiative that seeks to identify a pharmacological modulator of most human proteins by the year 2035. As part of an ongoing series of annual updates of this initiative, we summarise here the efforts of the EUbOPEN project whose objectives and results are making a strong contribution to the goals of Target 2035. EUbOPEN is a public-private partnership with four pillars of activity: (1) chemogenomic library collections, (2) chemical probe discovery and technology development for hit-to-lead chemistry, (3) profiling of bioactive compounds in patient-derived disease assays, and (4) collection, storage and dissemination of project-wide data and reagents. The substantial outputs of this programme include a chemogenomic compound library covering one third of the druggable proteome, as well as 100 chemical probes, both profiled in patient derived assays, as well as hundreds of data sets deposited in existing public data repositories and a project-specific data resource for exploring EUbOPEN outputs.
Hydrogen-deuterium eXchange mass spectrometry (HDX-MS) is increasingly used in drug development to locate binding sites and to identify allosteric effects in drug/target interactions. However, the potential of this technique to quantitatively analyze drug candidate libraries remains largely unexplored. Here, a collection of 13 WDR5-targeting small molecules with surface plasmon resonance (SPR) dissociation coefficients (KD) ranging from 20 nM to ∼116 μM were characterized using differential HDX-MS (ΔHDX-MS). Conventional qualitative analysis of the ΔHDX-MS data set revealed the binding interfaces for all compounds and allosteric effects where present. We then demonstrated that ΔHDX-MS signal-to-noise (S/N) not only can rank library-relative affinity but also can accurately predict KD from a calibration curve constructed from high-quality SPR data. Three methods for S/N calculation are explored, each suitable for libraries with different characteristics. Our results demonstrate the potential for ΔHDX-MS use in drug candidate library affinity validation and/or determination while simultaneously characterizing structure.