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.

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.

Immunomodulatory drug (IMiD) resistance is a key clinical challenge in myeloma treatment. Previous data suggest almost one-third of myeloma patients acquire genetic alteration of the key IMiD effector cereblon (CRBN) by the time they are pomalidomide refractory. Some events, including stop codons/frameshift mutations and copy loss, have clearly explicable effects on CRBN protein function. Missense mutations have also been reported throughout the length of CRBN but their functional impact has not been systematically studied. This study modeled selected missense mutations and examined their effect on CRBN function also analyzing whether any mutations deleterious to IMiD action could be overcome using the novel cereblon E3 ligase modulators (CELMoDs). Three patterns of response to missense mutations were apparent: mutations that led to complete loss of CRBN function for all agents, those that had no effect on CRBN function, and those with agent-dependent effect on CRBN function. The latter group of 4 mutations were profiled in more detail with confirmatory experiments demonstrating an ability of the more potent CELMoDs to lead to neosubstrate degradation and cell death even though IMiDs were not active. Dynamic modeling based on a newly generated crystal structure of the DDB1/CRBN/lenalidomide complex, with greater resolution than those published to date, helped to understand the impact of these mutations. These results have important implications for the interpretation of CRBN sequencing results from patients for future therapy decisions, particularly differentiating those who may, despite relapsing on IMiDs with CRBN mutations, have the potential to still benefit from the use of CELMoD agents.

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.