Turning high-throughput structural biology into predictive inhibitor design.

A common challenge in drug design pertains to finding chemical modifications to a ligand that increases its affinity to the target protein. An underutilized advance is the increase in structural biology throughput, which has progressed from an artisanal endeavor to a monthly throughput of hundreds of different ligands against a protein in modern synchrotrons. However, the missing piece is a framework that turns high-throughput crystallography data into predictive models for ligand design. Here, we designed a simple machine learning approach that predicts protein-ligand affinity from experimental structures of diverse ligands against a single protein paired with biochemical measurements. Our key insight is using physics-based energy descriptors to represent protein-ligand complexes and a learning-to-rank approach that infers the relevant differences between binding modes. We ran a high-throughput crystallography campaign against the SARS-CoV-2 main protease (MPro), obtaining parallel measurements of over 200 protein-ligand complexes and their binding activities. This allows us to design one-step library syntheses which improved the potency of two distinct micromolar hits by over 10-fold, arriving at a noncovalent and nonpeptidomimetic inhibitor with 120 nM antiviral efficacy. Crucially, our approach successfully extends ligands to unexplored regions of the binding pocket, executing large and fruitful moves in chemical space with simple chemistry.

Sulfamate Acetamides as Self-Immolative Electrophiles for Covalent Ligand-Directed Release Chemistry.

Electrophiles for covalent inhibitors that are suitable for in vivo administration are rare. While acrylamides are prevalent in FDA-approved covalent drugs, chloroacetamides are considered too reactive for such purposes. We report sulfamate-based electrophiles that maintain chloroacetamide-like geometry with tunable reactivity. In the context of the BTK inhibitor ibrutinib, sulfamate analogues showed low reactivity with comparable potency in protein labeling, in vitro, and cellular kinase activity assays and were effective in a mouse model of CLL. In a second example, we converted a chloroacetamide Pin1 inhibitor to a potent and selective sulfamate acetamide with improved buffer stability. Finally, we show that sulfamate acetamides can be used for covalent ligand-directed release (CoLDR) chemistry, both for the generation of "turn-on" probes as well as for traceless ligand-directed site-specific labeling of proteins. Taken together, this chemistry represents a promising addition to the list of electrophiles suitable for in vivo covalent targeting.

Small molecule inhibitor of Igf2bp1 represses Kras and a pro-oncogenic phenotype in cancer cells.

Igf2bp1 is an oncofetal RNA binding protein whose expression in numerous types of cancers is associated with upregulation of key pro-oncogenic RNAs, poor prognosis, and reduced survival. Importantly, Igf2bp1 synergizes with mutations in Kras to enhance signalling and oncogenic activity, suggesting that molecules inhibiting Igf2bp1 could have therapeutic potential. Here, we isolate a small molecule that interacts with a hydrophobic surface at the boundary of Igf2bp1 KH3 and KH4 domains, and inhibits binding to Kras RNA. In cells, the compound reduces the level of Kras and other Igf2bp1 mRNA targets, lowers Kras protein, and inhibits downstream signalling, wound healing, and growth in soft agar, all in the absence of any toxicity. This work presents an avenue for improving the prognosis of Igf2bp1-expressing tumours in lung, and potentially other, cancer(s).

Site-Specific Labeling of Endogenous Proteins Using CoLDR Chemistry.

Chemical modifications of native proteins can affect their stability, activity, interactions, localization, and more. However, there are few nongenetic methods for the installation of chemical modifications at a specific protein site in cells. Here we report a covalent ligand directed release (CoLDR) site-specific labeling strategy, which enables the installation of a variety of functional tags on a target protein while releasing the directing ligand. Using this approach, we were able to label various proteins such as BTK, K-RasG12C, and SARS-CoV-2 PLpro with different tags. For BTK we have shown selective labeling in cells of both alkyne and fluorophores tags. Protein labeling by traditional affinity methods often inhibits protein activity since the directing ligand permanently occupies the target binding pocket. We have shown that using CoLDR chemistry, modification of BTK by these probes in cells preserves its activity. We demonstrated several applications for this approach including determining the half-life of BTK in its native environment with minimal perturbation, as well as quantification of BTK degradation by a noncovalent proteolysis targeting chimera (PROTAC) by in-gel fluorescence. Using an environment-sensitive "turn-on" fluorescent probe, we were able to monitor ligand binding to the active site of BTK. Finally, we have demonstrated efficient CoLDR-based BTK PROTACs (DC50 < 100 nM), which installed a CRBN binder onto BTK. This approach joins very few available labeling strategies that maintain the target protein activity and thus makes an important addition to the toolbox of chemical biology.

Tunable Methacrylamides for Covalent Ligand Directed Release Chemistry.

Targeted covalent inhibitors are an important class of drugs and chemical probes. However, relatively few electrophiles meet the criteria for successful covalent inhibitor design. Here we describe α-substituted methacrylamides as a new class of electrophiles suitable for targeted covalent inhibitors. While typically α-substitutions inactivate acrylamides, we show that hetero α-substituted methacrylamides have higher thiol reactivity and undergo a conjugated addition-elimination reaction ultimately releasing the substituent. Their reactivity toward thiols is tunable and correlates with the pKa/pKb of the leaving group. In the context of the BTK inhibitor ibrutinib, these electrophiles showed lower intrinsic thiol reactivity than the unsubstituted ibrutinib acrylamide. This translated to comparable potency in protein labeling, in vitro kinase assays, and functional cellular assays, with improved selectivity. The conjugate addition-elimination reaction upon covalent binding to their target cysteine allows functionalizing α-substituted methacrylamides as turn-on probes. To demonstrate this, we prepared covalent ligand directed release (CoLDR) turn-on fluorescent probes for BTK, EGFR, and K-RasG12C. We further demonstrate a BTK CoLDR chemiluminescent probe that enabled a high-throughput screen for BTK inhibitors. Altogether we show that α-substituted methacrylamides represent a new and versatile addition to the toolbox of targeted covalent inhibitor design.

MitoTimer-based high-content screen identifies two chemically-related benzothiophene derivatives that enhance basal mitophagy.

Mitochondrial turnover is required for proper cellular function. Both mitochondrial biogenesis and mitophagy are impaired in several degenerative and age-related diseases. The search for mitophagy activators recently emerged as a new therapeutical approach; however, there is a lack in suitable tools to follow mitochondrial turnover in a high-throughput manner. We demonstrate that the fluorescent protein, MitoTimer, is a reliable and robust probe to follow mitochondrial turnover. The screening of 15 000 small molecules led us to two chemically-related benzothiophenes that stimulate basal mitophagy in the beta-cell line, INS1. Enhancing basal mitophagy was associated with improved mitochondrial function, higher Complex I activity and Complex II and III expressions in INS1 cells, as well as better insulin secretion performance in mouse islets. The possibility of further enhancing mitophagy in the absence of mitochondrial stressors points to the existence of a 'basal mitophagy spare capacity'. To this end, we found two small molecules that can be used as models to better understand the physiological regulation of mitophagy.

Dasatinib response in acute myeloid leukemia is correlated with FLT3/ITD, PTPN11 mutations and a unique gene expression signature.

Novel targeted therapies demonstrate improved survival in specific subgroups (defined by genetic variants) of acute myeloid leukemia (AML) patients, validating the paradigm of molecularly targeted therapy. However, identifying correlations between AML molecular attributes and effective therapies is challenging. Recent advances in high-throughput in vitro drug sensitivity screening applied to primary AML blasts were used to uncover such correlations; however, these methods cannot predict the response of leukemic stem cells (LSCs). Our study aimed to predict in vitro response to targeted therapies, based on molecular markers, with subsequent validation in LSCs. We performed ex vivo sensitivity screening to 46 drugs on 29 primary AML samples at diagnosis or relapse. Using unsupervised hierarchical clustering analysis we identified group with sensitivity to several tyrosine kinase inhibitors (TKIs), including the multi-TKI, dasatinib, and searched for correlations between dasatinib response, exome sequencing and gene expression from our dataset and from the Beat AML dataset. Unsupervised hierarchical clustering analysis of gene expression resulted in clustering of dasatinib responders and non-responders. In vitro response to dasatinib could be predicted based on gene expression (AUC=0.78). Furthermore, mutations in FLT3/ITD and PTPN11 were enriched in the dasatinib sensitive samples as opposed to mutations in TP53 which were enriched in resistant samples. Based on these results, we selected FLT3/ITD AML samples and injected them to NSG-SGM3 mice. Our results demonstrate that in a subgroup of FLT3/ITD AML (4 out of 9) dasatinib significantly inhibits LSC engraftment. In summary we show that dasatinib has an anti-leukemic effect both on bulk blasts and, more importantly, LSCs from a subset of AML patients that can be identified based on mutational and expression profiles. Our data provide a rational basis for clinical trials of dasatinib in a molecularly selected subset of AML patients.

Structural basis for producing selective MAP2K7 inhibitors.

Mitogen-activated protein kinase kinase 7 (MAP2K7) in the c-Jun N-terminal kinase signal cascade is an attractive drug target for a variety of diseases. The selectivity of MAP2K7 inhibitors against off-target kinases is a major barrier in drug development. We report a crystal structure of MAP2K7 complexed with a potent covalent inhibitor bearing an acrylamide moiety as an electrophile, which discloses a structural basis for producing selective and potent MAP2K7 inhibitors.

Rapid Covalent-Probe Discovery by Electrophile-Fragment Screening.

Covalent probes can display unmatched potency, selectivity, and duration of action; however, their discovery is challenging. In principle, fragments that can irreversibly bind their target can overcome the low affinity that limits reversible fragment screening, but such electrophilic fragments were considered nonselective and were rarely screened. We hypothesized that mild electrophiles might overcome the selectivity challenge and constructed a library of 993 mildly electrophilic fragments. We characterized this library by a new high-throughput thiol-reactivity assay and screened them against 10 cysteine-containing proteins. Highly reactive and promiscuous fragments were rare and could be easily eliminated. In contrast, we found hits for most targets. Combining our approach with high-throughput crystallography allowed rapid progression to potent and selective probes for two enzymes, the deubiquitinase OTUB2 and the pyrophosphatase NUDT7. No inhibitors were previously known for either. This study highlights the potential of electrophile-fragment screening as a practical and efficient tool for covalent-ligand discovery.

A generic, cost-effective, and scalable cell lineage analysis platform.

Advances in single-cell genomics enable commensurate improvements in methods for uncovering lineage relations among individual cells. Current sequencing-based methods for cell lineage analysis depend on low-resolution bulk analysis or rely on extensive single-cell sequencing, which is not scalable and could be biased by functional dependencies. Here we show an integrated biochemical-computational platform for generic single-cell lineage analysis that is retrospective, cost-effective, and scalable. It consists of a biochemical-computational pipeline that inputs individual cells, produces targeted single-cell sequencing data, and uses it to generate a lineage tree of the input cells. We validated the platform by applying it to cells sampled from an ex vivo grown tree and analyzed its feasibility landscape by computer simulations. We conclude that the platform may serve as a generic tool for lineage analysis and thus pave the way toward large-scale human cell lineage discovery.

High-throughput screening identifies broad-spectrum Coronavirus entry inhibitors.

The COVID-19 pandemic highlighted the need for antivirals against emerging coronaviruses (CoV). Inhibiting spike (S) glycoprotein-mediated viral entry is a promising strategy. To identify small molecule inhibitors that block entry downstream of receptor binding, we established a high-throughput screening (HTS) platform based on pseudoviruses. We employed a three-step process to screen nearly 200,000 small molecules. First, we identified hits that inhibit pseudoviruses bearing the SARS-CoV-2 S glycoprotein. Counter-screening against pseudoviruses with the vesicular stomatitis virus glycoprotein (VSV-G), yielded sixty-five SARS-CoV-2 S-specific inhibitors. These were further tested against pseudoviruses bearing the MERS-CoV S glycoprotein, which uses a different receptor. Out of these, five compounds, which included the known broad-spectrum inhibitor Nafamostat, were subjected to further validation and tested against pseudoviruses bearing the S glycoprotein of the Alpha, Delta, and Omicron variants as well as bona fide SARS-CoV-2. This rigorous approach revealed an unreported inhibitor and its derivative as potential broad-spectrum antivirals.

Targeting phenylalanine assemblies as a prospective disease-modifying therapy for phenylketonuria.

Phenylketonuria is characterized by the accumulation of phenylalanine, resulting in severe cognitive and neurological disorders if not treated by a remarkably strict diet. There are two approved drugs today, yet both provide only a partial solution. We have previously demonstrated the formation of amyloid-like toxic assemblies by aggregation of phenylalanine, suggesting a new therapeutic target to be further pursued. Moreover, we showed that compounds that halt the formation of these assemblies also prevent their resulting toxicity. Here, we performed high-throughput screening, searching for compounds with inhibitory effects on phenylalanine aggregation. Morin hydrate, one of the most promising hits revealed during the screen, was chosen to be tested in vivo using a phenylketonuria mouse model. Morin hydrate significantly improved cognitive and motor function with a reduction in the number of phenylalanine brain deposits. Moreover, while phenylalanine levels remained high, we observed a recovery in dopaminergic, adrenergic, and neuronal markers. To conclude, the ability of Morin hydrate to halt phenylalanine aggregation without reducing phenylalanine levels implies the toxic role of the phenylalanine assemblies in phenylketonuria and opens new avenues for disease-modifying treatment.

Targeting SRSF2 mutations in leukemia with RKI-1447: A strategy to impair cellular division and nuclear structure.

Spliceosome machinery mutations are common early mutations in myeloid malignancies; however, effective targeted therapies against them are still lacking. In the current study, we used an in vitro high-throughput drug screen among four different isogenic cell lines and identified RKI-1447, a Rho-associated protein kinase inhibitor, as selective cytotoxic effector of SRSF2 mutant cells. RKI-1447 targeted SRSF2 mutated primary human samples in xenografts models. RKI-1447 induced mitotic catastrophe and induced major reorganization of the microtubule system and severe nuclear deformation. Transmission electron microscopy and 3D light microscopy revealed that SRSF2 mutations induce deep nuclear indentation and segmentation that are apparently driven by microtubule-rich cytoplasmic intrusions, which are exacerbated by RKI-1447. The severe nuclear deformation in RKI-1447-treated SRSF2 mutant cells prevents cells from completing mitosis. These findings shed new light on the interplay between microtubules and the nucleus and offers new ways for targeting pre-leukemic SRSF2 mutant cells.

Identifying a selective inhibitor of autophagy that targets ATG12-ATG3 protein-protein interaction.

Macroautophagy/autophagy is a catabolic process by which cytosolic content is engulfed, degraded and recycled. It has been implicated as a critical pathway in advanced stages of cancer, as it maintains tumor cell homeostasis and continuous growth by nourishing hypoxic or nutrient-starved tumors. Autophagy also supports alternative cellular trafficking pathways, providing a mechanism of non-canonical secretion of inflammatory cytokines. This opens a significant therapeutic opportunity for using autophagy inhibitors in cancer and acute inflammatory responses. Here we developed a high throughput compound screen to identify inhibitors of protein-protein interaction (PPI) in autophagy, based on the protein-fragment complementation assay (PCA). We chose to target the ATG12-ATG3 PPI, as this interaction is indispensable for autophagosome formation, and the analyzed structure of the interaction interface predicts that it may be amenable to inhibition by small molecules. We screened 41,161 compounds yielding 17 compounds that effectively inhibit the ATG12-ATG3 interaction in the PCA platform, and which were subsequently filtered by their ability to inhibit autophagosome formation in viable cells. We describe a lead compound (#189) that inhibited GFP-fused MAP1LC3B/LC3B (microtubule associated protein 1 light chain 3 beta) puncta formation in cells with IC50 value corresponding to 9.3 µM. This compound displayed a selective inhibitory effect on the growth of autophagy addicted tumor cells and inhibited secretion of IL1B/IL-1ß (interleukin 1 beta) by macrophage-like cells. Compound 189 has the potential to be developed into a therapeutic drug and its discovery documents the power of targeting PPIs for acquiring specific and selective compound inhibitors of autophagy.Abbreviations: ANOVA: analysis of variance; ATG: autophagy related; CQ: chloroquine; GFP: green fluorescent protein; GLuc: Gaussia Luciferase; HEK: human embryonic kidney; IL1B: interleukin 1 beta; LPS: lipopolysaccharide; MAP1LC3B/LC3B: microtubule associated protein 1 light chain 3 beta; PCA: protein-fragment complementation assay; PDAC: pancreatic ductal adenocarcinoma; PMA: phorbol 12-myristate 13-acetate; PPI: protein-protein interaction. VCL: vinculin.

Azapeptide activity-based probes for the SARS-CoV-2 main protease enable visualization of inhibition in infected cells.

The COVID-19 pandemic has revealed the vulnerability of the modern, global society. With expected waves of future infections by SARS-CoV-2, treatment options for infected individuals will be crucial in order to decrease mortality and hospitalizations. The SARS-CoV-2 main protease is a validated drug target, for which the first inhibitor has been approved for use in patients. To facilitate future work on this drug target, we designed a solid-phase synthesis route towards azapeptide activity-based probes that are capped with a cysteine-reactive electrophile for covalent modification of the active site of Mpro. This design led to the most potent ABP for Mpro and one of the most potent inhibitors reported thus far. We demonstrate that this ABP can be used to visualize Mpro activity and target engagement by drugs in infected cells.

Specific quinone reductase 2 inhibitors reduce metabolic burden and reverse Alzheimer's disease phenotype in mice.

Biological aging can be described as accumulative, prolonged metabolic stress and is the major risk factor for cognitive decline and Alzheimer's disease (AD). Recently, we identified and described a quinone reductase 2 (QR2) pathway in the brain, in which QR2 acts as a removable memory constraint and metabolic buffer within neurons. QR2 becomes overexpressed with age, and it is possibly a novel contributing factor to age-related metabolic stress and cognitive deficit. We found that, in human cells, genetic removal of QR2 produced a shift in the proteome opposing that found in AD brains while simultaneously reducing oxidative stress. We therefore created highly specific QR2 inhibitors (QR2is) to enable evaluation of chronic QR2 inhibition as a means to reduce biological age-related metabolic stress and cognitive decline. QR2is replicated results obtained by genetic removal of QR2, while local QR2i microinjection improved hippocampal and cortical-dependent learning in rats and mice. Continuous consumption of QR2is in drinking water improved cognition and reduced pathology in the brains of AD-model mice (5xFAD), with a noticeable between-sex effect on treatment duration. These results demonstrate the importance of QR2 activity and pathway function in the healthy and neurodegenerative brain and what we believe to be the great therapeutic potential of QR2is as first-in-class drugs.

Open science discovery of potent noncovalent SARS-CoV-2 main protease inhibitors.

We report the results of the COVID Moonshot, a fully open-science, crowdsourced, and structure-enabled drug discovery campaign targeting the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease. We discovered a noncovalent, nonpeptidic inhibitor scaffold with lead-like properties that is differentiated from current main protease inhibitors. Our approach leveraged crowdsourcing, machine learning, exascale molecular simulations, and high-throughput structural biology and chemistry. We generated a detailed map of the structural plasticity of the SARS-CoV-2 main protease, extensive structure-activity relationships for multiple chemotypes, and a wealth of biochemical activity data. All compound designs (>18,000 designs), crystallographic data (>490 ligand-bound x-ray structures), assay data (>10,000 measurements), and synthesized molecules (>2400 compounds) for this campaign were shared rapidly and openly, creating a rich, open, and intellectual property-free knowledge base for future anticoronavirus drug discovery.

Optimization of Covalent MKK7 Inhibitors via Crude Nanomole-Scale Libraries.

High-throughput nanomole-scale synthesis allows for late-stage functionalization (LSF) of compounds in an efficient and economical manner. Here, we demonstrated that copper-catalyzed azide-alkyne cycloaddition could be used for the LSF of covalent kinase inhibitors at the nanoscale, enabling the synthesis of hundreds of compounds that did not require purification for biological assay screening, thus reducing experimental time drastically. We generated crude libraries of inhibitors for the kinase MKK7, derived from two different parental precursors, and analyzed them via the high-throughput In-Cell Western assay. Select inhibitors were resynthesized, validated via conventional biological and biochemical methods such as western blots and liquid chromatography-mass spectrometry (LC-MS) labeling, and successfully co-crystallized. Two of these compounds showed over 20-fold increased inhibitory activity compared to the parental compound. This study demonstrates that high-throughput LSF of covalent inhibitors at the nanomole-scale level can be an auspicious approach in improving the properties of lead chemical matter.

Non-quaternary oximes detoxify nerve agents and reactivate nerve agent-inhibited human butyrylcholinesterase.

Government-sanctioned use of nerve agents (NA) has escalated dramatically in recent years. Oxime reactivators of organophosphate (OP)-inhibited acetylcholinesterase (AChE) and butyrylcholinesterase (BChE) serve as antidotes toward poisoning by OPNAs. The oximes used as therapeutics are quaternary compounds that cannot penetrate the blood-brain barrier (BBB). There remains an urgent need for the development of next generation OPNA therapeutics. We have developed two high-throughput screening (HTS) assays using a fluorogenic NA surrogate, O-ethyl methylphosphonyl O-4-methyl-3-cyano-coumarin (EMP-MeCyC). EMP-MeCyC detoxification and EMP-BChE reactivation screening campaigns of ~155,000 small molecules resulted in the identification of 33 nucleophile candidates, including non-quaternary oximes. Four of the oximes were reactivators of both Sarin- and VX-inhibited BChE and directly detoxified Sarin. One oxime also detoxified VX. The novel reactivators included a non-quaternary pyridine amidoxime, benzamidoxime, benzaldoxime and a piperidyl-ketoxime. The VX-inhibited BChE reactivation reaction rates by these novel molecules were similar to those observed with known bis-quaternary reactivators and faster than mono-quaternary pyridinium oximes. Notably, we discovered the first ketoxime reactivator of OP-ChEs and detoxifier of OPNAs. Preliminary toxicological studies demonstrated that the newly discovered non-quaternary oximes were relatively non-toxic in mice. The discovery of unique non-quaternary oximes opens the door to the design of novel therapeutics and decontamination agents following OPNA exposure.

An automatic pipeline for the design of irreversible derivatives identifies a potent SARS-CoV-2 Mpro inhibitor.

Designing covalent inhibitors is increasingly important, although it remains challenging. Here, we present covalentizer, a computational pipeline for identifying irreversible inhibitors based on structures of targets with non-covalent binders. Through covalent docking of tailored focused libraries, we identify candidates that can bind covalently to a nearby cysteine while preserving the interactions of the original molecule. We found âˆ¼11,000 cysteines proximal to a ligand across 8,386 complexes in the PDB. Of these, the protocol identified 1,553 structures with covalent predictions. In a prospective evaluation, five out of nine predicted covalent kinase inhibitors showed half-maximal inhibitory concentration (IC50) values between 155 nM and 4.5 µM. Application against an existing SARS-CoV Mpro reversible inhibitor led to an acrylamide inhibitor series with low micromolar IC50 values against SARS-CoV-2 Mpro. The docking was validated by 12 co-crystal structures. Together these examples hint at the vast number of covalent inhibitors accessible through our protocol.