A road map to evaluate the proteome-wide selectivity of covalent kinase inhibitors.

Kinases are principal components of signal transduction pathways and the focus of intense basic and drug discovery research. Irreversible inhibitors that covalently modify non-catalytic cysteines in kinase active sites have emerged as valuable probes and approved drugs. Many protein classes, however, have functional cysteines, and therefore understanding the proteome-wide selectivity of covalent kinase inhibitors is imperative. Here, we accomplish this objective using activity-based protein profiling coupled with quantitative MS to globally map the targets, both specific and nonspecific, of covalent kinase inhibitors in human cells. Many of the specific off-targets represent nonkinase proteins that, notably, have conserved active site cysteines. We define windows of selectivity for covalent kinase inhibitors and show that, when these windows are exceeded, rampant proteome-wide reactivity and kinase target-independent cell death conjointly occur. Our findings, taken together, provide an experimental road map to illuminate opportunities and surmount challenges for the development of covalent kinase inhibitors.

Protein kinases are one of the largest and most diverse enzyme classes in Nature, with more than 500 members in the human proteome[1]. The size of the kinome, coupled with the central roles that kinases play in cell signaling, physiology, and disease, has inspired widespread effort to create selective kinase inhibitors as basic research probes and therapeutics[2,3].

Most kinase inhibitors described to date are reversible compounds that interact with the conserved ATP-binding pocket of kinases. A number of technology platforms have been developed to profile inhibitor selectivity across the kinome (and, in some cases, the broader ATP-binding proteome) in purified protein and native biological systems[4-6]. The implementation of these platforms, combined with structure-guided medicinal chemistry, has greatly increased the number of kinases for which selective inhibitors have been developed. Despite these advances, the vast majority of human kinases still lack selective inhibitors for investigating their functions in biological systems[3].

A second, emerging class of inhibitors block kinase activity irreversibly by forming covalent bonds with nucleophilic residues, most commonly cysteine, found in the ATP-binding pockets of a substantial fraction of kinases[3,7-9]. Irreversible inhibitors can exhibit advantages over reversible compounds[3,7-11] such as achieving more complete and sustained target engagement in living systems by being less susceptible to competition by high intracellular concentrations of ATP and requiring the physical turnover of kinase proteins to restore inhibited signaling pathways. The cysteine(s) targeted by covalent kinase inhibitors also offer a potential selectivity filter[12], as these residues are not uniformly conserved across the kinome[12,13]. Thus far, covalent inhibitors have been developed for several protein kinases[3], including oncogenic drug targets like the epidermal growth factor receptor (EGFR) and Bruton's tyrosine kinase (BTK), for which the corresponding irreversible inhibitors afatinib (or BIBW-2992)[14] and ibrutinib (or PCI-32675)[15] were recently approved to treat non-small cell lung cancer and chronic lymphocytic leukemia (CLL), respectively. Initial studies with fluorescent or clickable probes suggest that at least some covalent kinase inhibitors can exhibit good selectivity[16-18], although this critical parameter has not, to our knowledge, been thoroughly examined across the proteome for any covalent kinase inhibitor. It has been suggested that, as long as covalent kinase inhibitors are selective within the kinome, specificity across the greater proteome may be inferred[9]. Cysteine, however, is the most intrinsically nucleophilic proteinaceous residue, and many protein classes rely on cysteines for function[19-21]. These factors suggest that covalent kinase inhibitors have the potential to cross-react, either specifically or non-specifically, with proteins outside of the kinome. Such “off-target” activity complicates the assignment of biological functions to kinases in chemical biology experiments and could lead to unanticipated toxicities in drug development programs.

Here, we use activity-based protein profiling (ABPP[22]) combined with quantitative mass spectrometry (MS) to perform a global and in-depth analysis of proteins targeted by covalent kinase inhibitors in human cancer cells. We find that covalent kinase inhibitors, including approved drugs, have defined, but limited concentration windows across which selective target inhibition can be achieved. Once this selectivity range is breached, substantial off-target protein reactivity and kinase target-independent cytotoxicity is observed. Our results thus indicate that medicinal chemistry efforts aimed at optimizing the selectivity of covalent kinase inhibitors should account for their reactivity across the entire human proteome in order to ensure suitable windows of selectivity for basic pharmacology and drug development initiatives. Herein we present an experimental framework for achieving this objective.

Results

Proteomic profiling of covalent kinase inhibitors

For our studies of covalent kinase inhibitors, we selected representative agents that target EGFR (PF-6274484, 1)[23] and BTK (ibrutinib, 2)[18] (). Inhibitors 1 and 2 both use αβ-unsaturated amide electrophiles (hereafter referred to as Michael acceptors) to covalently react with active-site cysteines in the ATP-binding pockets of their target kinases. We synthesized alkynylated variants of inhibitors 1 and 2 (probes 3 and 4, respectively; ), taking care to install the alkyne at positions that would not perturb the intrinsic reactivity of the Michael acceptors nor interfere with target kinase interactions based on established structural-activity relationships[18,23].

Initial ABPP experiments were performed in A431 and Ramos cells, which possess high levels of the primary inhibitor targets EGFR and BTK, respectively. We first evaluated the in vitro and in situ reactivity of probes 3 and 4 by gel-based ABPP, where cell lysates or cells were treated with probes (0.001-10 μM, 1 hr), and then reacted with an azide-rhodamine (N3-Rh) reporter tag under copper-catalyzed azide-alkyne cycloaddition (CuACC or click chemistry) conditions[24] and probe-labeled proteins visualized by SDS-PAGE and in-gel fluorescence scanning[25]. Competitive ABPP experiments, where proteomes or cells were pre-treated with inhibitors 1 and 2 (10 μM, 30 min) facilitated detection of specific targets. In probe 3-treated A431 cells, we observed the concentration-dependent and inhibitor 1-competed labeling of an ~150 kDa membrane protein that presumably represented EGFR ( and Supplementary Results, Supplementary Fig. 1a). The concentration-dependent labeling of several additional proteins was also detected in A431 cells and most of these labeling events were not competed by inhibitor 1. In probe 4-treated Ramos cells, we detected the concentration-dependent and inhibitor 2-competed labeling of an ~70 kDa soluble protein that presumably represented BTK, along with a handful of additional inhibitor 2-competed labeling events visible at low concentrations of probe (≤ 100 nM) and many additional labeling events that appeared at higher concentrations of probe (1 μM or above) and were not competed by inhibitor 2 ( and Supplementary Fig. 1a, b).

Competitive ABPP experiments performed in vitro and in situ produced IC50 values for inhibitors 1 and 2 blockade of EGFR and BTK probe-labeling, respectively, that generally matched those reported for the same[18] or related[23] compounds using substrate assays (, Supplementary Fig. 1c, d, and Supplementary Table 1) and were also similar to estimated half-maximal labeling values for probes 3 and 4 (Supplementary Fig. 1e), indicating that structural modifications to install alkynes did not substantially alter inhibitor potency for the target kinases. Probe 3 could detect EGFR across a range of endogenous expression levels in human cancer cell lines (Supplementary Fig. 1f). Inhibitor 1 exhibited similar in vitro potencies against EGFR in all of the cancer cell lines examined, but consistently showed lower in situ IC50 values in cancer lines expressing mutant (L858R or exon 19 deletion) forms of EGFR compared to cancer lines expressing wild type (WT)-EGFR (Supplementary Fig. 1c and Supplementary Table 1). While the basis for this differential potency is not clear, it suggests that the activation status of EGFR can impact the extent of target engagement by inhibitors in cancer cells.

ABPP studies identified additional inhibitor-competed targets in A431 and Ramos cells, including EGFR as a potent target of inhibitor 2 (a known interaction[18]) (; also see Supplementary Fig. 1g). Other inhibitor-competed proteins were selectively labeled by probe 3 or 4 (Supplementary Fig. 1h; arrows), which mirrored their respective inhibition by 1 and 2 (), albeit with weaker potencies than the target kinases EGFR and BTK. Finally, it is interesting to note that certain inhibitor-competed proteins were detected in situ, but not in vitro (Supplementary Fig. 1g), underscoring the importance of performing competitive ABPP studies in living systems.

Taken together, these gel-based ABPP data indicated that clickable probes can detect “on”- and “off”-targets for covalent kinase inhibitors in cancer cells, where off-targets appear to include both specific (inhibitor-competed and/or probe-selective) and non-specific reactivity events. We next set out to identify the targets of covalent kinase inhibitors by combining ABPP with quantitative mass spectrometry (MS).

Protein targets of covalent kinase inhibitors

We enriched and identified proteins that reacted with probes 3 and 4 by combining ABPP with SILAC (stable isotope labeling by amino acids in cell culture[26]) MS analysis (see Supplementary Fig. 2 and refs. 27 and 28 for more details on the ABPP-SILAC methodology). In our initial ABPP-SILAC studies, we treated cancer cells grown in isotopically ‘light’ or ‘heavy’ media with probe (3 or 4; 1 μM, 1 hr) and DMSO, respectively, which enabled a full inventory of all probe-labeled proteins, defined as proteins that showed light:heavy ratios > 5. Combined analyses of both A431 and Ramos cells treated with probes 3 and 4 identified a total of 29 probe targets ( and Supplementary Tables 2 and 3), including EGFR and BTK, additional known kinase targets of the tested covalent inhibitors – ERBB2[18,23], BLK[18,29], and TEC[18] – and many other non-kinase proteins. Competitive ABPP-SILAC experiments revealed that probe labeling of a subset of these targets, including the aforementioned kinases and a handful of additional proteins, were blocked (> 3-fold) by pre-treatment with inhibitor (1 or 2, 10 μM, 1 hr) ( and Supplementary Tables 2 and 3). An overlapping, but distinct subset of targets showed preferential reactivity (> 5-fold) with probe 3 or 4 in ABPP-SILAC experiments that compared cells treated with equal concentrations of each probe (1 μM; 1 hr; and Supplementary Tables 2 and 3).

Combining data from our ABPP-SILAC experiments, we could categorize 24 of the 29 probe targets as shown in . We first note that the majority of the probe targets exhibited evidence of specificity in terms of being competed by parent inhibitors (“competed”; green, ), selectively reacting with probes 3 or 4 (“selective”; orange, ), or both (“competed and selective”; pink, ). Only a limited number of putatively non-specific targets (not-competed, not-selective; blue, ) were detected in cancer cells at the tested concentrations of covalent kinase inhibitors. Consistent with previous studies[18,23], EGFR and ERBB2 were identified as specific targets of both inhibitor 1/probe 3 and inhibitor 2/probe 4 ( and Supplementary Table 2 and 3). BTK, on the other hand, was primarily competed by inhibitor 2 ( and Supplementary Table 2). Additional protein kinase targets also distributed into the specific categories of competed (BLK) or competed and selective (MAP2K7, TEC, MLTK) ( and Supplementary Table 2 and 3). Most of these kinases are established off-targets that were identified in previous in vitro kinome selectivity screens[18]. Of particular interest, however, was the discovery of several specific nonkinase targets of probes 3 and 4 that included enzymes (e.g., DUS2L, PTGES2, ALDH1A1), receptors (e.g., AHR), and uncharacterized proteins (e.g., FAM213A) ( and Supplementary Table 2 and 3). The non-kinase targets did not share any obvious structural or functional relatedness to kinases or to one another. Intriguingly, however, most of these proteins possess conserved active site or functional cysteines (; Supplementary Table 2 and Supplementary Fig. 3), pointing to a potential unifying basis for their cross-reactivity with covalent kinase inhibitors.

We next recombinantly expressed a subset of off-targets in mammalian cells and assayed their probe reactivity by gel-based ABPP. Consistent with our ABPP-SILAC studies, recombinant DUS2L and FAM213A reacted selectively with probes 3 and 4, respectively (), and these reactions were competed by inhibitors 1 and 2 () with corresponding IC50 values of 1.6 and 1.0 μM (Supplementary Fig. 4a). These inhibition values were similar to those observed for endogenous proteins in cancer cells predicted by molecular mass to represent FAM213A and DUS2L (). Mutation of the conserved catalytic cysteine in DUS2L to alanine (C116A) modestly reduced protein expression and completely eliminated 3-labeling (). FAM213A is an uncharacterized protein that contains two conserved cysteines that form a CXXC motif as part of its predicted thioredoxin-like domain (Supplementary Fig. 3). Mutation of one of these cysteines to alanine (C85A) completely blocked labeling by probe 4 ().

Most of the kinase off-targets identified by ABPP-SILAC share the active-site cysteine residue targeted by covalent kinase inhibitors in EGFR and BTK[3]. One exception is MLTK (or ZAK), which was identified as a specific target of probe 4. We confirmed that probe 4, but not 3 reacted with recombinant MLTK and this reaction depended on a distinct active site-proximal cysteine (C22) and was competed by inhibitor 2 (but not 1) with an IC50 value of 0.7 μM (Supplementary Fig. 4b-d).

These data, taken together, indicate that off-targets of covalent kinase inhibitors originate from structurally and functionally diverse protein families and that these reactions depend on functional cysteines in off-target active sites.

Inhibitor modifications that alter selectivity

We next asked whether modifications to the Michael acceptor might alter covalent kinase inhibitor reactivity. Simple methylation of the terminal Michael acceptor of inhibitor 2/probe 4 produced an inhibitor/probe pair [inhibitor 5 (ref. 30) and probe 6, respectively; ] that still covalently inhibited BTK with good potency in cells (in situ IC50 = 80 nM for inhibitor 5, compared to 2 nM for 2) and exhibited dramatically lower background proteome reactivity as measured by both concentration-dependent ( and Supplementary Fig. 5a) and time-course (Supplementary Fig. 5b) experiments comparing probes 4 and 6. The reduced cross-reactivity observed with inhibitor 5/probe 6 included both specific (e.g., BLK, MAP2K7, MLTK, FAM213A; Supplementary Fig. 5c, d) and non-specific targets (Supplementary Fig. 5d), as determined by gel-based ABPP and ABPP-SILAC studies. Indeed, the only specific off-target still labeled well by 6 in RAMOS cells was TEC kinase (Supplementary Fig. 5d). A direct comparison of the in situ proteome-reactivity profiles of probes 4 and 6 enabled rigorous assignment of suitable concentration windows across which these probes could achieve near-complete inhibition of BTK without substantial non-specific proteome reactivity. For probe 4, this window extended from 10 to 300 nM, while for probe 6, the window spanned from 500 nM to ≥ 10 μM (). These data should help to guide ongoing and future studies that aim to study the biological effects of selective BTK inactivation (see ref. 31 and Discussion below).

Inhibitor modifications that alter in situ reactivity

One common substitution to the Michael acceptor of covalent kinase inhibitors is the dimethylaminomethyl (DMAM) group, which can be found on several covalent kinase inhibitors, including afatinib[3]. The DMAM group is thought to improve drug solubility and, at the same time, maintain potency, possibly by providing a proximal basic amine to activate the attacking cysteine thiol[32].

We synthesized DMAM-substituted analogues of inhibitor 1 and probe 3 (inhibitor 7 and probe 8, respectively; ) and found that 1 and 7 showed similar in vitro and in situ potencies for inhibiting EGFR in cancer cells, as measured by competitive ABPP (Supplementary Table 1). Probes 3 and 8 also showed similar potencies for labeling EGFR in cancer cell lysates, where 3 appeared to exhibit a somewhat greater degree of background proteome reactivity (Supplementary Fig. 6a). In striking contrast to these in vitro profiles, however, probe 8 exhibited much greater time- ( and Supplementary Fig. 6b) and concentration- (Supplementary Fig. 6c) dependent proteome reactivity in living cancer cells. Pre-treatment with inhibitor 7 failed to block most of the extensive proteome reactivity of probe 8 (with EGFR labeling being a notable exception; Supplementary Fig. 6d). ABPP-SILAC experiments confirmed that many more protein targets in cancer cells exhibited greater reactivity with probe 8 compared to probe 3 ( and Supplementary Table 3). The only exceptions were DUS2L and PTGES2, two proteins that we earlier identified as selective targets of probe 3 ( and ) and found here to preferentially react with 3 over 8 (). We also compared the protein labeling profiles of cells treated with 1 versus 10 μM of probe 8 and found that virtually all targets showed substantially higher labeling at 10 μM probe (EGFR and ERBB2 being notable exceptions; Supplementary Fig. 6e and Supplementary Table 3). Taken together, these findings indicate that introduction of a DMAM group onto the inhibitor 1/probe 3 scaffold dramatically increased protein reactivity in cancer cells. Examination of a DMAM analogue of probe 4 (probe 9) revealed that the DMAM modification did not confer higher cellular protein reactivity to the ibrutinib scaffold of covalent kinase inhibitors (Supplementary Fig. 7).

We found that probe 8 displayed much greater stability in cancer cells compared to probe 1 (), providing a plausible explanation for its enhanced in situ proteomic reactivity. Inhibitor 7 also showed more prolonged accumulation in cells compared to inhibitor 1, as did afatinib (10; ) compared to its unsubstituted Michael acceptor analogue (11) (Supplementary Fig. 8a, b), indicating that the DMAM group conferred a general increase in the cellular half-lives of EGFR inhibitors.

Finally, we evaluated the proteome reactivity of probes 3 and 8 in vivo by treating mice with 10 or 20 mg/kg of each compound for 1 hr, after which animals were sacrificed and tissues removed for ABPP. Both probes 3 and 8 labeled an ~150 kDa protein in mouse liver and these labeling events were competed by pre-treatment with afatinib (20 mg/kg, 1 hr) (). We presume that this protein represents EGFR, which is highly expressed in mouse liver (BioGPS http://biogps.org/). Probe 8 also labeled many additional proteins in mouse liver compared to probe 3 ( and Supplementary Fig. 8c). Similar in vivo studies were performed with the BTK probes 4 and 6, where we found that probe 4 showed higher proteomic reactivity than probe 6 in vivo (Supplementary Fig. 8d). These data indicate that the differences in proteomic reactivity displayed by covalent kinase inhibitors in cell culture are also observed in vivo. Importantly, both the EGFR and BTK probes were tested at pharmacologically relevant doses that others have used for structurally related inhibitors (e.g., afatinib and ibrutinib, respectively) to block cancer progression in mouse xenograft models[33,34].

Extensive proteomic reactivity promotes cytotoxicity

A critical question is whether off-target protein reactivity is pharmacologically relevant at concentrations of covalent inhibitors required to produce biological effects through blocking target kinases. The EGFR system provided an excellent opportunity to address this question, as afatinib has been approved for the treatment of exon 19-deletion and L858R EGFR mutant cancers[14] and remains under clinical development for the treatment of T790M EGFR mutant lung tumors (http://clinicaltrials.gov/ct2/show/NCT01647711?term=t790m&rank=1). Past preclinical studies have shown that afatinib displays much greater potency for killing exon 19-deletion and L858R-EGFR cancer cells compared to T790M-EGFR cancer cells (low-nM versus high-nM, respectively)[33,35]. We replicated these findings using the HCC827 (exon 19-del EGFR) and H1975 (L858R/T790M EGFR) cancer lines, where afatinib (and inhibitor 7 and probe 8) showed cytotoxicity IC50 values of < 1 nM and 300-400 nM, respectively ( and Supplementary Fig. 9). We also tested these compounds in three EGFR-independent cancer cell lines (MDA-MB-435S, NIH3T3, SW620), where they displayed cytotoxicity values that were within two to five-fold (0.7-2.2 μM) of those observed in H1975 cells ( and Supplementary Fig. 9). Interestingly, the unreactive, saturated analogues of the EGFR inhibitors (12-14; Supplementary Fig. 9) showed uniformly weaker cytotoxic activity across the EGFR-independent cell lines ( and Supplementary Fig. 9).

We next correlated the cytotoxicity profiles of EGFR inhibitors with their proteome reactivity as measured by gel-based ABPP. In the low-nM range required to kill HCC827 cells, probe 8 showed high selectivity for EGFR (Supplementary Fig. 6c). In contrast, at high nM-low μM concentrations (approximating the IC50 values for killing H1975 and EGFR-independent cancer cells), probe 8 showed extensive concentration-(Supplementary Fig. 6c) and time-dependent (, Supplementary Fig. 6b, and Supplementary Fig. 10) reactivity with the proteomes of all cell lines examined. Negligible time-dependent proteomic reactivity was observed with saturated analogue probe 13 ( and Supplementary Fig. 10). ABPP-SILAC experiments performed with 3 μM of 8 [corresponding to the IC90 value for killing of H1975 cells ()] confirmed the remarkably broad proteomic reactivity of this probe, which canvased > 80 and 400 proteins after 1 and 8 h probe treatments, respectively (Supplementary Fig. 11 and Supplementary Table 3). Pre-treatment with afatinib (3 μM) revealed that the 8-labeled proteins corresponded to a limited number of specific off-targets (e.g., afatinib-competed) and many additional non-specific (afatinib-insensitive) off-targets (Supplementary Fig. 11 and Supplementary Table 3).

These studies, taken together, indicate that covalent EGFR inhibitors bearing a DMAM modification cause cell death in EGFR-independent cancer cells through a mechanism that depends on the reactivity of the inhibitors. Additionally, the potency of this effect approaches the potency at which these compounds kill T790M-EGFR cancer cells (e.g., the IC50 for killing EGFR-independent cancer cells equals the IC70-80 for killing H1975 cells; and Supplementary Fig. 9).

Discussion

Understanding target engagement and off-target interactions is important for all small-molecule probes[36]. These parameters are, somewhat ironically, much easier to measure for irreversible inhibitors than reversibly acting agents, as the very feature that has historically concerned scientists about irreversible probes, namely their covalent binding to proteins[37], also provides a means to detect, enrich, and identify interacting proteins in a proteome-wide manner directly in living systems. We[27,38-40], and others[17,41,42], have used bioorthogonal probes to assess the proteome-wide specificity of irreversible inhibitors in cell and animal models. Nonetheless, these methods have only been applied to a handful of probes and the basis for and pharmacological relevance of off-target activity (specific or non-specific) has gone mostly uninvestigated. Recognizing a growing interest in the development of irreversible inhibitors in both academia and industry[3,7-11], including the recent approval of the first covalent kinase inhibitors afatinib and ibrutinib for treatment of non-small cell lung cancer[14] and CLL[15], respectively, we felt that establishing a more general and rigorous strategy for characterizing the targets of irreversible probes was warranted. Here, we have presented a chemoproteomic platform that meets this goal and applied it to determine the protein reactivity profiles of irreversible kinase inhibitors in both cultured cancer cells and in vivo.

Several important insights were gained from our analysis. First, for most of the covalent kinase inhibitors examined, the off-targets were principally specific in that they showed evidence of either competition by excess parent inhibitor or selective labeling by a subset of the probes under investigation. Specific off-targets were not, however, restricted to kinases, and included proteins from diverse mechanistic classes, including enzymes, receptors, and proteins of uncharacterized function. A unifying feature among these proteins was the presence of active site or functional cysteines, which we confirmed in several instances were required for inhibitor reactivity. These data indicate that proteins bearing functional cysteine residues represent a common source for off-target activity of covalent kinase inhibitors.

An exceptional probe was the DMAM agent 8, which exhibited extensive, time-dependent proteome reactivity in cancer cells. This outcome could be explained, at least in part, by the high and persistent cellular concentrations of probe 8, which converted a selective labeling profile for EGFR observed shortly after treatment to widespread proteome-wide reactivity at later time points (). That the concentrations at which DMAM-modified EGFR inhibitors started to exhibit substantial proteomic reactivity (≥ 0.1 μM) correlated with the emergence of cytotoxic activity in EGFR-independent cancer cells (IC50 values of 0.7-2.2 μM) leads us to conclude that off-target protein modification is a principal driver of the EGFR-independent cell death effects of these compounds. Also consistent with this interpretation is the uniformly greater cytotoxicity displayed by probe 8, inhibitor 2, and afatinib compared to their saturated analogues. We suspect that the cytotoxic effects of DMAM-modified EGFR inhibitors are caused by reactivity with many proteins in cells, rather than being due to the specific inactivation of a single protein for multiple reasons: 1) the tested DMAM-modified EGFR inhibitors killed all three EGFR-independent cancer cell lines with similar potencies; and 2) a very large number of time-dependent protein labeling events were observed (> 400) at inhibitor concentrations that killed these cancer cells.

While the DMAM-modified inhibitors killed exon 19-del-EGFR cancer cells at concentrations where selective EGFR labeling was observed, the concentration of these compounds required to kill T790M-EGFR cancer cells traversed into the range where extensive off-target reactivity and EGFR-independent cytotoxicity was observed ( and Supplementary Fig. 6, 9, and 10). Our results thus highlight potential challenges for current efforts to expand the clinical utility of afatinib beyond exon 19-deletion and L858R-EGFR tumors to include T790M-EGFR tumors. We should note that the mechanistic basis for the reduced activity of DMAM-modified inhibitors against T790MEGFR cancer cells remains unclear, with some[33,43], but not all studies[44] suggesting reduced potency for this mutant form of the kinase. Regardless, the effect appears to be intrinsic to the T790M-EGFR protein, since introduction of this mutant kinase into cancer cells is sufficient to reduce their sensitivity to afatinib[35].

The reduced proteomic reactivity of probe 6, which still maintained good activity against its target kinase BTK, indicates that it may be possible, for exceptionally potent covalent inhibitors, to trade some potency in exchange for improved selectivity. This could be particularly important for tool compounds, where a defined, but restricted window of selectivity, as we observed for inhibitor 2/probe 4, may complicate biological experiments. Indeed, we have found many papers in the literature where 2 (ibrutinib) was applied to cells at ≥ 5 μM concentrations and its effects interpreted to occur through BTK[45-48]. Based on our chemoproteomic profiles, it would be difficult to draw such conclusions with confidence given the large number of proteins that react with probe 4 in this concentration range. We should also note, however, that reductions in covalent reactivity may not correlate with improvements in selectivity for reversible binding to off-targets. As we have shown for EGFR inhibitors, unreactive, saturated analogues of Michael acceptor-containing kinase inhibitors can serve as useful control probes for discerning pharmacological effects that are mediated by covalent versus non-covalent mechanisms. Finally, from a technical perspective, our finding that modifications to the Michael acceptor can dramatically alter proteome reactivity indicates that efforts to create chemoproteomic probes that authentically report on inhibitor selectivity should avoid attaching reporter tags directly to the inhibitor's reactive group (as has been done previously for 2[18]).

Projecting forward, we believe that the chemoproteomic platform described herein should be applicable to any irreversible inhibitor, assuming that it can be modified by an alkyne without substantial alterations in structure or activity. We elected to characterize advanced kinase inhibitors in this study, but chemoproteomics can also facilitate the discovery of kinases that are targeted by covalent inhibitors in phenotypic screens[41]. As a greater number of covalent kinase inhibitors are examined using chemoproteomics, it will be interesting to determine whether a finite set of specific off-targets is identified. Considering the limited overlap in specific off-targets observed for 1 and 2 in this study, combined with the large and diverse number of proteins that require cysteine residues for function[19-21], we suspect that many additional cysteine-dependent proteins in the human proteome may prove sensitive to covalent kinase inhibitors. Finally, we are emboldened by the discovery that irreversible kinase inhibitors targeting cysteine can have, in general, only a limited number of off-target interactions in cells and these off-targets appear to be confined to proteins that themselves possess active-site or functional cysteines. That the off-targets also exhibit a clear structure-activity relationship indicates achieving selectivity for covalent kinase inhibitors, or for that matter, any compound that acts via cysteine adduction, should have a medicinal chemistry solution, as long as these efforts are integrated with chemoproteomics as a guiding platform.

Online Methods

Cell culture

Cells were grown at 37 °C under a humidified 5% CO2 atmosphere, in a culture medium consisting of high-glucose DMEM (Caisson Labs) for A431, HEK293T, NIH-3T3 and HeLa cells or RPMI (Caisson Labs) for Ramos, HCC827, and A549 cells. The SW620 and MDA-MB-435S cells were grown in L-15 media in a humidified atmosphere with no CO2 added. All media was supplemented with 10% FBS (Gemini), penicillin (50 IU/mL), streptomycin (50 μg/mL), and glutamine (2 mM) (Cellgro). For SILAC experiments, the culture medium was replaced with either SILAC DMEM or SILAC RPMI (Thermo) as appropriate, and the medium was supplemented with 10% dialyzed FBS, penicillin, streptomycin, and glutamine. For the isotopically heavy cell samples, 100 μg/mL of both [13C6,15N4]L-arginine-HCl and [13C6,15N2]L-lysine-HCl (Sigma-Aldrich) was also added to the culture medium. For the isotopically light cell samples, the culture medium was supplemented with 100 μg/mL of both L-arginine-HCl and L-lysine-HCl. Cells were passaged at least six times in isotope-containing medium before being used for ABPP-SILAC experiments.

In situ cell treatment, gel-based ABPP

Adherent cell lines were grown to ~ 90% confluence in 10 mL of growth media in 10 cm plates. The growth media was aspirated off, and the cells were washed 2X with DPBS (~ 10 mL), followed by the addition of fresh media (10 mL). The inhibitor or probe (10 μL of 1000X DMSO stock) was then added, the plate was swirled to disperse the compound and then incubated at 37 °C for 1 hr. For IC50 measurements, a 1 hr inhibitor pre-incubation was followed by a 1 hr probe chase. Following incubation, the media was removed and the cells were washed 3X with 10 mL DPBS. The cells were then harvested and stored as pellets at − 80 °C. For suspension cell lines (Ramos cells), the cells were grown to ~ 3 × 106 cells/mL in large culture flasks and then were aliquoted (10 mL) into 10 cm plates for treatment. The cells were then treated and stored as described above for the adherent cell lines.

In situ cell treatment, ABPP-SILAC

For adherent cell lines, the same procedure as above was followed, except that 10 cm plates were used and the treatment volume was increased to 10 mL. A total of 4 plates were treated for each biological replicate. For Ramos cells, the cells were grown as described above and then were aliquoted (15 mL) into 100 mL tissue culture flasks for treatment. The cells were then treated and harvested as described above.

Sample processing for in situ-labeled samples, gel-based ABPP

Cell pellets were thawed on ice and cold DPBS (200 μL) plus protease inhibitor cocktail (Roche) was added to each sample. The cells were lysed using a probe sonicator and the lysate was fractionated by centrifugation at 100,000 × g for 45 min. The supernatant (soluble proteome) was transferred to a separate microfuge tube and the pellet (membrane proteome) was resuspended in 200 μL cold DPBS (plus protease inhibitors) using a probe sonicator. Protein concentrations for each sample were determined using the DC protein assay (BioRad) and then nomalized to 1 mg/mL in a volume of 43 μL. Click chemistry was performed on each sample using final concentrations of 25 μM rhodamine azide, 1 mM Tris(2-carboxyethyl)phosphine (TCEP, Sigma Aldrich), 100 μM Tris[(1-benzyl-1H-1,2,3-triazol-4-yl)methyl]amine (TBTA, Sigma Aldrich) and 1 mM CuSO4 in a final volume of 50 μL. The samples were incubated at room temperature for 1 hr, followed by addition of 20 μL 4X SDS loading buffer. The samples were then loaded (30 μL) and resolved on a 10% SDS-PAGE gel and imaged using a Hitachi FMBIO-II flatbed fluorescence scanner. Fluorescence images are shown in grayscale.

In vitro cell lysate labeling, gel-based ABPP

Soluble and membrane proteome samples from untreated cells were obtained by lysis and ultracentrifuge fractionation as described above and were normalized to a concentration of 1 mg/mL in a volume of 49 μL. Samples were then treated with inhibitor or probe (1 uL 50X stock in DMSO), mixed by agitation, and incubated at 37 °C for 1 hr. For IC50 measurements, a 1 hr inhibitor treatment was followed by a 1 hr probe treatment. Click chemistry was then set up (reagents and concentrations same as described above) in a total volume of 57 μL and the samples were incubated at room temperature for 1 hr. The samples were treated with 20 μL 4X SDS-PAGE loading buffer and then loaded (30 μL) and resolved on a 10% SDS-PAGE gel. The gel was then imaged as described above.

Determination of IC50 values by competitive gel-based ABPP

Determination of IC50 values from gel-based competitive ABPP experiments was performed as previously described[49].

Recombinant expression and labeling of proteins in HEK293T cells

The cDNA constructs for DUS2L, MAP2K7, and FAM213A were obtained from Origene in the pCMV6-entry vector, with expression of C-terminal FLAG tags. The cDNA construct for ALDH1A1 was obtained from Thermo Fisher and cloned into the pFLAG-CMV2 vector, with expression of an N-terminal FLAG tag. For MLTK, a HEK293T cell line with stable overexpression of an MLTK construct containing an N-terminal FLAG tag was used. Overexpression of the other proteins in HEK293T cells was obtained by transient transfection using the FuGENE HD transfection system according to the manufacturers protocol. Forty-eight hours after transfection, in situ gel-based ABPP was then performed as detailed above. Western blotting was accomplished as previously described[50], using an anti-FLAG primary antibody, IRdye 680 (Licor Biosciences) secondary antibody and Licor Odyssey CLx near-infrared imager.

Sample processing for analysis by LC-MS/MS

Sample preparation for LC-MS/MS was performed as previously described[51], with minor modifications. In brief, light and heavy cell pellets were processed as described above for in-situ labeled gel-based ABPP samples to obtain the membrane and soluble proteomes of both light and heavy samples. All samples were normalized to 2 mg/mL in 430 μL of DPBS and click chemistry was performed with final concentrations of 200 μM biotin-azide (Chempep Inc.), 1 mM TCEP, 100 μM TBTA, and 1 mM CuSO4 in a final volume of 500 μL. The reaction was allowed to proceed at room temperature for 1 hr, and the excess reagents were then removed by CHCl3/MeOH precipitation with 2 mL MeOH, 1.5 mL H2O, and 0.5 mL CHCl3. Following isolation, the heavy and light proteome pellets were then combined into 600 μL of MeOH and the sample was precipitated again by addition of 600 μL of H2O and 150 μL of CHCl3. The pellet obtained was air dried and then resuspended in 500 μL of aqueous 6M urea, 25 mM NH4HCO3. The peptides were reduced (DTT) and alkylated at cysteine (iodoacetamide) and then diluted to 6 mL with DPBS and enriched over avidin agarose (Sigma; 100 μL slurry) for 4 hrs at room temperature. The avidin beads were washed and an overnight on-bead trypsin digestion was performed as previously described[51]. The final trypsin digest sample (300 μL) was acidified with 16 μL of formic acid and pressure loaded onto a biphasic (strong cation exchange/reverse phase) capillary column for analysis by two-dimensional LC-MS/MS.

Mass spectrometry and data analysis

Mass spectrometry was performed using a Thermo Orbitrap Velos mass spectrometer, using a previously described protocol[51,52]. Peptides were eluted using a five-step multidimensional LC-MS (MudPIT[53]) protocol (using 0%, 25%, 50%, 80% and 100% salt bumps of 500 mM aqueous ammonium acetate, followed by an increasing gradient of aqueous acetonitrile and 0.1% formic acid in each step), and data were collected in data-dependent acquisition mode (400–1800 mass to charge ratio (m/z) and 30 data-dependent fragmentation (MS2) scans) with dynamic exclusion enabled (repeat count of 1, exclusion duration of 20 s). Prolucid searches allowed for variable oxidation of methionine (+15.9949 m/z), static modification of cysteine residues (+57.0215 m/z; iodoacetamide alkylation) and one end specificity for trypsin. Each data set was independently searched with light and heavy parameter files; for the light search, all other amino acids were left at default masses; for the heavy search, static modifications on lysine (+8.0142 m/z) and arginine (+10.0082 m/z) were specified. The precursor-ion mass tolerance was set to 50 ppm and the fragment-ion mass tolerance was the default assignment of 600 ppm. The data were searched using a human reverse-concatenated nonredundant (gene-centric) FASTA database that was assembled from the Uniprot database (www.uniprot.org). The resulting matched MS2 spectra were assembled into protein identifications, then filtered using DTASelect (version 2.0.47), and only half-tryptic or fully-tryptic peptides were accepted for identification, and only fully-tryptic peptides were considered for quantification. Peptide false positive rates for all datasets were < 0.3%. Redundant peptide identifications common between multiple proteins were allowed, but the database was restricted to a single consensus splice variant. SILAC ratios were quantified using in-house software as described (CIMAGE[19]). Briefly, extracted MS1 ion chromatograms (± 10 ppm) from both 'light' and 'heavy' target peptide masses (m/z) are generated using a retention time window (± 10 min) centered on the time when the peptide ion was selected for MS/MS fragmentation, and subsequently identified. Next, the ratio of the peak areas under the light and heavy signals (signal-to-noise ratio > 2.5) are calculated. Computational filters used to ensure that the correct peak-pair is used for quantification include a co-elution correlation score filter (R2 ≥ 0.8), removing target peptides with bad co-elution profile, and an 'envelope correlation score' filter (R2 > 0.8) that eliminates target peptides whose predicted pattern of the isotopic envelope distribution does not match the experimentally observed high-resolution MS1 spectrum. Also, peptides detected as singletons, where only the heavy or light isotopically labeled peptide was detected and sequenced, but which passed all other filtering parameters, were given a standard ratio of 20, which is the maximum SILAC ratio reported here.

Metabolite analysis

A431 cells were seeded at 2 × 106 cells/60 mm dish and cultured in complete DMEM medium (4 mL) for 24 h. The media was removed and the cells were washed once with 3 mL of DPBS. Fresh media (3 mL) was added to the plates followed by treatment with probe 3 or probe 6 (1000X stock in DMSO, 3 μL). Plates were incubated at 37 °C until time point, after which the media was aspirated and the plates were washed 5 x 5 mL of DPBS and the cells were harvested and the cell pellet collected. H2O-soluble cellular metabolites were extracted using a MeOH/H2O extraction protocol, essentially as previously described[3]. In brief, cell pellets were resuspended in 200 μL of an 80:20 mixture of MeOH/H2O and an internal deuterated standard, 1mM of deuterated serine-d3 (Cambridge Isotope Laboratories), was added to the extraction solution for absolute quantification and sample normalization. The mixture was sonicated for 2 x 8 pulses followed by a 5 min centrifugation at 14,000g. The supernatant was collected and stored at −80 °C until LCMS analysis.

Samples of metabolites from both probe treatments (3 replicates/probe) were run sequentially. H2O-soluble cellular metabolites were separated by reverse phase chromatography (Phenomenex, C18) using elution conditions as described previously[1]. MS analysis was performed on an Agilent G6410B tandem mass spectrometer using ESI source setting as described previously[1]. The mass spectrometer was run in MRM mode, monitoring the transition of m/z from 454.1 to 343.1 for probe 6 , 397.1 to 343.1 for probe 3, and 109.1 to 63.2 for serine-d3. Metabolites were quantified by measuring the area under the peak in comparison to an external calibration curve. Ion abundances were normalized against serine-d3.

Cytotoxicity Assays

All cell lines were plated at 5,000 cells/well in 96-well plates in100 μL growth media (except SW620, which was plated at 10,000 cells/well). The cells were then treated with compounds or DMSO (1% DMSO final) and incubated at 37 °C for 72 hrs. Cell density was then measured using the Sulforhodamine B assay[54]. In brief, the cells were fixed by addition of 100 μL of cold 10% aqueous trichloroacetic acid and incubation for 1 hr at 4 °C. The TCA solution was decanted and the plates were washed gently with H2O (3X) and incubated at room temperature until dry. The plates were then stained by addition of 100 μL of 0.057% (wt/vol) Sulforhodamine B solution and incubation at room temperature for 30 min. The staining solution was decanted and the plates were washed gently with 1% acetic acid (3X) and then allowed to dry at room temperature. The plates were then treated with 200 μL 10 mM Tris (pH 10.5) solution and incubated at room temperature for 30 minutes followed by OD measurements at 535 nM using a Tecan Infinite F500 instrument.

In vivo studies

C57BL/6 mice were injected (i.p.) with compounds using a vehicle composed of 17:1:1:1 (v/v/v/v) solution of saline/EtOH/DMSO/Cremophor EL (10 μL g−1). Mice were treated with inhibitors (20 mg/kg) or vehicle for 1 hr, followed by treatment with probes (10 or 20 mg/kg) for 1 hr. Mice were then anesthetized with isoflurane and euthanized by cervical dislocation. Tissues were harvested and washed with cold DPBS (3X) followed by homogenization using a NextAdvance bullet blender and Zirconium Oxide beads. Samples were then centrifuged at 1,500 × g for 5 minutes and supernatant was removed and sonicated using a probe sonicator to ensure complete cell lysis. The samples were centrifuged at 100,000 × g for 45 min to provide the membrane proteome as a pellet and cytosolic proteome as supernatant. The pellet was washed with cold DPBS and resuspended by sonication. Protein concentrations were normalized to 2 mg/mL and analyzed by gel-based ABPP as described above. All mouse studies were performed following protocols that received approval from the TSRIIACUC office.