Computational analysis of phosphopeptide binding to the polo-box domain of the mitotic kinase PLK1 using molecular dynamics simulation.

The Polo-Like Kinase 1 (PLK1) acts as a central regulator of mitosis and is over-expressed in a wide range of human tumours where high levels of expression correlate with a poor prognosis. PLK1 comprises two structural elements, a kinase domain and a polo-box domain (PBD). The PBD binds phosphorylated substrates to control substrate phosphorylation by the kinase domain. Although the PBD preferentially binds to phosphopeptides, it has a relatively broad sequence specificity in comparison with other phosphopeptide binding domains. We analysed the molecular determinants of recognition by performing molecular dynamics simulations of the PBD with one of its natural substrates, CDC25c. Predicted binding free energies were calculated using a molecular mechanics, Poisson-Boltzmann surface area approach. We calculated the per-residue contributions to the binding free energy change, showing that the phosphothreonine residue and the mainchain account for the vast majority of the interaction energy. This explains the very broad sequence specificity with respect to other sidechain residues. Finally, we considered the key role of bridging water molecules at the binding interface. We employed inhomogeneous fluid solvation theory to consider the free energy of water molecules on the protein surface with respect to bulk water molecules. Such an analysis highlights binding hotspots created by elimination of water molecules from hydrophobic surfaces. It also predicts that a number of water molecules are stabilized by the presence of the charged phosphate group, and that this will have a significant effect on the binding affinity. Our findings suggest a molecular rationale for the promiscuous binding of the PBD and highlight a role for bridging water molecules at the interface. We expect that this method of analysis will be very useful for probing other protein surfaces to identify binding hotspots for natural binding partners and small molecule inhibitors.

A Role for Hydration in Interleukin-2 Inducible T Cell Kinase (Itk) Selectivity.

The selectivity profile of kinase inhibitors is commonly explained in terms of favourable or unfavourable interactions between the inhibitor and active site residues. In practice residue sequence differences have not always been sufficient to explain observed selectivity profiles. A series of interleukin-2 inducible T cell kinase (Itk, EC 2.7.10.2) inhibitors that achieve selectivity through the introduction of a single nitrogen atom in an aromatic ring has recently been discovered. Structures of these inhibitors bound to Itk showed this nitrogen to be solvent exposed and not involved in any direct interactions with the enzyme. By analysing active site hydration, using the molecular dynamics tool WaterMap, the observed selectivity profile can be explained in terms of the replacement of a thermodynamically unfavourable water molecule by the inhibitor and improved hydration of the bound ligand. The location of this hydration site was successfully used to enrich virtual screening results in their content of selective Itk inhibitors.

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.

Identification of ligand binding sites on proteins using a multi-scale approach.

Identification of a ligand binding site on a protein is pivotal to drug discovery. To date, no reliable and computationally feasible general approach to this problem has been published. Here we present an automated efficient method for determining binding sites on proteins for potential ligands without any a priori knowledge. Our method is based upon the multiscale concept where we deal with a hierarchy of models generated using a k-means clustering algorithm for the potential ligand. This is done in a simple approach whereby a potential ligand is represented by a growing number of feature points. At each increasing level of detail, a pruning of potential binding site is performed. A nonbonding energy function is used to score the interactions between molecules at each step. The technique was successfully employed to seven protein-ligand complexes. In the current paper we show that the algorithm considerably reduces the computational effort required to solve this problem. This approach offers real opportunities for exploiting the large number of structures that will evolve from structural genomics.