A New Paradigm for KIM-PTP Drug Discovery: Identification of Allosteric Sites with Potential for Selective Inhibition Using Virtual Screening and LEI Analysis

The kinase interaction motif protein tyrosine phosphatases (KIM-PTPs), HePTP, PTPSL and STEP, are involved in the negative regulation of mitogen-activated protein kinase (MAPK) signalling pathways and are important therapeutic targets for a number of diseases. We have used VSpipe, a virtual screening pipeline, to identify a ligand cluster distribution that is unique to this subfamily of PTPs. Several clusters map onto KIM-PTP specific sequence motifs in contrast to the cluster distribution obtained for PTP1B, a classic PTP that mapped to general PTP motifs. Importantly, the ligand clusters coincide with previously reported functional and substrate binding sites in KIM-PTPs. Assessment of the KIM-PTP specific clusters, using ligand efficiency index (LEI) plots generated by the VSpipe, ascertained that the binders in these clusters reside in a more drug-like chemical–biological space than those at the active site. LEI analysis showed differences between clusters across all KIM-PTPs, highlighting a distinct and specific profile for each phosphatase. The most druggable cluster sites are unexplored allosteric functional sites unique to each target. Exploiting these sites may facilitate the delivery of inhibitors with improved drug-like properties, with selectivity amongst the KIM-PTPs and over other classical PTPs.     

Communication pathways between the nucleotide pocket and distal regulatory sites in protein kinases

Protein kinases control many cellular processes via the ATP-dependent phosphorylation of specific amino acids on target proteins. Despite the availability of the three-dimensional structures of a variety of protein kinases, it has been particularly difficult to explain how noncatalytic domains removed from the active site regulate catalytic function. In this review, we describe how solution methodologies complement the available structural data and help explain how protein kinases may utilize medium-to-long-range effects to regulate substrate phosphorylation. For illustration, two protein kinases, cAMP-dependent protein kinase and the C-terminal Src kinase, are presented as paradigms for the serine/threonine- and tyrosine-specific families. While active-site residues provide an optimal environment for fast phosphoryl group transfer in these and other kinases, the overall rate of protein phosphorylation is limited by nucleotide binding and associated structural changes. Hydrogen-deuterium exchange studies reveal that nucleotide binding induces changes that radiate from a central structural assembly composed of the catalytic loop, glycine-rich loop, and helix alpha C to unique peripheral regions inside and outside the kinase core. This collection of conserved and unique elements delivers information from the active site to distal regions and possibly provides information flow back to the active site. This "push-pull" hypothesis offers a means for understanding how protein kinases can be regulated by protein-protein interactions far from the active site.     

Comparison of conservation within and between the Ser/Thr and Tyr protein kinase family: proposed model for the catalytic domain of the epidermal growth factor receptor

The protein kinase family can be subdivided into two main groupsbased on their ability to phosphorylate Ser/Thr or Tyr substrates.In order to understand the basis of this functional difference,we have carried out a comparative analysis of sequence conservationwithin and between the Ser/Thr and Tyr protein kinases. A multiplesequence alignment of 86 protein kinase sequences was generated.For each position in the alignment we have computed the conservationof residue type in the Ser/Thr, in the Tyr and in both of thekinase subfamilies. To understand the structural and/or functionalbasis for the conservation, we have mapped these conservationproperties onto the backbone of the recently determined structureof the cAMP–dependent Ser/Thr kinase. The results showthat the kinase structure can be roughly segregated, based uponconservation, into three zones. The inner zone contains residueshighly conserved in all the kinase family and describes thehydrophobic core of the enzyme together with residues essentialfor substrate and ATP binding and catalysis. The outer zonecontains residues highly variable in all kinases and representsthe solvent–exposed surface of the protein. The thirdzone is comprised of residues conserved in either the Ser/Thror Tyr kinases or in both, but which are not conserved betweenthem. These are sandwiched between the hydrophobic core andthe solvent-exposed surface. In addition to analyzing overallconservation hi the kinase family, we have also looked at conservationof its substrate and ATP binding sites. The ATP site is highlyconserved throughout the kinases, whereas the substrate bindingsite is more variable. The active site contains several positionswhich differ between the Ser/Thr and Tyr kinases and may beresponsible for discriminating between hydroxyl bearing sidechains. Using this information we propose a model for Tyr substratebinding to the catalytic domain of the epidermal growth factorreceptor (EGFR).     

Comparing Drug Images and Repurposing Drugs with BioGPS and FLAPdock: The Thymidylate Synthase Case

Repurposing and repositioning drugs has become a frequently pursued and successful strategy in the current era, as new chemical entities are increasingly difficult to find and get approved. Herein we report an integrated BioGPS/FLAPdock pipeline for rapid and effective off‐target identification and drug repurposing. Our method is based on the structural and chemical properties of protein binding sites, that is, the ligand image, encoded in the GRID molecular interaction fields (MIFs). Protein similarity is disclosed through the BioGPS algorithm by measuring the pockets’ overlap according to which pockets are clustered. Co‐crystallized and known ligands can be cross‐docked among similar targets, selected for subsequent in vitro binding experiments, and possibly improved for inhibitory potency. We used human thymidylate synthase (TS) as a test case and searched the entire RCSB Protein Data Bank (PDB) for similar target pockets. We chose casein kinase IIα as a control and tested a series of its inhibitors against the TS template. Ellagic acid and apigenin were identified as TS inhibitors, and various flavonoids were selected and synthesized in a second‐round selection. The compounds were demonstrated to be active in the low‐micromolar range.     

Scanning Protein Surfaces with DNA-Encoded Libraries

Understanding the ligandability of a target protein, defined as the capability of a protein to bind drug-like compounds on any site, can give important stimuli to drug-development projects. For instance, inhibition of protein–protein interactions usually depends on the identification of protein surface binders. DNA-encoded chemical libraries (DELs) allow scanning of protein surfaces with large chemical space. Encoded library selection screens uncovered several protein–protein interaction inhibitors and compounds binding to the surface of G protein-coupled receptors (GPCRs) and kinases. The protein surface-binding chemotypes from DELs are predominantly chemically modified and cyclized peptides, and functional small-molecule peptidomimetics. Peptoid libraries and structural peptidomimetics have been less studied in the DEL field, hinting at hitherto less populated chemical space and suggesting alternative library designs. Roughly a third of bioactive molecules evolved from smaller, target-focused libraries. They showcase the potential of encoded libraries to identify more potent molecules from weak, for example, fragment-like, starting points.     

Identification of putative binding sites of P-glycoprotein based on its homology model

A homology model of P-glycoprotein based on the crystal structure of the multidrug transporter Sav1866 is developed, incorporated into a membrane environment, and optimized. The resulting model is analyzed in relation to the functional state and potential binding sites. The comparison of modeled distances to distances reported in experimental studies between particular residues suggests that the model corresponds most closely to the first ATP hydrolysis step of the protein transport cycle. Comparison to the protein 3D structure confirms this suggestion. Using SiteID and Site Finder programs three membrane related binding regions are identified: a region at the interface between the membrane and cytosol and two regions located in the transmembrane domains. The regions contain binding pockets of different size, orientation, and amino acids. A binding pocket located inside the membrane cavity is also identified. The pockets are analyzed in relation to amino acids shown experimentally to influence the protein function. The results suggest that the protein has multiple binding sites and may bind and/or release substrates in multiple pathways.     

NMR backbone assignment of a protein kinase catalytic domain by a combination of several approaches: application to the catalytic subunit of cAMP-dependent protein kinase

Protein phosphorylation is one of the most important mechanisms used for intracellular regulation in eukaryotic cells. Currently, one of the best-characterized protein kinases is the catalytic subunit of cAMP-dependent protein kinase or protein kinase A (PKA). PKA has the typical bilobular structure of kinases, with the active site consisting of a cleft between the two structural lobes. For full kinase activity, the catalytic subunit has to be phosphorylated. The catalytic subunit of PKA has two main phosphorylation sites: Thr197 and Ser338. Binding of ATP or inhibitors to the ATP site induces large structural changes. Here we describe the partial backbone assignment of the PKA catalytic domain by NMR spectroscopy, which represents the first NMR assignment of any protein kinase catalytic domain. Backbone resonance assignment for the 42 kDa protein was accomplished by an approach employing 1) triply ((2)H,(13)C,(15)N) labeled protein and classical NMR assignment experiments, 2) back-calculation of chemical shifts from known X-ray structures, 3) use of paramagnetic adenosine derivatives as spin-labels, and 4) selective amino acid labeling. Interpretation of chemical-shift perturbations allowed mapping of the interaction surface with the protein kinase inhibitor H7. Furthermore, structural conformational changes were observed by comparison of backbone amide shifts obtained by 2D (1)H,(15)N TROSY of an inactive Thr197Ala mutant with the wild-type enzyme.     

Targeting Protein Kinases in Blood Cancer: Focusing on CK1α and CK2

Disturbance of protein kinase activity may result in dramatic consequences that often lead to cancer development and progression. In tumors of blood origin, both tyrosine kinases and serine/threonine kinases are altered by different types of mutations, critically regulating cancer hallmarks. CK1α and CK2 are highly conserved, ubiquitously expressed and constitutively active pleiotropic kinases, which participate in multiple biological processes. The involvement of these kinases in solid and blood cancers is well documented. CK1α and CK2 are overactive in multiple myeloma, leukemias and lymphomas. Intriguingly, they are not required to the same degree for the viability of normal cells, corroborating the idea of “druggable” kinases. Different to other kinases, mutations on the gene encoding CK1α and CK2 are rare or not reported. Actually, these two kinases are outside the paradigm of oncogene addiction, since cancer cells’ dependency on these proteins resembles the phenomenon of “non-oncogene” addiction. In this review, we will summarize the general features of CK1α and CK2 and the most relevant oncogenic and stress-related signaling nodes, regulated by kinase phosphorylation, that may lead to tumor progression. Finally, we will report the current data, which support the positioning of these two kinases in the therapeutic scene of hematological cancers.     

A Fragment‐Based Approach to Probing Adenosine Recognition Sites by Using Dynamic Combinatorial Chemistry

A new strategy that combines the concepts of fragment‐based drug design and dynamic combinatorial chemistry (DCC) for targeting adenosine recognition sites on enzymes is reported. We demonstrate the use of 5′‐deoxy‐5′‐thioadenosine as a noncovalent anchor fragment in dynamic combinatorial libraries templated by Mycobacterium tuberculosis pantothenate synthetase. A benzyl disulfide derivative was identified upon library analysis by HPLC. Structural and binding studies of protein–ligand complexes by X‐ray crystallography and isothermal titration calorimetry informed the subsequent optimisation of the DCC hit into a disulfide containing the novel meta‐nitrobenzyl fragment that targets the pantoate binding site of pantothenate synthetase. Given the prevalence of adenosine‐recognition motifs in enzymes, our results provide a proof‐of‐concept for using this strategy to probe adjacent pockets for a range of adenosine binding enzymes, including other related adenylate‐forming ligases, kinases, and ATPases, as well as NAD(P)(H), CoA and FAD(H₂) binding proteins.     

Experimental and computational active site mapping as a starting point to fragment-based lead discovery

Small highly soluble probe molecules such as aniline, urea, N‐methylurea, 2‐bromoacetate, 1,2‐propanediol, nitrous oxide, benzamidine, and phenol were soaked into crystals of various proteins to map their binding pockets and to detect hot spots of binding with respect to hydrophobic and hydrophilic properties. The selected probe molecules were first tested at the zinc protease thermolysin. They were then applied to a wider range of proteins such as protein kinase A, D ‐xylose isomerase, 4‐diphosphocytidyl‐2C‐methyl‐D ‐erythritol synthase, endothiapepsin, and secreted aspartic protease 2. The crystal structures obtained clearly show that the probe molecules populate the protein binding pockets in an ordered fashion. The thus characterized, experimentally observed hot spots of binding were subjected to computational active site mapping using HotspotsX. This approach uses knowledge‐based pair potentials to detect favorable binding positions for various atom types. Good agreement between the in silico hot spot predictions and the experimentally observed positions of the polar hydrogen bond forming functional groups and hydrophobic portions was obtained. Finally, we compared the observed poses of the small‐molecule probes with those of much larger structurally related ligands. They coincide remarkably well with the larger ligands, considering their spatial orientation and the experienced interaction patterns. This observation confirms the fundamental hypothesis of fragment‐based lead discovery: that binding poses, even of very small molecular probes, do not significantly deviate or move once a ligand is grown further into the binding site. This underscores the fact that these probes populate given hot spots and can be regarded as relevant seeds for further design.     

Understanding Kinase Selectivity Through Energetic Analysis of Binding Site Waters

Kinases remain an important drug target class within the pharmaceutical industry; however, the rational design of kinase inhibitors is plagued by the complexity of gaining selectivity for a small number of proteins within a family of more than 500 related enzymes. Herein we show how a computational method for identifying the location and thermodynamic properties of water molecules within a protein binding site can yield insight into previously inexplicable selectivity and structure–activity relationships. Four kinase systems (Src family, Abl/c‐Kit, Syk/ZAP‐70, and CDK2/4) were investigated, and differences in predicted water molecule locations and energetics were able to explain the experimentally observed binding selectivity profiles. The successful predictions across the range of kinases studied here suggest that this methodology could be generally applicable for predicting selectivity profiles in related targets.     

Mutational analysis supports a structural model for the cell cycle protein kinase p34

Structural models for the eukaryotic cell cycle control proteinp34 from human, S.pombe and S.cerevisiae have been derived fromthe crystallographic coordinates of the cAMP-dependent proteinkinase (cAPK) catalytic subunit (active conformation) and comparedwith the structure of Inactive CDK2 apoenzyme. Differences betweenthe p34 and cAPK catalytic sites provide a possible explanationfor their different substrate specificities. The p34 modelslocalize Tyrl5 and Thrl4 close to the sites of catalysis andsubstrate recognition where their phosphorylatlon could inhibitp34 kinase activity either by blocking MgATP or substrate binding.The conserved sequences PSTAIRE and LYLIFEFL are both closeto the catalytic site and accessible on the protein surfaceavailable to mediate interactions with other proteins. It ispredicted that p34 has an active-site cleft composed almostentirely of sequences common to all protein kinases and sequencesunique to the p34 protein family. Genetic and biochemical analysesof p34 have shown that it interacts extensively with a numberof other proteins. The model allows the relative dispositionof these sites of mutation to each other and to the sites ofcatalysis and substrate recognition to be appreciated. Surfaceregions on p34 that are important for function have been identified.These sites identify residues that may interact with p13^SUCL,cydin, plO7^WEEL and p80^cdc25     

Phosphocatalytic Kinome Activity Profiling of Apoptotic and Ferroptotic Agents in Multiple Myeloma Cells

Through phosphorylation of their substrate proteins, protein kinases are crucial for transducing cellular signals and orchestrating biological processes, including cell death and survival. Recent studies have revealed that kinases are involved in ferroptosis, an iron-dependent mode of cell death associated with toxic lipid peroxidation. Given that ferroptosis is being explored as an alternative strategy to eliminate apoptosis-resistant tumor cells, further characterization of ferroptosis-dependent kinase changes might aid in identifying novel druggable targets for protein kinase inhibitors in the context of cancer treatment. To this end, we performed a phosphopeptidome based kinase activity profiling of glucocorticoid-resistant multiple myeloma cells treated with either the apoptosis inducer staurosporine (STS) or ferroptosis inducer RSL3 and compared their kinome activity signatures. Our data demonstrate that both cell death mechanisms inhibit the activity of kinases classified into the CMGC and AGC families, with STS showing a broader spectrum of serine/threonine kinase inhibition. In contrast, RSL3 targets a significant number of tyrosine kinases, including key players of the B-cell receptor signaling pathway. Remarkably, additional kinase profiling of the anti-cancer agent withaferin A revealed considerable overlap with ferroptosis and apoptosis kinome activity, explaining why withaferin A can induce mixed ferroptotic and apoptotic cell death features. Altogether, we show that apoptotic and ferroptotic cell death induce different kinase signaling changes and that kinome profiling might become a valid approach to identify cell death chemosensitization modalities of novel anti-cancer agents.     

Architectural Repertoire of Ligand‐Binding Pockets on Protein Surfaces

Knowledge of the three‐dimensional structure of ligand binding sites in proteins provides valuable information for computer‐assisted drug design. We present a method for the automated extraction and classification of ligand binding site topologies, in which protein surface cavities are represented as branched frameworks. The procedure employs a growing neural gas approach for pocket topology assignment and pocket framework generation. We assessed the structural diversity of 623 known ligand binding site topologies based on framework cluster analysis. At a resolution of 5 Å only 23 structurally distinct topology groups were formed; this suggests an overall limited structural diversity of ligand‐accommodating protein cavities. Higher resolution allowed for identification of protein‐family specific pocket features. Pocket frameworks highlight potentially preferred modes of ligand–receptor interactions and will help facilitate the identification of druggable subpockets suitable for ligand affinity and selectivity optimization.     

Mechanism and substrate specificity of tRNA-guanine transglycosylases (TGTs): tRNA-modifying enzymes from the three different kingdoms of life share a common catalytic mechanism

Transfer RNA-guanine transglycosylases (TGTs) are evolutionarily ancient enzymes, present in all kingdoms of life, catalyzing guanine exchange within their cognate tRNAs by modified 7-deazaguanine bases. Although distinct bases are incorporated into tRNA at different positions in a kingdom-specific manner, the catalytic subunits of TGTs are structurally well conserved. This review provides insight into the sequential steps along the reaction pathway, substrate specificity, and conformational adaptions of the binding pockets by comparison of TGT crystal structures in complex with RNA substrates of a eubacterial and an archaebacterial species. Substrate-binding modes indicate an evolutionarily conserved base-exchange mechanism with a conserved aspartate serving as a nucleophile through covalent binding to C1' of the guanosine ribose moiety in an intermediate state. A second conserved aspartate seems to control the spatial rearrangement of the ribose ring along the reaction pathway and supposedly operates as a general acid/base. Water molecules inside the binding pocket accommodating interaction sites subsequently occupied by polar atoms of substrates help to elucidate substrate-recognition and substrate-specificity features. This emphasizes the role of water molecules as general probes to map binding-site properties for structure-based drug design. Additionally, substrate-bound crystal structures allow the extraction of valuable information about the classification of the TGT superfamily into a subdivision of presumably homologous superfamilies adopting the triose-phosphate isomerase type barrel fold with a standard phosphate-binding motif.     

Probing novel 1-aza-9-oxafluorenes as selective GSK-3beta inhibitors

Within the histopathology of Alzheimer's disease (AD) certain hallmarks are beeing observed. The occurance of protein deposits belong to such characteristic features. Such deposits can be found extracellular as beta-amyloid (Abeta) plaques and intracellular as neurofibrillary tangles (NFTs). In the search for novel AD therapeutics it became of great interest to investigate the formation of NFTs and their contribution to the AD symptomatic. NFTs consist of hyperphosphorylated tau protein. Within the phosphorylation process of tau protein two kinases are of great importance: cyclin dependent kinase 5 (cdk5) and its truncated regulatory subunit p25 and glycogen synthase kinase 3beta (GSK-3beta). The role of both kinases within the NFT formation process is still under debate. To better understand the pathophysiological process highly selective inhibitors of both kinases are of value. Known inhibitors lack the necessary selectivity. We developed novel 1-aza-9-oxafluo-renes as selective GSK-3beta inhibitors. Structure-activity relationships of a series of 4-phenyl substituted derivatives are discussed. Variation of the 3-side chain led to selective carbonyl amide derivatives with selectivity factors of more than 100 at the tested ATP competitor concentrations. Such selectivities permit specific investigation of the role of GSK-3beta within the NFT formation processes.     

Dynamic substrate enhancement for the identification of specific, second-site-binding fragments targeting a set of protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) are key regulators in living systems and thus are attractive drug targets. The development of potent, selective PTP inhibitors has been a difficult challenge mainly due to the high homology of the phosphotyrosine substrate pockets. Here, a strategy of dynamic substrate enhancement is described targeting the secondary binding sites of PTPs. By screening four different PTPs from bacterial (MptpA) and human origin (PTP1B, HePtp, Shp2) with this assay, specific fragments were identified. One highly specific fragment that binds to the secondary site of Mycobacterium tuberculosis protein tyrosine phosphatase A (MptpA) was characterized in order to validate the assay concept. Finally by covalently linking the secondary fragment to a phosphotyrosine mimetic, a moderately active but highly specific inhibitor of MptpA was obtained.     

Probes for Photoaffinity Labelling of Kinases

Protein kinases, one of the largest enzyme superfamilies, regulate many physiological and pathological processes. They are drug targets for multiple human diseases, including various cancer types. Probes for the photoaffinity labelling of kinases are important research tools for the study of members of this enzyme superfamily. In this review, we discuss the design principles of these probes, which are mainly derived from inhibitors targeting the ATP pocket. Overall, insights from crystal structures guide the placement of photoreactive groups and detection tags. This has resulted in a wide variety of probes, of which we provide a comprehensive overview. We also discuss several areas of application of these probes, including the identification of targets and off-targets of kinase inhibitors, mapping of their binding sites, the development of inhibitor screening assays, the imaging of kinases, and identification of protein binding partners.