Introduction

Ibrutinib (IMBRUVICA®), has revolutionized the treatment of the B cell malignancy, CLL [1]. Ibrutinib is a covalent active site inhibitor of the multi-domain B-cell kinase, BTK (Fig. 1a), that was first approved by the FDA in 2013 [2]. Inhibition of BTK disrupts signaling downstream of the B-cell receptor (BCR), a pathway on which the survival of CLL cells are dependent [3, 4]. The success of Ibrutinib has spurred the development of other BTKi including the now clinically approved BTK inhibitors: Acalabrutinib (CALQUENCE®), Zanubrutinib (BRUKINSA®), Tirabrutinib/ONO-4059 (VELEXBRU®), Pirtobrutinib/LOXO-405 (JAYPIRCA®) and Orelabrutinib (HIBRUKA®) [59]. All of the clinically approved BTKi (with the exception of Pirtobrutinib) are covalent active site inhibitors that bind to BTK C481 residue within the kinase active site. Pirtobrutinib is currently the first and only clinically approved non-covalent BTK active site inhibitor [7]. Importantly, all clinically approved BTK inhibitors in this study, covalent and non-covalent, completely block nucleotide binding.

(a) Domain organization of full-length (FL) BTK and the BTK linker-kinase domain (LKD) fragment used in this study: PHTH, Pleckstrin homology-Tec homology domain; PRR, proline-rich region; SH3, Src homology 3 domain; SH2, Src homology 2 domain, SH2-kinase linker (L) and the catalytic kinase domain. Key residues are indicated above each domain. (b) Autoinhibited conformation of FL BTK based on the crystal structure of FL BTK [40]. The PHTH domain (purple) is dynamic, in transient contact with several regions on the core SH3-SH2-kinase domain and is not visible in the crystal structure of full-length BTK [40]. Dynamics of the PHTH domain is represented by the multiple poses of the PHTH domain and the double headed arrow. (c) Co-crystal structure of BTK LKD (light cyan cartoon) bound to Ibrutinib (PDB ID: 5P9J) showing the location of C481 (yellow spheres), T474 and L528 (red spheres) within the kinase active site (broken oval). (c) Pie charts showing the prevalence of the BTK resistance mutations in CLL patients treated with various BTK inhibitors. The total number of patients with mutations in BTK are indicated below each chart. See Supp. Table 1 for additional details.

Half-life of clinically approved BTK inhibitors: Ibrutinib[64], Acalabrutinib [65], Zanubrutinib [66], Tirabrutinib [67] and Pirtobrutinib [68].

While BTKi are highly effective in the treatment of CLL, they are also being used to treat other B-cell malignancies such as Mantle cell lymphoma (MCL), Waldenström’s macroglobulinemia, Marginal zone lymphoma (MZL) and are being evaluated in the treatment of multiple sclerosis and rheumatoid arthritis [5, 10, 11]. The plethora of clinically approved BTKi along with several other promising candidates currently in clinical trials pose a new challenge with respect to patient treatment. We now need a way to assess the suitability of a given BTKi for a given patient or disease state. Indeed, recent clinical trial data comparing one BTKi to another indicate that BTKi differ in their effectiveness in treating various conditions [1215]. While toxicity, specificity and other criteria are used in these clinical comparisons, a molecular level characterization of the interaction of the inhibitor with full-length BTK is lacking and could be key to understanding these differences. In our previous work, using high resolution structural biology techniques such as Nuclear Magnetic Resonance (NMR) and Hydrogen Deuterium Exchange Mass Spectrometry (HDX-MS), we have shown that some inhibitors (Ibrutinib) exert unanticipated allosteric effects upon binding to BTK while others (Fenebrutinib) do not [16]. Ibrutinib binding to the BTK kinase active site leads to a shift in the conformational ensemble of full-length BTK towards an open conformation where the regulatory SRC homology (SH3 and SH2) domains are released from the ‘back’ surface of the kinase domain, disrupting the closed autoinhibited conformation of the full-length protein (Fig. 1b, [16]). Allosteric conformational changes upon drug binding have been shown to alter protein-ligand interactions in other systems with functional consequences in vivo [1719]. The full impact of the panel of currently approved BTKi on the conformation of full-length BTK is not known.

Despite the success of BTKi in CLL, a recurring problem is the development of BTKi resistance. Resistance typically develops in patients ∼2 years after the start of BTK inhibitor treatment [20]. Analysis of CLL patients that stop responding to BTKi has revealed that most develop mutations within Btk or in the substrate of BTK: Phospholipase C gamma 2 (PLCγ2) [2123]. Interestingly, the specific resistance mutations that develop in Btk seem to be dependent on the specific BTKi used [2329]. While data available so far for the BTKi inhibitors Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib are low in number, there is an emerging trend (Fig. 1c,d). Perhaps not surprisingly, patients treated with the reversible inhibitor Pirtobrutinib do not develop mutations in C481 (Fig. 1d). Instead, they develop other kinase active site mutations including BTK T474I and L528W (Fig. 1c,d) [26]. In contrast, 90% or more of CLL patients treated with Ibrutinib and Acalabrutinib develop mutations in BTK C481 (Fig. 1d, Supp. Table 1 [21, 23]). Substitution of C481 with serine is the most common resistance mutation found in these patients [30]. While mutation of BTK C481 to residues other than serine have also been reported (C481F/Y/R), they occur at a much lower frequency [31]. Additionally, other sites within BTK such as T474 and L528 are rarely (if ever) mutated in Ibrutinib and Acalabrutinib-treated CLL patients (Fig. 1d, Supp. Table 1, [2124, 29, 30, 3234]). Surprisingly, BTK T474 and L528 are frequently mutated in CLL patients treated with the covalent inhibitors Zanubrutinib and Tirabrutinib (Fig. 1d, [24, 27, 28]). BTK T474I and L528W are found at almost equal or higher frequency compared to C481S in these patients (Fig. 1d, Supp. Table 1, [24, 27, 28]). Moreover, in Zanubrutinib-treated patients, the L528W mutation was often present together with C481S mutation (on different alleles), and with the L528W mutation present at a higher allelic frequency compared to C481S [28]. This mutational bias with respect to covalent BTKi is unexpected given the shared mode of action across the panel of covalent inhibitors, and the reasons for these differences are unclear. Additionally, the mechanism of action of the BTK T474I and L528W mutations is not known. Given the development of T474I and L528W mutations in multiple CLL patients treated with both covalent (Zanubrutinib and Tirabrutinib) and non-covalent (Pirtobrutinib) BTK inhibitors, we focused on these BTK mutations and explore their mechanisms of action.

Here, we probe the impact of four clinically-approved BTKi: Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib on the conformation of full-length BTK. We find that each of these BTKi brings about a unique combination of changes in full-length BTK. Acalabrutinib, Zanubrutinib and Tirabrutinib disrupted the autoinhibited conformation of full-length BTK similarly to Ibrutinib but did so to varying degrees. Interestingly, Acalabrutinib and Tirabrutinib altered the dynamics of the kinase G-helix, a region that has been previously characterized as the PLCγ substrate docking site [35]. Pirtobrutinib on the other hand stabilized the compact autoinhibited conformation of full-length BTK and is the first BTK inhibitor observed to do so. Additionally, we probed the mechanism of action of the BTK resistance mutations T474I and L528W. We show that the catalytically inactive BTK L528W mutant activated the SRC kinase HCK and that this activation is dependent on the proline-rich region within BTK. The BTK T474I mutation disrupted binding to Zanubrutinib, Tirabrutinib and Pirtobrutinib and likely evades the action of these drugs due to reduced binding to these inhibitors. Furthermore, we provide an explanation for the mutational bias observed in patients treated with different covalent BTKi. The development of the C481S resistance mutation is dependent on the half-life of the inhibitor which likely explains the low prevalence of the C481S resistance mutation in patients treated with the covalent inhibitors Tirabrutinib and Zanubrutinib (which have a long half-life) as compared to Ibrutinib and Acalabrutinib (which have a short half-life). Characterization of the interaction of BTKi with full-length BTK allows us to better interpret clinical trial results and will help guide the choice of BTKi to be used for treatment. Furthermore, understanding the mechanism of action of resistance mutations allows us to develop treatment strategies that either prevent or delay development of resistance mutations and ways to treat them when they arise.

Results

Assessing the impact of inhibitor binding on the isolated BTK linker kinase domain by NMR

To evaluate the impact of inhibitor binding on BTK, we first monitored its effect on the conformation of the linker kinase domain (LKD) fragment of BTK (Fig. 1a). The catalytic kinase domain of BTK can adopt an active or inactive conformation in solution and previous studies have shown that inhibitor binding can stabilize one or more of these conformations [16]. The switch from an inactive to an active kinase conformation involves changes in key structural elements within the kinase domain. These changes include the inward movement of the C-helix from an inactive ‘αC out’ to the active ‘αC in’ position, changes in the side chain rotamer conformation of BTK W395, formation of a conserved BTK K430:E445 salt bridge and unfurling of the collapsed activation loop, leading to the exposure of the conserved Y551 on the activation loop for phosphorylation (Fig. 2a, [16]).

The BTK kinase domain can interconvert between active and inactive conformations. (a) Superposition of the structure of BTK linker-kinase domain (LKD) bound to Dasatinib (PDB ID: 3K54) in the active kinase conformation (grey cartoon) with the Ibrutinib bound structure (PDB ID: 5P9J) in the inactive conformation (red cartoon). The expanded inset shows the inward movement of the αC-helix, the change in W395 rotamer conformation and the K430/E445 salt bridge formation that accompanies kinase activation. (b-f) Co-crystal structures of BTK LKD (light cyan cartoon) bound to Ibrutinib (PDB ID: 5P9J), Acalabrutinib (PDB ID: 8FD9), Zanubrutinib (PDB ID: 6J6M), Tirabrutinib (PDB ID: 5P9M) and Pirtobrutinib (PDB ID: 8FLL) in the inactive kinase conformation. The inhibitors are shown as dark blue sticks, the kinase activation loop is purple and C481, Y551 and W395 residues are shown as sticks with transparent spheres. Electron density for part of the activation loop is missing in the Ibrutinib co-crystal structure and is indicated as dotted lines (b). In the Acalabrutinib structure, the activation loop has several mutations [36] and the SH2-kinase linker (including W395) is absent (c). Electron density for W395 sidechain is missing in the BTK:Pirtobrutinib co-crystal structure (f). (g) Overlay of the BTK:Ibrutinib, Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib co-crystal structures. With the exception of the Tirabrutinib co-crystal structure (e), no major structural variation is observed in the kinase domains. The activation loop in the Tirabrutinib bound structure adopts a different conformation compared to the other co-crystal structures.

Crystal structures of the BTK kinase domain in complex with Acalabrutinib, Zanubrutinib, Tirabrutinib or Pirtobrutinib are available and show that each of these BTK inhibitors stabilize the kinase domain in an inactive conformation similar to that of Ibrutinib-bound BTK (Fig. 2b-f, [16, 3638]). In all the inhibitor bound structures the C-helix is in the ‘αC out’ position and the activation loop is collapsed into the kinase active site, burying the conserved BTK Y551 (Fig. 2b-f). Moreover, superposition of these BTK/inhibitor complexes shows that the structures of the kinase domain remains largely the same regardless of which inhibitor is bound (Fig. 2g). The structure of the Tirabrutinib bound BTK is the only structure that shows a difference; the activation loop in that complex adopts an alternative inactive loop conformation compared to the other inhibitor bound structures (Fig. 2e,g).

Previous NMR and HDX-MS analysis of the interaction of BTK with a different panel of active site inhibitors has shown that there are differences between the solution and crystal behavior [16]. To test whether the BTK inhibitors Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib can stabilize an inactive BTK kinase conformation in solution as predicted by the crystal structures, we evaluated the BTK inhibitor bound complexes by NMR. The LKD fragment of BTK was isotopically labeled with 15N and a 1H-15N TROSY-HSQC spectrum was obtained in the presence or absence of the inhibitor (Fig. 3). We have previously shown that BTK W395 within the linker region (L) that precedes the kinase domain provides a useful probe to monitor the conformational state adopted by the kinase domain in solution [16, 39]. Assignments for the apo BTK linker kinase domain, which adopts the active (αC in) conformation, show that the W395 indole 1H resonates at 10.21 ppm in the 1H-15N TROSY-HSQC spectrum (Fig. 3, top panel, black spectrum). An upfield shift in the BTK W395 side chain indole NH resonance is observed upon Ibrutinib binding and is consistent with the outward movement of the C-helix (αC out), and stabilization of the inactive kinase domain conformation by Ibrutinib (Fig.3, top panel, cyan spectrum [16]).

BTK inhibitors stabilize the inactive kinase conformation in solution. The tryptophan side chain region of the 1H-15N TROSY HSQC spectra of 15N-labelled apo BTK linker-kinase domain (black spectrum) overlaid with that of the inhibitor bound spectrum (cyan spectrum). Here and in subsequent figures, the broken black and grey lines indicate the position of the BTK W395 resonance in the active (αC-in) and inactive (αC-out) states respectively as shown in Figure 2a. The shift in the BTK W395 resonance upon inhibitor binding is indicated by an arrow in each spectrum. The structures of each inhibitor are shown on the right. The BTK W395 indole NH resonance is in the inactive (αC-out) position in the Ibrutinib (published earlier [16]), Acalabrutinib, Zanubrutinib and Tirabrutinib bound BTK LKD samples. Multiple peaks corresponding to W395 are seen in the Tirabrutinib bound spectrum suggesting that the kinase adopts multiple conformations in solution. The downfield shift observed in W395 in the Pirtobrutinib bound structure is likely due to local changes in the chemical environment due to the distinct chemical structure of Pirtobrutinib. W395 assignments in the inhibitor bound spectra were confirmed by acquiring inhibitor bound spectra with the BTK LKD W395A mutant (see Supp. Fig. S1).

Comparing the tryptophan indole region of the 1H-15N TROSY-HSQC spectra for the BTK LKD inhibitor bound to Acalabrutinib, Zanubrutinib or Tirabrutinib with that of the apo BTK LKD protein shows that W395 undergoes an upfield shift in the presence of inhibitor (Fig. 3), suggesting that, like Ibrutinib [16] and consistent with the crystal structures, all of these inhibitors stabilize the inactive kinase domain conformation in solution. Interestingly, the Tirabrutinib bound spectrum shows several peaks corresponding to W395, suggesting that the Tirabrutinib-bound BTK kinase domain is likely adopting multiple kinase conformations in solution. Additionally, the magnitude of the upfield shift in W395 resonance is smaller in both the Tirabrutinib and Zanubrutinib spectra (compared to Ibrutinib and Acalabrutinib), which may reflect a relatively larger active state population in these samples.

In contrast to the covalent inhibitors, Pirtobrutinib causes a downfield shift in W395, which might suggest that this inhibitor is stabilizing an active conformation of BTK (Fig. 3, bottom panel). However, the fluorinated benzene ring in Pirtobrutinib, which is adjacent to the C-helix and W395 in the Pirtobrutinib bound structure, may cause changes in the local environment of W395 and thereby give rise to the unusual chemical shift change. Indeed, HDX-MS data (see below) confirms the stabilization of the inactive kinase conformation by Pirtobrutinib. All W395 assignments in the WT BTK inhibitor bound spectra were confirmed by comparison to the corresponding 1H-15N TROSY-HSQC spectra of BTK LKD W395A mutant bound to the inhibitors (Supp. Fig. S1).

Probing the effects of inhibitor binding on full-length BTK by HDX-MS

We next evaluated the impact of inhibitor binding on the conformation of full-length BTK by HDX-MS. We have previously shown that HDX-MS can be used to probe conformational changes in full-length BTK that are brought about by inhibitor binding [16]. The structure of the full-length BTK protein in the autoinhibited conformation was recently solved [40]. However, electron density for both the N-terminal PHTH domain and proline rich region is missing in this structure, suggesting that the N-terminal region of BTK is dynamic. Indeed, CryoEM, SAXS and solution data all indicate that the PHTH domain is highly dynamic and likely transiently contacts multiple sites on the core SH3-SH2-kinase domain (Fig. 1b, [40]). Changes in deuterium incorporation must therefore be mapped on the crystal structure of the SH3-SH2-kinase fragment of BTK, which adopts a compact inactive autoinhibited conformation (Fig. 1b, [41]). Additionally, our previous work [16] has shown that for the majority of the inhibitor complexes studied changes in deuterium incorporation are rarely observed in the PHTH domain suggesting that the PHTH domain dynamics are not affected by active site occupancy. Inhibitor binding typically leads to decreased deuterium incorporation in peptides derived from the kinase active site due to stabilization of this region in the presence of the bound inhibitor [16]. Additionally, inhibitors such as Ibrutinib that have allosteric effects lead to an increase in deuterium incorporation in peptides derived from the BTK SH3 and SH2 domains as well as the SH2-kinase linker, indicating a shift away from the autoinhibited conformation [16]. Importantly, HDX-MS provides information about the conformational average adopted by the ensemble and in most cases cannot distinguish between individual intermediate conformations.

Full-length BTK was mixed with Acalabrutinib, Zanubrutinib, Tirabrutinib or Pirtobrutinib and subjected to HDX-MS analysis. Peptides that could be followed in each of the six experimental conditions (apo, and inhibitor bound) were used for comparison (see Supplemental Datafile). Intact mass analysis supports a single binding site on BTK for the covalent inhibitors similar to what has been observed previously with Ibrutinib (Fig. 4a, [16]).

Assessing the impact of BTK inhibitors on full-length BTK by HDX-MS. (a) Intact mass analysis of wild-type FL BTK before (bottom spectrum, black) and after one hour incubation with a 2-fold molar excess of covalent BTK inhibitors: Ibrutinib (red), Acalabrutinib (purple), Zanubrutinib (blue) and Tirabrutinib (orange) show a mass increase of one inhibitor molecule. (b) Clinically approved BTK inhibitors induce allosteric changes in full-length BTK. Relative deuterium level of peptides in apo BTK was subtracted from the deuterium level of the corresponding peptide from each drug-bound form of BTK (DWT drug-bound-DWT apo) and the differences colored according to the scale shown. In this and subsequent figures, peptic peptides are shown from N-to C-terminus, top to bottom, and the amount of time in deuterium is shown left to right. The relative difference data shown here represents a curated set of peptides that are coincident across all 6 states (apo and five drug-bound BTK forms). The identification of these chosen peptides, the relative difference values, and the complete data set for each state can be found in the Supplemental Datafile. The approximate position of the domains of BTK, as described in Figure 1a, is shown at the left. Deuterium incorporation curves of selected peptides (indicated with a gray box in panel b and labelled i-vi) from various regions of the protein are shown below. Data for Ibrutinib has been previously published [16].

Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib show decreased deuterium incorporation in peptides derived from the kinase N-lobe and kinase activation loop similar to that observed previously with Ibrutinib (Fig. 4b,c, Fig. 5a-e, [16]). This is consistent with the binding of these inhibitors in the kinase active site and stabilization of the inactive kinase conformation observed by NMR as well as in the crystal structures. Furthermore, we note that each of these inhibitors induce allosteric effects. Binding of Acalabrutinib, Zanubrutinib and Tirabrutinib to BTK causes increased deuterium incorporation in peptides derived from the SH3 domain suggesting that the autoinhibited conformation of full-length BTK is destabilized (Fig. 4b,c, Fig. 5b-d). This is similar to what has been previously observed with Ibrutinib [16]. Additionally, Acalabrutinib and Tirabrutinib showed increased deuterium incorporation in peptides derived from the G-helix of the kinase domain, a region that has been previously identified as the PLCγ substrate docking site (Fig. 4b,c, Fig. 5b-d, [35]). Taken together, Acalabrutinib, Zanubrutinib and Tirabrutinib binding to BTK leads to hybrid conformations of full-length BTK, where the kinase domain is stabilized in an inactive conformation and the regulatory domains are disrupted from their autoinhibitory conformation. In stark contrast, Pirtobrutinib binding to BTK shows decreased deuterium incorporation in peptides derived from the SH3 and SH2 domains, suggesting that Pirtobrutinib stabilizes the compact, autoinhibited conformation of full-length BTK (Fig. 4b,c, Fig. 5e). To date, Pirtobrutinib is the only inhibitor that we have tested that uniformly stabilizes both the kinase domain as well as the regulatory domains in the inactive autoinhibited conformation. Thus, each BTK inhibitor causes unique changes in the overall conformation of full-length BTK that is not readily predicted from crystal structures. We next turn our attention to the distinct resistance mutations that arise upon treatment with the different BTK inhibitors.

(a-e) Mapping the HDX-MS changes induced by each BTK inhibitor on the structure of the BTK SH3-SH2-kinase fragment (PDB ID: 4XI2). Changes mapped onto the structure represent the maximal HDX-MS change that occurred at any time point. Major difference greater than 1.0 Da are shown as dark blue (decrease) or dark green (increase); modest differences between 0.5 Da and 1.0 Da are shown as light blue (decrease) and light green (increase). Localization of the changes in deuterium incorporation was accomplished using overlapping peptides included in the complete peptide data set provided in the Supplemental Datafile. The location of peptides i – vi from Figure 4 are indicated in panel a. Data corresponding to Ibrutinib has been previously published [16].

Probing the intrinsic effects of the BTK T474I and L528W mutations

BTK resistance mutations can potentially confer a selective advantage to cells in several ways: by increasing the activity of the kinase, by changing the conformation/stability of the kinase which in turn can alter protein-protein interactions, or by disrupting drug binding. To probe the mechanism by which the BTK T474I and L528W mutations confer resistance to Zanubrutinib, Tirabrutinib and Pirtobrutinib we first tested the impact of each mutation on the catalytic activity of the kinase and then investigated the impact on the overall conformation of the protein by HDX-MS and NMR.

To test the catalytic activity of the BTK mutants we set up in vitro kinase assays using purified full-length WT and mutant BTK proteins. Kinase activity was monitored by following phosphorylation on the activation loop Y551 by western blotting. The BTK T474I mutant shows phosphorylation on Y551, however it is lower than that of the WT protein (Fig.6a,b) suggesting that the T474I mutation reduces the activity of the kinase. In contrast, the BTK L528W mutant is completely inactive with no detectable phosphorylation on Y551 throughout the time course (Fig.6c,d). These results are consistent with previous activity reports on both the BTK T474I and L528W mutants [26, 28, 42, 43]. Given that neither the BTK T474I or the L528W mutant had increased catalytic activity, augmenting kinase activity is not the mechanism by which these mutations confer resistance.

Resistance mutations can change the conformation of the protein which in turn can alter protein-protein interactions. To test whether the BTK T474I and L528W mutations alter the overall conformation of the protein, we carried out HDX-MS analysis on the apo mutant proteins. Comparing the BTK T474I mutant to WT BTK, the BTK T474I mutant shows increased deuterium incorporation within the kinase domain N-lobe β2-β3 strands and the activation loop, but no changes elsewhere on the protein (Fig.6e,f,h). These increased dynamics within the N-lobe of the BTK T474I mutant could potentially alter drug binding to the active site. The BTK L528W mutant on the other hand shows a decrease in deuterium incorporation in the kinase domain and the SH3 domain suggesting that the L528W mutation has a stabilizing effect on the full-length autoinhibited conformation of BTK (Fig.6e,g,h).

Additionally, we tested the BTK T474I and L528W mutants by NMR. An overlay of the 1H-15N TROSY-HSQC spectra of apo BTK LKD T474I with that of apo WT BTK LKD shows some changes in the mutant spectrum (Fig.6i). The resonance corresponding to BTK A524 located close (∼18) to T474 in the kinase active site shows line broadening in the BTK T474I mutant spectrum suggesting increased dynamics in the mutant. This is consistent with the increased deuterium incorporation observed in the kinase domain for this mutant by HDX-MS. The BTK W395 resonance in the BTK T474I mutant spectrum overlaps with that of the WT suggesting that the mutant kinase domain adopts an active conformation similar to that of the WT protein. The BTK L528W mutant on the other hand shows a small upfield shift in the W395 resonance, suggesting that the BTK L528W mutant kinase domain is shifted towards an inactive conformation compared to the WT protein. This change in the BTK L528W mutant is consistent with the HDX-MS changes for this protein which suggest that the mutation has a stabilizing effect on the compact autoinhibited conformation. Additionally, a new peak is observed as expected in the BTK L528W mutant spectrum consistent with the introduction of an additional tryptophan due to the mutation. Overall, as the BTK T474I and L528W mutations by themselves cause only minor conformational changes in the nucleotide free apo form, such changes are unlikely to constitute the mechanism by which these mutations confer resistance.

BTK L528W mutant activates the SRC family kinase HCK

Previous studies have shown that the BTK L528W mutant propagates BCR signaling despite being catalytically dead [26, 43]. PLCγ phosphorylation and calcium signaling are maintained in cells carrying this BTK mutation [26, 43]. This suggests that the BTK L528W mutant might recruit other kinases to compensate for its lack of activity. The mechanism of action of another catalytically inactive BTK resistant mutation: BTK C481F/Y, which arises in Ibrutinib treated CLL patients, has been recently reported [44]. In that work the BTK C481F/Y resistance mutant was shown to bind and activate the SRC family kinase HCK, thereby propagating signals from the BCR signaling pathway [44]. This activation of HCK by BTK C481F/Y requires phosphorylation on the BTK kinase activation loop, Y551. The phosphorylated Y551 is suggested to bind to the HCK SH2 domain, displacing the autoinhibited conformation of HCK, thereby activating the HCK catalytic function. To determine if a similar mechanism is at work for the catalytically dead BTK L528W mutant, we tested the ability of the BTK L528W mutant to activate HCK in a western blot assay by monitoring the phosphorylation levels on two substrates: PLCγ1 (pY783 phosphorylation) and HCK itself (pY levels). In addition, we also included BTK WT that had been pre-incubated with Zanubrutinib as a control. As shown in Fig.7a,b, phosphorylation increased on both PLCγ and HCK in the presence of full-length BTK L528W, suggesting that the catalytically dead BTK mutant is able to activate HCK. Interestingly, Zanubrutinib bound WT BTK also promoted increased phosphorylation levels on both PLCγ1 and HCK suggesting that drug-inactivated BTK also activates HCK. This is distinct from what has been previously observed for the BTK C481Y/F resistant mutation where, compared to WT BTK, the C481Y/F mutant BTK preferentially engages HCK. Nevertheless, since the BTK L528W mutant activates HCK, we next designed experiments to probe the requirements for this activation.

Activation of HCK by BTK C481F/Y requires phosphorylation on the BTK kinase activation loop, Y551. Therefore, we first tested whether phosphorylation on Y551 is required for activation of HCK by the BTK L528W mutant. While BTK cannot autophosphorylate itself on BTK Y551 in this assay (BTK FL WT is inhibited by Zanubrutinib and BTK FL L528W is catalytically dead), both proteins are able to be phosphorylated by HCK on BTK Y551. Interestingly, phosphorylation on BTK Y551 is lower in the BTK L528W mutant compared to the WT BTK. This is likely due to the BTK kinase activation loop Y551 being less accessible to phosphorylation by HCK in the mutant (which adopts the autoinhibited conformation) compared to WT BTK (Fig. 7a). Despite the lower pY551 level on BTK FL L528W mutant compared to BTK FL WT, the BTK FL L528W mutant is able to activate HCK similar to the BTK WT protein, suggesting that phosphorylation levels on BTK Y551 do not impact HCK activation. Additionally, we tested whether the BTK L528W/Y551F double mutant is able to activate HCK. The BTK L528W/Y551F double mutant is able to activate HCK similar to BTK L528W, suggesting that phosphorylation on Y551 is not required for HCK activation by BTK L528W (Supp. Fig. S2a). Furthermore, BTK FL L528W increased phosphorylation on HCK and PLCγ1 even in the absence of Zanubrutinib (Fig. S2b) showing that drug interactions are not responsible for the increased phosphorylation on HCK and PLCγ1.

Probing the impact of the BTK resistance mutations T474I and L528W on BTK. (a-d) Western blot comparing the kinase activity of full-length (FL) BTK wild-type (WT), T474I and L528W mutants. BTK autophosphorylation was monitored using the BTK pY551 antibody and the total protein levels monitored using the Anti-His antibody. (b, d) Histogram quantifying the western blots shown in (a and c). The blots were quantified and normalized as described in the Materials and Methods. Data shown are the average of three independent experiments. (e) HDX difference data for the BTK T474I and L528W mutants (DMutant apo -DWT apo). Color scale and peptide/time course arrangement are the same as in Figure 4. See the Supplemental Datafile for additional information, including all peptide identifications and deuterium values. (f,g) Mapping the mutational induced HDX-MS changes on the structure of the BTK SH3-SH2-kinase fragment. Changes mapped onto the structure represent the maximal HDX-MS change that occurred at any time point. (h) Deuterium incorporation curves of selected peptides (indicated with a gray box in panel e and labelled i-iv) from various regions of the protein are shown. (i) The tryptophan side chain region of the 1H-15N TROSY HSQC spectra of 15N-labelled apo WT BTK linker-kinase domain (black spectrum) overlaid with that of the apo mutant kinase spectrum (cyan spectrum). The boxed region in the T474I spectral overlay is expanded above.

The BTK L528W mutant can activate HCK. (a-f) Kinase activity of HCK in the presence or absence of full-length BTK L528W mutant was compared in a western blot assay by monitoring PLCγ1 phosphorylation (pY783 antibody) and HCK autophosphorylation (pY antibody). BTK phosphorylation on Y551 was monitored using the Anti-pY551 antibody. Total protein levels monitored using the Anti-His antibody. Full-length WT BTK preincubated with Zanubrutinib was used as a control. (b, d and f) Histogram quantifying the western blots shown in (a, c and e). The blots were quantified and normalized as described in the Materials and Methods. Data shown are the average of three independent experiments. (b and c) Kinase activity of HCK in the presence or absence of the isolated linker kinase domain (LKD) fragment of the BTK L528W mutant (b) or the full-length proline mutant of BTK L528W (BTK FL L528W/Pro: BTK L528W/ P189A/P192A/P203A/P206A, (c)) was compared as in (a). (g) Thermal stability analysis of BTK FL WT and BTK FL L528W. Data shown are the average of three independent experiments.

We next tested whether the regulatory domains of BTK were required for HCK activation by the BTK L528W mutant by using the isolated linker kinase domain of BTK (a construct in which the BTK regulatory domains have been deleted). As shown in Fig. 7c,d, the isolated linker kinase domain of BTK L528W does not activate HCK. This suggests that the regulatory domain/s of BTK are required for the activation of HCK and that the isolated linker kinase domain of BTK alone is insufficient for HCK activation. This again is different from what has been previously reported; the BTK C481Y/F resistant mutant requires the kinase domain for HCK activation [44]. We therefore turned our focus onto the BTK regulatory domains to isolate the specific region/s required for HCK activation by the BTK L528W resistance mutation.

Displacement of the SH3 domain from its autoinhibitory conformation upon binding of a proline-rich ligand is a classic mechanism by which SRC family kinases such as HCK are activated [45]. BTK contains such a proline-rich ligand sequence within its regulatory region and so we tested HCK activation by a BTK mutant lacking the proline-rich sequence (BTK L528W/Pro: BTK L528W/ P189A/P192A/P203A/P206A). As shown in Fig. 7e,f, the BTK L528W proline-depleted mutant does not activate HCK. This data again emphasizes the difference between the mechanism of action of the BTK L528W mutant from the previously characterized BTK C481Y/F mutant. Additionally, this data is consistent with the observation that WT BTK can also activate HCK, as the proline-rich region is present in both WT and mutant L528W BTK. The ability of both WT and L528W mutant BTK to activate HCK, however, suggests that there must be additional reasons for the selection of the BTK L528W resistance mutation in patients.

BTK L528W is more stable than WT BTK

Changes in protein levels have been shown to cause resistance in several cancers [46]. Cellular protein levels can be altered due to changes in protein expression, stability or degradation, or changes at the RNA level [46]. An increase in BTK protein levels in vivo could drive protein-protein interactions typically not observed in the WT background. Interestingly, we note that protein yield after bacterial expression and purification of the full-length BTK L528W mutant is higher than that of the WT protein, suggesting that the L528W mutant is more stable than the WT protein. To test whether the BTK L528W mutation has a stabilizing effect on BTK, we measured the melting temperature (Tm) of the BTK L528W mutant and compared it to that of the WT protein in a thermal shift assay. As shown in Fig.7g, the BTK L528W mutant is more stable than the WT protein (a one degree increase in Tm compared to the WT). The increased stability of the L528W mutant, albeit small, could result in increased protein levels, an increase in the lifetime of the protein or altered protein trafficking in vivo, all of which would drive the interaction of the mutant protein with HCK unlike the WT protein.

HDX-MS reveals that the BTK T474I and L528W mutants show reduced binding to BTK inhibitors

Disruption of inhibitor binding due to mutations is a classic mechanism by which resistance arises [47, 48]. While loss of inhibitor binding is likely irrelevant to the catalytically dead BTK L528W mutant, we nevertheless evaluated it along with the BTK T474I mutant using HDX-MS and NMR. HDX-MS data show that Tirabrutinib, Zanubrutinib and Pirtobrutinib induce overall changes in deuterium incorporation for both WT BTK and the T474I mutant (Fig. 8a,b,c) These data suggest that the inhibitors are still capable of binding to the T474I mutant. However, in the presence of Tirabrutinib, Zanubrutinib or Pirtobrutinib, the magnitude of the protection in the kinase domain (N-lobe β2-β3 strands and activation loop) is consistently reduced in drug bound T474I mutant as compared to the WT protein (Fig. 8a,b,c,e). The differences in deuterium exchange for drug binding to WT and mutant BTK suggest that the T474I mutation causes a reduction in inhibitor binding. Specifically, the peptides derived from the kinase activation loop in particular do not show as much protection in the mutant relative to the WT suggesting that the inhibitors do not stabilize the activation loop in the inactive conformation upon binding to the T474I mutant. Interestingly, the changes observed in the BTK regulatory domains (SH3 and SH2 domains) upon Zanubrutinib, Tirabrutinib or Pirtobrutinib binding to WT BTK as well as increased deuterium incorporation in the G-helix of the kinase domain upon Tirabrutinib binding to WT BTK are maintained or reduced in the drug bound BTK T474I mutant. In stark contrast to the BTK T474I mutant, the BTK L528W mutant does not show any change in deuterium incorporation in the presence of Zanubrutinib, Tirabrutinib or Pirtobrutinib, providing strong evidence that the BTK L528W mutant does not bind the inhibitors (Fig.8d). Taken together, both the BTK T474I and the BTK L528W mutation impact inhibitor binding. Indeed, these results are consistent with previous studies that show that the BTK T474I and the BTK L528W mutations disrupt binding to Zanubrutinib and Pirtobrutinib [26]. We next tested the impact of the resistance mutations on inhibitor binding by NMR.

HDX-MS analysis of BTK inhibitor binding to BTK T474I and L528W mutants. (a-d) Relative deuterium level of peptides in apo mutant BTK was subtracted from the deuterium level of the corresponding peptide from each inhibitor-bound form of BTK (Ddrug-bound-Dapo) and compared to the changes in the WT protein. The differences are colored according to the scale shown. (e) Deuterium incorporation curves of selected peptides (indicated with a gray box in panels a-c, and labelled i-iv) are shown.

Probing the effects of BTK T474I and L528W mutations on inhibitor binding by NMR

Comparison of the Zanubrutinib bound WT BTK LKD spectrum to that of BTK T474I shows that the BTK T474I mutant undergoes minor changes in the presence of Zanubrutinib, suggesting that the mutant is binding weakly to the inhibitor (Fig. 9a, left and middle panels). The W395 resonance shows little to no chemical shift change upon addition of Zanubrutinib to the BTK T474I mutant, suggesting that this drug does not stabilize the inactive conformation of the mutant kinase (Fig. 9a, middle panel). These results are consistent with HDX-MS changes in the T474I mutant upon Zanubrutinib binding (Fig. 8a).

NMR analysis of BTK inhibitor binding to the BTK T474I and L528W mutants. Th tryptophan side chain region of the 1H-15N TROSY HSQC spectra of 15N-labelled apo BTK linker-kinase domain (black spectrum) overlaid with that of the inhibitor bound spectrum (cyan spectrum). The broken black and grey lines indicate the position of the BTK W395 resonanc and have been described earlier in Figure 3. The shift in the BTK W395 indole NH resonance upon inhibitor binding is indicated by an arrow in each spectrum. The red asterisks indicate the presence of unbound kinase domain in the inhibitor bound NMR sample.

Spectral comparison of Tirabrutinib bound to WT and BTK T474I (Fig. 9b, left and middle panels) shows drug induced chemical shift changes for both complexes, supporting the conclusion from the HDX-MS data (Fig. 8b) that the drug binds to both WT and mutant BTK. However, the spectral changes observed in the spectrum of Tirabrutinib bound to BTK T474I are different from that of the WT/Tirabrutinib bound spectrum (Fig. 9b). Although BTK W395 gives rise to multiple peaks in the spectra of Tirabrutinib bound to both WT and T474I, the W395 resonances in the BTK T474I/Tirabrutinib complex are shifted upfield to a lesser extent compared to the same resonance in the WT/Tirabrutinib spectrum. This suggests that the BTK T474I mutation reduces binding to Tirabrutinib and the inactive kinase domain conformation (αC out) is less populated in the drug bound mutant protein compared to WT BTK bound to Tirabrutinib. These results are consistent with the HDX-MS results that show reduced protection in peptides corresponding to the kinase activation loop in the mutant compared to WT due to the lack of stabilization of the kinase inactive conformation upon Tirabrutinib binding.

Spectral overlay of the Pirtobrutinib bound BTK T474I with the apo T474I protein shows pronounced chemical shift changes in the presence of the inhibitor, providing evidence consistent with the HDX-MS data (Fig. 8c), that Pirtobrutinib binds the BTK T474I mutant (Fig. 9c, middle). However, unlike Pirtobrutinib bound to WT BTK, multiple W395 peaks are visible in the spectrum of Pirtobrutinib bound to BTK T474I. The resonance frequencies of the additional peaks are shifted toward that of the apo protein suggesting the possibility of fast exchange between inhibitor bound and unbound states. These observations suggest that the T474I mutation reduces affinity towards Pirtobrutinib and are consistent with previously published SPR binding data [26].

BTK inhibitors that retain binding to the BTK T474I mutant could serve as potential alternate treatment option for patients with this mutation. Additionally, the lack of or low frequency of the BTK T474I mutation in patients treated with Ibrutinib and Acalabrutinib could be due to the inability of this mutation to disrupt inhibitor binding. We therefore tested the ability of the BTK T474I mutant to bind Ibrutinib, Acalabrutinib and Fenebrutinib (GDC-0853). The NMR data suggest that the BTK T474I mutation does not impact Ibrutinib or Fenebrutinib binding (Fig. 9d,f), but Acalabrutinib binding is reduced (Fig. 9e). Taken together, the T474I mutation significantly reduces binding to Zanubrutinib, Tirabrutinib and Pirtobrutinib, retains binding to Ibrutinib and Fenebrutinib, and may have an adverse effect on Acalabrutinib binding. Given that the BTK T474I mutant retains activity (albeit reduced activity), disruption of drug binding is a possible mechanism by which this mutation escapes inhibition.

Spectral overlays of the BTK L528W mutant with and without Zanubrutinib show no chemical shift changes (Fig. 9a, right panel) suggesting that the mutation completely disrupts inhibitor binding in complete agreement with the HDX-MS data (Fig. 8d). Tirabrutinib does show chemical shift changes, but the changes are markedly different from that observed in the WT protein (Fig. 9b, right). In the L528W spectrum in the presence of Tirabrutinib, W395 shifts in the opposite downfield direction compared to the WT spectrum suggesting that Tirabrutinib could be binding to the BTK L528W in a different orientation. Alternatively, the distinct chemical shift change could be due to the mutated L528W residue altering the local chemical environment within the kinase.

The Pirtobrutinib-bound BTK L528W spectrum (Fig. 9c, Supp. Fig. S1c) shows two resonance positions, one of which overlaps with the resonance in the apo protein and the other that corresponds to that of the mutant protein bound to Pirtobrutinib. This data suggests a mixture of inhibitor bound and unbound BTK kinase domain in solution, likely due to a reduction in Pirtobrutinib affinity caused by the L528W mutation. Although the L528W mutation reduces binding to both Tirabrutinib and Pirtobrutinib, the NMR data suggests that it retains partial binding unlike the HDX-MS data that suggests complete disruption of binding. The higher inhibitor concentrations used in the NMR experiments compared to the HDX-MS experiments likely explain this discrepancy. Interestingly, we note similarities between the BTK L528W Pirtobrutinib bound spectrum and that of Tirabrutinib bound spectrum, suggesting that the BTK kinase domain adopts similar conformations in solution when bound to these different drugs. Additionally, the BTK L528W mutant retains binding to Fenebrutinib, however Ibrutinib and Acalabrutinib binding are disrupted (Fig. 9d-f, right). Taken together, the BTK L528W mutation significantly disrupts Zanubrutinib, Tirabrutinib and Pirtobrutinib binding. However, drugs based on Fenebrutinib could be developed to treat patients carrying this mutation (see discussion below).

Discussion

Binding of BTK active site inhibitors can have long-range effects on the protein. Here, we build on earlier work to show that three of the clinically approved BTKi (Acalabrutinib, Zanubrutinib and Tirabrutinib) shift the conformational ensemble of full-length BTK, destabilizing the autoinhibited conformation of the SH3 and SH2 domains to varying degrees (Fig.4 and 5). In marked contrast, Pirtobrutinib led to the stabilization of the compact autoinhibited conformation of full-length BTK. The exposure or stabilization of the regulatory domains of BTK by active site BTK inhibitors can alter the interaction of BTK with ligands in vivo. Indeed, BTK inhibitors have been shown to differ in their effectiveness at terminating signals downstream of the Fc receptor versus the B cell receptor [49, 50]. As more BTKi become available, a molecular-level understanding of the interaction of the inhibitor with full-length BTK will aid the interpretation of the efficacy of different BTK inhibitors in treating disease states driven by different signaling pathways.

The development of specific resistance mutations in patients treated with different BTKi is intriguing. Covalent BTKi rely on BTK C481 within the kinase active site. The C481S mutation prevents the covalent attachment of these inhibitors to BTK and converts the mode of binding of these drugs to that of a reversible inhibitor. Importantly, we and others have shown that the BTK C481S mutation does not prevent binding of Ibrutinib to the BTK kinase domain [16, 26]. In fact, the BTK C481S mutant binds to Ibrutinib just as well as the WT protein in vitro under equilibrium conditions where drug concentrations do not vary over time [16]. However, in patients, drug concentrations change with time; inhibitor concentrations peak rapidly after intake followed by a decrease as the drug is cleared (represented by the half-life of the drug). The occupancy of BTK C481S protein by these covalent inhibitors in patients is therefore dictated by the availability of these covalent inhibitors over time. Interestingly, Ibrutinib and Acalabrutinib have short half-lives compared to Zanubrutinib and Tirabrutinib (Table 1). The short half-lives of Ibrutinib and Acalabrutinib suggests that the BTK C481S mutant is likely to be unoccupied when inhibitor concentrations fall, thereby allowing the BTK C481S mutant to escape from inhibition. Conversely, the longer half-life of Zanubrutinib and Tirabrutinib in patients suggests that the BTK C481S mutant is likely to remain occupied (inhibited) for a longer period of time. The C481S mutation therefore does not offer a selective advantage as a resistance mutation in Zanubrutinib and Tirabrutinib treated patients and likely explains the lower frequency of occurrence of this mutation in these patients compared to Ibrutinib and Acalabrutinib treated patients. Taken together, the BTK C481S mutant can escape covalent inhibition at low drug concentrations in vivo and is therefore predicted to arise in patients treated with BTK covalent inhibitors that have short half-lives.

The T474I mutation arises frequently in Tirabrutinib and Pirtobrutinib treated CLL patients (Fig. 1d). Our HDX-MS and NMR binding studies show that the T474I mutation disrupts binding to both inhibitors. Reduced binding of BTK T474I mutant to Tirabrutinib and Pirtobrutinib, along with the partial catalytic activity retained by this mutant, would allow BTK T474I mutant to escape inhibition and maintain BCR signaling. This gives a selective advantage to the cells that carry this mutation. Additionally, our NMR binding studies show that the BTK T474I mutation does not disrupt binding of Ibrutinib. This is consistent with the low occurrence of this mutation in Ibrutinib-treated CLL patients. While the BTK T474I mutation does seem to disrupt binding to Acalabrutinib, the lower activity of the T474I mutant (compared to the BTK C481S mutant) may explain the higher prevalence of the BTK C481S mutation (compared to T474I) in Acalabrutinib treated patients.

Analysis of the crystal structures of BTK bound to inhibitors provide clues as to why some BTK inhibitors are more susceptible to loss of binding due to the BTK T474I and L528W mutations. All BTK inhibitors evaluated in this study occupy the kinase active site, contact the kinase hinge region (BTK residues E475-M477) and (with the exception of Fenebrutinib) extend into the back pocket of the kinase (Fig. 10). BTK T474 is adjacent to the kinase hinge region and BTK L528 occupies the ‘base’ of the kinase active site. Given the central location of both these BTK residues, introduction of the bulky T474I and L528W mutations are likely to impact multiple inhibitor interactions. Zanubrutinib and Pirtobrutinib in particular, make extensive contacts with the kinase hinge region with three hydrogen bonds to BTK E475 and M477 backbone atoms. In contrast, Ibrutinib makes only a single hydrogen bond with BTK E475 backbone. This likely explains why the BTK T474I and L528W mutations disrupt Zanubrutinib and Pirtobrutinib binding compared to Ibrutinib. The BTK inhibitor Fenebrutinib also seems to be less impacted by the BTK T474I and L528W mutations. The BTK T474I and L528W mutations are likely accommodated as Fenebrutinib does not occupy the ‘back pocket’ of the kinase active site. While these BTK WT inhibitor bound structures provide clues as to why the efficacy of certain BTK inhibitors are more impacted by the BTK T474I and L528W mutations, some questions remain. Tirabrutinib for example makes a single hydrogen bond with the hinge region, however, it is severely impacted by the BTK T474I and L528W mutations. These mutations likely disrupt additional interactions that Tirabrutinib makes with the kinase such as the glycine loop of the kinase. Crystal structures of inhibitors bound to the BTK mutants are needed to conclusively evaluate their effects.

Close up of BTK inhibitors bound to the BTK kinase active site. The covalent BTK inhibitors (Zanubrutinib, Tirabrutinib, Acalabrutinib and Ibrutinib, represented as sticks) are bound to BTK C481. All BTK inhibitors, covalent and non-covalent (Pirtobrutinib and Fenebrutinib) interact with the hinge region of the kinase and with the exception of Fenebrutinib, extend into the back pocket of the kinase. Hydrogen bonds between the inhibitor and hinge region are shown as dotted lines. The binding of these BTK inhibitors are incompatible with nucleotide binding. The location of BTK T474 and L528, residues that are frequently mutated in CLL patients are shown as sticks and spheres. BTK T474 is adjacent to the hinge region and BTK L528 is at the ‘base’ of the kinase active site.

Ideas and Speculation

Successful treatment of CLL has relied on the inhibition of BTK catalytic activity to curb BCR signaling. Inhibitor-bound, catalytically-inactive BTK is incapable of propagating BCR signaling. Paradoxically, the catalytically inactive BTK L528W mutant is able to propagate BCR signals [26, 43]. A similar catalytically inactive BTK mutant, the BTK C481Y/F resistance mutation also arises in Ibrutinib and Acalabrutinib treated patients, albeit at a low frequency [44]. Since these BTK mutants are catalytically dead, other kinases must compensate for the absence of BTK activity to account for the intact BCR signaling in these patients. Exogenous kinases do not appear to compensate for WT BTK that is inactivated by BTK inhibitors, and so there must be additional differences between the catalytically dead (BTK L528W and BTK C481Y/F) mutant kinases and inhibitor bound WT BTK that is rendered catalytically inactive by virtue of drug binding. Our HDX-MS and NMR data show that BTK L528W mutation does not drastically change the conformation of BTK. However, this mutation does increase the stability of BTK compared to the WT protein. This increased stability could alter protein levels or change the trafficking/localization of the mutant protein in vivo which in turn could alter protein-protein interactions that drive activation of an alternate kinase such as HCK in vivo. Comparison of BTK protein levels in patients pre- and post-development of resistance will be needed to test this hypothesis. Alternatively, the BTK L528W mutation may require an additional change/s in the cell that allow for this dead kinase to recruit other kinases such as HCK to propagate BCR signaling. This requirement for additional changes may explain the low frequency of the L528W mutation in Ibrutinib and Acalabrutinib treated patients; the single amino acid change, C481S, combined with the rapid clearance of drug seems sufficient to drive resistance. Finally, the absence of the BTK L528W resistance mutation in Ibrutinib treated CLL patients could be due to the promiscuity of Ibrutinib; candidate compensatory kinases may be inactivated by the drug [51]. The increased specificity of Zanubrutinib, Tirabrutinib and Pirtobrutinib on the other hand could allow for the utilization of compensatory kinases such as HCK by BTK L528W in order to propagate BCR signaling.

Understanding the mechanism/s by which resistance mutations evade inhibition allows for counter strategies to be devised [5254]. The reversible BTK inhibitor Pirtobrutinib is able to inhibit the BTK C481S mutant due to its long serum half-life and has been suggested as a treatment option for patients that develop this mutation [55]. Alternatively, drug occupancy in the context of the BTK C481S resistance mutant could be increased by administering a twice daily dose of the inhibitor as opposed to the single daily dose possibly circumventing their short half-lives [56]. This increased exposure time to the covalent inhibitor could cause increased side effects and so such risks should be considered before altering the dosage regimen. Alternatively, Tirabrutinib, which exhibits the longest half-life among the available covalent BTK inhibitors, could be used to slow progression of disease in patients that develop the C481S resistance mutation.

The BTK T474I mutation disrupts binding to Zanubrutinib, Tirabrutinib, Pirtobrutinib and Acalabrutinib (Fig.8, 9) but retains binding to Ibrutinib and Fenebrutinib creating possible options for patients that develop the BTK T474I resistance mutation. The BTK L528W mutant is catalytically inactive, hence treatment with any BTK active site inhibitor is futile. Previous studies have shown that the similarly inactive BTK C481Y/F mutant can be suppressed by the use of PROTACS to induce degradation of the full-length protein [44]. Importantly, our data indicate that although the BTK L528W mutation disrupted binding to most BTK inhibitors it retains binding to Fenebrutinib (Fig.9), suggesting that PROTACS based on Fenebrutinib or similar backbones could be developed to treat this resistance mutation [57]. Alternate treatment approaches such as targeting HCK or other BCR signaling proteins such as BCL-2 [58] could also be used to counter this resistance mutation. As additional patient data becomes available, the patterns of resistance mutations for different inhibitors will become clearer and strategies to circumvent resistance using existing inhibitors should improve. As for many areas of medicine, the treatment landscape should benefit from the ongoing era of personalized medicine.

Materials and methods

Constructs and reagents

The bacterial expression constructs for murine BTK linker kinase domain (LKD) and full-length (FL) have been described previously [39]. All BTK constructs carry the solubilizing Y617P mutation for bacterial expression [39]. All mutations were made using the site directed mutagenesis kit (Agilent), and the sequences of all constructs were confirmed by sequencing at the Iowa State University DNA facility. Acalabrutinib, Zanubrutinib, Tirabrutinib and Pirtobrutinib were purchased from MedChem Express. Ibrutinib was purchased from Selleckchem. The pCDF-1 Duet HCK SH3-SH2-Kinase domain/YopH construct with a ‘YEEI’ tail, a kind gift from Dr. Tom Smithgall, was mutated to convert the ‘YEEI’ tail residues to that of the WT protein (YQQQ). The pGEX PLCγ1 cSH2-linker Y771F/Y775F protein substrate construct has been described previously [59].

Protein expression and purification

Expression and purification of the PLCγ1 cSH2-linker Y771F/Y775F protein substrate has been described previously [59]. HCK, Full-length BTK and BTK linker kinase domain (WT and mutants) were produced by co-expressing with YopH in BL21(DE3) (Millipore Sigma) or BL21-Gold(DE3) cells (Agilent Technologies) as described previously [39]. Briefly, the culture was grown at 37 °C to an O.D. 600 nm of 0.6 to 0.8. The temperature of the culture was lowered to 18°C and then induced with either 1.0 mM IPTG for BTK LKD construct or 0.1 mM IPTG for BTK full-length and HCK. The culture was harvested 24 hours after induction and the pellets were resuspended in lysis buffer (50 mM KH2PO4, pH 8.0, 150 mM NaCl, 20 mM imidazole and 0.5 mg/ml lysozyme) and stored at −80°C. Cells were lysed by thawing and the action of lysozyme, and 3000 U DNAse I (Sigma) and 1 mM PMSF were added to the lysate, incubated at RT for 20 minutes and then spun at 16,000 rpm for one hour at 4°C. Glycerol was added to the supernatant to a final concentration of 10 % and was then incubated with Ni-NTA resin (QIAGEN) for two hours, washed with Tris pH 8.0, 75 mM NaCl, 40 mM imidazole, and eluted in 20 mM Tris pH 8.0, 150 mM NaCl, 250 mM Imidazole, and 10% glycerol. Eluted protein was flash frozen in liquid nitrogen and stored at −80°C. The proteins were concentrated and further purified by size exclusion chromatography (Hiload Superdex 26/60 200 pg or Hiload Superdex 26/60 75 pg, GE Healthcare). The fractions containing pure protein were pooled, concentrated, snap frozen and stored at −80° C. The final buffer consists of 20 mM Tris pH 8.0, 150 mM Sodium chloride, 0.02 % Sodium azide and 10% glycerol. Initial phosphorylation levels of all purified BTK and HCK proteins used in this study is below western immuno-detection.

NMR

Uniformly 15N labeled BTK samples were produced as described earlier by growth in modified M9 minimal media containing 15N ammonium chloride (1g/L, Cambridge Isotope Laboratories, Inc.) as the sole source of nitrogen [39]. The final NMR sample buffer consists of 20 mM Tris, 150 mM Sodium chloride, 10% glycerol, and 0.02 % Sodium azide at pH 8.0. All NMR spectra were collected at 298 K on a Bruker AVIII HD 800 spectrometer equipped with a 5 mm HCN z-gradient cryoprobe operating at a 1H frequency of 800.37 NMR samples with inhibitors consisted of 150 µM 15N labeled BTK, mixed with 200 µM inhibitor in 2% DMSO. All data were analyzed using NMRViewJ [60].

HDX-MS

General procedures for HDX-MS of BTK have been described in detail previously [39]. Details specific to experiments conducted here are provided in the Supplemental Datafile in the format recommended [61] for HDX-MS experimental descriptions. All HDX-MS data have been deposited to the ProteomeXchange Consortium via the PRIDE [62] partner repository with the dataset identifier PXD047865. Briefly, prior to continuous labeling HDX experiments, purified BTK full-length wild-type, T474I or L528W (20 µM) and inhibitor (40 µM), (20 mM Tris pH 8.0, 150mM NaCl, 10% glycerol, 2% DMSO) were allowed to interact at 21 °C for 1 hour. After the binding reactions, both the free kinase and kinase bound to inhibitor were placed on ice prior to deuterium labeling. Deuterium labeling proceeded for the times described using labeling buffer, and labeling was stopped with an equal volume of quench buffer at 0°C (details in Supplemental Datafile). Quenched samples were immediately analyzed using a Waters nanoACQUITY with HDX technology using online pepsin digestion with a Waters Enzymate immobilized pepsin column and UPLC separation of the resulting peptic peptides. Mass spectra were acquired using a Waters Synapt HDMSE mass spectrometer. Peptides generated from online pepsin digestion were identified with Waters Protein Lynx Global Server 3.0 using separate unlabeled protein that was prepared in the same manner as protein labeled with deuterium. Deuterium incorporation was quantified using Waters DynamX 3.0. Deuterium levels for each peptide were calculated by subtracting the average mass of the undeuterated control sample from that of the deuterium labeled sample; the data were not corrected for back exchange and are therefore reported as relative [63]. Vertical difference maps in Figures 4, 6 and 8 do not represent a linear sequence of non-overlapping peptides. All coincident and overlapping peptides for comparisons in each figure are provided in figure identified tabs of the Supplemental Datafile.

Activity assays

In vitro kinase assays were performed as described previously [39]. Briefly, 1μM BTK FL, BTK FL T474I or BTK FL L528W proteins were incubated in a kinase assay buffer (50 mM Hepes pH 7.0, 10 mM MgCl2, 1 mM DTT, 5 % glycerol, 1 mM Pefabloc, and 200 μM ATP) at room temperature for varying time. The reactions were stopped by the addition of SDS-PAGE loading buffer and the samples were boiled, separated by SDS−PAGE, and Western blotted with the anti-BTK pY551 antibody (BD PharmingenTM) or anti-His antibody (EMD Millipore) as described previously [39]. The bands were quantified using the ChemiDocTM (Biorad) gel imaging system. The phosphorylation levels (Anti-BTK pY551 blot) were normalized to the total protein level (Anti-His blot). The BTK FL value at 40 minutes was set to 1 and compared to BTK FL T474I or BTK FL L528W. The HCK activation experiments were performed in the absence of Zanubrutinib or by preincubating 0.2 µM of BTK WT, L528W or L528W/Y551F proteins with 0.5 µM Zanubrutinib for 15 min at RT. The reaction was initiated with the addition of 0.2 µM HCK in a kinase assay buffer (50 mM Hepes pH 7.0, 10 mM MgCl2, 1 mM DTT, 5 % glycerol, 1 mM Pefabloc, 30 µM PLCγ1 cSH2 linker Y771F/Y775F substrate and 200 μM ATP). Time points were taken at 5, 10 and 15 minutes. The reactions were stopped by the addition of SDS-PAGE loading buffer and the samples were boiled, separated by SDS−PAGE, and Western blotted with the anti-PLCγ1-pY783 antibody (Cell Signaling) or anti-pY antibody (4G10, Millipore Sigma) or anti-His antibody (EMD Millipore) as described previously [39]. The bands were quantified using the ChemiDocTM (Biorad) gel imaging system. The HCK alone value was set to 1 and compared to HCK + BTK WT/L528W. Initial phosphorylation levels of BTK, prior to the start of the activity assay were undetectable by western immuno-detection.

Thermal shift assays

BTK FL WT or mutant at 5 µM was incubated with 10X final concentration of SYPRO Orange dye (Thermo Fisher Scientific) in a total volume of 20 µL Tris buffer (20 mM Tris, 150 mM Sodium chloride, 10% glycerol, pH 8.0). Thermal shift assays were performed on the StepOnePlus RT-PCR machine in the Iowa State University DNA Facility. The temperature was raised from 25 °C to 95 °C and measurements recorded with every 0.5 °C increment in temperature. The derivative plot of the data was analyzed for the melting temperature (Tm). Assays were performed thrice and the average Tm was calculated.

Abbreviations

  • BTKi: BTK inhibitor

  • LKD: linker-kinase domain

  • CLL: chronic lymphocytic leukemia

  • MCL: Mantle cell lymphoma

  • MZL: Marginal zone lymphoma

  • PRR: proline-rich region

  • HDX-MS: hydrogen/deuterium exchange-mass spectrometry

  • NMR: Nuclear Magnetic Resonance

  • WT: wild-type.

Acknowledgements

This work is supported by a grant from the National Institutes of Health (National Institute of Allergy and Infectious Diseases, AI43957) to A.H.A., J.R.E and T.E.W., and by funds from the Scott-Waudby Trust, the Hope Against Cancer charity, Cancer Research UK in conjunction with the UK Department of Health on an Experimental Cancer Medicine Centre grant [C10604/A25151] to S.J, R.G.B and M.J.SD. Research at the University of Leicester was carried out at the National Institute for Health and Care Research (NIHR) Leicester Biomedical Research Centre (BRC). The authors also thank the Roy J. Carver Charitable Trust, Muscatine, Iowa for ongoing research support.

NMR analysis of inhibitor binding. (a) Cartoon of the BTK linker kinase domain (LKD) fragment used for NMR analysis. (b) Assignment of W395 in inhibitor bound spectra of BTK LKD. The tryptophan side chain region of the 1H-15N TROSY HSQC spectra of 15N-labelled inhibitor bound BTK LKD WT (black spectrum) overlaid with that of the inhibitor bound BTK LKD W395A spectrum (cyan spectrum). The boxed peak indicates the W395 resonance in each of the WT inhibitor bound spectrum. (c) Overlay of the 1H-15N TROSY HSQC spectra of 15N-labelled apo BTK LKD L528W (black spectrum) with that of the Pirtobrutinib bound BTK LKD L528W spectrum (cyan spectrum). The boxed resonances show resonance doubling in the Pirtobrutinib bound spectrum, one of which overlaps with the apo BTK LKD L528W spectrum, indicating a mixture of bound and unbound protein.

BTK FL L528W is able to activate HCK. (a) Phosphorylation on BTK Y551 is not required for HCK activation by BTK FL L528W. Kinase activity of HCK in the presence or absence of full-length BTK L528W single or L528W/Y551F double mutant was compared in a western blot assay by monitoring PLCγ1 phosphorylation (pY783 antibody) and HCK autophosphorylation (pY antibody). Total protein levels monitored using the Anti-His antibody. (b) Zanubrutinib is not required activation of HCK by BTK FL L528W. Kinase activity of HCK in the presence or absence of full-length BTK L528W mutant, without Zanubrutinib wa compared in a western blot assay by monitoring PLCγ1 phosphorylation and HCK autophosphorylation (pY antibody). Total protein levels monitored using the Anti-His antibody.

BTK mutations detected in CLL patients treated with BTK inhibitors. In addition to mutations in BTK C481, T474 and L528, other BTK mutations that have been detected include: R28S, E108K, G164D, V416L, A428D, M437R, R490H, Q516K, V537I and T316A [2124, 2630, 3234].