Structural Characterization of Proline-rich Tyrosine Kinase 2 (PYK2) Reveals a Unique (DFG-out) Conformation and Enables Inhibitor Design

Proline-rich tyrosine kinase 2 (PYK2) is a cytoplasmic, non-receptor tyrosine kinase implicated in multiple signaling pathways. It is a negative regulator of osteogenesis and considered a viable drug target for osteoporosis treatment. The high-resolution structures of the human PYK2 kinase domain with different inhibitor complexes establish the conventional bilobal kinase architecture and show the conformational variability of the DFG loop. The basis for the lack of selectivity for the classical kinase inhibitor, PF-431396, within the FAK family is explained by our structural analyses. Importantly, the novel DFG-out conformation with two diarylurea inhibitors (BIRB796, PF-4618433) reveals a distinct subclass of non-receptor tyrosine kinases identifiable by the gatekeeper Met-502 and the unique hinge loop conformation of Leu-504. This is the first example of a leucine residue in the hinge loop that blocks the ATP binding site in the DFG-out conformation. Our structural, biophysical, and pharmacological studies suggest that the unique features of the DFG motif, including Leu-504 hinge-loop variability, can be exploited for the development of selective protein kinase inhibitors.

Proline-rich tyrosine kinase 2 (PYK2) is a cytoplasmic, nonreceptor tyrosine kinase implicated in multiple signaling pathways. It is a negative regulator of osteogenesis and considered a viable drug target for osteoporosis treatment. The high-resolution structures of the human PYK2 kinase domain with different inhibitor complexes establish the conventional bilobal kinase architecture and show the conformational variability of the DFG loop. The basis for the lack of selectivity for the classical kinase inhibitor, PF-431396, within the FAK family is explained by our structural analyses. Importantly, the novel DFG-out conformation with two diarylurea inhibitors (BIRB796, PF-4618433) reveals a distinct subclass of non-receptor tyrosine kinases identifiable by the gatekeeper Met-502 and the unique hinge loop conformation of Leu-504. This is the first example of a leucine residue in the hinge loop that blocks the ATP binding site in the DFG-out conformation. Our structural, biophysical, and pharmacological studies suggest that the unique features of the DFG motif, including Leu-504 hinge-loop variability, can be exploited for the development of selective protein kinase inhibitors.
Proline-rich tyrosine kinase 2 (PYK2) 2 and focal adhesion kinase (FAK) comprise the focal adhesion kinase subfamily of non-receptor tyrosine kinases. PYK2 and FAK are large multidomain proteins containing an N-terminal FERM domain, a central catalytic domain, and a C-terminal segment containing dual proline rich (PR) subdomains and a focal adhesion targeting (FAT) region (1,2). While FAK is widely expressed, PYK2 expression is relatively restricted with highest levels in brain and the hematopoeitic system. Unlike FAK, optimal PYK2 activation is dependent on Ca 2ϩ mobilization. PYK2 (Ϫ/Ϫ) animals have been described previously, and develop normally (3,4). Characterization of the immune system of PYK2(Ϫ/Ϫ) animals revealed the absence of marginal zone B-cells along with abnormal T-cell independent type II responses (4), and altered macrophage morphology, migration and signaling in response to cell attachment or chemokine treatment (3). These studies strengthen the link between PYK2 and signaling through chemokine and integrin receptors. In addition, PYK2(Ϫ/Ϫ) mice were shown to have increased susceptibility to diet-induced obesity and diabetes (5).
Recently, the characterization of PYK2(Ϫ/Ϫ) mice showed a high bone mass phenotype resulting from increased osteogenesis and osteoblast activity. Using PYK2(Ϫ/Ϫ) mouse bone marrow cultures and hMSCs expressing a PYK2 shRNA, elimination or reduction of PYK2 protein levels resulted in significantly enhanced osteogeogenesis. Importantly, the daily administration of a pyrimidine-based PYK2 inhibitor, PF-431396, increased bone formation, and protected against bone loss in ovariectomized rats (6). PYK2(Ϫ/Ϫ) mice showed mild osteopetrosis which was attributed to the impairment in osteoclast function (7). Therefore, the high bone mass phenotype may result from both enhanced osteoblast and impaired osteoclast elements.
PYK2 is one member of a family of over 500 evolutionarily conserved enzymes with high amino acid and structural conservation within the catalytic ATP binding pocket. Classical kinase inhibitors bind to the ATP site and compete for substrate binding. Thus, while classical inhibitors based on ATP binding analogs have been readily identified, the inherent promiscuity of action for this class has presented significant challenges to drug design (8). With the exception of cancer therapeutics, where additional therapeutic benefits may be gained by the inhibition of multiple kinase targets (e.g. Sutent, Sorafenib), minimizing off-target activity is most often desired. Therefore, there is great interest in identifying unique allosteric regulatory domains for specific kinase targets. Despite intense effort, small molecule inhibitors exploiting extra-catalytic allosteric sites have been limited to a few examples including IKK (9) and MEK (10). Alternatively, bipartite inhibitors have been developed that stabilize an inactive conformation of the protein kinase, the prototypical example being BIRB796 binding to p38 and Gleevec binding to Abl. Such compounds make contact with both the conserved ATP site and less conserved regions of the activation loop, thus offering the potential for improved selectivity (11). The N terminus of the activation loop contains an invariant Asp-Phe-Gly (DFG) motif, and is an important determinant of enzyme activity. In the active or "DFG-in" conformation, these amino acids are involved in the coordination of ATP. Conversely, the "DFG-out" state does not bind ATP and the kinase is inactive. While a handful of kinases are known to adopt a DFG-out conformation (e.g. p38, Abl, etc), it remains to be determined how general this strategy might be in the design of selective kinase inhibitors.
To help elucidate the molecular mechanism of PYK2 and its substrate specificity, we used biophysical methods and determined multiple x-ray structures of the PYK2 kinase domain. High-resolution structures of apo and ATP␥S-bound forms were obtained as well as a complex with PF-431396, a "classical" kinase inhibitor. Empirical screening identified BIRB796 as a weak PYK2 kinase inhibitor. Surface plasmon resonance (SPR) and NMR studies indicated that PYK2 could adopt a "DFG-out" conformation. Despite the low affinity, a 1.75-Å co-crystal structure was obtained with BIRB796 revealing a novel binding mode. Our biophysical and structural results provide insight into the enzyme-substrate complex and allowed us to advance the rational design of a selective DFG-out inhibitor with improved PYK2 selectivity and potency. The compound, PF-4618433, showed robust osteogenic activity in hMSC cultures.

EXPERIMENTAL PROCEDURES
Protein Cloning, Expression, and Purification-A modified pFastBac vector was constructed to have an N-terminal His 6 tag followed by a thrombin cleavage site. PfuTurbo (from Stratagene) was used to amplify the catalytic domain of PYK2 from a full-length clone. The 854-bp DNA fragment was digested with BamHI and XhoI and cloned into the modified pFastBac vector at these sites to yield an N-terminal His 6 -tagged thrombin cleavable fusion with PYK2 catalytic domain residues Pro-416 to Glu-692 based on accession U33284. Standard expression conditions using the Bac-to-Bac Baculovirus Expression System from Invitrogen and Sf9 cells were employed for large-scale expression in shake flask and WAVE BioReactors.
Baculovirus cell paste containing the overexpressed His-PYK2 catalytic domain (416 -692) recombinant protein was resuspended in 3 volumes of buffer A (100 mM KH 2 PO 4 / K 2 HPO 4 ; pH 8.0, 200 mM NaCl, 1 mM TCEP, Complete TM protease inhibitor (1 tablet/50 ml), Roche, 2 mM phenylmethylsulfonyl fluoride, 10% glycerol, 0.5 l of Benzonase (Sigma) per 50 ml volume). All reagents and cells were kept on ice. Cells were lysed with 2 passes on a microfluidizer (pressure at 18 K). The lysate was clarified by centrifugation at 4°C for 35 min at 30,000 ϫ g (16,000 rpm) in a Sorval SL-50T (or appropriate) rotor. The supernatant was extracted and dispensed into a tube containing 20 ml of TALON Superflow Metal Affinity Resin (BD-Clonetech), pre-equilibrated with buffer B (100 mM KH 2 PO 4 /K 2 HPO 4 ; pH 8.0, 200 mM NaCl, 1 mM TCEP, 2 mM phenylmethylsulfonyl fluoride, 5 mM imidazole; pH 7.5). The suspension was incubated at 4°C for 3 h, with gentle rocking, then centrifuged at 700 ϫ g in a swinging bucket rotor for 2-4 min, and the supernatant was extracted. The supernatant was kept, and the resin resuspended in 20.0 ml of buffer B (this wash step was repeated four more times). The resin was then transferred to an XK-26 column (Amersham Biosciences) connected to a FPLC TM and washed with buffer B until the OD reached 0.05 units. Peak elution was achieved with a gradient using buffer C (buffer B with 250 mM imidazole). Eluted fractions were pooled and buffer exchanged into buffer D (20 mM Tris; pH 7.5, 150 mM NaCl, 1 mM TCEP) using multiple runs on a HiPrep 26/10 desalt column. Pooled fractions from the desalt step were concentrated to no more than 10 mg/ml and 1-ml aliquots loaded onto a Superdex 200 16/30 column, equilibrated with buffer E (25 ml HEPES; pH7.6, 100 mM NaCl, 1 mM TCEP). Peak fractions were pooled from multiple runs and concentrated to 10 -15 mg/ml using a Vivaspin 6 ml concentrator (10 K cutoff), aliquoted into 110-l volumes and flash-frozen at Ϫ80°C for use in crystallization trials.
NMR and SPR Studies-The NMR spectra of 15 N-PHE-labeled PYK2 were collected at 600 MHz on a Bruker DRX600 spectrometer equipped with cryoprobe. The concentration of the PYK2 protein was 100 M in 50 mM deuterated Tris buffer, 100 mM NaCl, and pH 7.5. Ligands were added from deuterated DMSO stock solution to yield final ligand concentration of 150 M. To assign the Phe of the DFG loop, a mutated sample of [ 15 N]Phe-labeled PYK2 with the Phe of the DFG loop mutated to Trp was made, and TROSY spectrum collected. To confirm the assignment, a sample of 13 C-Asp-15 N-Phe labeled PYK2 was also made, and high resolution TROSY spectra of the protein with and without 13 C decoupling were collected (data not shown). Because the DFG loop contains the only Asp-Phe pattern in the PYK2 kinase sequence, the Phe of the DFG loop was identified by the change of 13 C-coupling pattern in the TROSY spectra.
SPR kinetic measurements were performed on a Biacore 3000 biosensor (GE Healthcare, Piscataway, NJ). PYK2 catalytic domain was biotinylated with maleimide PEO biotin (Pierce Chemical) and immobilized on a neutravidin surface prepared on a CM4 chip by standard direct amine coupling procedures. Kinetic analyses were performed at 25°C in 50 mM Tris/HCl pH 7.5, 150 mM NaCl, 10 mM MgCl 2 , 1 mM MnCl 2 , 0.05% Tween 20, 3% DMSO. Compounds were injected over sample and reference flow cells at 50 or 100 l/min. Surfaces were regenerated by prolonged buffer flow. Data were processed using Scrubber 2 software (BioLogic Software) to y-zero, x-align, and double reference against buffer injections and a neutravidin surface lacking any immobilized ligand. Kinetic analyses were performed using BiaEval software (GE Healthcare).
Crystallization-Crystallization of apo-PYK2, with the BIRB796, with PF-431396, and with the PF-4618433 inhibitor was achieved using a 5 mg/ml protein stock. For co-crystallization experiments, inhibitor (from a DMSO stock of 30 mM) was used to give a final concentration of between 0.5-1 mM. Hanging drops of 2 l ϩ 2 l were setup over 750 l of well solution in standard Linbro plates and incubated at 22°C. Initial screening using Hampton screens produced a condition, 0.1 M bis-Tris: pH 6.0 -7.0, 0.2 M MgCl 2 , 20 -27% PEG 3350, 1 mM TCEP, which gave needle crystals that were optimized by refinement of the conditions as well as multiple rounds of streak and microseeding. For the inhibitor co-complex crystallization, the optimized crystals were further soaked overnight in 3 mM inhibitor to increase the occupancy of the inhibitor.
Crystallization of PYK2 with ATP␥S was achieved using a protein solution at 15 mg/ml with 1-2 mM ATP␥S added, 2 ϩ 2 l hanging drops were setup over 750 l of well solutions con-taining 0.1 M citrate; pH 5.0 -6.0, 0.5-1.5 M lithium sulfate, and 0.2-0.4 M ammonium sulfate. Plates were incubated at 22°C, and crystals resembling plate clusters were obtained after 3-5 days.
Structure Determination-Diffraction data were collected from flash-frozen crystals at 100 K at beamline 17-ID at Advanced Photon Source of Argonne National Laboratory. Data were processed using the HKL2000 suite of software (12). Data collection statistics are summarized in Table 1.
The structure of PYK2-apo was solved by molecular replacement methods with the CCP4 version of PHASER (13), using the FAK structure (PDB code 1MP8) as a search model. (During the manuscript preparation, a PYK2-apo structure was deposited with PDB code of 3CC6). After molecular replacement, maximum likelihood-based refinement of the atomic position and temperature factors were performed with REFMAC (14), and the atomic model was built with the program O (15). The refined PYK2-apo was then used as a starting model for all other complexes. The stereochemical quality of the final model was assessed by PROCHECK (16). Crystallographic statistics for the final models are shown in Table 1. Structural alignments with PYK2-apo and other complexes were performed using the lsq options in O (15). Figures were prepared with PYMOL (Delano Scientific).
Osteogenesis Assays-hMSC cultures, alkaline phosphatase assays, and mineralization analysis were performed as described previously (6), except that for mineralization analysis blue calcein was added on day 7, cells fixed on day 14, and analyzed for total intensity.
Image Analysis of Calcein-labeled Bone Nodules-For hMSC cultures, osteoblast nodule formation was quantified by application of the Cellomics Spot Detector Bio Application. Images were taken at a fixed exposure and the blue fluorescence of calcein-labeled nodules was captured with a 12-bit camera and a 10ϫ objective. One hundred fields were collected per well. If more than 10 fields were assayed with 5 or less objects per field, the scan for that well was aborted. This criterion was set so as not to completely scan wells devoid of nodules. Objects with fluorescence above a minimum threshold were evaluated for size, shape, and average/total intensity. The object selection criteria were applied uniformly to all wells, and designed to eliminate fluorescent objects that did not resemble osteoblast nodules. All analyzed fields were also manually examined for image processing artifacts (bubble, flare, or dust), and all fields containing artifacts were excluded from the analysis.
Kinase Assays and Selectivity Profiling-Kinase assays for PYK2 and FAK were as described (6) except that unactivated full-length, C-terminal his-tagged human proteins were used in place of the Src-activated kinase domain protein and that test compounds were preincubated with enzyme for 1 h prior to the addition of 500 M ATP. IC 50 determinations were generated from 8 point dose response curves. Selectivity profiling against a diverse panel of kinases was performed by a commercial vendor (Invitrogen) and reported as average percent inhibition of duplicate tests.

RESULTS AND DISCUSSION
Overall Structures of PYK2-Five crystal structures were determined for unphosphorylated human PYK2 kinase domain, all of which were refined to acceptable R-values and geometry ( Table 1). The overall structure is very similar to other kinase catalytic domains, displaying a bilobal architecture. The smaller N-terminal lobe contains a five-stranded antiparallel ␤-sheet and a single ␣-helix (␣C). The larger C-terminal lobe is mostly ␣ helical and contains the activation loop involved in polypeptide substrate binding. N-and C-terminal lobes are connected by the so-called hinge sequence, which partially defines the binding site for ATP and ATP-competitive kinase inhibitors (Fig. 1). Two structures contain ATP-mimetic inhibitors (ATP␥S and PF-431396) and adopt an active, closed conformation in the DFG motif. In contrast, the structures in the presence of diaryl urea inhibitors (BIRB796 and PF-4618433) reveal an allosteric binding pocket distinct from the ATP pocket. The activation loop residues beyond DFG motif are partially disordered in all five structures and are not visible in the electron density maps; residues 577-586 (Apo),   D. ϭ 2.3 Å). These substantial struc-tural differences are in part due to the unusual feature of an intramolecular disulfide bond observed in the FAK crystal structure. In the FAK-ADP complex, Cys-456 and Cys-459, which are widely conserved in vertebrate FAK sequences, form a disulfide bond on the surface of the N-terminal lobe adjacent to the ␣C helix (17). In all five PYK2 structures, however, this disulfide bond is absent and the corresponding structural region is directly connected to an ␣C helix (Fig. 1).

Structures of ATP-␥S and PF-431396 Complexes-
The ATP binding sites of protein kinases are the most common targets for the design of small molecule inhibitors. Discovery and optimization of ATP-competitive inhibitors had been perceived as a difficult obstacle to overcome due to the highly conserved B, chemical structures of three PYK2 inhibitors used for co-crystallization. C, structure-based sequence alignment of the PYK2 with FAK kinase domains. Secondary structure elements of the PYK2 are represented by noodles (␣-helices) and arrows (␤-strands). They are colored cyan for helices and orange for strands. Identical residues are highlighted in magenta. The DFG motif is boxed in red. The conserved residues involved in forming the back pocket are shown as F below the sequences. The gatekeeper residue is shown as an asterisk below sequences. The Leu-504, which undergoes significant conformational changes with BIRB796, is shown as OE. The two conserved cysteine residues that form a disulfide bond in FAK are shown as f. nature of the ATP binding site among kinase domains. However, small structural differences and plasticity between the ATP binding sites of even closely homologous kinases has been successfully exploited to achieve selectivity and potency (18,19). In the PYK2-ATP␥s complex, the adenine base inserts into a hydrophobic pocket formed by the side chains of Leu-431, Ala-455, Val-487, Leu-504, and Leu-556. The N1 nitrogen and the amino group of the adenine ring hydrogen bond to the main chain of Tyr-505 and Glu-503, respectively. This H-bond pattern is at the heart of the interaction of ATP with kinases, and has been used extensively to design and optimize ATP competitive inhibitors (Fig. 2A).
Inhibitors of the trifluoromethyl pyrimidine series (PF-431396) occupy the nucleotide binding pocket, with the pyrimidine ring located in the adenine pocket (Fig. 2B). Both the N1 nitrogen of pyrimidine ring and amide nitrogen of the 5-amino oxindole form hydrogen-bonding interactions with the carbonyl and amide backbone of Tyr-505. The oxindole ring is further stabilized by water-mediated hydrogen bonding interaction with the side chain of Glu-509. The trifluoromethyl group is involved in hydrophobic interactions with Val-487, Met-502, and Ala-455. The hydrophobic phenyl ring occupies the ribose binding pocket and is perpendicular to the pyrimi-dine ring. The amine group connecting pyrimidine and phenyl ring forms a water-mediated hydrogen bond interaction with Asp-567 of the DFG motif. The methyl group attached to the sulfonamide is wedged into the glycine-rich loop by hydrophobic interactions with Leu-431, Gly-432, and Val-439. PF-431396 exhibits strong dual inhibitory activity for both PYK2 and FAK with an IC50 of 11 and 2 nM, respectively. Molecular modeling of PF-431396 predicts the same binding mode in FAK with small conformational changes in the glycinerich loop.
The structural analysis of PF-431396 offered an explanation for the lack of compound selectivity between these two FAK family enzymes. In fact, more than half of the kinases tested with PF-431396 in a diverse subset panel at 10 M showed Ͼ50% inhibition (Fig. 3), and greater than 90% inhibition was observed with this compound for JAK3, TrkA, and Aur2 even at 1 M (data not shown). Thus, while PF-431396 shows an overall preference for the FAK family, it is not absolutely selective. Off-target cross activity may be related to the modest activity of bone efficacy in a cell culture model (see Fig. 7, below). Therefore, we sought to determine whether the PYK2 kinase might be a candidate for a DFG-out inhibitor approach, a strategy that has led to compounds with improved selectivity against other kinases. The prototype DFG-out p38 kinase inhibitor, BIRB796, was profiled in a PYK2 enzyme assay and was active with an IC 50 of 1.5 M (Fig. 3). Surface Plasmon Resonance (SPR)-Classical kinase inhibitors bind catalytic sites at rates that approach diffusion limits (20,21,22), whereas DFG-out inhibitors of p38 MAP kinase have been shown to bind two to three orders of magnitude slower (23,24). Results of SPR kinetic analysis show that this is also the case for PYK2 (Fig. 4). As expected, the classical inhibitor, PF-431396 binds with an on-rate constant of 10 6 /Ms, whereas nonclassical inhibitors, BIRB796 and PF-4618433 bind with on-rate con-  stants of 10 3 /Ms. It is of interest to note that BIRB796 binds to both activated and non-activated PYK2 catalytic domains with similar on-rate constants (Fig. 4) as well as the activated and non-activated full-length enzyme. Similar effects were also observed with DFG-out inhibitors of the p38 MAP kinase catalytic domain. 3 These data may have implications about the solution structure of activated forms and the mechanism of DFG-out binding. It is often assumed, based on crystal structures that the activated kinase resides primarily in the DFG-in conformation. If this were indeed the case, one might expect that the on-constants for the active kinase would be much slower than for non-activated kinase as the latter exhibits more conformational heterogeneity based on NMR. The fact that these constants are essentially the same suggests that either equivalent sampling of the appropriate DFG-out conformation occurs or that "DFG-out compounds" can induce the conformational change regardless of activation state.
Another point of interest is that despite considerable efforts to develop more potent DFG-out inhibitors of PYK2, significantly faster on-rate constants were not achieved. This may imply that PYK2 only infrequently samples the particular DFGout conformation required for binding resulting in only a small fraction of the total population as competent for binding. This upper limit on the rate of association may represent an insurmountable energy barrier and has implications for the successful development of DFG-out PYK2 inhibitors, necessitating a significant decrease in off-rate as observed in the p38-BIRB796 interaction.
NMR Spectroscopy-It has been demonstrated for p38 kinase that binding of DFG-out inhibitors to P38 kinase causes sharpening of TROSY spectra due to their ability to stabilize the DFG-out conformation (25).
Using the same approach, BIRB796 was assessed as a DFGout inhibitor for PYK2. All the phenylalanine residues of the PYK2 kinase domain were selectively labeled with 15 N, and TROSY spectra were collected. In the apo state, the TROSY 3 Non-published results for BIRB796 and Pfizer inhibitors. spectrum of 15 N-Phe labeled PYK2 showed only a few phenylalanine peaks, with the rest of the peaks missing (Fig. 5). The missing peaks are the result of extreme line broadening caused by conformational exchange in the intermediate exchange regime of the NMR time scale (millisecond). Tight binding of a classical ATP site inhibitor of the trifluoromethyl pyrimidine series does not change the characteristic of the TROSY spectra, indicating that the DFG-loop is still flexible. However, the addi-tion of BIRB796 caused a dramatic change in the TROSY spectra. The linewidth of all the peaks became uniform with the peaks in the rigid region of the structure, and the missing Phe in the DFG loop was clearly visible. These data indicate that BIRB796 was able to freeze the DFG loop motion, consistent with it being a DFG-out binder as seen in P38 kinase. The millisecond exchange rate observed by NMR line broadening is consistent with the slow 10 3 /Ms on-rate of DFG out inhibitors as measured by SPR. This suggests that the slow on-rate of DFG-out inhibitors is limited by the intrinsic rate of opening of the DFG-out conformation. Therefore, further optimization of DFG-out inhibitors will likely result from improvement in the off-rates instead of the on-rates.
DFG-out Conformation of PYK2 upon BIRB796 Binding-In the PYK2-BIRB796 complex, the activation loop forms an inactive conformation with the DFG motif in a catalytically incompetent DFG-out state. The Phe-568 side chain in the DFG motif moves by ϳ11 Å and sandwiches the naphthyl moiety of the BIRB796 with the gatekeeper residue, Met-502 (Fig. 6A). Interestingly, the naphthyl moiety of BIRB796 in the PYK2 structure is perpendicular to that of BIRB796 in the P38␣ structure (23) (Fig. 6B). Similar to the P38-BIRB796 structures, several distinct inhibitors such as imatinib, sorafenib and AAL993, which have been crystallographically proven to stabilize the DFG-out conformation (26,27,28), show significantly different orientation of the phenyl moieties in the crystal structures of ABL, B-RAF, and VEGFR (Fig. 6C). The different binding mode of the naphthyl moiety seen in PYK2 compared with the corresponding moieties in P38␣, ABL, B-RAF, and VEGFR is in part due to the difference in the gatekeeper residue. In the PYK2-BIRB796 complex, a conserved gatekeeper residue, Met502, connecting  the N-terminal domain and the hinge loop (Glu-503-Gly-508) plays an important structural role and is involved in forming an additional lipophilic pocket with Met-478, Leu-489, and Ile-500. In contrast, the crystal structures of P38␣-BIRB796, ABLimanitib, B-RAF-sorafenib, and VEGFR-AAL993 all contain smaller threonine or valine side chains as the gatekeeper residue. As a result, the rotation of the naphthyl ring to 90°will lead to steric clashes with Met-506 in PYK2.
The tolyl group on the pyrazole ring interacts with the hydrophobic portion of the side chain of the conserved Glu-474, which makes hydrogen bonds with both NH groups of the urea moiety. The t-butyl group is buried in a hydrophobic pocket composed of Leu-481, Ile-486, Leu-540, and Leu-565. The DFG-out conformation is stabilized by H-bond interactions of the urea carbonyl group with the amide backbone of Asp-567 of the DFG-motif and the amide backbone of Phe-568 with the Asn-554 side chain. The Phe side chain is further stabilized by hydrophobic interactions with Leu-556, Val-487 and the naphthyl group of BIRB796.
Importantly, the morpholino substituent in PYK2-BIRB796 complex does not form a hydrogen bond at the hinge loop. In the P38-BIRB796 structure, the ethoxymorpholino substituent on the naphthyl ring enhanced the affinity by ϳ11 fold (23) probably due to a hydrogen bond with the backbone amide of residue109, which is equivalent to the H-bond made by the N1 atom of the adenine ring of ATP␥S. Upon BIRB796 binding in PYK2, the ␤-turn in the hinge loop formed by the carbonyl and amide backbones of Tyr-505 and Gly-508 vanishes due to hinge movement of the N-terminal lobe and the Leu-504 side chain moves ϳ2.5 Å to block the ATP-binding pocket to form a hydrophobic pocket with Leu-431, Ala-455, Val-487, Met-502, Leu-556, and Phe-568 (Fig. 6, A and D). This movement of the leucine side chain prevents hydrogen bond interactions with backbone residues in the hinge loop and changes the shape and nature of the ATP binding pocket. Hence, the DFG-out conformation with Leu-504 movement is incompatible with ATP binding because of steric overlap of both Phe-568 and Leu-504 with the phosphate and adenine groups of ATP. To our knowledge, this is the first example that a leucine residue in the hinge loop blocks the ATP binding site in the DFG-out conformation.
Conformation of PYK2 in Complex with PF-4618433 and Its Pharmacological Activity-With a definitive crystal structure in hand demonstrating a DFG-out binding mode for PYK2, we undertook efforts to modify BIRB796 to improve potency and selectivity for PYK2. Initial efforts were focused on modification of the tolyl-tert-butylpyrazole moiety. Modification of this group generally decreased potency, which is surprising given that this group is frequently viewed as a selectivity element. Conversely, modification of the naphthyl group led to modest improvements in potency. Through a combination of library and singleton chemistry efforts, we identified PF-4618433 wherein the N-morpholinylethoxynaphthalene of BIRB796 has been replaced with a pyridinyloxymethylpyrazole. This modification gave rise to a PYK2 inhibitor with an IC 50 of 637 nM.
In contrast to the increased osteogenesis observed with shPYK2 or expression of DN PYK2 in hMSC cultures previously reported (6), PF-431396 treatment led to a modest albeit significant increase in day 7 alkaline phosphatase activity, but only when dosing was restricted to days 4 -7 (Fig. 7A). Offtarget activity may be obscuring the osteogenic action of PF-431396, as the compound-treated cells appeared to be stressed and the cultures did not mineralize (data not shown). While kinase selectivity profiling indicated that the two DFGout compounds have a more restricted spectrum of activity, PF-4618433 showed improved PYK2 potency, reduced p38 activity, and superior overall selectivity relative to the prototype BIRB796 (Fig. 3). Thus, PF-4618433 was tested for pharmacological activity in the functional osteogenesis assay. Cultures treated with increasing concentrations of PF-4618433 starting at day 1 and maintained throughout showed a dose-dependent increase in both alkaline phosphatase activity (Fig. 7B) and mineralization (Fig. 7, C and D). The improved osteogenic activity of the DFG-out PYK2 inhibitor PF-4618433 was revealed across a series of DFG-out analogs (LB, PB, and AB, data not shown), , and alkaline phosphatase activity determined on day 7 (A and B, respectively). hMSC cultures treated with PF-4618433 were analyzed for mineral by incorporation of blue calcein and the well imaged by fluorescent micrograph (see "Experimental Procedures" for details). The individual field corresponding to the median total fluorescence intensity from each treatment group is shown in C and total fluorescence quantitated in D. All data are presented as means Ϯ S.E. *, p Ͻ 0.01. Significance was determined by unpaired t test to compare differences with the DMSO vehicle control. and correlates with the enhanced kinase selectivity profile when compared with classical inhibitors such as PF-431396.
In the crystal structure of PYK2 complexed with PF-4618433, the pyrazole ring rotates ϳ40°relative to the naphthyl ring in BIRB796 to form a hydrogen-bond interaction with the catalytic Lys-457. Unlike the Leu-504 blocking the adenine binding site in PYK2-BIRB796, the Leu-504 remains in a similar conformation to that seen in the apo structure (Fig. 8). As a result, the pyridyl nitrogen forms a hydrogen-bond interaction with the backbone amide of Tyr-505 in the hinge loop. Notably, three additional residues in the activation loop, Leu-570 -Arg-572 following the DFG motif, are visible in the co-crystal structure and the guanidium group of Arg-572 is stabilized by a stacking interaction with the tolyl group of the PF-4618433. Additionally, Asp-576 and Glu-577 from either the same molecule or a neighboring packing molecule are visible and shield the entire tolyl group of the PF-4618433.
In conclusion, the crystal structures of PYK2 with three different inhibitor complexes show the conformational variability of the DFG motif. The DFG-out conformation in two diarylurea inhibitors reveals a distinct subclass of non-receptor tyrosine kinases identifiable in part by gatekeeper Met-502 and the unique hinge loop conformation of Leu-504. Our structural and solution studies suggest that the important features of the DFG motif including Leu-504 variability in the hinge loop can be utilized in the design and development of new selective inhibitors against protein kinases.