Structural evolution of the protein kinase-like superfamily

The protein kinase family is large and important, but it is only one family in a larger superfamily of homologous kinases that phosphorylate a variety of substrates and play important roles in all three superkingdoms of life. We used a carefully constructed structural alignment of selected kinases as the basis for a study of the structural evolution of the protein kinase-like superfamily. The comparison of structures revealed a "universal core" domain consisting only of regions required for ATP binding and the phosphotransfer reaction. Remarkably, even within the universal core some kinase structures display notable changes, while still retaining essential activity. Hence, the protein kinase-like superfamily has undergone substantial structural and sequence revision over long evolutionary timescales. We constructed a phylogenetic tree for the superfamily using a novel approach that allowed for the combination of sequence and structure information into a unified quantitative analysis. When considered against the backdrop of species distribution and other metrics, our tree provides a compelling scenario for the development of the various kinase families from a shared common ancestor. We propose that most of the so-called "atypical kinases" are not intermittently derived from protein kinases, but rather diverged early in evolution to form a distinct phyletic group. Within the atypical kinases, the aminoglycoside and choline kinase families appear to share the closest relationship. These two families in turn appear to be the most closely related to the protein kinase family. In addition, our analysis suggests that the actin-fragmin kinase, an atypical protein kinase, is more closely related to the phosphoinositide-3 kinase family than to the protein kinase family. The two most divergent families, alpha-kinases and phosphatidylinositol phosphate kinases (PIPKs), appear to have distinct evolutionary histories. While the PIPKs probably have an evolutionary relationship with the rest of the kinase superfamily, the relationship appears to be very distant (and perhaps indirect). Conversely, the alpha-kinases appear to be an exception to the scenario of early divergence for the atypical kinases: they apparently arose relatively recently in eukaryotes. We present possible scenarios for the derivation of the alpha-kinases from an extant kinase fold.     

Microbial Protein-tyrosine Kinases

Microbial ester kinases identified in the past 3 decades came as a surprise, as protein phosphorylation on Ser, Thr, and Tyr amino acids was thought to be unique to eukaryotes. Current analysis of available microbial genomes reveals that “eukaryote-like” protein kinases are prevalent in prokaryotes and can converge in the same signaling pathway with the classical microbial “two-component” systems. Most microbial tyrosine kinases lack the “eukaryotic” Hanks domain signature and are designated tyrosine kinases based upon their biochemical activity. These include the tyrosine kinases termed bacterial tyrosine kinases (BY-kinases), which are responsible for the majority of known bacterial tyrosine phosphorylation events. Although termed generally as bacterial tyrosine kinases, BY-kinases can be considered as one family belonging to the superfamily of prokaryotic protein-tyrosine kinases in bacteria. Other members of this superfamily include atypical “odd” tyrosine kinases with diverse mechanisms of protein phosphorylation and the “eukaryote-like” Hanks-type tyrosine kinases. Here, we discuss the distribution, phylogeny, and function of the various prokaryotic protein-tyrosine kinases, focusing on the recently discovered Mycobacterium tuberculosis PtkA and its relationship with other members of this diverse family of proteins.     

The Structure of an NDR/LATS Kinase–Mob Complex Reveals a Novel Kinase–Coactivator System and Substrate Docking Mechanism

Eukaryotic cells commonly use protein kinases in signaling systems that relay information and control a wide range of processes. These enzymes have a fundamentally similar structure, but achieve functional diversity through variable regions that determine how the catalytic core is activated and recruited to phosphorylation targets. “Hippo” pathways are ancient protein kinase signaling systems that control cell proliferation and morphogenesis; the NDR/LATS family protein kinases, which associate with “Mob” coactivator proteins, are central but incompletely understood components of these pathways. Here we describe the crystal structure of budding yeast Cbk1–Mob2, to our knowledge the first of an NDR/LATS kinase–Mob complex. It shows a novel coactivator-organized activation region that may be unique to NDR/LATS kinases, in which a key regulatory motif apparently shifts from an inactive binding mode to an active one upon phosphorylation. We also provide a structural basis for a substrate docking mechanism previously unknown in AGC family kinases, and show that docking interaction provides robustness to Cbk1’s regulation of its two known in vivo substrates. Co-evolution of docking motifs and phosphorylation consensus sites strongly indicates that a protein is an in vivo regulatory target of this hippo pathway, and predicts a new group of high-confidence Cbk1 substrates that function at sites of cytokinesis and cell growth. Moreover, docking peptides arise in unstructured regions of proteins that are probably already kinase substrates, suggesting a broad sequential model for adaptive acquisition of kinase docking in rapidly evolving intrinsically disordered polypeptides.     

Phylogenomics of phosphoinositide lipid kinases: perspectives on the evolution of second messenger signaling and drug discovery

Background  

Phosphoinositide lipid kinases (PIKs) generate specific phosphorylated variants of phosatidylinositols (PtdIns) that are critical for second messenger signaling and cellular membrane remodeling. Mammals have 19 PIK isoforms spread across three major families: the PtIns 3-kinases (PI3Ks), PtdIns 4-kinases (PI4Ks), and PtdIns-P (PIP) kinases (PIPKs). Other eukaryotes have fewer yet varying PIK complements. PIKs are also an important, emerging class of drug targets for many therapeutic areas including cancer, inflammatory and metabolic diseases and host-pathogen interactions. Here, we report the genomic occurrences and evolutionary relationships or phylogenomics of all three PIK families across major eukaryotic groups and suggest potential ramifications for drug discovery.     

Protein proximity networks and functional evaluation of the casein kinase 1 gamma family reveal unique roles for CK1γ3 in WNT signaling

Aberrant activation or suppression of WNT/β-catenin signaling contributes to cancer initiation and progression, neurodegeneration, and bone disease. However, despite great need and more than 40 years of research, targeted therapies for the WNT pathway have yet to be fully realized. Kinases are considered exceptionally druggable and occupy key nodes within the WNT signaling network, but several pathway-relevant kinases remain understudied and “dark.” Here, we studied the function of the casein kinase 1γ (CSNK1γ) subfamily of human kinases and their roles in WNT signaling. miniTurbo-based proximity biotinylation and mass spectrometry analysis of CSNK1γ1, CSNK1γ2, and CSNK1γ3 revealed numerous components of the β-catenin–dependent and β-catenin–independent WNT pathways. In gain-of-function experiments, we found that CSNK1γ3 but not CSNK1γ1 or CSNK1γ2 activated β-catenin–dependent WNT signaling, with minimal effect on other signaling pathways. We also show that within the family, CSNK1γ3 expression uniquely induced low-density lipoprotein receptor–related protein 6 phosphorylation, which mediates downstream WNT signaling transduction. Conversely, siRNA-mediated silencing of CSNK1γ3 alone had no impact on WNT signaling, though cosilencing of all three family members decreased WNT pathway activity. Finally, we characterized two moderately selective and potent small-molecule inhibitors of the CSNK1γ family. We show that these inhibitors and a CSNK1γ3 kinase–dead mutant suppressed but did not eliminate WNT-driven low-density lipoprotein receptor–related protein 6 phosphorylation and β-catenin stabilization. Our data suggest that while CSNK1γ3 expression uniquely drives pathway activity, potential functional redundancy within the family necessitates loss of all three family members to suppress the WNT signaling pathway.     

Oncogenic Mutations of p110α Isoform of PI 3-Kinase Upregulate Its Protein Kinase Activity

In addition to lipid kinase activity, the class-I PI 3-kinases also function as protein kinases targeting regulatory autophosphorylation sites and exogenous substrates. The latter include a recently identified regulatory phosphorylation of the GM-CSF/IL-3 βc receptor contributing to survival of acute myeloid leukaemia cells. Previous studies suggested differences in the protein kinase activity of the 4 isoforms of class-I PI 3-kinase so we compared the ability of all class-I PI 3-kinases and 2 common oncogenic mutants to autophosphorylate, and to phosphorylate an intracellular fragment of the GM-CSF/IL-3 βc receptor (βic). We find p110α, p110β and p110γ all phosphorylate βic but p110δ is much less effective. The two most common oncogenic mutants of p110α, H1047R and E545K have stronger protein kinase activity than wildtype p110α, both in terms of autophosphorylation and towards βic. Importantly, the lipid kinase activity of the oncogenic mutants is still inhibited by autophosphorylation to a similar extent as wildtype p110α. Previous evidence indicates the protein kinase activity of p110α is Mn²⁺ dependent, casting doubt over its role in vivo. However, we show that the oncogenic mutants of p110α plus p110β and p110γ all display significant activity in the presence of Mg²⁺. Furthermore we demonstrate that some small molecule inhibitors of p110α lipid kinase activity (PIK-75 and A66) are equally effective against the protein kinase activity, but other inhibitors (e.g. wortmannin and TGX221) show different patterns of inhibition against the lipid and protein kinases activities. These findings have implications for the function of PI 3-kinase, especially in tumours carrying p110α mutations.     

A human kinase yeast array for the identification of kinases modulating phosphorylation‐dependent protein–protein interactions

Protein kinases play an important role in cellular signaling pathways and their dysregulation leads to multiple diseases, making kinases prime drug targets. While more than 500 human protein kinases are known to collectively mediate phosphorylation of over 290,000 S/T/Y sites, the activities have been characterized only for a minor, intensively studied subset. To systematically address this discrepancy, we developed a human kinase array in Saccharomyces cerevisiae as a simple readout tool to systematically assess kinase activities. For this array, we expressed 266 human kinases in four different S. cerevisiae strains and profiled ectopic growth as a proxy for kinase activity across 33 conditions. More than half of the kinases showed an activity‐dependent phenotype across many conditions and in more than one strain. We then employed the kinase array to identify the kinase(s) that can modulate protein–protein interactions (PPIs). Two characterized, phosphorylation‐dependent PPIs with unknown kinase–substrate relationships were analyzed in a phospho‐yeast two‐hybrid assay. CK2α1 and SGK2 kinases can abrogate the interaction between the spliceosomal proteins AAR2 and PRPF8, and NEK6 kinase was found to mediate the estrogen receptor (ERα) interaction with 14‐3‐3 proteins. The human kinase yeast array can thus be used for a variety of kinase activity‐dependent readouts.     

A Conserved Loop in the Catalytic Domain of Eukaryotic Elongation Factor 2 Kinase Plays a Key Role in Its Substrate Specificity

Eukaryotic elongation factor 2 kinase (eEF2K) is the best-characterized member of the α-kinase family. Within this group, only eEF2K and myosin heavy chain kinases (MHCKs) have known substrates. Here we have studied the roles of specific residues, selected on the basis of structural data for MHCK A and TRPM7, in the function of eEF2K. Our data provide the first information regarding the basis of the substrate specificity of α-kinases, in particular the roles of residues in the so-called N/D loop, which appears to occupy a position in the structure of α-kinases similar to that of the activation loop in other kinases. Several mutations in the EEF2K gene occur in tumors, one of which (Arg303Cys) is at a highly conserved residue in the N/D loop. This mutation greatly enhances eEF2K activity and may be cytoprotective. Our data support the concept that the major autophosphorylation site (Thr348 in eEF2K) docks into a binding pocket to help create the kinase-competent conformation. This is similar to the situation for MHCK A and is consistent with this being a common feature of α-kinases.     

STE20 phosphorylation of AMPK-related kinases revealed by biochemical purifications combined with genetics

We have recently purified mammalian sterile 20 (STE20)–like kinase 3 (MST3) as a kinase for the multifunctional kinases, AMP-activated protein kinase–related kinases (ARKs). However, unresolved questions from this study, such as remaining phosphorylation activities following deletion of the Mst3 gene from human embryonic kidney cells and mice, led us to conclude that there were additional kinases for ARKs. Further purification recovered Ca²⁺/calmodulin-dependent protein kinase kinases 1 and 2 (CaMKK1 and 2), and a third round of purification revealed mitogen-activated protein kinase kinase kinase kinase 5 (MAP4K5) as potential kinases of ARKs. We then demonstrated that MST3 and MAP4K5, both belonging to the STE20-like kinase family, could phosphorylate all 14 ARKs both in vivo and in vitro. Further examination of all 28 STE20 kinases detected variable phosphorylation activity on AMP-activated protein kinase (AMPK) and the salt-inducible kinase 3 (SIK3). Taken together, our results have revealed novel relationships between STE20 kinases and ARKs, with potential physiological and pathological implications.     

Diversity, classification and function of the plant protein kinase superfamily

Eukaryotic protein kinases belong to a large superfamily with hundreds to thousands of copies and are components of essentially all cellular functions. The goals of this study are to classify protein kinases from 25 plant species and to assess their evolutionary history in conjunction with consideration of their molecular functions. The protein kinase superfamily has expanded in the flowering plant lineage, in part through recent duplications. As a result, the flowering plant protein kinase repertoire, or kinome, is in general significantly larger than other eukaryotes, ranging in size from 600 to 2500 members. This large variation in kinome size is mainly due to the expansion and contraction of a few families, particularly the receptor-like kinase/Pelle family. A number of protein kinases reside in highly conserved, low copy number families and often play broadly conserved regulatory roles in metabolism and cell division, although functions of plant homologues have often diverged from their metazoan counterparts. Members of expanded plant kinase families often have roles in plant-specific processes and some may have contributed to adaptive evolution. Nonetheless, non-adaptive explanations, such as kinase duplicate subfunctionalization and insufficient time for pseudogenization, may also contribute to the large number of seemingly functional protein kinases in plants.     

UPF201 Archaeal Specific Family Members Reveal Structural Similarity to RNA-Binding Proteins but Low Likelihood for RNA-Binding Function

We have determined X-ray crystal structures of four members of an archaeal specific family of proteins of unknown function (UPF0201; Pfam classification: DUF54) to advance our understanding of the genetic repertoire of archaea. Despite low pairwise amino acid sequence identities (10–40%) and the absence of conserved sequence motifs, the three-dimensional structures of these proteins are remarkably similar to one another. Their common polypeptide chain fold, encompassing a five-stranded antiparallel β-sheet and five α-helices, proved to be quite unexpectedly similar to that of the RRM-type RNA-binding domain of the ribosomal L5 protein, which is responsible for binding the 5S- rRNA. Structure-based sequence alignments enabled construction of a phylogenetic tree relating UPF0201 family members to L5 ribosomal proteins and other structurally similar RNA binding proteins, thereby expanding our understanding of the evolutionary purview of the RRM superfamily. Analyses of the surfaces of these newly determined UPF0201 structures suggest that they probably do not function as RNA binding proteins, and that this domain specific family of proteins has acquired a novel function in archaebacteria, which awaits experimental elucidation.     

Structure-Based Network Analysis of Activation Mechanisms in the ErbB Family of Receptor Tyrosine Kinases: The Regulatory Spine Residues Are Global Mediators of Structural Stability and Allosteric Interactions

The ErbB protein tyrosine kinases are among the most important cell signaling families and mutation-induced modulation of their activity is associated with diverse functions in biological networks and human disease. We have combined molecular dynamics simulations of the ErbB kinases with the protein structure network modeling to characterize the reorganization of the residue interaction networks during conformational equilibrium changes in the normal and oncogenic forms. Structural stability and network analyses have identified local communities integrated around high centrality sites that correspond to the regulatory spine residues. This analysis has provided a quantitative insight to the mechanism of mutation-induced “superacceptor” activity in oncogenic EGFR dimers. We have found that kinase activation may be determined by allosteric interactions between modules of structurally stable residues that synchronize the dynamics in the nucleotide binding site and the αC-helix with the collective motions of the integrating αF-helix and the substrate binding site. The results of this study have pointed to a central role of the conserved His-Arg-Asp (HRD) motif in the catalytic loop and the Asp-Phe-Gly (DFG) motif as key mediators of structural stability and allosteric communications in the ErbB kinases. We have determined that residues that are indispensable for kinase regulation and catalysis often corresponded to the high centrality nodes within the protein structure network and could be distinguished by their unique network signatures. The optimal communication pathways are also controlled by these nodes and may ensure efficient allosteric signaling in the functional kinase state. Structure-based network analysis has quantified subtle effects of ATP binding on conformational dynamics and stability of the EGFR structures. Consistent with the NMR studies, we have found that nucleotide-induced modulation of the residue interaction networks is not limited to the ATP site, and may enhance allosteric cooperativity with the substrate binding region by increasing communication capabilities of mediating residues.     

赤散囊菌MAPKs超家族的鉴定及分析

本研究以赤散囊菌Eurotium rubrum全基因组序列为对象,利用HMMER软件构建隐马尔可夫模型（hidden markov models,HMM）结合BLAST的方法鉴定了促分裂原活化蛋白激酶（mitogen-activated protein kinase,MAPK）超家族。通过构建系统发育树对鉴定蛋白进行分析,并利用MEME软件进行了保守性基序的预测及活性位点注释。分析结果表明,赤散囊菌基因组包含了4个MAPK蛋白,分别属于Hog1-type、MpkC-type、Slt2-type和Fus3/Kss1-type类型;3个MAPK kinase（MAPKK）蛋白,分别属于MKK1-type、Pbs2-type和Ste7-type类型;3个MAPK kinase kinase（MAPKKK）蛋白,分别属于BCK1-type、Ste11-type和Ssk22-type类型。保守性基序分析及注释结果表明,MAPKs超家族蛋白都包含了蛋白激酶活性位点“-D[L/I/V]K-”以及保守性的ATP-binding标签序列。MAPK与MAPKK蛋白分别包含了“-TxY-”和“-SD[I/V]WS-”磷酸化位点,且MAPK蛋白还包含一个保守性的common docking基序（CD motif）,而MAPKKK蛋白则包含了一个功能不明的保守性基序,其一致性序列为“-GTPYWMAPEV-”。研究结果为揭示MAPKs信号途径在赤散囊菌中参与调控的生物学过程奠定了基础。     

IRBIT Interacts with the Catalytic Core of Phosphatidylinositol Phosphate Kinase Type Iα and IIα through Conserved Catalytic Aspartate Residues

Phosphatidylinositol phosphate kinases (PIPKs) are lipid kinases that generate phosphatidylinositol 4,5-bisphosphate (PI(4,5)P₂), a critical lipid signaling molecule that regulates diverse cellular functions, including the activities of membrane channels and transporters. IRBIT (IP₃ R-binding protein released with inositol 1,4,5-trisphosphate) is a multifunctional protein that regulates diverse target proteins. Here, we report that IRBIT forms signaling complexes with members of the PIPK family. IRBIT bound to all PIPK isoforms in heterologous expression systems and specifically interacted with PIPK type Iα (PIPKIα) and type IIα (PIPKIIα) in mouse cerebellum. Site-directed mutagenesis revealed that two conserved catalytic aspartate residues of PIPKIα and PIPKIIα are involved in the interaction with IRBIT. Furthermore, phosphatidylinositol 4-phosphate, Mg²⁺, and/or ATP interfered with the interaction, suggesting that IRBIT interacts with catalytic cores of PIPKs. Mutations of phosphorylation sites in the serine-rich region of IRBIT affected the selectivity of its interaction with PIPKIα and PIPKIIα. The structural flexibility of the serine-rich region, located in the intrinsically disordered protein region, is assumed to underlie the mechanism of this interaction. Furthermore, in vitro binding experiments and immunocytochemistry suggest that IRBIT and PIPKIα interact with the Na⁺/HCO₃ ⁻ cotransporter NBCe1-B. These results suggest that IRBIT forms signaling complexes with PIPKIα and NBCe1-B, whose activity is regulated by PI(4,5)P₂.     

Co-Conserved MAPK Features Couple D-Domain Docking Groove to Distal Allosteric Sites via the C-Terminal Flanking Tail

Mitogen activated protein kinases (MAPKs) form a closely related family of kinases that control critical pathways associated with cell growth and survival. Although MAPKs have been extensively characterized at the biochemical, cellular, and structural level, an integrated evolutionary understanding of how MAPKs differ from other closely related protein kinases is currently lacking. Here, we perform statistical sequence comparisons of MAPKs and related protein kinases to identify sequence and structural features associated with MAPK functional divergence. We show, for the first time, that virtually all MAPK-distinguishing sequence features, including an unappreciated short insert segment in the β4-β5 loop, physically couple distal functional sites in the kinase domain to the D-domain peptide docking groove via the C-terminal flanking tail (C-tail). The coupling mediated by MAPK-specific residues confers an allosteric regulatory mechanism unique to MAPKs. In particular, the regulatory αC-helix conformation is controlled by a MAPK-conserved salt bridge interaction between an arginine in the αC-helix and an acidic residue in the C-tail. The salt-bridge interaction is modulated in unique ways in individual sub-families to achieve regulatory specificity. Our study is consistent with a model in which the C-tail co-evolved with the D-domain docking site to allosterically control MAPK activity. Our study provides testable mechanistic hypotheses for biochemical characterization of MAPK-conserved residues and new avenues for the design of allosteric MAPK inhibitors.     

ATP and MO25α Regulate the Conformational State of the STRADα Pseudokinase and Activation of the LKB1 Tumour Suppressor

Pseudokinases lack essential residues for kinase activity, yet are emerging as important regulators of signal transduction networks. The pseudokinase STRAD activates the LKB1 tumour suppressor by forming a heterotrimeric complex with LKB1 and the scaffolding protein MO25. Here, we describe the structure of STRADα in complex with MO25α. The structure reveals an intricate web of interactions between STRADα and MO25α involving the αC-helix of STRADα, reminiscent of the mechanism by which CDK2 interacts with cyclin A. Surprisingly, STRADα binds ATP and displays a closed conformation and an ordered activation loop, typical of active protein kinases. Inactivity is accounted for by nonconservative substitution of almost all essential catalytic residues. We demonstrate that binding of ATP enhances the affinity of STRADα for MO25α, and conversely, binding of MO25α promotes interaction of STRADα with ATP. Mutagenesis studies reveal that association of STRADα with either ATP or MO25α is essential for LKB1 activation. We conclude that ATP and MO25α cooperate to maintain STRADα in an “active” closed conformation required for LKB1 activation. It has recently been demonstrated that a mutation in human STRADα that truncates a C-terminal region of the pseudokinase domain leads to the polyhydramnios, megalencephaly, symptomatic epilepsy (PMSE) syndrome. We demonstrate this mutation destabilizes STRADα and prevents association with LKB1. In summary, our findings describe one of the first structures of a genuinely inactive pseudokinase. The ability of STRADα to activate LKB1 is dependent on a closed “active” conformation, aided by ATP and MO25α binding. Thus, the function of STRADα is mediated through an active kinase conformation rather than kinase activity. It is possible that other pseudokinases exert their function through nucleotide binding and active conformations.     

ATP and MO25α Regulate the Conformational State of the STRADα Pseudokinase and Activation of the LKB1 Tumour Suppressor

Pseudokinases lack essential residues for kinase activity, yet are emerging as important regulators of signal transduction networks. The pseudokinase STRAD activates the LKB1 tumour suppressor by forming a heterotrimeric complex with LKB1 and the scaffolding protein MO25. Here, we describe the structure of STRADα in complex with MO25α. The structure reveals an intricate web of interactions between STRADα and MO25α involving the αC-helix of STRADα, reminiscent of the mechanism by which CDK2 interacts with cyclin A. Surprisingly, STRADα binds ATP and displays a closed conformation and an ordered activation loop, typical of active protein kinases. Inactivity is accounted for by nonconservative substitution of almost all essential catalytic residues. We demonstrate that binding of ATP enhances the affinity of STRADα for MO25α, and conversely, binding of MO25α promotes interaction of STRADα with ATP. Mutagenesis studies reveal that association of STRADα with either ATP or MO25α is essential for LKB1 activation. We conclude that ATP and MO25α cooperate to maintain STRADα in an “active” closed conformation required for LKB1 activation. It has recently been demonstrated that a mutation in human STRADα that truncates a C-terminal region of the pseudokinase domain leads to the polyhydramnios, megalencephaly, symptomatic epilepsy (PMSE) syndrome. We demonstrate this mutation destabilizes STRADα and prevents association with LKB1. In summary, our findings describe one of the first structures of a genuinely inactive pseudokinase. The ability of STRADα to activate LKB1 is dependent on a closed “active” conformation, aided by ATP and MO25α binding. Thus, the function of STRADα is mediated through an active kinase conformation rather than kinase activity. It is possible that other pseudokinases exert their function through nucleotide binding and active conformations.     

Conservation of structural fluctuations in homologous protein kinases and its implications on functional sites

Our aim is to explore the similarities in structural fluctuations of homologous kinases. Gaussian Network Model based Normal Mode Analysis was performed on 73 active conformation structures in Ser/Thr/Tyr kinase superfamily. Categories of kinases with progressive evolutionary divergence, viz. (i) Same kinase with many crystal structures, (ii) Within‐Subfamily, (iii) Within‐Family, (iv) Within‐Group, and (v) Across‐Group, were analyzed. We identified a flexibility signature conserved in all kinases involving residues in and around the catalytic loop with consistent low‐magnitude fluctuations. However, the overall structural fluctuation profiles are conserved better in closely related kinases (Within‐Subfamily and Within‐family) than in distant ones (Within‐Group and Across‐Group). A substantial 65.4% of variation in flexibility was not accounted by variation in sequences or structures. Interestingly, we identified substructural residue‐wise fluctuation patterns characteristic of kinases of different categories. Specifically, we recognized statistically significant fluctuations unique to families of protein kinase A, cyclin‐dependent kinases, and nonreceptor tyrosine kinases. These fluctuation signatures localized to sites known to participate in protein‐protein interactions typical of these kinase families. We report for the first time that residues characterized by fluctuations unique to the group/family are involved in interactions specific to the group/family. As highlighted for Src family, local regions with differential fluctuations are proposed as attractive targets for drug design. Overall, our study underscores the importance of consideration of fluctuations, over and above sequence and structural features, in understanding the roles of sites characteristic of kinases. Proteins 2016; 84:957–978. © 2016 Wiley Periodicals, Inc.     

Evolutionary Paths of the cAMP-Dependent Protein Kinase (PKA) Catalytic Subunits

3′,5′-cyclic adenosine monophosphate (cAMP) dependent protein kinase or protein kinase A (PKA) has served as a prototype for the large family of protein kinases that are crucially important for signal transduction in eukaryotic cells. The PKA catalytic subunits Cα and Cβ, encoded by the two genes PRKACA and PRKACB, respectively, are among the best understood and characterized human kinases. Here we have studied the evolution of this gene family in chordates, arthropods, mollusks and other animals employing probabilistic methods and show that Cα and Cβ arose by duplication of an ancestral PKA catalytic subunit in a common ancestor of vertebrates. The two genes have subsequently been duplicated in teleost fishes. The evolution of the PRKACG retroposon in simians was also investigated. Although the degree of sequence conservation in the PKA Cα/Cβ kinase family is exceptionally high, a small set of signature residues defining Cα and Cβ subfamilies were identified. These conserved residues might be important for functions that are unique to the Cα or Cβ clades. This study also provides a good example of a seemingly simple phylogenetic problem which, due to a very high degree of sequence conservation and corresponding weak phylogenetic signals, combined with problematic nonphylogenetic signals, is nontrivial for state-of-the-art probabilistic phylogenetic methods.