Erratum: iTAG an optimized IMiD-induced degron for targeted protein degradation in human and murine cells.

[This corrects the article DOI: 10.1016/j.isci.2023.107059.].

A Degron Blocking Strategy Towards Improved CRL4CRBN Recruiting PROTAC Selectivity.

Small molecules inducing protein degradation are important pharmacological tools to interrogate complex biology and are rapidly translating into clinical agents. However, to fully realise the potential of these molecules, selectivity remains a limiting challenge. Herein, we addressed the issue of selectivity in the design of CRL4CRBN recruiting PROteolysis TArgeting Chimeras (PROTACs). Thalidomide derivatives used to generate CRL4CRBN recruiting PROTACs have well described intrinsic monovalent degradation profiles by inducing the recruitment of neo-substrates, such as GSPT1, Ikaros and Aiolos. We leveraged structural insights from known CRL4CRBN neo-substrates to attenuate and indeed remove this monovalent degradation function in well-known CRL4CRBN molecular glues degraders, namely CC-885 and Pomalidomide. We then applied these design principles on a previously published BRD9 PROTAC (dBRD9-A) and generated an analogue with improved selectivity profile. Finally, we implemented a computational modelling pipeline to show that our degron blocking design does not impact PROTAC-induced ternary complex formation. We believe that the tools and principles presented in this work will be valuable to support the development of targeted protein degradation.

iTAG an optimized IMiD-induced degron for targeted protein degradation in human and murine cells.

To address the limitation associated with degron based systems, we have developed iTAG, a synthetic tag based on IMiDs/CELMoDs mechanism of action that improves and addresses the limitations of both PROTAC and previous IMiDs/CeLMoDs based tags. Using structural and sequence analysis, we systematically explored native and chimeric degron containing domains (DCDs) and evaluated their ability to induce degradation. We identified the optimal chimeric iTAG(DCD23 60aa) that elicits robust degradation of targets across cell types and subcellular localizations without exhibiting the well documented "hook effect" of PROTAC-based systems. We showed that iTAG can also induce target degradation by murine CRBN and enabled the exploration of natural neo-substrates that can be degraded by murine CRBN. Hence, the iTAG system constitutes a versatile tool to degrade targets across the human and murine proteome.

Binding Mode Prediction and Virtual Screening Applications by Covalent Docking.

The mechanism of action of covalent drugs involves the formation of a bond between their electrophilic warhead group and a nucleophilic residue of the protein target. The recent advances in covalent drug discovery have accelerated the development of computational tools for the design and characterization of covalent binders. Covalent docking algorithms can predict the binding mode of covalent ligands by modeling the bonds and interactions formed at the reaction site. Their scoring functions can estimate the relative binding affinity of ligands towards the target of interest, thus allowing virtual screening of compound libraries. However, most of the scoring schemes have no specific terms for the bond formation, and therefore it prevents the direct comparison of warheads with different intrinsic reactivity. Herein, we describe a protocol for the binding mode prediction of covalent ligands, a typical virtual screening of compound sets with a single warhead chemistry, and an alternative approach to screen libraries that include various warhead types, as applied in recently validated studies.

WIDOCK: a reactive docking protocol for virtual screening of covalent inhibitors.

Here we present WIDOCK, a virtual screening protocol that supports the selection of diverse electrophiles as covalent inhibitors by incorporating ligand reactivity towards cysteine residues into AutoDock4. WIDOCK applies the reactive docking method (Backus et al. in Nature 534:570-574, 2016) and extends it into a virtual screening tool by introducing facile experimental or computational parametrization and a ligand focused evaluation scheme together with a retrospective and prospective validation against various therapeutically relevant targets. Parameters accounting for ligand reactivity are derived from experimental reaction kinetic data or alternatively from computed reaction barriers. The performance of this docking protocol was first evaluated by investigating compound series with diverse warhead chemotypes against KRASG12C, MurA and cathepsin B. In addition, WIDOCK was challenged on larger electrophilic libraries screened against OTUB2 and NUDT7. These retrospective analyses showed high sensitivity in retrieving experimental actives, by also leading to superior ROC curves, AUC values and better enrichments than the standard covalent docking tool available in AutoDock4 when compound collections with diverse warheads were investigated. Finally, we applied WIDOCK for the prospective identification of covalent human MAO-A inhibitors acting via a new mechanism by binding to Cys323. The inhibitory activity of several predicted compounds was experimentally confirmed and the labelling of Cys323 was proved by subsequent MS/MS measurements. These findings demonstrate the usefulness of WIDOCK as a warhead-sensitive, covalent virtual screening protocol.

Covalent Docking in Drug Discovery: Scope and Limitations.

Drug discovery efforts for new covalent inhibitors have drastically increased in the last few years. The binding mechanism of covalent compounds entails the formation of a chemical bond between their electrophilic warhead group and the protein of interest. The use of moderately reactive warheads targeting nonconserved nucleophilic residues can improve the affinity and selectivity profiles of covalent binders as compared to their non-covalent analogs. Recent advances have also enabled their use as chemical probes to disclose novel and also less tractable targets. Increasing interest in covalent drug discovery prompted the development of new computational tools, including covalent docking methods, that are available to predict the binding mode and affinity of covalent ligands. These tools integrate conventional non-covalent docking and scoring schemes by modeling the newly formed covalent bond and the interactions occurring at the reaction site. In this review, we provide a thorough analysis of state-of-the-art covalent docking programs by highlighting their main features and current limitations. Focusing on the implemented algorithms, we show the differences in handling the formation of the new covalent bond and their relative impact on the prediction. This analysis provides a comprehensive overview of the current technology and suggests future improvements in computer-aided covalent drug design. Finally, discussing successful retrospective and prospective covalent docking-based virtual screening applications, we intend to identify best practices for the drug discovery community.

Discovery of a novel kinase hinge binder fragment by dynamic undocking.

One of the key motifs of type I kinase inhibitors is their interactions with the hinge region of ATP binding sites. These interactions contribute significantly to the potency of the inhibitors; however, only a tiny fraction of the available chemical space has been explored with kinase inhibitors reported in the last twenty years. This paper describes a workflow utilizing docking with rDock and dynamic undocking (DUck) for the virtual screening of fragment libraries in order to identify fragments that bind to the kinase hinge region. We have identified 8-amino-2H-isoquinolin-1-one (MR1), a novel and potent hinge binding fragment, which was experimentally tested on a diverse set of kinases, and is hereby suggested for future fragment growing or merging efforts against various kinases, particularly MELK. Direct binding of MR1 to MELK was confirmed by STD-NMR, and its binding to the ATP-pocket was confirmed by a new competitive binding assay based on microscale thermophoresis.

Discovery of Immunoproteasome Inhibitors Using Large-Scale Covalent Virtual Screening.

Large-scale virtual screening of boronic acid derivatives was performed to identify nonpeptidic covalent inhibitors of the ß5i subunit of the immunoproteasome. A hierarchical virtual screening cascade including noncovalent and covalent docking steps was applied to a virtual library of over 104,000 compounds. Then, 32 virtual hits were selected, out of which five were experimentally confirmed. Biophysical and biochemical tests showed micromolar binding affinity and time-dependent inhibitory potency for two compounds. These results validate the computational protocol that allows the screening of large compound collections. One of the lead-like boronic acid derivatives identified as a covalent immunoproteasome inhibitor is a suitable starting point for chemical optimization.

DUckCov: a Dynamic Undocking-Based Virtual Screening Protocol for Covalent Binders.

Thanks to recent guidelines, the design of safe and effective covalent drugs has gained significant interest. Other than targeting non-conserved nucleophilic residues, optimizing the noncovalent binding framework is important to improve potency and selectivity of covalent binders toward the desired target. Significant efforts have been made in extending the computational toolkits to include a covalent mechanism of protein targeting, like in the development of covalent docking methods for binding mode prediction. To highlight the value of the noncovalent complex in the covalent binding process, here we describe a new protocol using tethered and constrained docking in combination with Dynamic Undocking (DUck) as a tool to privilege strong protein binders for the identification of novel covalent inhibitors. At the end of the protocol, dedicated covalent docking methods were used to rank and select the virtual hits based on the predicted binding mode. By validating the method on JAK3 and KRas, we demonstrate how this fast iterative protocol can be applied to explore a wide chemical space and identify potent targeted covalent inhibitors.

A road map for prioritizing warheads for cysteine targeting covalent inhibitors.

Targeted covalent inhibitors have become an integral part of a number of therapeutic protocols and are the subject of intense research. The mechanism of action of these compounds involves the formation of a covalent bond with protein nucleophiles, mostly cysteines. Given the abundance of cysteines in the proteome, the specificity of the covalent inhibitors is of utmost importance and requires careful optimization of the applied warheads. In most of the cysteine targeting covalent inhibitor programs the design strategy involves incorporating Michael acceptors into a ligand that is already known to bind non-covalently. In contrast, we suggest that the reactive warhead itself should be tailored to the reactivity of the specific cysteine being targeted, and we describe a strategy to achieve this goal. Here, we have extended and systematically explored the available organic chemistry toolbox and characterized a large number of warheads representing different chemistries. We demonstrate that in addition to the common Michael addition, there are other nucleophilic addition, addition-elimination, nucleophilic substitution and oxidation reactions suitable for specific covalent protein modification. Importantly, we reveal that warheads for these chemistries impact the reactivity and specificity of covalent fragments at both protein and proteome levels. By integrating surrogate reactivity and selectivity models and subsequent protein assays, we define a road map to help enable new or largely unexplored covalent chemistries for the optimization of cysteine targeting inhibitors.

Comparative Evaluation of Covalent Docking Tools.

Increased interest in covalent drug discovery led to the development of computer programs predicting binding mode and affinity of covalent inhibitors. Here we compare the performance of six covalent docking tools, AutoDock4, CovDock, FITTED, GOLD, ICM-Pro, and MOE, for reproducing experimental binding modes in an unprecedently large and diverse set of covalent complexes. It was found that 40-60% of the top scoring ligand poses are within 2.0 Å RMSD from the experimental binding mode. This rate showed program dependent increase and achieved 50-90% when the best RMSD among the top ten scoring poses was considered. This performance is comparable to that of noncovalent docking tools and therefore suggests that anchoring the ligand does not necessarily improve the accuracy of the prediction. The effect of various ligand and protein features on the docking performance was investigated. At the level of warhead chemistry, higher success rate was found for Michael additions, nucleophilic additions and nucleophilic substitutions than for ring opening reactions and disulfide formation. Increasing ligand size and flexibility generally affects pose predictions unfavorably, although AutoDock4, FITTED, and ICM-Pro were found to be less sensitive up to 35 heavy atoms. Increasing the accessibility of the target cysteine tends to result in improved binding mode predictions. Docking programs show protein dependent performance suggesting a target-dependent choice of the optimal docking tool. It was found that noncovalent docking into Cys/Ala mutated proteins by ICM-Pro and Glide reproduced experimental binding modes with only slightly lower performance and at a significantly lower computational expense than covalent docking did. Overall, our results highlight the key factors influencing the docking performance of the investigated tools and they give guidelines for selecting the optimal combination of warheads, ligands, and tools for the system investigated. Results also identify the most important aspects to be considered for developing improved protocols for docking and virtual screening of covalent ligands.