Identifying Human Kinase-Specific Protein Phosphorylation Sites by Integrating Heterogeneous Information from Various Sources

Phosphorylation is an important type of protein post-translational modification. Identification of possible phosphorylation sites of a protein is important for understanding its functions. Unbiased screening for phosphorylation sites by in vitro or in vivo experiments is time consuming and expensive; in silico prediction can provide functional candidates and help narrow down the experimental efforts. Most of the existing prediction algorithms take only the polypeptide sequence around the phosphorylation sites into consideration. However, protein phosphorylation is a very complex biological process in vivo. The polypeptide sequences around the potential sites are not sufficient to determine the phosphorylation status of those residues. In the current work, we integrated various data sources such as protein functional domains, protein subcellular location and protein-protein interactions, along with the polypeptide sequences to predict protein phosphorylation sites. The heterogeneous information significantly boosted the prediction accuracy for some kinase families. To demonstrate potential application of our method, we scanned a set of human proteins and predicted putative phosphorylation sites for Cyclin-dependent kinases, Casein kinase 2, Glycogen synthase kinase 3, Mitogen-activated protein kinases, protein kinase A, and protein kinase C families (avaiable at http://cmbi.bjmu.edu.cn/huphospho). The predicted phosphorylation sites can serve as candidates for further experimental validation. Our strategy may also be applicable for the in silico identification of other post-translational modification substrates

Role of structural dynamics at the receptor G protein interface for signal transduction

GPCRs catalyze GDP/GTP exchange in the α-subunit of heterotrimeric G proteins (Gαßγ) through displacement of the Gα C-terminal α5 helix, which directly connects the interface of the active receptor (R*) to the nucleotide binding pocket of G. Hydrogen-deuterium exchange mass spectrometry and kinetic analysis of R* catalysed G protein activation have suggested that displacement of α5 starts from an intermediate GDP bound complex (R*•GGDP). To elucidate the structural basis of receptor-catalysed displacement of α5, we modelled the structure of R*•GGDP. A flexible docking protocol yielded an intermediate R*•GGDP complex, with a similar overall arrangement as in the X-ray structure of the nucleotide free complex (R*•Gempty), however with the α5 C-terminus (GαCT) forming different polar contacts with R*. Starting molecular dynamics simulations of GαCT bound to R* in the intermediate position, we observe a screw-like motion, which restores the specific interactions of α5 with R* in R*•Gempty. The observed rotation of α5 by 60° is in line with experimental data. Reformation of hydrogen bonds, water expulsion and formation of hydrophobic interactions are driving forces of the α5 displacement. We conclude that the identified interactions between R* and G protein define a structural framework in which the α5 displacement promotes direct transmission of the signal from R* to the GDP binding pocket