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

Perspectives on Computational Enzyme Modeling: From Mechanisms to Design and Drug Development

Understanding enzyme mechanisms is essential for unraveling the complex molecular machinery of life. In this review, we survey the field of computational enzymology, highlighting key principles governing enzyme mechanisms and discussing ongoing challenges and promising advances. Over the years, computer simulations have become indispensable in the study of enzyme mechanisms, with the integration of experimental and computational exploration now established as a holistic approach to gain deep insights into enzymatic catalysis. Numerous studies have demonstrated the power of computer simulations in characterizing reaction pathways, transition states, substrate selectivity, product distribution, and dynamic conformational changes for various enzymes. Nevertheless, significant challenges remain in investigating the mechanisms of complex multistep reactions, large-scale conformational changes, and allosteric regulation. Beyond mechanistic studies, computational enzyme modeling has emerged as an essential tool for computer-aided enzyme design and the rational discovery of covalent drugs for targeted therapies. Overall, enzyme design/engineering and covalent drug development can greatly benefit from our understanding of the detailed mechanisms of enzymes, such as protein dynamics, entropy contributions, and allostery, as revealed by computational studies. Such a convergence of different research approaches is expected to continue, creating synergies in enzyme research. This review, by outlining the ever-expanding field of enzyme research, aims to provide guidance for future research directions and facilitate new developments in this important and evolving field.

Two-Step Covalent Docking with Attracting Cavities

Due to their various advantages, interest in the development of covalent drugs has been renewed in the past few years. It is therefore important to accurately describe and predict their interactions with biological targets by computer-aided drug design tools such as docking algorithms. Here, we report a covalent docking procedure for our in-house docking code Attracting Cavities (AC), which mimics the two-step mechanism of covalent ligand binding. Ligand binding to the protein cavity is driven by nonbonded interactions, followed by the formation of a covalent bond between the ligand and the protein through a chemical reaction. To test the performance of this method, we developed a diverse, high-quality, openly accessible re-docking benchmark set of 95 covalent complexes bound by 8 chemical reactions to 5 different reactive amino acids. Combination with structures from previous studies resulted in a set of 304 complexes, on which AC obtained a success rate (rmsd ≤ 2 Å) of 78%, outperforming two state-of-the-art covalent docking codes, genetic optimization for ligand docking (GOLD (66%)) and AutoDock (AD (35%)). Using a more stringent success criterion (rmsd ≤ 1.5 Å), AC reached a success rate of 71 vs 55% for GOLD and 26% for AD. We additionally assessed the cross-docking performance of AC on a set of 76 covalent complexes of the SARS-CoV-2 main protease. On this challenging test set of mainly small and highly solvent-exposed ligands, AC yielded success rates of 58 and 28% for re-docking and cross-docking, respectively, compared to 45 and 17% for GOLD.

Landscape of In Silico Tools for Modeling Covalent Modification of Proteins: A Review on Computational Covalent Drug Discovery

Covalent drug discovery has been a challenging research area given the struggle of finding a sweet balance between selectivity and reactivity for these drugs, the lack of which often leads to off-target activities and hence undesirable side effects. However, there has been a resurgence in covalent drug design following the success of several covalent drugs such as boceprevir (2011), ibrutinib (2013), neratinib (2017), dacomitinib (2018), zanubrutinib (2019), and many others. Design of covalent drugs includes many crucial factors, where "evaluation of the binding affinity" and "a detailed mechanistic understanding on covalent inhibition" are at the top of the list. Well-defined experimental techniques are available to elucidate these factors; however, often they are expensive and/or time-consuming and hence not suitable for high throughput screens. Recent developments in in silico methods provide promise in this direction. In this report, we review a set of recent publications that focused on developing and/or implementing novel in silico techniques in "Computational Covalent Drug Discovery (CCDD)". We also discuss the advantages and disadvantages of these approaches along with what improvements are required to make it a great tool in medicinal chemistry in the near future.

Computational and biophysical methods for the discovery and optimization of covalent drugs

Drugs that act by covalently attaching to their targets have been used to treat human diseases for over a hundred years. However, the deliberate design of covalent drugs was discouraged due to concerns of toxicity and off-target effects. Recent successes in covalent drug discovery have sparked fresh interest in this field. New screening and testing methods aimed at covalent inhibitors can play pivotal roles in facilitating the discovery process. This feature article focuses on computational and biophysical advances originating from our labs over the past decade and how these approaches have contributed to the design of prolyl oligopeptidase (POP) and SARS-CoV-2 3CLpro covalent inhibitors.

DOX_BDW: Incorporating Solvation and Desolvation Effects of Cavity Water into Nonfitting Protein-Ligand Binding Affinity Prediction

Accurate prediction of the protein-ligand binding affinity (PLBA) with an affordable cost is one of the ultimate goals in the field of structure-based drug design (SBDD), as well as a great challenge in the computational and theoretical chemistry. Herein, we have systematically addressed the complicated solvation and desolvation effects on the PLBA brought by the difference of the explicit water in the protein cavity before and after ligands bind to the protein-binding site. Based on the new solvation model, a nonfitting method at the first-principles level for the PLBA prediction was developed by taking the bridging and displaced water (BDW) molecules into account simultaneously. The newly developed method, DOX_BDW, was validated against a total of 765 noncovalent and covalent protein-ligand binding pairs, including the CASF2016 core set, Cov_2022 covalent binding testing set, and six testing sets for the hit and lead compound optimization (HLO) simulation. In all of the testing sets, the DOX_BDW method was able to produce PLBA predictions that were strongly correlated with the corresponding experimental data (R = 0.66-0.85). The overall performance of DOX_BDW is better than the current empirical scoring functions that are heavily parameterized. DOX_BDW is particularly outstanding for the covalent binding situation, implying the need for considering an electronic structure in covalent drug design. Furthermore, the method is especially recommended to be used in the HLO scenario of SBDD, where hundreds of similar derivatives need to be screened and refined. The computational cost of DOX_BDW is affordable, and its accuracy is remarkable.

Improvement of binding pose prediction of the MR1 covalent ligands by inclusion of simple pharmacophore constraints and structural waters in the docking process

The major histocompatibility complex (MHC) class I-related molecule, MR1, is a key component of the immune system, presenting antigens to T-cell receptors (TCRs) and modulating the immune response against various antigens. MR1 possesses a compact ligand-binding pocket despite its ability to interact with ligands that can have either agonistic or antagonistic effects on the immune system. Agonistic ligands can stimulate the immune response, while antagonistic ligands do not elicit an immune response. In most cases, ligand binding to MR1 is mediated through a covalent bond with Lys43. However, recent studies have suggested that a variety of small molecules can interact with the MR1-binding site. In this study, we have used several approaches to improve the binding pose prediction of covalent ligands to MR1, including docking in mutated receptors, and imposing simple pharmacophore constraints and structural water molecules. The careful assignment of pharmacophore constraints and inclusion of structural water molecules in the challenging docking process of covalent docking improved the binding pose prediction and virtual screening performance. In a retrospective virtual screening, the proposed approach exhibited EF1% and EF2% values of 7.4 and 5.5, respectively. Conversely, when using the mutated receptor, both EF1% and EF2% were recorded as 0 for the conventional docking method. The performance of the pharmacophore constraints was also evaluated on other covalent docking cases, and compared to previously reported results for common covalent docking methods. The proposed approach achieved an average RMSD of 2.55, while AutoDock4, CovDock, FITTED, GOLD, ICM-Pro, and MOE exhibited average RMSD values of 3.0, 2.93, 3.04, 4.93, 2.44, and 3.36, respectively. Our results demonstrate that the inclusion of simple pharmacophore constraints and structural waters can improve the prediction of binding poses of covalent ligands to MR1, which can aid in the discovery of novel immunotherapeutic agents.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-023-03694-w.

Docking covalent targets for drug discovery: stimulating the computer-aided drug design community of possible pitfalls and erroneous practices

The continuous approval of covalent drugs in recent years for the treatment of diseases has led to an increased search for covalent agents by medicinal chemists and computational scientists worldwide. In the computational parlance, molecular docking which is a popular tool to investigate the interaction of a ligand and a protein target, does not account for the formation of covalent bond, and the increasing application of these conventional programs to covalent targets in early drug discovery practice is a matter of utmost concern. Thus, in this comprehensive review, we sought to educate the docking community about the realization of covalent docking and the existence of suitable programs to make their future virtual-screening events on covalent targets worthwhile and scientifically rational. More interestingly, we went beyond the classical description of the functionality of covalent-docking programs down to selecting the 'best' program to consult with during a virtual-screening campaign based on receptor class and covalent warhead chemistry. In addition, we made a highlight on how covalent docking could be achieved using random conventional docking software. And lastly, we raised an alert on the growing erroneous molecular docking practices with covalent targets. Our aim is to guide scientists in the rational docking pursuit when dealing with covalent targets, as this will reduce false-positive results and also increase the reliability of their work for translational research.

SCARdock: A Web Server and Manually Curated Resource for Discovering Covalent Ligands

BACKGROUND

Covalent drugs have been intentionally discarded historically due to the concern of off-target side effects, but the past decade has seen a fast resurgence of the discovery of covalent drugs. Compared to noncovalent ligands, covalent ligands might have better biochemical efficiency, lower patient burden, less dosing frequency, less drug resistance, and improved target specificity.

RESULTS

Previously, we proposed the steric-clashes alleviating receptor (SCAR) strategy for screening and repurposing covalent inhibitors. To help the discovery of covalent ligands targeting protein targets, we have developed a web server dedicated to providing the SCARdock protocol to general users. Along with this server, we presented three high-quality data sets for the discovery of covalent ligands: a manually curated data set containing 954 high-quality complex structures of covalent ligands and proteins, a manually curated data set of 68 experimentally confirmed covalent warheads targeting 11 different residues, and a prefiltered, classified, and ready-to-use data set of 690,018 entries of purchasable virtual compounds containing these experimentally verified warheads.

CONCLUSIONS

The SCARdock server and the accompanied data sets would be of great value to the discovery of covalent ligands and are available freely at http://www.liugroup.site/scardock/ or https://scardock.com.

HCovDock: an efficient docking method for modeling covalent protein-ligand interactions

Covalent inhibitors have received extensive attentions in the past few decades because of their long residence time, high binding efficiency and strong selectivity. Therefore, it is valuable to develop computational tools like molecular docking for modeling of covalent protein-ligand interactions or screening of potential covalent drugs. Meeting the needs, we have proposed HCovDock, an efficient docking algorithm for covalent protein-ligand interactions by integrating a ligand sampling method of incremental construction and a scoring function with covalent bond-based energy. Tested on a benchmark containing 207 diverse protein-ligand complexes, HCovDock exhibits a significantly better performance than seven other state-of-the-art covalent docking programs (AutoDock, Cov_DOX, CovDock, FITTED, GOLD, ICM-Pro and MOE). With the criterion of ligand root-mean-squared distance < 2.0 Å, HCovDock obtains a high success rate of 70.5% and 93.2% in reproducing experimentally observed structures for top 1 and top 10 predictions. In addition, HCovDock is also validated in virtual screening against 10 receptors of three proteins. HCovDock is computationally efficient and the average running time for docking a ligand is only 5 min with as fast as 1 sec for ligands with one rotatable bond and about 18 min for ligands with 23 rotational bonds. HCovDock can be freely assessed at http://huanglab.phys.hust.edu.cn/hcovdock/.

BIreactive: Expanding the Scope of Reactivity Predictions to Propynamides

We present the first comprehensive study on the prediction of reactivity for propynamides. Covalent inhibitors like propynamides often show improved potency, selectivity, and unique pharmacologic properties compared to their non-covalent counterparts. In order to achieve this, it is essential to tune the reactivity of the warhead. This study shows how three different in silico methods can predict the in vitro properties of propynamides, a covalent warhead class integrated into approved drugs on the market. Whereas the electrophilicity index is only applicable to individual subclasses of substitutions, adduct formation and transition state energies have a good predictability for the in vitro reactivity with glutathione (GSH). In summary, the reported methods are well suited to estimate the reactivity of propynamides. With this knowledge, the fine tuning of the reactivity is possible which leads to a speed up of the design process of covalent drugs.

Current Challenges in Small Molecule Proximity-Inducing Compound Development for Targeted Protein Degradation Using the Ubiquitin Proteasomal System

Bivalent proximity-inducing compounds represent a novel class of small molecule therapeutics with exciting potential and new challenges. The most prominent examples of such compounds are utilized in targeted protein degradation where E3 ligases are hijacked to recruit a substrate protein to the proteasome via ubiquitination. In this review we provide an overview of the current state of E3 ligases used in targeted protein degradation, their respective ligands as well as challenges and opportunities that present themselves with these compounds.

Modeling receptor flexibility in the structure-based design of KRASG12C inhibitors

KRAS has long been referred to as an ‘undruggable’ target due to its high affinity for its cognate ligands (GDP and GTP) and its lack of readily exploited allosteric binding pockets. Recent progress in the development of covalent inhibitors of KRAS^G12C has revealed that occupancy of an allosteric binding site located between the α3-helix and switch-II loop of KRAS^G12C—sometimes referred to as the ‘switch-II pocket’—holds great potential in the design of direct inhibitors of KRAS^G12C. In studying diverse switch-II pocket binders during the development of sotorasib (AMG 510), the first FDA-approved inhibitor of KRAS^G12C, we found the dramatic conformational flexibility of the switch-II pocket posing significant challenges toward the structure-based design of inhibitors. Here, we present our computational approaches for dealing with receptor flexibility in the prediction of ligand binding pose and binding affinity. For binding pose prediction, we modified the covalent docking program CovDock to allow for protein conformational mobility. This new docking approach, termed as FlexCovDock, improves success rates from 55 to 89% for binding pose prediction on a dataset of 10 cross-docking cases and has been prospectively validated across diverse ligand chemotypes. For binding affinity prediction, we found standard free energy perturbation (FEP) methods could not adequately handle the significant conformational change of the switch-II loop. We developed a new computational strategy to accelerate conformational transitions through the use of targeted protein mutations. Using this methodology, the mean unsigned error (MUE) of binding affinity prediction were reduced from 1.44 to 0.89 kcal/mol on a set of 14 compounds. These approaches were of significant use in facilitating the structure-based design of KRAS^G12C inhibitors and are anticipated to be of further use in the design of covalent (and noncovalent) inhibitors of other conformationally labile protein targets.

Advances in UDP-N-Acetylglucosamine Enolpyruvyl Transferase (MurA) Covalent Inhibition

Peptidoglycan is a cross-linked polymer responsible for maintaining the bacterial cell wall integrity and morphology in Gram-negative and Gram-positive bacteria. The peptidoglycan pathway consists of the enzymatic reactions held in three steps: cytoplasmic, membrane-associated, and periplasmic. The Mur enzymes (MurA-MurF) are involved in a cytoplasmic stage. The UDP-N-acetylglucosamine enolpyruvyl transferase (MurA) enzyme is responsible for transferring the enolpyruvate group from phosphoenolpyruvate (PEP) to UDP-N-acetylglucosamine (UNAG) to form UDP-N-acetylglucosamine enolpyruvate (EP-UNAG). Fosfomycin is a natural product analogous to PEP that acts on the MurA target enzyme via binding covalently to the key cysteine residue in the active site. Similar to fosfomycin, other MurA covalent inhibitors have been described with a warhead in their structure that forms a covalent bond with the molecular target. In MurA, the nucleophilic thiolate of Cys115 is pointed as the main group involved in the warhead binding. Thus, in this minireview, we briefly describe the main recent advances in the design of MurA covalent inhibitors.

Cov_DOX: A Method for Structure Prediction of Covalent Protein-Ligand Bindings

A handful of molecular docking tools have been extended to enable a covalent docking. However, all of them face the challenge brought by the covalent bond between proteins and ligands. Many covalent drug design scenarios still heavily rely on demanding crystallographic experiments for accurate binding structures. Aiming at filling the gap between covalent dockings and crystallographic experiments, we develop and validate a hybrid method, dubbed as Cov_DOX, in this work. Cov_DOX achieves an overall success rate of 81% with RMSD < 2 Å for the Top 1 pose prediction in the validation against a test set including 405 crystal structures for covalent protein-ligand complexes, covering various types of the warhead chemistry and receptors. Such accuracy is not far from the much more demanding crystallographic experiments, in sharp contrast to the performance of the covalent docking front runners (success rate: 40-60%).

Docking Ligands into Flexible and Solvated Macromolecules. 8. Forming New Bonds─Challenges and Opportunities

Over the years, structure-based design programs and specifically docking small molecules to proteins have become prominent in drug discovery. However, many of these computational tools have been developed to primarily dock enzyme inhibitors (and ligands to other protein classes) relying heavily on hydrogen bonds and electrostatic and hydrophobic interactions. In reality, many drug targets either feature metal ions, can be targeted covalently, or are simply not even proteins (e.g., nucleic acids). Herein, we describe several new features that we have implemented into Fitted to broaden its applicability to a wide range of covalent enzyme inhibitors and to metalloenzymes, where metal coordination is essential for drug binding. This updated version of our docking program was tested for its ability to predict the correct binding mode of drug-sized molecules in a large variety of proteins. We also report new datasets that were essential to demonstrate areas of success and those where additional efforts are required. This resource could be used by other program developers to assess their own software.

Chemical reactivity prediction: current methods and different application areas

The ability to predict chemical reactivity of a molecule is highly desirable in drug discovery, both ex vivo (synthetic route planning, formulation, stability) and in vivo: metabolic reactions determine pharmacodynamics, pharmacokinetics and potential toxic effects, and early assessment of liabilities is vital to reduce attrition rates in later stages of development. Quantum mechanics offer a precise description of the interactions between electrons and orbitals in the breaking and forming of new bonds. Modern algorithms and faster computers have allowed the study of more complex systems in a punctual and accurate fashion, and answers for chemical questions around stability and reactivity can now be provided. Through machine learning, predictive models can be built out of descriptors derived from quantum mechanics and cheminformatics, even in the absence of experimental data to train on. In this article, current progress on computational reactivity prediction is reviewed: applications to problems in drug design, such as modelling of metabolism and covalent inhibition, are highlighted and unmet challenges are posed.

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.