News center
Boasting modern processing equipment

Discovery of Small Molecules that Inhibit the Disordered Protein, p27Kip1

Dec 16, 2023

Scientific Reports volume 5, Article number: 15686 (2015) Cite this article

6163 Accesses

63 Citations

87 Altmetric

Metrics details

Disordered proteins are highly prevalent in biological systems, they control myriad signaling and regulatory processes and their levels and/or cellular localization are often altered in human disease. In contrast to folded proteins, disordered proteins, due to conformational heterogeneity and dynamics, are not considered viable drug targets. We challenged this paradigm by identifying through NMR-based screening small molecules that bound specifically, albeit weakly, to the disordered cell cycle regulator, p27Kip1 (p27). Two groups of molecules bound to sites created by transient clusters of aromatic residues within p27. Conserved chemical features within these two groups of small molecules exhibited complementarity to their binding sites within p27, establishing structure-activity relationships for small molecule:disordered protein interactions. Finally, one compound counteracted the Cdk2/cyclin A inhibitory function of p27 in vitro, providing proof-of-principle that small molecules can inhibit the function of a disordered protein (p27) through sequestration in a conformation incapable of folding and binding to a natural regulatory target (Cdk2/cyclin A).

Intrinsically disordered proteins (IDPs) are highly prevalent in eukaryotic cells and often perform regulatory and signaling functions1. As a structural class, IDPs are more tightly regulated than folded proteins from yeast to humans2 and, upon overexpression, are more prone to alter cell phenotype than their folded counterparts3. Thus, maintaining normal disordered protein homeostasis is critical for cellular behavior. The phenotypes associated with overexpression or hyperactivity of proteins with folded domains can be counteracted by small molecule inhibitors of, for example, enzyme function (e.g., Gleevec which inhibits BCR-Abl in chronic myelogenous leukemia4) or protein-protein interactions (e.g., ABT-263 and ABT-199 which inhibit BCL-2 and BCL-xL in hematological and lymphoid malignancies5,6). In contrast, disordered proteins represent challenging targets for inhibition by small molecules due to their dynamic and heterogeneous conformations. Nonetheless, progress has been made. For example, an inhibitor of the phosphatase PTP1B (MSI-1436), with known roles in diabetes, obesity and breast cancer, was recently shown through structural, biochemical and functional assays to act through an allosteric mechanism by binding to a disordered regulatory region of the enzyme7. Also, a small molecule (YK-4-279) that binds to the disordered EWS-FLI1 fusion oncoprotein associated with Ewing's sarcoma family tumors inhibited direct binding to RNA helicase A (RHA)8, a functional partner of EWS-FLI1 in tumor cells8 and altered RHA-dependent protein interactions and RNA splicing9. Both of these compounds (MSI-1436 and YK-4-279) exhibited on-target effects in cellular assays8,9. Other studies have identified small molecules that target disordered cMyc10,11,12, α-synuclein13 and Alzheimer β-amyloid peptide14. A recent computational study15 showed that a small molecule (10074-A4) previously reported to modulate cMyc function10 bound in different ways to different cMyc molecules within an ensemble of many disordered conformations, leading the authors to suggest the concept of "ligand clouds binding to protein clouds". In the studies discussed above, small molecules that target disordered proteins were discovered using a variety of approaches, including functional screens, in vitro binding screens and/or computational screens. Nuclear magnetic resonance (NMR)-based screening of low molecular weight small molecules termed fragments (reviewed in16) binding to folded protein targets is a well-established method for identifying initial "hits" in the process of drug discovery17,18. However, NMR-based fragment screening has not, to our knowledge, been applied to identify small molecules that bind to a disordered protein target. Here, we utilized NMR-based screening to identify fragment molecules that bind to and modulate the function of the prototypical disordered protein, p27Kip1 (p27; also known as CDKN1B), a regulator of the cyclin-dependent kinases that control eukaryotic cell division19.

The motivation for targeting p27 was two-fold. First, the structural and functional features of p27 are well understood20,21,22,23, providing an ideal model system for studying small molecule:disordered protein interactions. Second, the ability to chemically modulate p27 function would be beneficial in several biological settings. For example, p27 is inappropriately phosphorylated in breast cancer on threonine 157, which is associated with abnormal cytoplasmic localization and up-regulation of cell migration24,25,26,27. The availability of a small molecule inhibitor of p27 would be beneficial to prevent abnormal migration of breast cancer cells. Alternatively, in both sensory and non-sensory epithelial cells of the inner ear, p27 maintains cell cycle exit and terminal differentiation28 and its inhibition resulted in their cell cycle reentry and regeneration for hearing restoration29,30. While small molecules that inhibit the transcription of p27 have been reported31, here we developed approaches to identify small molecules that bind directly to p27 and have potential to alter its function in the two cellular settings discussed above.

The target of our studies was the N-terminal, kinase inhibitory domain of p27 (p27-KID), which binds to and regulates the catalytic activity of nuclear cyclin-dependent kinase (Cdk)/cyclin complexes that control eukaryotic cell division32. p27-KID, which is highly disordered in isolation20,21, adopts an extended conformation upon binding to Cdk2/cyclin A (Fig. 1a) that can be subdivided into three, functionally distinct sub-domains. Sub-domain D1 binds to a conserved pocket on cyclin A and blocks substrate recruitment33; sub-domain D2 forms intra- and inter-molecular (between p27 and Cdk2) β-strands upon binding to Cdk2 and also inserts a turn of helix into its ATP binding pocket, inhibiting kinase activity34; and sub-domain LH forms an α-helix that connects sub-domains D1 and D2. We hypothesized that, if small molecules that bind to p27-KID could be identified, they may induce the disordered polypeptide to adopt conformations that are incompetent for binding to Cdk/cyclin complexes. We tested this hypothesis by screening a library of fragment molecules for binding to p27-KID using NMR spectroscopy. We identified two subsets of fragment molecules (36 in total) that differentially bound weakly but specifically to two partially overlapped regions of p27-KID. From these subsets, we then generated pharmacophore models that allowed identification of additional small molecules that bound to p27-KID and further clarification of structure-activity relationships. A variety of assays, including fluorescence anisotropy, NMR spectroscopy and a Cdk2 kinase activity assay, were used to demonstrate that one of the identified small molecules displaced the kinase binding region of p27 from Cdk2 and partially restored catalytic activity of Cdk2. In addition, molecular dynamics computations provided insights into the dynamic "structure" of the region of p27 targeted by small molecules. Our results provide insights into the nature of interactions between small molecules and a disordered protein and demonstrate that such interactions can alter disordered protein regulatory function.

Identification of two groups of small molecules that bind specifically to distinct but overlapping regions of disordered p27-KID.

(a) Structure of p27-KID bound to Cdk2/Cyclin A (PDB ID 1JSU); the sub-domains of p27-KID, including D1, LH, D2.1, D2.2 and D2.3, are indicated. (b,c) Chemical structures of two small molecules, SJ572710 (SJ710) (b) and SJ572403 (SJ403) (c), that bound to p27-KID and are members of Group 1 and Group 2, respectively (see text for Group definitions). (d,e) Histograms showing individual 1H chemical shift perturbation values (through analysis of 2D 1H-15N "in-phase" HSQC spectra) from titrations of the compounds in (b) and (c), respectively, into 15N-p27-KID. The threshold for identifying specific interactions with p27-KID residues was defined as two standard deviations above the average of the perturbation values (represented by a dotted black line in the graph). The experimental spectral resolution in the 1H dimension (2.4 Hz) is represented by the dotted magenta line. Chemical shift perturbations for molar ratios of 15N-p27-KID (100 μM) to inhibitors of 1:30 are shown. (f,g) Amide proton chemical shift perturbations plotted versus the concentration of SJ710 (f) and SJ403 (g), respectively. Binding isotherms of select residues show specific binding between p27-KID and fragment hits. The trajectories of the chemical shifts (solid black lines) report global dissociation constants of 4.8 ± 1.3 and 2.2 ± 0.3 mM for the interactions of p27-KID with SJ710 (f) and SJ403 (g), respectively.

We used one-dimensional (1D) 1H WaterLOGSY35 and STD36 NMR methods to identify "fragment-like" small molecules37,38,39 that bound to p27-KID. Fragment molecules were selected from either a commercial library (1,100 compounds from the Ro3 collection, Maybridge/Thermo Fisher Scientific) or an in-house library of 1,222 compounds based on the "Rule of Three"40 and other criteria (see Methods). Two and seven molecules each were identified from the respective libraries to bind p27-KID (termed "hits"; representative 1D 1H NMR data is shown in Suppl. Fig. 1a,b and all preliminary hits are presented in Suppl. Table 1a). Binding sites for these compounds were identified by titration into 15N-p27-KID and analysis of two-dimensional (2D) 1H-15N HSQC spectra. Significant chemical shift perturbations (CSPs), which were largest for amide 1H resonances (see Methods), were only observed for amide groups of residues within the D2 sub-domain of p27-KID (p27-D2); eight hits caused CSPs within a short region with the sequence F87YY89 [(F, phenylalanine; and Y, tyrosine), termed sub-region D2.3; Fig. 1a,d] and one hit caused CSPs within the same region as well as within two other regions near residues W60N61 and E75WQ77 [(W, tryptophan; N, asparagine; E, glutamic acid; and Q, glutamine), termed sub-domains D2.1 and D2.2, respectively; Fig. 1a,e]. We termed these molecules Groups 1 and 2, respectively (Fig. 1b,c); representative 1H CSP histograms are shown in Fig. 1d,e and 15N CSPs are presented in Suppl. Fig. 2. Data for all other small molecules are shown in Suppl. Fig. 3.

We analyzed the Group 1 and 2 molecules using computational modeling and identified additional candidate p27-KID binding molecules in the two fragment libraries that were not detected by the original 1D NMR screens. 1D and 2D NMR analysis of these molecules led to the identification of 15 additional p27-D2-binding compounds (six with Group 1-like binding features and nine with Group 2-like features; Suppl. Fig. 3b and Suppl. Table 1b). CSPs were observed at high compound concentrations, consistent with relatively weak binding to p27-D2, but were specific to the noted regions within p27-D2 and rigorously reproducible. Sixteen-point 2D 1H-15N "in phase" HSQC NMR titrations of a Group 1 (SJ572710, hereafter termed SJ710) and 2 (SJ572403, hereafter termed SJ403) hits, respectively, into a constant concentration of 15N-p27-KID (100 μM), provided interpretable binding isotherms for specific resonances of p27-KID; the determined dissociation constant (Kd) values were 4.8 ± 1.3 and 2.2 ± 0.3 mM, respectively (Fig. 1f,g). The overlaid 2D 1H-15N HSQC NMR spectra are displayed in Fig. 2.

Fragment hits show specific interactions with specific regions within the D2 subdomain of p27-KID.

Overlaid 2D 1H-15N "in-phase" HSQC NMR spectra showing chemical shift perturbations of specific residues within the D2 subdomain of p27-KID after titration of fragment hits SJ710 (a, Group 1) and SJ403 (d, Group 2). For each hit, a total of sixteen spectra were recorded. The concentration of 15N-p27-KID (100 μM) was kept constant throughout the titration and the concentration of the inhibitor was varied from 0 (red) to 3 mM (blue). All experiments were performed on a 800 MHz spectrometer and utilized acquisition parameters that provided a spectral resolution of 2.4 Hz and 1.7 Hz in the 1H and 15N dimensions, respectively. The insets in the lower left show resonances for the tryptophan indole NH moieties (s-W60, s-W76). (b,e) Expanded regions (blue box) showing a subset of the p27-KID residues involved in interaction with SJ710 (b) and SJ403 (e), respectively. (c,f) Expanded regions (black box) showing two of the p27-KID residues that do not show any perturbations upon interaction with SJ710 (c) and SJ403 (f), respectively.

Three-dimensional (3D) molecular interaction field analysis41 (see Methods section) identified common chemical features of the two groups of p27-D2 binding molecules (Fig. 3a,b). This analysis revealed interaction "field points" around the small molecules that can favorably participate in electrostatic, van der Waals and hydrophobic interactions. These field points are used to calculate molecular similarity and to align molecules, even those with different two-dimensional topologies, to define a common pharmacophore.

Consensus field maps for the two groups of small molecules that bound to p27-KID.

Blue and red field points indicate where electropositive and electronegative groups are favored, respectively, in the protein binding partner; gold field points indicate regions where hydrophobic interactions are favored; and yellow field points indicate regions where favorable van der Waals contacts are possible. The size of the field point increases with the magnitude of the favorable interaction energy. Two representative molecules from Group 1 (a) and Group 2 (b) are shown aligned to their respective consensus field maps.

The molecules in Groups 1 and 2 have two or three heterocyclic aromatic rings but significantly differ in the distribution of favorable field points. Group 1 is mainly defined by a large hydrophobic core (multiple gold polygons in close proximity; Fig. 3a), one large electropositive interaction region (two cyan polygons in close proximity) and two smaller electropositive (cyan) and electronegative (red) field points at the opposite end of the hydrophobic core. In contrast, the field map for Group 2 (Fig. 3b) exhibited a smaller hydrophobic core relative to Group 1 and two equally-sized regions of favorable electropositive interaction. To expand the diversity of our chemical screen, we used the consensus field maps for Group 1 and 2 molecules to identify 184 additional possible p27-D2-binding molecules within a library of 10,455 commercially-available fragment-like molecules. 1D and 2D NMR analysis of these identified 12 additional p27-D2-binding compounds (8 with Group 1-like binding features and 4 with Group 2-like features; Suppl. Fig. 3c and Suppl. Table 1c). These additional molecules had CSP profiles comparable to the previously identified hits. To further test the specificity of small molecule:disordered protein interactions, we determined whether the simple amino acids tryptophan and tyrosine bound to p27-KID. However, even when titrated to a 30-fold molar excess, neither aromatic amino acid cause chemical shift perturbations in 2D HSQC spectra upon titration into p27-KID (Suppl. Fig. 4). This is most likely because they lack the specific chemical features embodied in the pharmacophore models for the two groups of fragment hits.

We used the Group 2 molecule, SJ403 (Suppl. Table 1b), to test our hypothesis that molecules that bind to p27 can alter its Cdk regulatory function. For these experiments, because SJ403 binds weakly to p27-KID, we studied its ability to modulate the binding of p27-D2 to Cdk2/cyclin A (Kd value binding to Cdk2/cyclin A, 73 ± 8 nM; Fig. 4a versus 5 nM for p27-KID20). We first used fluorescence anisotropy (FA) to monitor displacement of a single-cysteine (Cys) mutant of p27-D2, with arginine 93 mutated to Cys (R93C) and labeled with Alexa Fluor 488 (p27-D2-FL), from Cdk2/cyclin A by SJ403. Titration of SJ403 caused concentration-dependent reduction of FA of p27-D2-FL, with an IC50 value of 475 ± 67 μM (Fig. 4b), suggesting that SJ403 displaced p27-D2 from Cdk2/cyclin A.

Fluorescence anisotropy analysis of displacement of p27-D2-FL from Cdk2/cyclin A by SJ403.

(a) Fluorescence anisotropy analysis of Cdk2/cyclin A binding to p27-D2-FL. The starting condition (20 nM p27-D2-FL, 300 nM Cdk2/cyclin A) for assessing the effect of fragment hit SJ403 on p27 function is represented as the red circle on the binding isotherm. (b) Titration of SJ403 (0–3 mM) caused the displacement of p27-D2 from Cdk2/cyclin A, resulting in a decrease in fluorescence anisotropy values. Experiments were performed in triplicate and average values and the standard deviations of the means are shown.

We next used NMR spectroscopy to monitor the displacement of 2H/13C/15N-labeled p27-D2 from Cdk2/cyclin A by SJ403; the complex of p27-D2 with Cdk2/cyclin A (100 μM) was prepared with a slight excess of 2H/13C/15N-p27-D2 (mole ratio, 1.1:1.0 p27-D2:Cdk2/cyclin A). Peaks for unbound p27-D2 dramatically increased in intensity in the presence of SJ403, while the resonances for p27-D2 bound to Cdk2/cyclin A were reduced in intensity (Fig. 5 and Suppl. Fig. 7a–c), consistent with partial displacement of p27-D2 from Cdk2/cyclin A. A ratiometric method was used to analyze the 2D TROSY-HSQC NMR spectra. In this method, the relative population of the free state resonances (pf; for the resonance of each p27-D2 residue) was determined as a fraction of the total intensity of both free and bound resonances for a given residue. These values were compared by forming the ratio of the two relative free state populations (pfw/SJ403/pfw/o SJ403); values greater than 1 indicated compound-dependent displacement. The relative populations for the bound state (pb) were also determined for samples without (pbw/o SJ403) and with SJ403 (pbw/SJ403). In this case, ratios less than 1 indicated displacement of p27-D2 from Cdk2/cyclin A. Furthermore, the chemical shift values of the free state resonances reflected binding of p27-D2 to SJ403 (Suppl. Fig. 7d,e), further supporting the conclusion that SJ403 shifts the binding equilibrium between p27-D2 and Cdk2/cyclin A to increase the population of unbound p27-D2 through protein:small molecule interactions.

2D 1H-15N TROSY-HSQC NMR analysis of displacement of 15N-p27-D2 from Cdk2/cyclin A by SJ403.

2D 1H-15N TROSY-HSQC NMR spectra of p27-D2/Cdk2/cyclin A in the absence (control) and presence of SJ403 (3 mM), respectively, were recorded. For the residues reported, both free and bound p27-D2 resonances were observed in the TROSY-HSQC spectra. An increase in the relative intensity of free p27-D2 resonances was observed in the TROSY-HSQC spectrum of p27-D2/Cdk2/cyclin A in the presence of SJ403, with respect to the control spectrum (a), while the resonances of p27-D2 bound to Cdk2/cyclin A showed a decrease in relative intensity (b). This permits the calculation of populations of free (f) and bound (b) p27-D2, in the absence (w/o SJ403) and presence of SJ403 (w/SJ403)). The pfw/o SJ403 values represent the relative populations of resonances for free p27-D2 before addition of SJ403 and the pfw/SJ403 values are the relative populations of free state p27-D2 resonances after SJ403 addition. Correspondingly, the pbw/o SJ403 values represent the relative populations of resonances for p27-D2 bound to Cdk2/cyclin A before addition of SJ403 and the pbw/SJ403 values are the relative populations of bound state p27-D2 resonances after SJ403 addition. The mean values and standard deviations of the means for triplicate measurements are shown; one set of triplicate 2D 1H-15N TROSY-HSQC spectra are illustrated in Suppl. Fig. 7.

The FA and NMR results indicated that SJ403 partially displaced p27-D2 from Cdk2/cyclin A; therefore, we next investigated whether displacement modulated the kinase activity of Cdk2. p27-D2, while lacking the D1 and LH sub-domains found in p27-KID, is still a moderately potent inhibitor of Cdk2/cyclin A (IC50 value, 152 ± 39 nM; Suppl. Fig. 9a). Under conditions wherein Cdk2/cyclin A was almost completely inhibited by p27-D2 (~13% of full kinase activity was maintained), titration of SJ403 over the concentration range shown in FA and NMR experiments to cause p27-D2 displacement led to increased kinase activity (from 13% to ~20%, a >50% increase; Fig. 6a and Suppl. Fig. 9b), consistent with partial displacement of p27-D2 from Cdk2/cyclin A by SJ403. Furthermore, in the absence of p27-D2, SJ403 substantially inhibited Cdk2/cyclin A activity at concentrations that, in the presence of p27-D2, were associated with p27-D2 displacement and increased kinase activity (Fig. 6b and Suppl. 9c). Thus, we conclude that the primary effect of SJ403 is displacement of p27-D2 from Cdk2/cyclin A (through the binding of SJ403 to p27-D2) and partial restoration of kinase activity, even as a secondary effect of SJ403 is to inhibit kinase activity (through the binding of SJ403 to Cdk2/cyclin A). These results provide proof-of-principle that a small molecule (SJ403) inhibits the function of a disordered protein (p27-D2) through sequestration in a conformation incapable of binding and inhibiting Cdk2/cyclin A.

The small molecule SJ403 displaces p27-D2 from Cdk2/cyclin A and partially restores Cdk2 kinase activity.

Analysis of Cdk2 kinase activity [within the Cdk2/cyclin A complex (200 pM) using Histone H1 as the substrate] in the presence of p27-D2 (200 nM) and varied concentrations of SJ403 (from 10 μM to 3 mM) (a) and varied concentrations of SJ403 alone (b). Cdk2 is partially active under the starting conditions in (a) and activity is enhanced through addition of SJ403. In the absence of p27-D2 (b), addition of SJ403 is associated with Cdk2 inhibition. All experiments were performed in triplicate and average values and the standard deviations of the means are shown.

We previously showed based on NMR and molecular dynamics (MD) computational data that p27-KID exhibits different types of partially populated secondary structure in the free state, including helical structure within the LH sub-domain, a β-hairpin conformation within sub-region D2.1 and helical structure within sub-region D2.3 of sub-domain D2 (ref. 21 ; see Fig. 1a for sub-domain/sub-region nomenclature). Interestingly, the Group 1 and 2 molecules variably interacted with two of these partially structured regions of p27-KID (sub-regions D2.1 and D2.3), but also with sub-region D2.2, which is highly dynamic and does not exhibit persistent secondary structure21 (Fig. 7a–c); these interaction sites are all within the Cdk2-binding D2 sub-domain. Notably, the p27 binding molecules did not interact with sub-domains D1 or LH. Importantly, each of the regions within p27-KID that interacted with small molecules contained several aromatic amino acids (Fig. 7a). In fact, sub-domain D2 contains eight of the nine aromatic amino acids found within p27; a ninth, F33, is found within sub-domain D1 but did bind to small molecules. Group 1 molecules caused the largest CSPs for residues F87, Y88 and Y89 within sub-region D2.3 and Group 2 molecules perturbed resonances for these residues as well as for residues W60, N61 (sub-region D2.1) and E75, W76 and Q77 (sub-region D2.2). Phenylalanine (F) or tyrosine (Y) residues flank the tryptophan (W) residues, but exhibited smaller CSP values in the presence of Group 2 compounds than did the W residues, suggesting that the small molecules preferentially perturbed the electronic environment of the indole rings of these residues. We tested this hypothesis through mutagenesis of W60 or W76, or both, within p27-KID to either F or alanine (A) and used 2D NMR to map binding to SJ403. The singly mutant p27-KID variants exhibited patterns of CSPs similar to those observed for wild-type p27-KID except that perturbations near the mutated W residues were absent (Suppl. Fig. 5a,b,d,e). However, the two doubly mutated p27-KID variants exhibited only very weak CSPs within the F87YY89 aromatic cluster (Suppl. Fig. 5c,f). Together, these results indicated that each of the two W residues, W60 and W76, contributed substantially to binding to SJ403. Mutations of F87, Y88 or Y89 to A were associated with substantially reduced CSPs within the Group 1 compound binding region of p27-KID (Suppl. Fig. 5j–l), whereas mutation of R90 to A (Suppl. Fig. 5m) had little effect on the binding, implying that the clustering of the aromatic side chains is critical to the small molecule:protein interaction. With the Y88 to A p27-KID mutant, the Group 2 compound, SJ403, caused CSPs near the two W residues (W60 or W76) but not within the F87A88Y89 region of p27-KID (Suppl. Fig. 5h). However, with the F87 to A and Y89 to A mutants, SJ403 caused CSP patterns similar to that observed with wild-type p27-KID (Suppl. Fig. 5g,i), suggesting that Y88 contributes to a greater extent to interactions with SJ403 (and possibly other Group 2 compounds) than does Y89.

Small molecules bind to clusters of aromatic residues with the D2 sub-domain of p27-KID.

(a) The amino acid sequence of p27-D2 showing sub-domain D2 and the individual small molecule binding regions (labeled D2.1, D2.2 and D2.3). The aromatic amino acids within these regions are illustrated in bold. These represent eight of the nine aromatic residues within p27-KID. (b,c) Average 1H chemical shift perturbation values for all Group 1 (b) and Group 2 (c) compounds illustrating preferences for interactions with Y (b) and W and Y (c) residues, respectively.

The two W residues and Y88 of p27 are conserved in the related disordered cell cycle regulatory protein, p21Waf1/Cip1 42 (p21; Fig. 8a), allowing their dominance in interactions with Group 1 & 2 compounds to be further tested (Fig. 8b,c). The Group 2 compound, SJ403, interacted with residues near W49 and W65 of p21 (within the p21 kinase inhibitory domain, p21-KID) but not with Y77 (homologous to Y88 of p27; Fig. 8c), supporting the importance of W residues for interactions with this small molecule but also suggesting that clustering of at least two Y and F residues is required for additional interactions. The two leucine (L) residues flanking Y77 in p21 (L76 and L78) do not substitute for F87 and Y89 in p27. Similarly, the L76YL78 region of p21 does not support binding to a Group 1 compound (SJ319843; Fig. 8b). These findings with p21-KID further strengthen our explanation of the molecular basis for interactions of Group 1 and 2 compounds with p27.

Interaction of p21-KID with representative p27 fragment hits.

(a) Comparison of amino acid sequences of the D2 subdomains of p27-KID and p21-KID; (b,c) Individual 1H (ΔδHN) chemical shift perturbations (through analysis of 2D 1H-15N HSQC NMR spectra) caused by addition of SJ319843 (b, Group 1) and SJ403 (c, Group 2), respectively, to 2H/15N-p21-KID (20 μM). The threshold for identifying specific binding interactions with the p27-KID residues was defined as two standard deviations above the average of chemical shift perturbations (Δδave+2σ, represented as the dotted black line) and considered if perturbations were larger than the experimental 1H digital resolution (shown as dotted magenta line).

To gain insight into the structural features of the small molecule binding sites within p27, we performed molecular dynamics (MD) simulations (over 400 ns) with p27-D2. The results recapitulated past MD findings (over 100 ns) for the longer p27-KID construct21 which revealed a transient β-hairpin involving residues W60-F64 and H67-L70 (within sub-region D2.1) and two α-helical turns involving residues E80-G82 and F87-Y89 (within sub-region D2.3). In the new MD trajectory, these secondary structures were stable on short nanosecond time-scales but unfolded and refolded over the longer time periods sampled in this experiment. These longer time-scale transitions were coupled with transient formation of hydrophobic clusters involving residues within sub-regions D2.1 and 2.2 (containing W60 and W76) and sub-region D2.3 (containing Y88; Fig. 9a) which gave rise to extended and compact conformations (labeled E and C in Fig. 9b). Representative extended and compact conformers are presented in Fig. 9c and d, respectively and all others in Suppl. Fig. 10a,b. (Note that conformers are defined as compact if two of the three aromatic residues of p27 critical for binding to small molecules, W60, W76 and Y88, are less than 20 Å from each other.) We speculate that the extended and compact conformations create binding sites for Group 1 and Group 2 compounds, respectively. Our interpretation of these results is that Group 1 compounds bind to the F87YY89 motif within sub-region D2.3 in p27-D2 molecules in which Y88 is far away from W60 and W76 (>20 Å) and that Group 2 molecules bind to compact conformations when at least two of the three critical aromatic residues within the different sub-regions (sub-domains D2.1, D2.2 and D2.3) are clustered. Interestingly, analysis of the MD trajectory showed that Y88 and either W60 or W76 were frequently in close contact but that all three residues were rarely in close proximity (Suppl. Fig. 10c). This suggested that there are several different conformations with clustered aromatic residues (in particular, Y88 and either W60 or W76) capable of binding to Group 2 compounds, consistent with mutagenesis results showing that either W60 or W76, but not both, are dispensable for Group 2 compound binding (Suppl. Fig. 5a–f). In summary, the new MD results for p27-D2 suggest strongly that transient conformational fluctuations that create and disrupt clusters of aromatic residues modulate the binding of Group 1 and Group 2 small molecules to p27-D2. These results are consistent with the identification of W60, W76 and Y88 by NMR as the principal sites for compound binding and with results showing that binding is altered through mutation of these residues.

Distance profiles of W60, W76 and Y88 Cβ atoms monitored as a function of time from MD simulations reveal extended and compact states in the p27-D2 conformational ensemble.

(a) Distance between W60-W76 (left), W60-Y88 (center) and W76-Y88 (right) Cβ atoms monitored over time from MD with gray lines representing the instantaneous distance (every 2 ps) and red lines representing the time average (every one ns). The times at which representative extended (E) and compact (C) conformers were extracted for display in Suppl. Fig. 10 are indicated. (b) Illustration of the distributions of the distances in (a) shown as histograms (red lines). The distributions show the existence of two distinct states for W60-Y88 (center) and W76-Y88 (right) indicating compact (C) and extended (E) conformations (see Suppl. Fig. 10). (c,d) Representative extended (c) and compact (d) conformations, respectively, showing relative positions of residues of p27-D2 involved in interaction with fragment hits.

Drug discovery against folded proteins often involves identification of compounds that bind to sites that are naturally utilized for interactions with small molecule or macromolecular ligands. These types of binding sites are temporally stable and enable specific and tight interactions with chemically complementary small molecules. In contrast, disordered proteins (or disordered protein regions) exhibit dynamic and heterogeneous conformations that do not display similar, temporally stable small molecule binding sites. However, despite the lack of temporally stable feature, many disordered proteins/regions interact with macromolecular partners through the process of folding upon binding. We hypothesized that the ability of a disordered protein to bind other proteins would also confer the ability to bind small molecules and tested this idea through the studies of p27 described herein.

The entire p27 protein is disordered and its N-terminal domain (p27-KID) becomes ordered upon binding to Cdk2/cyclin A. A short linear motif within the D1 sub-domain of p27-KID (with the sequence R30NLFG34; L, Leucine; Suppl. Fig. 11a,b) binds to a pocket on the surface of cyclin A. Sub-domain LH, which links sub-domains D1 and D2, contacts the surfaces of cyclin A and Cdk2 but contributes little to the binding energy43. Sub-domain D2 of p27 (p27-D2; 34 amino acids in length) adopts extensive secondary structure and makes extensive hydrophobic interactions with Cdk2 upon binding (Suppl. Fig. 11a,c and d). In addition, numerous hydrogen bonds form between residues within the D2 sub-domain and Cdk2. While p27-D2 exhibits 11 hydrophobic and aromatic residues (including I, L, V, F, W and Y residues), many of which contribute to interactions with Cdk2, this sub-domain does not independently form a stable hydrophobic core.

Despite extensive disorder in the unbound state, we identified two groups of small molecules that bound with exquisite specificity albeit with low affinity to two overlapping regions within p27. These molecules bound to regions containing aromatic rings, with preference for W and Y residues. Molecules in Group 1 bound to a localized region, F87YY89, while molecules in Group 2 bound to this region as well as two others containing two W residues (W60 and W76). Strikingly, molecules within Groups 1 (25 molecules) and 2 (14 molecules), respectively, were chemically similar, demonstrating chemical structure-binding activity relationships for these two types of small molecule: protein interactions. We refer to this phenomenon as "fuzzy structure activity relationships (SAR)" due to the "fuzzy" character of the respective interaction potential maps (Fig. 3a,b).

These chemical features of the p27 binding molecules can be rationalized based upon the features of residues within and flanking the NMR-identified binding sites within p27-D2. Group 1 and 2 molecules exhibited two or three heterocyclic aromatic rings and exhibited the potential for these rings to participate in hydrophobic interactions (Fig. 3a,b; gold polygons). Molecules in both groups also displayed the potential to bind electropositive moieties (Fig. 3a,b; cyan polygons), consistent with the F87YY89 region in the Group 1/2 binding sites being flanked on the C-terminal end by R90PPR93 and W60 in the Group 2 binding site being flanked on the N-terminal end by R58K59. Residue W76 within the Group 2 binding site is flanked by amino acids with both electronegative and electropositive, as well as polar, features (E71GKYEW76QEVEK81) although the potential for interactions with electronegative moieties was only weakly represented (Fig. 3b; red polygons). In addition to electrostatic features, the small molecules, which exhibit hydrogen bond donors and acceptors, may achieve specificity through transient hydrogen bonds with complementary groups within p27-D2, as observed for small molecules binding to Myc44. We argue that these and other currently unappreciated features of the p27-D2 polypeptide chain create the potential to specifically bind Group 1 and 2 small molecules.

What are the conformational features of small molecule binding sites within p27? Our MD computations with p27-D2 revealed dynamic fluctuations of pairwise distances between the aromatic residues (specifically, W60, W76 and Y88) within the Group 2 binding sites (Fig. 9). While W60 and W76 sampled a relatively narrow distance range (16.5 ± 2.3 Å), the distances between W60 and Y88 and W76 and Y88, fluctuated over a wider range (20.5 ± 5.5 Å and 19.6 ± 4 Å, respectively). Furthermore, these latter residue pair distances each exhibited two discrete populations that we term compact (C) and extended (E; Fig. 9b, middle and right panels). We propose that the compact conformations create binding sites for Group 2 compounds, consistent with NMR CSP patterns (Fig. 7c). In addition, we propose that extended conformations favor the binding of Group 1 compounds to F87YY89 region, also consistent with NMR CSP patterns (Fig. 7b). W60 is flanked by F62 and F64 and W76 by Y74, but mutagenesis results (Suppl. Fig. 5a–f) suggested that the W residues are the principal determinants of Group 2 compound binding. This is probably due to the potential of side chains of W and Y residues (but not of F residues) to form hydrogen bonds with small molecules, in addition to participating in hydrophobic and π-stacking interactions. Interestingly, distance correlation analysis (Suppl. Fig. 10c) indicated that, in the compact conformations, Y88 is most frequently close to either W60 or W76 but rarely close to both W residues. This suggested that close proximity of Y88 and one or the other of the two W residues created binding sites within p27-D2 for Group 2 compounds. This observation is also consistent with mutagenesis results, which showed that mutation of both W residues, but not of individual W residues, abrogated binding to Group 2 compounds (Suppl. Fig. 5c, f versus a, b, d, e). In summary, p27-D2 exhibited transient close contacts between W60 or W76 and Y88, which we propose creates binding sites for Group 2 compounds. Additionally, when neither W residue and Y88 are close, the spatial proximity of the three aromatic residues in the F87YY89 region created binding sites for Group 1 compounds. It is interesting that, in the Cdk2/cyclin A-bound structure of p27-KID34, the five aromatic residues within sub-regions D2.1 and D2.2 (that bind to Group 2 compounds) are in close proximity (but separated from the F87YY89 region; Suppl. Fig. 11a,c and d). The MD results show that, in the absence of Cdk2/cyclin A, subsets of these eight aromatic residues transiently interact, sometimes creating binding sites for different types of heterocyclic aromatic small molecules (Group 1 or 2).

The compound, SJ403, was demonstrated to sequester p27-D2 away from Cdk2/cyclin A and activate Cdk2 catalytic activity, effectively fulfilling our goal of inhibiting p27's cell cycle inhibitory function. While the affinity of Group 1 and 2 molecules is low, they do exhibit high specificity for particular regions of p27. As discussed above, residues within these regions otherwise engage Cdk2. Remarkably, the >2,300 compounds that were screened failed to bind other regions of p27 (sub-domains D1 and LH), suggesting that these other regions lack a sufficient density of aromatic residues (specifically, W and Y residues) to specifically recognize small heterocyclic aromatic molecules. Sub-domain LH, in isolation, does not bind to Cdk2/cyclin A but sub-domain D1 binds cyclin A with high affinity (Kd, 42 nM)45. We speculate that the RxLFG motif within this latter region, due to its limited length, cannot adopt conformations that create binding pockets for small molecules, as is possible for the much longer D2 sub-domain. However, the low affinity of the Group 1 and 2 compounds for p27 limits their applicability toward our broader goal of modulating p27 function in cells and, ultimately, humans. How can the affinity of small molecules for p27 be increased? We propose that the Group 1 and 2 molecules cause a degree of conformational restriction within p27-D2 and that molecules that enhance this restriction will exhibit higher affinity. We envision that small molecules with greater "three-dimensionality", that present chemically diverse and complex features, will be better templates for binding and sequestering p27. Efforts based on two strategies are underway to optimize our fragment hits using synthetic chemistry. First, we are "growing" the Group 1 and 2 scaffolds by introducing diverse chemical moieties at various positions on the heterocyclic ring systems to enable additional interactions with residues near W60, W76 and Y88 within p27-D2. Second, when the growing experiments are complete, we will synthetically "link" the optimal Group 1 and Group 2 molecules to further enhance binding to p27-D2. The results of these future experiments will indicate whether synthetic strategies for compound optimization that have emerged from structure-based drug discovery can be applied to a disordered protein. In conclusion, we have discovered small molecules with "fuzzy SAR" that mediate specific binding to and inhibition of p27, demonstrating the potential to rationally "drug" disordered protein targets in the future.

The p27 constructs were expressed in E. Coli with an N-terminal 6xHis affinity tag after sub-cloning into pET28a (Novagen) using established procedures20. This included p27-KID (residues 22-105 of human p27) and the following mutants: W60A, W60F, W76A, W76F, W60A-W76A, W60F-W76F, F87A, Y88A, Y89A and R90A. p27-D2 (residues 58-105 of human p27) and the mutant, C99S-R93C, were expressed similarly. Isotope-labeled proteins (15N, 13C/15N and 2H/13C/15N) were expressed in a MOPS-based minimal media using established procedures22. All p27 constructs were purified by nickel affinity chromatography, digested with thrombin to remove the 6xHis tag and further purified using reverse-phase high performance liquid chromatography (HPLC) using a C4 column (Vydac) and 0.1% trifluoroacetic acid-containing water/acetonitrile solvent system. Protein concentrations were determined by UV absorbance at 280 nm using a molar extinction coefficient of 15,470 M−1 cm−1 for p27-KID, p27-KID-F87A, p27-KID-R90A, p27-D2 and p27-D2-C99S-R93C; 13980 M−1cm−1 for p27-KID-Y88A and p27-KID-Y89A; 9,970 M−1 cm−1 for p27-KID variants with a single tryptophan residue; and 4,470 M−1 cm−1 for p27-KID variants without a tryptophan residue. Full length human Cdk2, active Cdk2 (phosphorylated at threonine 160), truncated human cyclin A (residues 173–432) and p21-KID were expressed and purified using established protocols20,46.

The Cdk2/cyclin A complex was prepared by equilibrating Cdk2 and cyclin A (1:1 mole ratio) for 1 hour at 4 °C followed by purification using size exclusion chromatography (S75 resin, GE Healthcare, Piscataway, NJ) in buffer containing 20 mM Tris, pH 7.5, 200 mM NaCl, 5 mM TCEP). The ternary complex was prepared by equilibrating 2H/13C/15N-p27-D2 with Cdk2/cyclin A(1.1:1 mole ratio) at 4 °C overnight followed by purification using size exclusion chromatography (S200 resin, GE Healthcare, Piscataway, NJ) in the same buffer as the Cdk2/cyclin A complex. For NMR experiments, the 2H/13C/15N-p27-D2/Cdk2/cyclin A complex was buffer exchanged into 20 mM Na phosphate, pH 6.5, 200 mM NaCl, 10 mM DTT-D10 and 10% 2H2O.

The fragment library used to screen for binding to p27-KID was composed of two collections that were designed according to different criteria: (a) 1,100 fragments were purchased from the Maybridge Ro3 Diversity Fragment Library (‘Maybridge’)47 and (b) 1,222 fragments were purchased from Enamine (N = 250) and Life Chemicals (N = 972) using an in-house algorithm (‘In-House’).

The Maybridge fragment collection was designed to provide broad coverage of chemical space for fragment-based drug discovery. Each fragment satisfies Congreve's Rule of Three40: (a) molecular weight < 300; (b) number of hydrogen bond donors ≤ 3; (c) number of hydrogen bond acceptors ≤ 3; and (d) clogP ≤ 3. All compounds have experimentally determined equilibrium solubility ≥200 mM (DMSO) and ≥1 mM (PBS) and confirmed purity ≥95% (based upon analysis by liquid chromatography/mass spectrometry) and are free of reactive or toxic functional groups. The 1,100 fragments constitute a ‘core’ set that encompasses the chemical diversity of the entire collection (1,823 fragments).

The In-House fragment collection consists of 1,222 commercially-available compounds selected using a custom algorithm designed to identify structurally complex, low molecular weight molecules with scaffolds that were well-sampled within the separate St. Jude high-throughput screening library (HTS library, >500,000 compounds; see below). First, commercial fragment collections (subsets of larger diversity collections filtered for ‘fragment-like’ characteristics) were filtered to remove molecules containing inorganic atoms, isotopes, or invalid structures and to remove molecules that were not available in sufficient quantity (<50 mg). Passing molecules were abstracted to Murcko scaffolds using Pipeline Pilot (‘Generate Fragments’ component in Accelrys v. 8.5 with alpha atoms preserved, see ref. 48 for the general method). These scaffolds were further filtered according to the following rules: number of reactive substructures = 0 (‘REOS’ filters49,50,51, number of rotatable bonds <= 3, number of heavy atoms >= 10, number of rings >= 1 and number of ring substitutions >1 for single ring systems and number of molecules present in the St. Jude HTS library containing the scaffold >= 8. Molecules containing these scaffolds were identified in the commercial fragment libraries and then prioritized for purchase according to highest Oprea complexity52. This library has the following average calculated physicochemical properties: MW = 246 ± 39 Da, number of atoms = 17 ± 3, log P = 1.7 ± 1.0, polar surface area = 63 ± 19 Å2, number of H-bond acceptors = 4.3 ± 1.4 and number of H-bond donors = 1.3 ± 0.9. The distributions of these and other chemical features of the two fragment libraries are summarized in Suppl. Fig. 12a.

Consensus field maps for Group 1 and 2 molecules were defined using the FieldTemplater module in Forge (v10, Cresset, REF= The FieldTemplater module takes as input a set of reference molecules (the active compounds for Groups 1 and 2 identified in the first two rounds of screening; Suppl. Table 1A, B) and seeks to align them in order to maximize overlap between their interaction fields. Reference molecules are first aligned using field points, followed by slower, more accurate optimization of the alignment using the full interaction field. Prior to alignment, conformations were generated using the ‘Very Accurate and Slow’ conformation hunt option in Forge and default settings. The similarity scoring function used in the alignment was based on 100% field similarity and 0% shape similarity in order to maximize the topological diversity of molecules retrieved using the model in the next round of screening.

The consensus field maps for Group 1 and 2 molecules were used to query a database of 10,455 fragment-like molecules from ChemDiv ( using the FIeldScreen module in Forge. Of the 215 molecules with the highest field similarity scores, 184 compounds were purchased. However, 106 of these were poorly soluble under our assay conditions and not screened. From the remaining 78 compounds, 12 additional hits were identified using 2D 1H-15N HSQC NMR (see below).

Screening and validation NMR experiments were performed at 298 K (25 °C) using either a Varian Inova 600 MHz spectrometer equipped with a triple resonance (HCN) room temperature gradient probe or a Bruker Avance 600 MHz spectrometer equipped with TCI cryogenic gradient probe and a SampleJet sample changer. Fragment molecules were initially analyzed as pools of five molecules dissolved at 10 mM each in DMSO-D6. The fragment pools contained in 96-well plates were mixed using a Gilson 215 liquid handler with buffer (20 mM Na phosphate, pH 6.5, 20 mM NaCl, 10% 2H2O, 5 mM DTT-D10) or buffer containing the p27-KID protein (10 μM p27-KID) to give final compound concentrations of 200 μM each. For initial fragment screening experiments, one-dimensional (1D) 1H- and WaterLOGSY35 NMR spectra were recorded for compound pools without and with protein. Pools exhibiting hits were deconvoluted by analyzing pure compounds using 1D 1H and validated by two-dimensional (2D) heteronuclear NMR experiments (2D 1H-15N HSQC titrations) using the Bruker Avance 600 MHz spectrometer. NMR samples used for 2D 1H-15N HSQC titrations contained 100 μM 15N-labeled p27 protein (p27-KID, p27-KID mutants, or p27-D2) in 20 mM Na phosphate, pH 6.5, 200 mM NaCl, 10% 2H2O, 5 mM DTT-D10; compounds dissolved in DMSO-D6 were titrated to the desired concentrations. DMSO-D6 was added to maintain a constant concentration (2% vol/vol). A spectral resolution of 3.5 and 5.7 Hz was achieved in the 1H and 15N dimensions, respectively. Tryptophan and tyrosine, respectively, were titrated into 15N-p27-KID (100 μM) up to 3 mM. Three-dimensional (3D) backbone triple-resonance experiments to establish resonance assignments for the p27 constructs were performed using the Bruker Avance 600 MHz spectrometer. Assignments for p27-KID and p27-D2 are illustrated in Suppl. Fig. 12b,c and d,e, respectively. 2D 1H-15N TROSY-HSQC NMR experiments with 2H/13C/15N-p27-D2/Cdk2/cyclin A and SJ403 have been recorded at 308 K (35 °C). 2D 1H-15N HSQC experiments of representative fragment hits (1 mM of SJ319843, group 1 and SJ403, group 2, respectively) with 15N-p21-KID (20 μM) were recorded at 298 K using a Bruker Avance 600 MHz spectrometer. NMR spectra were processed using Bruker Topspin software and analyzed using computer-aided resonance assignment (CARA) software53.

Chemical shift perturbations are generally quantified as combined 1H and 15N chemical shift values. However, the analysis of primary data for small molecules binding to p27 showed that the largest CSPs were for the 1H dimension, which were statistically significant and that the corresponding 15N CSP values were often small and not significant. Thus, the use of a combined chemical shift values would mask the effects of compound binding. We also rigorously considered the magnitude of CSPs relative to the experimental digital resolution of the 1H and 15N dimensions of HSQC spectra and a threshold defined as the average CSP value plus two times the standard deviation of the mean (Δδave + 2σ). Thus, the assessment of statistical significance was based upon whether a particular CSP value was larger than, i) the experimental spectral resolution in the given dimension (1H or 15N) and ii) the quantity, Δδave + 2σ, for that dimension.

Sixteen point titrations of representative fragment hits from Group 1 and Group 2 (SJ572710 and SJ572403, respectively) into 15N-p27-KID (100 μM) were recorded using an adopted "in-phase" 2D 1H-15N HISQC with proton decoupling during 15N chemical shift labeling achieved with a WALTZ16 composite pulse with an amplitude of 7.1 kHz54,55. All experiments were collected at a Bruker Avance 800 MHz spectrometer equipped with TCI cryogenic gradient probe. The following molar ratios of 15N-p27-KID (100 μM) to inhibitor were used: 1:0, 1:0.1, 1:1.5, 1:3, 1:4.5, 1:6, 1:7.5, 1:9, 1:10.5, 1:12, 1:15, 1:18, 1:21, 1:24, 1:27 and 1:30. Each spectrum was recorded with 256 (t1,max = 31.0 ms) and 1024 (t2,max = 106.5 ms) complex points in the 15N and 1H dimensions, respectively, with eight transients collected per point. The spectral resolution of the 1H and 15N dimensions was 2.4 and 1.7 Hz, respectively. Chemical shift perturbations throughout the entire titration had to be greater than this resolution threshold to be considered for further analysis. The data was processed using the NMRPipe package56 and analyzed using in-house scripts written in Python using the Scipy computing libraries and Mathematica (Wolfram Research). In order to alleviate human bias in the peak position determination the automatic peak picking function in NMRpipe was utilized in which a spectral window was assigned for each resonance over the entire trajectory of the titration for a given chemical shift. The error in the peak position for a given resonance was taken as the 15N or 1HN line width divided by the signal-to-noise ratio. All resonances that exhibited chemical shift perturbations greater than the spectral resolutions and Δδave + 2σ were subsequently grouped and fitted globally for their respective maximum chemical shift difference and a global dissociation constant (Kd). The error in the Kd values was determined using a Monte-Carlo approach in which an error of 10% was imposed on the ligand concentration and the error in the chemical shift's peak position was considered.

The p27-D2-C99S-R93C mutant was conjugated with Alexa Fluor488-C5-Maleimide (Life Technologies, p27-D2-FL) in buffer containing 20 mM Na phosphate, pH7.3, 20 mM NaCl according to the manufacturer's protocol. The conjugated protein was further purified by reverse-phase high performance liquid chromatography (HPLC) using a C4 column (Vydac) and 0.1% trifluoroacetic acid-containing water/acetonitrile solvent system. Lyophilized HPLC fractions were resuspended in buffer containing 20 mM Na phosphate, pH 7, 200 mM NaCl, 3.5 mM TCEP. Fluorescence anisotropy measurements were performed at 25 °C on a Horiba Fluorolog 3 spectrofluorometer. Briefly, p27-D2-FL (20 nM) was mixed with Cdk2/Cyclin A (300 nM) and added to the required amount of lyophilized compound (SJ403). All samples have been incubated overnight in the dark at 4 °C prior to fluorescence measurements. Fluorescence anisotropy binding data was analyzed using KaleidaGraph. Curve fitting was performed on the average of three independent experiments by a non-linear regression binding model.

Cdk2/Cyclin A (200 pM) was mixed with Histone H1 (15 μM; EMD Millipore), p27-D2 (200 nM) and varied amounts of lyophilized SJ403 and incubated overnight at 4 °C. The effect of the SJ403 compound on Cdk2 activity was determined by performing similar experiments in the absence of p27-D2. Subsequently, ATP (50 μM total concentration, of which 6 μCi γ 32P-ATP (PerkinElmer, Inc)) was added to each reaction and further incubated for 35 minutes at 30 °C. Each reaction had a total volume of 20 μL. The sample buffer contained 20 mM HEPES pH 7.3, 25 mM sodium β-glycerolphosphate, 15 mM MgCl2, 16 mM EGTA, 0.5 mM Na3VO4 and 10 mM DTT. The reactions were quenched by addition of SDS-gel loading buffer (5 μL) and then analyzed by SDS-PAGE (10 μL). The gels were dried at 70 °C under vacuum and a phosphoimager (GE Healthcare, Piscataway, NJ) was used to quantify the 32P-Histone H1 bands. IC50 values were determined by curve fitting using KaleidaGraph software. Experiments were performed in triplicate and mean IC50 and standard deviations of the mean values are reported.

All-atom MD simulations using graphics processing unit (GPU)-optimized AMBER 12 software57 were used to explore the conformational landscape of the p27-D2 domain. The conformation of p27-D2 within the p27-KID/Cdk2/cyclin A structure (PDB ID, 1JSU) was used as the starting structure from MD computations with amino acid protonation states modified to reflect pH 7.0. The structure was placed in a rectangular box of TIP3P water that was 15 Å larger on all sides than the p27-D2 molecule. In addition to neutralizing the system by adding counter ions, 20 mM NaCl was added to mimic the experimental conditions in our simulations.

The system was equilibrated and stabilized using multi-step energy minimization at 298 K, as described previously58. All production runs were performed using the constant number of particles, volume and energy (NVE) ensemble with periodic boundary conditions. The particle-mesh Ewald (PME) method was used for electrostatic interactions and a 10 Å cut-off was used for Lennard-Jones interactions. The SHAKE algorithm was used to restrict the motions of all covalently bonded hydrogen atoms. Simulations were performed at 298 K and 1 atm pressure. The total time-scale for our simulation was 0.4 microseconds with snapshots being stored every 2 ps, resulting in a total of 200,000 snapshots from the trajectory.

How to cite this article: Iconaru, L. I. et al. Discovery of Small Molecules that Inhibit the Disordered Protein, p27Kip1. Sci. Rep. 5, 15686; doi: 10.1038/srep15686 (2015).

van der Lee, R. et al. Classification of intrinsically disordered regions and proteins. Chem Rev 114, 6589–6631 (2014). Pmc4095912.

Article CAS Google Scholar

Gsponer, J., Futschik, M. E., Teichmann, S. A. & Babu, M. M. Tight regulation of unstructured proteins: from transcript synthesis to protein degradation. Science 322, 1365–1368 (2008).

Article CAS ADS Google Scholar

Vavouri, T., Semple, J. I., Garcia-Verdugo, R. & Lehner, B. Intrinsic protein disorder and interaction promiscuity are widely associated with dosage sensitivity. Cell 138, 198–208 (2009).

Article CAS Google Scholar

Druker, B. J. et al. Effects of a selective inhibitor of the Abl tyrosine kinase on the growth of Bcr-Abl positive cells. Nature medicine 2, 561–566 (1996).

Article CAS Google Scholar

Park, C. M. et al. Discovery of an orally bioavailable small molecule inhibitor of prosurvival B-cell lymphoma 2 proteins. J Med Chem 51, 6902–6915 (2008).

Article CAS Google Scholar

Souers, A. J. et al. ABT-199, a potent and selective BCL-2 inhibitor, achieves antitumor activity while sparing platelets. Nature medicine 19, 202–208 (2013).

Article CAS Google Scholar

Krishnan, N. et al. Targeting the disordered C terminus of PTP1B with an allosteric inhibitor. Nat Chem Biol 10, 558–566 (2014). Pmc4062594.

Article CAS Google Scholar

Erkizan, H. V. et al. A small molecule blocking oncogenic protein EWS-FLI1 interaction with RNA helicase A inhibits growth of Ewing's sarcoma. Nat Med 15, 750–756 (2009). Pmc2777681.

Article CAS Google Scholar

Selvanathan, S. P. et al. Oncogenic fusion protein EWS-FLI1 is a network hub that regulates alternative splicing. Proc Natl Acad Sci USA 112, E1307–1316 (2015).

Article CAS Google Scholar

Yin, X., Giap, C., Lazo, J. S. & Prochownik, E. V. Low molecular weight inhibitors of Myc-Max interaction and function. Oncogene 22, 6151–6159 (2003).

Article CAS Google Scholar

Follis, A. V., Hammoudeh, D. I., Wang, H., Prochownik, E. & Metallo, S. J. Structural rationale for the coupled binding and unfolding of the c-Myc oncoprotein by small molecules. Chemistry & Biology 15, 1149–1155 (2008).

Article CAS Google Scholar

Hammoudeh, D. I., Follis, A. V., Prochownik, E. V. & Metallo, S. J. Multiple independent binding sites for small-molecule inhibitors on the oncoprotein c-Myc. J Am Chem Soc 131, 7390–7401 (2009).

Article CAS Google Scholar

Toth, G. et al. Targeting the intrinsically disordered structural ensemble of alpha-synuclein by small molecules as a potential therapeutic strategy for Parkinson's disease. PLoS One 9, e87133 (2014). Pmc3925190.

Article ADS Google Scholar

Convertino, M., Vitalis, A. & Caflisch, A. Disordered binding of small molecules to Abeta (12–28). J Biol Chem 286, 41578–41588 (2011). Pmc3308868.

Article CAS Google Scholar

Jin, F., Yu, C., Lai, L. & Liu, Z. Ligand clouds around protein clouds: a scenario of ligand binding with intrinsically disordered proteins. PLoS Comput Biol 9, e1003249 (2013). Pmc3789766.

Article ADS Google Scholar

Hajduk, P. J. & Greer, J. A decade of fragment-based drug design: strategic advances and lessons learned. Nat Rev Drug Discov 6, 211–219 (2007).

Article CAS Google Scholar

Shuker, S. B., Hajduk, P. J., Meadows, R. P. & Fesik, S. W. Discovering high-affinity ligands for proteins: SAR by NMR. Science 274, 1531–1534 (1996).

Article CAS ADS Google Scholar

Dalvit, C. NMR methods in fragment screening: theory and a comparison with other biophysical techniques. Drug Discov Today 14, 1051–1057 (2009).

Article CAS Google Scholar

Mitrea, D. M. & Kriwacki, R. W. Regulated unfolding of proteins in signaling. FEBS Lett 587, 1081–1088 (2013).

Article CAS Google Scholar

Lacy, E. R. et al. p27 binds cyclin-CDK complexes through a sequential mechanism involving binding-induced protein folding. Nat Struct Mol Biol 11, 358–364 (2004).

Article CAS Google Scholar

Sivakolundu, S. G., Bashford, D. & Kriwacki, R. W. Disordered p27Kip1 exhibits intrinsic structure resembling the Cdk2/cyclin A-bound conformation. J Mol Biol 353, 1118–1128 (2005).

Article CAS Google Scholar

Grimmler, M. et al. Cdk-inhibitory activity and stability of p27Kip1 are directly regulated by oncogenic tyrosine kinases. Cell 128, 269–280 (2007).

Article CAS Google Scholar

Galea, C. A. et al. Role of intrinsic flexibility in signal transduction mediated by the cell cycle regulator, p27 Kip1. J Mol Biol 376, 827–838 (2008). 2350195.

Article CAS Google Scholar

Liang, J. et al. PKB/Akt phosphorylates p27, impairs nuclear import of p27 and opposes p27-mediated G1 arrest. Nat. Med. 8, 1153–1160 (2002).

Article CAS Google Scholar

Shin, I. et al. PKB/Akt mediates cell-cycle progression by phosphorylation of p27(Kip1) at threonine 157 and modulation of its cellular localization. Nat Med 8, 1145–1152 (2002).

Article CAS Google Scholar

Viglietto, G. et al. Cytoplasmic relocalization and inhibition of the cyclin-dependent kinase inhibitor p27(Kip1) by PKB/Akt-mediated phosphorylation in breast cancer. Nat. Med. 8, 1136–1144 (2002).

Article CAS Google Scholar

Zhao, H. et al. Mislocalization of p27 to the cytoplasm of breast cancer cells confers resistance to anti-HER2 targeted therapy. Oncotarget 5, 12704–12714 (2014).

PubMed Google Scholar

Chen, P. & Segil, N. p27(Kip1) links cell proliferation to morphogenesis in the developing organ of Corti. Development 126, 1581–1590 (1999).

CAS PubMed Google Scholar

White, P. M., Doetzlhofer, A., Lee, Y. S., Groves, A. K. & Segil, N. Mammalian cochlear supporting cells can divide and trans-differentiate into hair cells. Nature 441, 984–987 (2006).

Article CAS ADS Google Scholar

Walters, B. J. et al. Auditory Hair Cell-Specific Deletion of p27Kip1 in Postnatal Mice Promotes Cell-Autonomous Generation of New Hair Cells and Normal Hearing. J Neurosci 34, 15751–15763 (2014). Pmc4236404.

Article Google Scholar

Walters, B. J. et al. High-throughput screening reveals alsterpaullone, 2-cyanoethyl as a potent p27Kip1 transcriptional inhibitor. PLoS One 9, e91173 (2014). Pmc3960108.

Article ADS Google Scholar

Sherr, C. J. & Roberts, J. M. CDK inhibitors: positive and negative regulators of G1-phase progression. Genes & development 13, 1501–1512 (1999).

Article CAS Google Scholar

Schulman, B. A., Lindstrom, D. L. & Harlow, E. Substrate recruitment to cyclin-dependent kinase 2 by a multipurpose docking site on cyclin A. Proc Natl Acad Sci USA 95, 10453–10458 (1998). Pmc27915.

Article CAS ADS Google Scholar

Russo, A. A., Jeffrey, P. D., Patten, A. K., Massague, J. & Pavletich, N. P. Crystal structure of the p27Kip1 cyclin-dependent-kinase inhibitor bound to the cyclin A-Cdk2 complex. Nature 382, 325–331 (1996).

Article CAS ADS Google Scholar

Dalvit, C., Fogliatto, G., Stewart, A., Veronesi, M. & Stockman, B. WaterLOGSY as a method for primary NMR screening: practical aspects and range of applicability. Journal of biomolecular NMR 21, 349–359 (2001).

Article CAS Google Scholar

Mayer, M. & Meyer, B. Characterization of Ligand Binding by Saturation Transfer Difference NMR Spectroscopy. Angew. Chem. Int. Ed. 38, 1784–1788 (1999).

3.0.CO;2-Q" data-track-action="article reference" href="" aria-label="Article reference 36" data-doi="10.1002/(SICI)1521-3773(19990614)38:123.0.CO;2-Q">Article CAS Google Scholar

Surade, S. & Blundell, T. L. Structural biology and drug discovery of difficult targets: the limits of ligandability. Chem Biol 19, 42–50 (2012).

Article CAS Google Scholar

Scott, D. E., Coyne, A. G., Hudson, S. A. & Abell, C. Fragment-based approaches in drug discovery and chemical biology. Biochemistry 51, 4990–5003 (2012).

Article CAS Google Scholar

Murray, C. W., Verdonk, M. L. & Rees, D. C. Experiences in fragment-based drug discovery. Trends Pharmacol Sci 33, 224–232 (2012).

Article CAS Google Scholar

Congreve, M., Carr, R., Murray, C. & Jhoti, H. A ‘rule of three’ for fragment-based lead discovery? Drug Discov Today 8, 876–877 (2003).

Article Google Scholar

Cheeseright, T., Mackey, M., Rose, S. & Vinter, A. Molecular field extrema as descriptors of biological activity: definition and validation. Journal of chemical information and modeling 46, 665–676 (2006).

Article CAS Google Scholar

Kriwacki, R. W., Hengst, L., Tennant, L., Reed, S. I. & Wright, P. E. Structural studies of p21(waf1/cip1/sdi1) in the free and Cdk2-bound state: Conformational disorder mediates binding diversity. Proc. Natl. Acad. Sci. USA 93, 11504–11509 (1996).

Article CAS ADS Google Scholar

Otieno, S. & Kriwacki, R. Probing the role of nascent helicity in p27 function as a cell cycle regulator. PLoS One 7, e47177 (2012). Pmc3470550.

Article CAS ADS Google Scholar

Follis, A. V., Hammoudeh, D. I., Wang, H., Prochownik, E. V. & Metallo, S. J. Structural rationale for the coupled binding and unfolding of the c-Myc oncoprotein by small molecules. Chem Biol 15, 1149–1155 (2008).

Article CAS Google Scholar

Otieno, S., Grace, C. R. & Kriwacki, R. W. The role of the LH subdomain in the function of the Cip/Kip cyclin-dependent kinase regulators. Biophys J 100, 2486–2494 (2011). 3093569.

Article CAS Google Scholar

Bowman, P., Galea, C. A., Lacy, E. & Kriwacki, R. W. Thermodynamic characterization of interactions between p27(Kip1) and activated and non-activated Cdk2: intrinsically unstructured proteins as thermodynamic tethers. Biochim Biophys Acta 1764, 182–189 (2006).

Article CAS Google Scholar

Maybridge, T. F. S. I. Maybridge Ro3 Diversity Fragment Library, (2010).

Shelat, A. A. & Guy, R. K. Scaffold composition and biological relevance of screening libraries. Nat Chem Biol 3, 442–446 (2007).

Article CAS Google Scholar

Rishton, G. M. Reactive compounds and in vitro false positives in HTS. Drug Discov Today 2, 382–384 (1997).

Article CAS Google Scholar

Hann, M. et al. Strategic pooling of compounds for high-throughput screening. J Chem Inf Comput Sci 39, 897–902 (1999).

Article CAS Google Scholar

Walters, W. P. & Murcko, M. A. Prediction of ‘drug-likeness’. Adv Drug Deliv Rev 54, 255–271 (2002).

Article CAS Google Scholar

Allu, T. K. & Oprea, T. I. Rapid evaluation of synthetic and molecular complexity for in silico chemistry. Journal of chemical information and modeling 45, 1237–1243 (2005).

Article CAS Google Scholar

Keller, R. The Computer Aided Resonance Assignment Tutorial (CANTINA Verlag, 2004).

Yuwen, T. & Skrynnikov, N. R. CP-HISQC: a better version of HSQC experiment for intrinsically disordered proteins under physiological conditions. J Biomol NMR 58, 175–192 (2014).

Article CAS Google Scholar

Iwahara, J., Jung, Y. S. & Clore, G. M. Heteronuclear NMR spectroscopy for lysine NH(3) groups in proteins: unique effect of water exchange on (15)N transverse relaxation. J Am Chem Soc 129, 2971–2980 (2007).

Article CAS Google Scholar

Delaglio, F. et al. NMR Pipe: A multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277–293 (1995).

Article CAS Google Scholar

Salomon-Ferrer, R., Götz, A. W., Poole, D., Le Grand, S. & Walker, R. C. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. J. Chem. Theory Comput. 9, 3878–3888 (2013).

Article CAS Google Scholar

Ramanathan, A. & Agarwal, P. K. Computational identification of slow conformational fluctuations in proteins. J Phys Chem B 113, 16669–16680 (2009). Pmc2872677.

Article CAS Google Scholar

Download references

The authors thank Christy R. Grace for assistance with NMR experiments, Heather Ross for compound library management and Ariele Viacava Follis for comments on the manuscript. This work was supported by US National Institutes of Health (NIH) grants R01CA082491 and 1R01GM083159 (to R.W.K.), 2R01DC006471 and 1R21DC013879-01 (to J.Z.), Office of Naval Research Grants N000140911014, N000141210191 and N000141210775 (to J.Z.), a US National Cancer Institute Cancer Center Support Grant P30CA21765 (at St. Jude Children's Research Hospital) and ALSAC. L.I.I. was the recipient of the Garwood Foundation Fellowship from St. Jude Children's Research Hospital. D.B. would like to acknowledge support from the US National Institute of General Medical Science (F32GM113290).

Kavitha Bharatham

Present address: Center for Chemical Biology and Therapeutics, Institute for Stem Cell Biology and Regenerative Medicine, GKVK, Bellary Road, Bangalore, 560065, India

Department of Structural Biology,, Memphis, TN, 38105

Luigi I. Iconaru, David Ban, Weixing Zhang & Richard W. Kriwacki

Department of Developmental Neurobiology,, Memphis, TN, 38105

Luigi I. Iconaru & Jian Zuo

Department of Chemical Biology and Therapeutics, St. Jude Children's Research Hospital, Memphis, TN, 38105

Kavitha Bharatham & Anang A. Shelat

Computational Science and Engineering Division, Health Data Sciences Institute, Oak Ridge National Laboratory, Oak Ridge, TN, 37830

Arvind Ramanathan

Department of Microbiology, Immunology and Biochemistry, The University of Tennessee Health Science Center, Memphis, TN, 38163

Richard W. Kriwacki

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

L.I.I., D.B., A.R., W.Z., A.S., J.Z. and R.W.K. designed research; L.I.I., D. B., K.B. and A.R. performed experiments; L.I.I., D.B., K.B., A.R., A.S. and R.W.K. analyzed data; and L.I.I., D.B., A.R., A.S., J.Z. and R.W.K. wrote the manuscript.

The authors declare no competing financial interests.

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit

Reprints and Permissions

Iconaru, L., Ban, D., Bharatham, K. et al. Discovery of Small Molecules that Inhibit the Disordered Protein, p27Kip1. Sci Rep 5, 15686 (2015).

Download citation

Received: 14 July 2015

Accepted: 25 September 2015

Published: 28 October 2015


Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Nature Communications (2022)

Nature Reviews Cancer (2022)

Oncogene (2022)

Nature Reviews Drug Discovery (2022)

Signal Transduction and Targeted Therapy (2022)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.