Protein kinase C Related Kinase 1 (PRK1) has been shown to be involved in the regulation of androgen receptor signaling and has been identified as a novel potential drug target for prostate cancer therapy. Since there is no PRK1 crystal structure available to date, multiple PRK1 homology models were generated in order to address the protein flexibility. An in-house library of compounds tested on PRK1 was docked into the ATP binding site of the generated models. In most cases a correct pose of the inhibitors could be identified by ensemble docking, while there is still a challenge of finding a reasonable scoring function that is able to rank compounds according to their biological activity. We estimated the binding free energy for our data set of structurally diverse PRK1 inhibitors using the MM-PB(GB)SA and QM/MM-GBSA methods. The obtained results demonstrate that a correlation between calculated binding free energies and experimental IC50 values was found to be usually higher than using docking scores. Furthermore, the developed approach was tested on a set of diverse PRK1 inhibitors taken from literature, which resulted in a significant correlation. The developed method is computationally inexpensive and can be applied as a postdocking filter in virtual screening as well as for optimization of PRK1 inhibitors in order to prioritize compounds for further biological characterization.

Identification of compounds that can bind to a target protein with high affinity is a nontrivial task in structure-based drug design. Several approaches ranging from simple scoring methods to more computationally demanding methods are usually applied for this purpose. In the current work, we used ligand docking in combination with QM/MM-GBSA, MM-GBSA, and MM-PBSA rescoring to discriminate between active and inactive Myt1 kinase inhibitors. Results show that QM/MM-GBSA rescoring performs better than normal docking scores or MM-GBSA rescoring in classifying active and inactive inhibitors. We also applied QM/MM-GBSA rescoring to estimate the binding affinities of compounds from different virtual screening runs. To prove our approach and to confirm its predictive power, a few compounds which were predicted to be active were purchased and experimentally tested. Among the five selected compounds, three showed significant inhibition of recombinant Myt1. PD-173952, which yielded a favorable QM/MM-GBSA binding free energy, showed a K(i) value of 8.1 nM. In addition, two compounds, PD-180970 and saracatinib, showed inhibition at the low micromolar level. Thus, the developed protocol might be useful for further virtual screening experiments to better discriminate between active and inactive compounds and to further optimize the identified hits.

BackgroundEpigenetics is defined as heritable changes in gene expression that are not based on changes in the DNA sequence. Posttranslational modification of histone proteins is a major mechanism of epigenetic regulation. The kinase PRK1 (protein kinase C related kinase 1, also known as PKN1) phosphorylates histone H3 at threonine 11 and is involved in the regulation of androgen receptor signalling. Thus, it has been identified as a novel drug target but little is known about PRK1 inhibitors and consequences of its inhibition.Methodology/Principal FindingUsing a focused library screening approach, we identified the clinical candidate lestaurtinib (also known as CEP-701) as a new inhibitor of PRK1. Based on a generated 3D model of the PRK1 kinase using the homolog PKC-theta (protein kinase c theta) protein as a template, the key interaction of lestaurtinib with PRK1 was analyzed by means of molecular docking studies. Furthermore, the effects on histone H3 threonine phosphorylation and androgen-dependent gene expression was evaluated in prostate cancer cells.Conclusions/SignificanceLestaurtinib inhibits PRK1 very potently in vitro and in vivo. Applied to cell culture it inhibits histone H3 threonine phosphorylation and androgen-dependent gene expression, a feature that has not been known yet. Thus our findings have implication both for understanding of the clinical activity of lestaurtinib as well as for future PRK1 inhibitors.

High affinity ligands for a given target tend to share key molecular interactions with important anchoring amino acids and therefore often present quite conserved interaction patterns. This simple concept was formalized in a topological knowledge-based scoring function (GRIM) for selecting the most appropriate docking poses from previously X-rayed interaction patterns. GRIM first converts protein-ligand atomic coordinates (docking poses) into a simple 3D graph describing the corresponding interaction pattern. In a second step, proposed graphs are compared to that found from template structures in the Protein Data Bank. Last, all docking poses are rescored according to an empirical score (GRIMscore) accounting for overlap of maximum common subgraphs. Taking the opportunity of the public D3R Grand Challenge 2015, GRIM was used to rescore docking poses for 36 ligands (6 HSP90α inhibitors, 30 MAP4K4 inhibitors) prior to the release of the corresponding protein-ligand X-ray structures. When applied to the HSP90α dataset, for which many protein-ligand X-ray structures are already available, GRIM provided very high quality solutions (mean rmsd = 1.06 Å, n = 6) as top-ranked poses, and significantly outperformed a state-of-the-art scoring function. In the case of MAP4K4 inhibitors, for which preexisting 3D knowledge is scarce and chemical diversity is much larger, the accuracy of GRIM poses decays (mean rmsd = 3.18 Å, n = 30) although GRIM still outperforms an energy-based scoring function. GRIM rescoring appears to be quite robust with comparison to the other approaches competing for the same challenge (42 submissions for the HSP90 dataset, 27 for the MAP4K4 dataset) as it ranked 3rd and 2nd respectively, for the two investigated datasets. The rescoring method is quite simple to implement, independent on a docking engine, and applicable to any target for which at least one holo X-ray structure is available.