Assessing Lysine and Cysteine Reactivities for Designing Targeted Covalent Kinase Inhibitors

Fragment-based covalent ligand discovery

Targeted covalent inhibitors have regained widespread attention in drug discovery and have emerged as powerful tools for basic biomedical research. Fueled by considerable improvements in mass spectrometry sensitivity and sample processing, chemoproteomic strategies have revealed thousands of proteins that can be covalently modified by reactive small molecules. Fragment-based drug discovery, which has traditionally been used in a target-centric fashion, is now being deployed on a proteome-wide scale thereby expanding its utility to both the discovery of novel covalent ligands and their cognate protein targets. This powerful approach is allowing ‘high-throughput’ serendipitous discovery of cryptic pockets leading to the identification of pharmacological modulators of proteins previously viewed as “undruggable”. The reactive fragment toolkit has been enabled by recent advances in the development of new chemistries that target residues other than cysteine including lysine and tyrosine. Here, we review the emerging area of covalent fragment-based ligand discovery, which integrates the benefits of covalent targeting and fragment-based medicinal chemistry. We discuss how the two strategies synergize to facilitate the efficient discovery of new pharmacological modulators of established and new therapeutic target proteins.

Covalent fragment-based ligand discovery greatly facilitates the discovery of useful fragments for drug discovery and helps unveil chemical-tractable biological targets in native biological systems.

Druggable Transient Pockets in Protein Kinases

Drug discovery using small molecule inhibitors is reaching a stalemate due to low selectivity, adverse off-target effects and inevitable failures in clinical trials. Conventional chemical screening methods may miss potent small molecules because of their use of simple but outdated kits composed of recombinant enzyme proteins. Non-canonical inhibitors targeting a hidden pocket in a protein have received considerable research attention. Kii and colleagues identified an inhibitor targeting a transient pocket in the kinase DYRK1A during its folding process and termed it FINDY. FINDY exhibits a unique inhibitory profile; that is, FINDY does not inhibit the fully folded form of DYRK1A, indicating that the FINDY-binding pocket is hidden in the folded form. This intriguing pocket opens during the folding process and then closes upon completion of folding. In this review, we discuss previously established kinase inhibitors and their inhibitory mechanisms in comparison with FINDY. We also compare the inhibitory mechanisms with the growing concept of “cryptic inhibitor-binding sites.” These sites are buried on the inhibitor-unbound surface but become apparent when the inhibitor is bound. In addition, an alternative method based on cell-free protein synthesis of protein kinases may allow the discovery of small molecules that occupy these mysterious binding sites. Transitional folding intermediates would become alternative targets in drug discovery, enabling the efficient development of potent kinase inhibitors.

Proton-Coupled Conformational Activation of SARS Coronavirus Main Proteases and Opportunity for Designing Small-Molecule Broad-Spectrum Targeted Covalent Inhibitors

The SARS coronavirus 2 (SARS-CoV-2) main protease (Mpro) is an attractive broad-spectrum antiviral drug target. Despite the enormous progress in structure elucidation, the Mpro's structure-function relationship remains poorly understood. Recently, a peptidomimetic inhibitor has entered clinical trial; however, small-molecule orally available antiviral drugs have yet to be developed. Intrigued by a long-standing controversy regarding the existence of an inactive state, we explored the proton-coupled dynamics of the Mpros of SARS-CoV-2 and the closely related SARS-CoV using a newly developed continuous constant pH molecular dynamics (MD) method and microsecond fixed-charge all-atom MD simulations. Our data supports a general base mechanism for Mpro's proteolytic function. The simulations revealed that protonation of His172 alters a conserved interaction network that upholds the oxyanion loop, leading to a partial collapse of the conserved S1 pocket, consistent with the first and controversial crystal structure of SARS-CoV Mpro determined at pH 6. Interestingly, a natural flavonoid binds SARS-CoV-2 Mpro in the close proximity to a conserved cysteine (Cys44), which is hyper-reactive according to the CpHMD titration. This finding offers an exciting new opportunity for small-molecule targeted covalent inhibitor design. Our work represents a first step toward the mechanistic understanding of the proton-coupled structure-dynamics-function relationship of CoV Mpros; the proposed strategy of designing small-molecule covalent inhibitors may help accelerate the development of orally available broad-spectrum antiviral drugs to stop the current pandemic and prevent future outbreaks.

Stability and Cell Permeability of Sulfonyl Fluorides in the Design of Lys-Covalent Antagonists of Protein-Protein Interactions

Recently we reported on aryl-fluorosulfates as possible stable and effective electrophiles for the design of lysine covalent, cell permeable antagonists of protein-protein interactions (PPIs). Here we revisit the use of aryl-sulfonyl fluorides as Lys-targeting moieties, incorporating these electrophiles in XIAP (X-linked inhibitor of apoptosis protein) targeting agents. We evaluated stability in buffer and reactivity with Lys311 of XIAP of various aryl-sulfonyl fluorides using biochemical and biophysical approaches, including displacement assays, mass spectrometry, SDS gel electrophoresis, and denaturation thermal shift measurements. To assess whether these modified electrophilic "warheads" can also react with Tyr, we repeated these evaluations with a Lys311Tyr XIAP mutant. Using a direct cellular assay, we could demonstrate that selected agents are cell permeable and interact covalently with their intended target in cell. These results suggest that certain substituted aryl-sulfonyl fluorides can be useful Lys- or Tyr-targeting electrophiles for the design of covalent pharmacological tools or even future therapeutics targeting protein-protein interactions.

Assessment of Tractable Cysteines for Covalent Targeting by Screening Covalent Fragments

Targeted covalent inhibition and the use of irreversible chemical probes are important strategies in chemical biology and drug discovery. To date, the availability and reactivity of cysteine residues amenable for covalent targeting have been evaluated by proteomic and computational tools. Herein, we present a toolbox of fragments containing a 3,5-bis(trifluoromethyl)phenyl core that was equipped with chemically diverse electrophilic warheads showing a range of reactivities. We characterized the library members for their reactivity, aqueous stability and specificity for nucleophilic amino acids. By screening this library against a set of enzymes amenable for covalent inhibition, we showed that this approach experimentally characterized the accessibility and reactivity of targeted cysteines. Interesting covalent fragment hits were obtained for all investigated cysteine-containing enzymes.

Assessment of proton-coupled conformational dynamics of SARS and MERS coronavirus papain-like proteases: Implication for designing broad-spectrum antiviral inhibitors

Broad-spectrum antiviral drugs are urgently needed to stop the Coronavirus Disease 2019 pandemic and prevent future ones. The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is related to the SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), which have caused the previous outbreaks. The papain-like protease (PLpro) is an attractive drug target due to its essential roles in the viral life cycle. As a cysteine protease, PLpro is rich in cysteines and histidines, and their protonation/deprotonation modulates catalysis and conformational plasticity. Here, we report the pK_a calculations and assessment of the proton-coupled conformational dynamics of SARS-CoV-2 in comparison to SARS-CoV and MERS-CoV PLpros using the recently developed graphical processing unit (GPU)-accelerated implicit-solvent continuous constant pH molecular dynamics method with a new asynchronous replica-exchange scheme, which allows computation on a single GPU card. The calculated pK_a's support the catalytic roles of the Cys-His-Asp triad. We also found that several residues can switch protonation states at physiological pH among which is C270/271 located on the flexible blocking loop 2 (BL2) of SARS-CoV-2/CoV PLpro. Simulations revealed that the BL2 can open and close depending on the protonation state of C271/270, consistent with the most recent crystal structure evidence. Interestingly, despite the lack of an analogous cysteine, BL2 in MERS-CoV PLpro is also very flexible, challenging a current hypothesis. These findings are supported by the all-atom fixed-charge simulations and provide a starting point for more detailed studies to assist the structure-based design of broad-spectrum inhibitors against CoV PLpros.

An electrophilic warhead library for mapping the reactivity and accessibility of tractable cysteines in protein kinases

Targeted covalent inhibitors represent a viable strategy to block protein kinases involved in different disease pathologies. Although a number of computational protocols have been published for identifying druggable cysteines, experimental approaches are limited for mapping the reactivity and accessibility of these residues. Here, we present a ligand based approach using a toolbox of fragment-sized molecules with identical scaffold but equipped with diverse covalent warheads. Our library represents a unique opportunity for the efficient integration of warhead-optimization and target-validation into the covalent drug development process. Screening this probe kit against multiple kinases could experimentally characterize the accessibility and reactivity of the targeted cysteines and helped to identify suitable warheads for designed covalent inhibitors. The usefulness of this approach has been confirmed retrospectively on Janus kinase 3 (JAK3). Furthermore, representing a prospective validation, we identified Maternal embryonic leucine zipper kinase (MELK), as a tractable covalent target. Covalently labelling and biochemical inhibition of MELK would suggest an alternative covalent strategy for MELK inhibitor programs.

Avoiding or Co-Opting ATP Inhibition: Overview of Type III, IV, V, and VI Kinase Inhibitors

As described in the previous chapter, most kinase inhibitors that have been developed for use in the clinic act by blocking ATP binding; however, there is growing interest in identifying compounds that target kinase activities and functions without interfering with the conserved features of the ATP-binding site. This chapter will highlight alternative approaches that exploit unique kinase structural features that are being targeted to identify more selective and potent inhibitors. The figure below, adapted from (Sammons et al., Molecular Carcinogenesis 58:1551–1570, 2019), provides a graphical description of the various approaches to manipulate kinase activity. In addition to the type I and II inhibitors, type III kinase inhibitors have been identified to target sites adjacent to the ATP-binding site in the catalytic domain. New information on kinase structure and substrate-binding sites has enabled the identification of type IV kinase inhibitor compounds that target regions outside the catalytic domain. The combination of targeting unique allosteric sites outside the catalytic domain with ATP-targeted compounds has yielded a number of novel bivalent type V kinase inhibitors. Finally, emerging interest in the development of irreversible compounds that form selective covalent interactions with key amino acids involved in kinase functions comprise the class of type VI kinase inhibitors.

Assessment of Proton-Coupled Conformational Dynamics of SARS and MERS Coronavirus Papain-like Proteases: Implication for Designing Broad-Spectrum Antiviral Inhibitors

Broad-spectrum antiviral drugs are urgently needed to stop the COVID-19 pandemic and prevent future ones. The novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is related to SARS-CoV and Middle East respiratory syndrome coronavirus (MERS-CoV), which have caused the previous outbreaks. The papain-like protease (PLpro) is an attractive drug target due to its essential roles in the viral life cycle. As a cysteine protease, PLpro is rich in cysteines and histidines and their protonation/deprotonation modulates catalysis and conformational plasticity. Here we report the pKa calculations and assessment of the proton-coupled conformational dynamics of SARS-CoV-2 in comparison to SARS-CoV and MERS-CoV PLpros using a newly developed GPU-accelerated implicit-solvent continuous constant pH molecular dynamics method with an asynchronous replica-exchange scheme. The calculated pKa's support the catalytic roles of the Cys-His-Asp triad. We also found that several residues can switch protonation states at physiological pH, among which is C270/271 located on the flexible blocking loop 2 (BL2) of SARS-CoV-2/CoV PLpro. Simulations revealed that the BL2 conformational dynamics is coupled to the titration of C271/270, in agreement with the crystal structures of SARS-CoV-2 PLpro. Simulations also revealed that BL2 in MERS-CoV PLpro is very flexible, sampling both open and closed states despite the lack of an analogous cysteine. Our work provides a starting point for more detailed mechanistic studies to assist structure-based design of broad-spectrum inhibitors against CoV PLpros.

Predicting Reactive Cysteines with Implicit-Solvent-Based Continuous Constant pH Molecular Dynamics in Amber

Cysteines existing in the deprotonated thiolate form or having a tendency to become deprotonated are important players in enzymatic and cellular redox functions and frequently exploited in covalent drug design; however, most computational studies assume cysteines as protonated. Thus, developing an efficient tool that can make accurate and reliable predictions of cysteine protonation states is timely needed. We recently implemented a generalized Born (GB) based continuous constant pH molecular dynamics (CpHMD) method in Amber for protein pKa calculations on CPUs and GPUs. Here we benchmark the performance of GB-CpHMD for predictions of cysteine pKa's and reactivities using a data set of 24 proteins with both down- and upshifted cysteine pKa's. We found that 10 ns single-pH or 4 ns replica-exchange CpHMD titrations gave root-mean-square errors of 1.2-1.3 and correlation coefficients of 0.8-0.9 with respect to experiment. The accuracy of predicting thiolates or reactive cysteines at physiological pH with single-pH titrations is 86 or 81% with a precision of 100 or 90%, respectively. This performance well surpasses the traditional structure-based methods, particularly a widely used empirical pKa tool that gives an accuracy less than 50%. We discuss simulation convergence, dependence on starting structures, common determinants of the pKa downshifts and upshifts, and the origin of the discrepancies from the structure-based calculations. Our work suggests that CpHMD titrations can be performed on a desktop computer equipped with a single GPU card to predict cysteine protonation states for a variety of applications, from understanding biological functions to covalent drug design.

Catalytic Site pKa Values of Aspartic, Cysteine, and Serine Proteases: Constant pH MD Simulations

Enzymatic function and activity of proteases is closely controlled by the pH value. The protonation states of titratable residues in the active site react to changes in the pH value, according to their pK_a, and thereby determine the functionality of the enzyme. Knowledge of the titration behavior of these residues is crucial for the development of drugs targeting the active site residues. However, experimental pK_a data are scarce, since the systems’ size and complexity make determination of these pK_a values inherently difficult. In this study, we use single pH constant pH MD simulations as a fast and robust tool to estimate the active site pK_a values of a set of aspartic, cysteine, and serine proteases. We capture characteristic pK_a shifts of the active site residues, which dictate the experimentally determined activity profiles of the respective protease family. We find clear differences of active site pK_a values within the respective families, which closely match the experimentally determined pH preferences of the respective proteases. These shifts are caused by a distinct network of electrostatic interactions characteristic for each protease family. While we find convincing agreement with experimental data for serine and aspartic proteases, we observe clear deficiencies in the description of the titration behavior of cysteines within the constant pH MD framework and highlight opportunities for improvement. Consequently, with this work, we provide a concise set of active site pK_a values of aspartic and serine proteases, which could serve as reference for future theoretical as well as experimental studies.

An Unbiased Immunoaffinity-Based Strategy for Profiling Covalent Drug Targets In Vivo

Activity-based chemical proteomics approaches used for identifying cellular targets of drugs are mainly dependent on the availability of probes derived from drugs. However, all chemical probes are structurally different from the drugs themselves and cannot fully mimic the real actions of drugs in cells. Here we present a concise and unbiased immunoaffinity-based strategy for identifying covalent drug targets in vivo. By using the specific antibody, we not only confirm the well-known ibrutinib-binding target BTK, but also identify some previously undescribed strongly binding proteins, such as CKAP4 in human cell lines and TAP1 in mouse organs. The observed target profiles between species may partially explain why certain drug candidates are very effective in mice but not in humans. This approach avoids the chemical modification of drugs, eliminates the nonspecific bindings of chemical probes, and allows to unbiasedly decode the underlying mechanisms of action of covalent drugs.

Investigation of Crystal Structures in Structure-Based Virtual Screening for Protein Kinase Inhibitors

Protein kinases are important drug targets in several therapeutic areas ,and structure-based virtual screening (SBVS) is an important strategy in discovering lead compounds for kinase targets. However, there are multiple crystal structures available for each target, and determining which one is the most favorable is a key step in molecular docking for SBVS due to the ligand induce-fit effect. This work aimed to find the most desirable crystal structures for molecular docking by a comprehensive analysis of the protein kinase database which covers 190 different kinases from all eight main kinase families. Through an integrated self-docking and cross-docking evaluation, 86 targets were eventually evaluated on a total of 2608 crystal structures. Results showed that molecular docking has great capability in reproducing conformation of crystallized ligands and for each target, the most favorable crystal structure was selected, and the AGC family outperformed the other family targets based on RMSD comparison. In addition, RMSD values, GlideScore, and corresponding bioactivity data were compared and demonstrated certain relationships. This work provides great convenience for researchers to directly select the optimal crystal structure in SBVS-based kinase drug design and further validates the effectiveness of molecular docking in drug discovery.

GPU-Accelerated Implementation of Continuous Constant pH Molecular Dynamics in Amber: pKa Predictions with Single-pH Simulations

We present a GPU implementation of the continuous constant pH molecular dynamics (CpHMD) based on the most recent generalized Born implicit-solvent model in the pmemd engine of the Amber molecular dynamics package. To test the accuracy of the tool for rapid pK_a predictions, a series of 2 ns single-pH simulations were performed for over 120 titratable residues in 10 benchmark proteins that were previously used to test the various continuous CpHMD methods. The calculated pK_a's showed a root-mean-square deviation of 0.80 and correlation coefficient of 0.83 with respect to experiment. Also, 90% of the pK_a's were converged with estimated errors below 0.1 pH units. Surprisingly, this level of accuracy is similar to our previous replica-exchange simulations with 2 ns per replica and an exchange attempt frequency of 2 ps^-1 (Huang, Harris, and Shen J. Chem. Inf. Model. 2018 , 58 , 1372 - 1383 ). Interestingly, for the linked titration sites in two enzymes, although residue-specific protonation state sampling in the single-pH simulations was not converged within 2 ns, the protonation fraction of the linked residues appeared to be largely converged, and the experimental macroscopic pK_a values were reproduced to within 1 pH unit. Comparison with replica-exchange simulations with different exchange attempt frequencies showed that the splitting between the two macroscopic pK_a's is underestimated with frequent exchange attempts such as 2 ps^-1, while single-pH simulations overestimate the splitting. The same trend is seen for the single-pH vs replica-exchange simulations of a hydrogen-bonded aspartyl dyad in a much larger protein. A 2 ns single-pH simulation of a 400-residue protein takes about 1 h on a single NVIDIA GeForce RTX 2080 graphics card, which is over 1000 times faster than a CpHMD run on a single CPU core of a high-performance computing cluster node. Thus, we envision that GPU-accelerated continuous CpHMD may be used in routine pK_a predictions for a variety of applications, from assisting MD simulations with protonation state assignment to offering pH-dependent corrections of binding free energies and identifying reactive hot spots for covalent drug design.

Small-Molecule Inhibitors of Necroptosis: Current Status and Perspectives

Necroptosis, an important form of programmed cell death (PCD), is a highly regulated caspase-independent type of cell death that plays a critical role in the pathophysiology of various inflammatory, infectious, and degenerative diseases. Currently, receptor-interacting protein kinase 1 (RIPK1), RIPK3, and mixed lineage kinase domain-like protein (MLKL) have been widely recognized as critical therapeutic targets of the necroptotic machinery. Targeting RIPK1, RIPK3, and/or MLKL is a promising strategy for necroptosis-related diseases. Following the identification of the first RIPK1 inhibitor Nec-1 in 2005, the antinecroptosis field is attracting increasing research interest from multiple disciplines, including the biological and medicinal chemistry communities. Herein, we will review the functions of necroptosis in human diseases, as well as the related targets and representative small-molecule inhibitors, mainly focusing on research articles published during the past 10 years. Outlooks and perspectives on the associated challenges are also discussed.

Aryl-fluorosulfate-based Lysine Covalent Pan-Inhibitors of Apoptosis Protein (IAP) Antagonists with Cellular Efficacy

We have recently investigated the reactivity of aryl-fluorosulfates as warheads to form covalent adducts with Lys, Tyr, and His residues. However, the rate of reaction of aryl-fluorosulfates seemed relatively slow, putting into question their effectiveness to form covalent adducts in cell. Unlike the previously reported agents that targeted a relatively remote Lys residue with respect to the target's binding site, the current agents were designed to more directly juxtapose an aryl-fluorosulfate with a Lys residue that is located within the binding pocket of the BIR3 domain of X-linked inhibitor of apoptosis protein (XIAP). We found that such new agents can effectively and rapidly form a covalent adduct with XIAP-BIR3 in vitro and in cell, approaching the rate of reaction, cellular permeability, and stability that are similar to what attained by acrylamides when targeting Cys residues. Our studies further validate aryl-fluorosulfates as valuable Lys-targeting electrophiles, for the design of inhibitors of both enzymes and protein-protein interactions.

How Electrostatic Coupling Enables Conformational Plasticity in a Tyrosine Kinase

Protein kinases are important cellular signaling molecules involved in cancer and a multitude of other diseases. It is well-known that inactive kinases display a remarkable conformational plasticity; however, the molecular mechanisms remain poorly understood. Conformational heterogeneity presents an opportunity but also a challenge in kinase drug discovery. The ability to predictively model various conformational states could accelerate selective inhibitor design. Here we performed a proton-coupled molecular dynamics study to explore the conformational landscape of a c-Src kinase. Starting from a completely inactive structure, the simulations captured all major types of conformational states without the use of a target structure, mutation, or bias. The simulations allowed us to test the experimental hypotheses regarding the mechanism of DFG flip, its coupling to the αC-helix movement, and the formation of regulatory spine. Perhaps the most significant finding is how key titratable residues, such as DFG-Asp, αC-Glu, and HRD-Asp, change protonation states dependent on the DFG, αC, and activation loop conformations. Our data offer direct evidence to support a long-standing hypothesis that protonation of Asp favors the DFG-out state and explain why DFG flip is also possible in simulations with deprotonated Asp. The simulations also revealed intermediate states, among which a unique DFG-out/α-C state formed as DFG-Asp is moved into a back pocket forming a salt bridge with catalytic Lys, which can be tested in selective inhibitor design. Our finding of how proton coupling enables the remarkable conformational plasticity may shift the paradigm of computational studies of kinases which assume fixed protonation states. Understanding proton-coupled conformational dynamics may hold a key to further innovation in kinase drug discovery.