High-Throughput Interrogation of Ligand Binding Mode Using a Fluorescence-Based Assay

Probing the pocket: A high-throughput fluorescence-based thermal shift (FTS) assay utilized different forms of a protein (in gray) to establish the binding mode of a ligand (see picture). The assay serves in the rapid evaluation of structure-activity binding-mode relationships for a series of ligands of Plk1, an important target of anticancer therapy.

Paweł Ś ledź, Steffen Lang, Christopher J. Stubbs, and Chris Abell* Structure-based ligand discovery has benefitted from extensive developments in recent years. [1][2][3][4] In particular, the application of high-throughput crystallography has vastly accelerated these efforts. [5] This approach is based on the rapid generation of protein-ligand complexes through the soaking of protein crystals with ligands and subsequent automated X-ray data collection and structure determination. However, suitable soakable crystals are not available for many proteins, especially for those complexes in which protein-protein interactions or structural changes play a role. Consequently, there is a need for more highthroughput and crystallographically independent methods for assessing the binding mode of molecules to select candidates for further analysis.
This has become particularly important to us while developing more potent ligands against the polo-box domain (PBD) of polo-like kinase 1 (Plk1). The PBD domain is responsible for the proper cellular localization of Plk1 through an array of phosphorylation-dependent proteinprotein interactions mediated by the phosphopeptide-binding groove. [6,7] Disrupting these interactions has shown potential as a strategy for anticancer therapy. [8] The discovery of an auxiliary flexible pocket on the PBD surface, which is involved in binding a peptide derived from polo-box interacting protein 1 (PBIP1), has been recently reported. [9,10] This flexible hydrophobic pocket recognizes its ligand through hydrophobic interactions of four aromatic residues (Y417, Y421, Y481, and F482) with F71 PBIP1 and an additional hydrogen bond between the hydroxy group of Y417 and the backbone nitrogen of D72 PBIP1 (Figure 1 a). [9] Our objective was to prepare more potent ligands of the PBD by improving binding to this auxiliary pocket by optimizing the hydrophobic interactions whilst retaining the key hydrogen bond. To efficiently distinguish specific recognition by the pocket from nonspecific hydrophobic interactions, we adapted a fluorescence-based thermal shift (FTS) assay. [11] In this assay the ability of a molecule to stabilize the protein during its thermal unfolding is quantified by its thermal shift (DT m ): the difference in the protein unfolding temperature in the presence and absence of a ligand. FTS has been successfully applied in fragment and high-throughput screening campaigns; [12,13] however, it generally provides very limited information on binding and it is not always appropriate for weakly binding ligands.
We designed a procedure through which information on the binding mode can be obtained using FTS. When a modified protein that cannot fully engage a ligand in the selected binding mode is used, ligands adopting such a binding mode are expected to stabilize the mutant protein against thermal unfolding to a lesser extent than the wild-type protein; ligands not engaging in the selected binding mode are expected to bind similarly to both mutant and wild-type protein (Figure 1 b). The difference in thermal shift of a ligand for the wild-type and mutant protein (DDT m ) is therefore a measure of how beneficial the pocket binding is for a particular ligand. Such an experimental setup eliminates errors arising from the concentration dependence of the assay, [14] as a single stock solution of the ligand can be used. Also noteworthy is that since in this particular case we were working with elaborated ligands with significant DT m values, the sensitivity of the assay (DT m > 1.0 K) was not a problem and identification of ligands with a significant DDT m value (higher than 2.0 K) was straightforward.
For the proof-of-principle study, we first obtained three different mutants of the PBD to assess their ability to distinguish binding in the hydrophobic pocket (Figure 1 a). The Y417A mutation was expected to prevent the formation of the key hydrogen bond with the ligand. A double mutant with an additional Y421A mutation was also prepared to further reduce the possible scope of hydrophobic interactions. [9] To assess the specificity of these effects, we prepared a Y481K mutant, which changed the distal end of the pocket, where interaction with the ligand is limited. These mutants were subsequently tested against two PBIP1-derived ligands-FDPPLHSpTA (1 a), a peptide that binds in the pocket, and ADPPLHSpTA (1 b), which does not. Both the Y417A and Y417A/Y421A mutants exhibited detectable and significant DDT m values for 1 a (3.7 and 3.2 K, respectively), while DDT m values for 1 b were within the experimental error (1.2 and 0.1 K). This was further confirmed by the study with the Y481K mutant, which gave negligible DDT m values for 1 a (À0.8 K) and 1 b (À1.0 K).
We used the double-mutant Y417A/Y421A to further probe the molecular recognition of the pocket (Table 1). A series of analogues of PBIP1-derived peptide 1 a with phenylalanine replaced by other acetylated hydrophobic and aromatic amino acids (2 a-f) were prepared and tested in the FTS binding mode assay. Interestingly, the pocket turned out to be very specific for hydrophobic aromatic residues. Tryptophan-bearing 2 e exhibited a significant DDT m (3.3 K), indicating hydrophobic pocket binding, in contrast to histidine-bearing 2 f (DDT m = 0.5 K). None of tested aliphatic residues (2 a-d) showed DDT m values higher than 1.0 K.
To further study the molecular recognition of the pocket, we expanded the series of PBIP1-derived peptides to include non-amino acid N-terminal residues, focusing on, but not limited to, hydrophobic aromatic ring systems ( Table 1). As expected, neither of the tested aliphatic ring compounds 3 a,b showed binding in the pocket. Only two out of the four simple phenylalanine analogues 3 c-f gave DDT m values exceeding 2.0 K. Peptide 3 d with a phenethyl substituent showed slightly better stabilization than the original PBIP1 peptide 1 a. Two simple analogues of 3 d, however, showed no indication of hydrophobic pocket binding. Ligand 3 c, lacking one methylene group, was too short to comfortably insert the phenyl ring to the bottom of the pocket, while the p-methoxy substituent in analogue 3 e appeared to cause a steric clash at the bottom of the pocket. Interestingly, 3,4-dichloro-substituted 3 f was found to bind well in the hydrophobic pocket. The two ortho-chloro substituents in the ring mimic the presence of an additional benzene ring because of their size and hydrophobicity. Three peptides incorporating heterocyclic single-ring systems 3 g-j were also prepared and tested in the assay. Interestingly, the pocket showed significantly higher selectivity for thiophene-derived 3 h over pyridine-based 3 g or furan-based 3 j, even though for 3 h the observed DDT m value of 2.5 K was not as high as observed for 1 a. In addition, the two double-ring aromatic systems 3 j,k were tested, both giving high DDT m values indicative of hydrophobic pocket binding. The structures of these compounds suggest that they are able to insert their benzene rings into the pocket, reinforcing the suggestion of its selectivity for aromatic moieties.
Analysis of peptide binding to the protein by isothermal titration calorimetry (ITC; Table 1 and Table S1 in the Supporting Information) revealed a clear correlation between Table 1: The results of the assay for the series of modified peptides (R-DPPLHSpTA-NH 2 ). DT m values are shown for the wild-type (WT) and double-mutant Y417A/Y421A (DM) protein; the difference between them (DDT m ) and the K D value against the wild-type protein is also shown. The color code corresponds to that used in Figure 2  the DT m for wild-type protein and binding affinity (pK D ) (Figure 2 a). The peptides are clustered into three groups (designated weak, medium, and tight binders), with the initial pocketbinding peptide 1 a being in the middle group. The correlation between DDT m and K D furnishes additional insight into the binding mode (Figure 2 b). Peptides identified in the assay as pocket binders have in general shown higher affinity than those not utilizing the pocket, and all the most potent ligands (3 d, 3 f, 3 j, and 3 k, colored purple) also bind to the pocket exhibiting the highest DDT m > 4.0 K. None of the peptides designated as weak binders (in blue) showed pocket binding. The seven compounds clustered together as medium binders in the center of the plot of pK D against DT m (Figure 2 a) are resolved in the plot of pK D against DDT m (Figure 2 b). Peptides 3 a-c, e (in red) gain additional binding affinity apparently from pocket-independent hydrophobic interactions, in contrast to similarly potent 1 a, 2 b, and 3 h (in green), shown to be pocket binders. These two different effects responsible for the increased affinity of our ligands can be readily distinguished by DDT m values, allowing for informed binder optimization.
To confirm the binding modes, we co-crystallized two representative examples with the PBD, the leucine-based peptide 2 a (DDT m = 0.3 K) and the most potent ligand 3 j (DDT m = 4.1 K). As predicted, peptide 2 a turned away from the pocket, corroborating the result of the assay (Figure 2 c). In contrast, 3 j bound in a way similar to 1 a, retaining the key hydrogen bond and the benzene ring of thiophene was inserted into the pocket, mimicking the phenylalanine of 1 a (Figure 2 d).
Our FTS-based assay allowed us to rapidly interrogate the binding mode of a series of peptides binding to the PBD of Plk1. The approach is high throughput, as all the FTS experiments described in Table 1 could be performed in duplicate on a single 96-well plate in less than two hours. The results of this assay allowed us to readily observe both the correlation between binding mode and ligand affinity as well as the selectivity of the pocket for particular aromatic moieties. They validate DDT m as a valuable parameter in the structure-based ligand-discovery process, allowing for much deeper insight into structure-activity binding-mode relationships than a purely activity-oriented screen. [15] As our experimental setup is generalizable, we believe that its adaptation for other targets may be useful, especially for challenging systems like protein interfaces. Interpretation and validation of the assay data. a) A plot of pK D vs. DT m shows a correlation between thermal shift and binding affinity; b) a plot of pK D vs. DDT m(WT) distinguishes two binding modes leading to increased affinity of modified peptides-pocket-dependent and independent; colors of the data points in (a, b) reflect classification of binders into one of four groups: weak nonpocket (blue), medium nonpocket (red), medium pocket (green), and tight pocket (magenta) and are as in Table 1; c) the structure of 2 a bound to the protein shows it does not bind to the pocket, which is occupied by PEG300; d) the structure of 3 j shows that it binds in the pocket, as indicated by the FTS assay.