Evolution of domain combinations in protein kinases and its implications for functional diversity

Proteome-Wide Detection and Annotation of Receptor Tyrosine Kinases (RTKs): RTK-PRED and the TyReK Database

Receptor tyrosine kinases (RTKs) form a highly important group of protein receptors of the eukaryotic cell membrane. They control many vital cellular functions and are involved in the regulation of complex signaling networks. Mutations in RTKs have been associated with different types of cancers and other diseases. Although they are very important for proper cell function, they have been experimentally studied in a limited range of eukaryotic species. Currently, there is no available database for RTKs providing information about their function, expression, and interactions. Therefore, the identification of RTKs in multiple organisms, the documentation of their characteristics, and the collection of related information would be very useful. In this paper, we present a novel RTK detection pipeline (RTK-PRED) and the Receptor Tyrosine Kinases Database (TyReK-DB). RTK-PRED combines profile HMMs with transmembrane topology prediction to identify and classify potential RTKs. Proteins of all eukaryotic reference proteomes of the UniProt database were used as input in RTK-PRED leading to a filtered dataset of 20,478 RTKs. Based on the information collected for these RTKs from multiple databases, the relational TyReK database was created.

Approved Small-Molecule ATP-Competitive Kinases Drugs Containing Indole/Azaindole/Oxindole Scaffolds: R&D and Binding Patterns Profiling

Kinases are among the most important families of biomolecules and play an essential role in the regulation of cell proliferation, apoptosis, metabolism, and other critical physiological processes. The dysregulation and gene mutation of kinases are linked to the occurrence and development of various human diseases, especially cancer. As a result, a growing number of small-molecule drugs based on kinase targets are being successfully developed and approved for the treatment of many diseases. The indole/azaindole/oxindole moieties are important key pharmacophores of many bioactive compounds and are generally used as excellent scaffolds for drug discovery in medicinal chemistry. To date, 30 ATP-competitive kinase inhibitors bearing the indole/azaindole/oxindole scaffold have been approved for the treatment of diseases. Herein, we summarize their research and development (R&D) process and describe their binding models to the ATP-binding sites of the target kinases. Moreover, we discuss the significant role of the indole/azaindole/oxindole skeletons in the interaction of their parent drug and target kinases, providing new medicinal chemistry inspiration and ideas for the subsequent development and optimization of kinase inhibitors.

Proteomics of Salt Gland-Secreted Sap Indicates a Pivotal Role for Vesicle Transport and Energy Metabolism in Plant Salt Secretion

Soil salinization is one of the major factors restricting crop growth and agricultural production worldwide. Recretohalophytes have developed unique epidermal structures in their aboveground tissues, such as salt glands or salt bladders, to secrete excess salt out of the plant body as a protective mechanism from ion damage. Three hypotheses were proposed to explain how salt glands secrete salts: the osmotic hypothesis, a hypothesis similar to animal fluid transport, and vesicle-mediated exocytosis. However, there is no direct evidence to show whether the salt gland-secreted liquid contains landmark proteins or peptides which would elucidate the salt secretion mechanism. In this study, we collected the secreted liquid of salt glands from Limonium bicolor, followed by extraction and identification of its constituent proteins and peptides by SDS-PAGE and mass spectrometry. We detected 214 proteins and 440 polypeptides in the salt gland-secreted droplets of plants grown under control conditions. Unexpectedly, the proportion of energy metabolism-related proteins increased significantly though only 16 proteins and 35 polypeptides in the droplets of salt-treated plants were detected. In addition, vesicle transport proteins such as the Golgi marker enzyme glycosyltransferase were present in the secreted sap of salt glands from both control and salt-treated plants. These results suggest that trans-Golgi network-mediated vesicular transport and energy production contributes to salt secretion in salt glands.

Update and elucidation of Plasmodium kinomes: Prioritization of kinases as potential drug targets for malaria

Malaria is a tropical disease caused by Plasmodium spp. and transmitted by the bite of infected Anopheles mosquitoes. Protein kinases (PKs) play key roles in the life cycle of the etiological agent of malaria, turning these proteins attractive targets for antimalarial drug discovery campaigns. As part of an effort to understand parasite signaling functions, we report the results of a bioinformatics pipeline analysis of PKs of eight Plasmodium species. To date, no P. malariae and P. ovale kinome assemble has been conducted. We classified, curated and annotated predicted kinases to update P. falciparum, P. vivax, P. yoelii, P. berghei, P. chabaudi, and P. knowlesi kinomes published to date, as well as report for the first time the kinomes of P. malariae and P. ovale. Overall, from 76 to 97 PKs were identified among all Plasmodium spp. kinomes. Most of the kinases were assigned to seven of nine major kinase groups: AGC, CAMK, CMGC, CK1, STE, TKL, OTHER; and the Plasmodium-specific group FIKK. About 30% of kinases have been deeply classified into group, family and subfamily levels and only about 10% remained unclassified. Furthermore, updating and comparing the kinomes of P. vivax and P. falciparum allowed for the prioritization and selection of kinases as potential drug targets that could be explored for discovering new drugs against malaria. This integrated approach resulted in the selection of 37 protein kinases as potential targets and the identification of investigational compounds with moderate in vitro activity against asexual P. falciparum (3D7 and Dd2 strains) stages that could serve as starting points for the search of potent antimalarial leads in the future.

Targeting the Non-Catalytic Functions: a New Paradigm for Kinase Drug Discovery?

Protein kinases have been highly fruitful targets for cancer drug discovery in the past two decades, while most of these drugs bind to the "adenosine triphosphate (ATP)-site" and inhibit kinase catalytic activity. Recently, accumulated evidence suggests that kinases possess functions beyond catalysis through their scaffolds, and the scaffolding functions could play critical roles in multiple cellular signaling and cell fate controls. Small molecules modulating the noncatalytic functions of kinases are rarely reported but emerge as new promising therapeutic strategies for various diseases. Herein, we summarize the characterized noncatalytic functions of kinases, and highlight the recent progress on developing small-molecule modulators of the noncatalytic functions of kinases. Mechanisms and characteristics of different kinds of modulators are also discussed. It is also speculated that targeting the noncatalytic functions would represent a new direction for kinase-based drug discovery.

Chemogenomics and bioinformatics approaches for prioritizing kinases as drug targets for neglected tropical diseases

Neglected tropical diseases (NTDs) are a group of twenty-one diseases classified by the World Health Organization that prevail in regions with tropical and subtropical climate and affect more than one billion people. There is an urgent need to develop new and safer drugs for these diseases. Protein kinases are a potential class of targets for developing new drugs against NTDs, since they play crucial role in many biological processes, such as signaling pathways, regulating cellular communication, division, metabolism and death. Bioinformatics is a field that aims to organize large amounts of biological data as well as develop and use tools for understanding and analyze them in order to produce meaningful information in a biological manner. In combination with chemogenomics, which analyzes chemical-biological interactions to screen ligands against selected targets families, these approaches can be used to stablish a rational strategy for prioritizing new drug targets for NTDs. Here, we describe how bioinformatics and chemogenomics tools can help to identify protein kinases and their potential inhibitors for the development of new drugs for NTDs. We present a review of bioinformatics tools and techniques that can be used to define an organisms kinome for drug prioritization, drug and target repurposing, multi-quinase inhibition approachs and selectivity profiling. We also present some successful examples of the application of such approaches in recent case studies.

Phosphoglycerate kinase: structural aspects and functions, with special emphasis on the enzyme from Kinetoplastea

Phosphoglycerate kinase (PGK) is a glycolytic enzyme that is well conserved among the three domains of life. PGK is usually a monomeric enzyme of about 45 kDa that catalyses one of the two ATP-producing reactions in the glycolytic pathway, through the conversion of 1,3-bisphosphoglycerate (1,3BPGA) to 3-phosphoglycerate (3PGA). It also participates in gluconeogenesis, catalysing the opposite reaction to produce 1,3BPGA and ADP. Like most other glycolytic enzymes, PGK has also been catalogued as a moonlighting protein, due to its involvement in different functions not associated with energy metabolism, which include pathogenesis, interaction with nucleic acids, tumorigenesis progression, cell death and viral replication. In this review, we have highlighted the overall aspects of this enzyme, such as its structure, reaction kinetics, activity regulation and possible moonlighting functions in different protistan organisms, especially both free-living and parasitic Kinetoplastea. Our analysis of the genomes of different kinetoplastids revealed the presence of open-reading frames (ORFs) for multiple PGK isoforms in several species. Some of these ORFs code for unusually large PGKs. The products appear to contain additional structural domains fused to the PGK domain. A striking aspect is that some of these PGK isoforms are predicted to be catalytically inactive enzymes or ‘dead’ enzymes. The roles of PGKs in kinetoplastid parasites are analysed, and the apparent significance of the PGK gene duplication that gave rise to the different isoforms and their expression in Trypanosoma cruzi is discussed.

Polycystic Kidney Disease Ryanodine Receptor Domain (PKDRR) Proteins in Oomycetes

In eukaryotes, two sources of Ca²⁺ are accessed to allow rapid changes in the cytosolic levels of this second messenger: the extracellular medium and intracellular Ca²⁺ stores, such as the endoplasmic reticulum. One class of channel that permits Ca²⁺ entry is the transient receptor potential (TRP) superfamily, including the polycystic kidney disease (PKD) proteins, or polycystins. Channels that release Ca²⁺ from intracellular stores include the inositol 1,4,5-trisphosphate/ryanodine receptor (ITPR/RyR) superfamily. Here, we characterise a family of proteins that are only encoded by oomycete genomes, that we have named PKDRR, since they share domains with both PKD and RyR channels. We provide evidence that these proteins belong to the TRP superfamily and are distinct from the ITPR/RyR superfamily in terms of their evolutionary relationships, protein domain architectures and predicted ion channel structures. We also demonstrate that a hypothetical PKDRR protein from Phytophthora infestans is produced by this organism, is located in the cell-surface membrane and forms multimeric protein complexes. Efforts to functionally characterise this protein in a heterologous expression system were unsuccessful but support a cell-surface localisation. These PKDRR proteins represent potential targets for the development of new "fungicides", since they are of a distinctive structure that is only found in oomycetes and not in any other cellular organisms.

Unveiling the Kinomes of Leishmania infantum and L. braziliensis Empowers the Discovery of New Kinase Targets and Antileishmanial Compounds

Leishmaniasis is a neglected tropical disease caused by parasites of the genus Leishmania (NTD) endemic in 98 countries. Although some drugs are available, current treatments deal with issues such as toxicity, low efficacy, and emergence of resistance. Therefore, there is an urgent need to identify new targets for the development of new antileishmanial drugs. Protein kinases (PKs), which play an essential role in many biological processes, have become potential drug targets for many parasitic diseases. A refined bioinformatics pipeline was applied in order to define and compare the kinomes of L. infantum and L. braziliensis, species that cause cutaneous and visceral manifestations of leishmaniasis in the Americas, the latter being potentially fatal if untreated. Respectively, 224 and 221 PKs were identified in L. infantum and L. braziliensis overall. Almost all unclassified eukaryotic PKs were assigned to six of nine major kinase groups and, consequently, most have been classified into family and subfamily. Furthermore, revealing the kinomes for both Leishmania species allowed for the prioritization of potential drug targets that could be explored for discovering new drugs against leishmaniasis. Finally, we used a drug repurposing approach and prioritized seven approved drugs and investigational compounds to be experimentally tested against Leishmania. Trametinib and NMS-1286937 inhibited the growth of L. infantum and L. braziliensis promastigotes and amastigotes and therefore might be good candidates for the drug repurposing pipeline.

Activation loop targeting strategy for design of receptor-interacting protein kinase 2 (RIPK2) inhibitors

Development of selective kinase inhibitors remains a challenge due to considerable amino acid sequence similarity among family members particularly in the ATP binding site. Targeting the activation loop might offer improved inhibitor selectivity since this region of kinases is less conserved. However, the strategy presents difficulties due to activation loop flexibility. Herein, we report the design of receptor-interacting protein kinase 2 (RIPK2) inhibitors based on pan-kinase inhibitor regorafenib that aim to engage basic activation loop residues Lys169 or Arg171. We report development of CSR35 that displayed >10-fold selective inhibition of RIPK2 versus VEGFR2, the target of regorafenib. A co-crystal structure of CSR35 with RIPK2 revealed a resolved activation loop with an ionic interaction between the carboxylic acid installed in the inhibitor and the side-chain of Lys169. Our data provides principle feasibility of developing activation loop targeting type II inhibitors as a complementary strategy for achieving improved selectivity.

Classification of Protein Kinases Influenced by Conservation of Substrate Binding Residues

With the advent of genome sequencing projects in the recent past, several kinases have come to light as regulating different signaling pathways. These kinases are generally classified into different subfamilies based on their sequence similarity with members of known subfamilies of kinases. A functional association is then defined to the kinase based on the subfamily to which it has been characterized. However, one of the key factors that give identity to a kinase in a subfamily is its ability to phosphorylate a given set of substrates. Substrate specificity of a kinase is largely determined by the residues at the substrate binding site. Though in general the sequence similarity based measure for classification more or less gives the preliminary idea on subfamily, understanding the molecular basis of kinase substrate recognition could further refine the classification scheme for kinases and render a better understanding of their functional role. In this analysis we emphasize on the possibility of using putative substrate binding information in the classification of a given kinase into a particular subfamily.

Dynamics of genomic innovation in the unicellular ancestry of animals

Which genomic innovations underpinned the origin of multicellular animals is still an open debate. Here, we investigate this question by reconstructing the genome architecture and gene family diversity of ancestral premetazoans, aiming to date the emergence of animal-like traits. Our comparative analysis involves genomes from animals and their closest unicellular relatives (the Holozoa), including four new genomes: three Ichthyosporea and Corallochytrium limacisporum. Here, we show that the earliest animals were shaped by dynamic changes in genome architecture before the emergence of multicellularity: an early burst of gene diversity in the ancestor of Holozoa, enriched in transcription factors and cell adhesion machinery, was followed by multiple and differently-timed episodes of synteny disruption, intron gain and genome expansions. Thus, the foundations of animal genome architecture were laid before the origin of complex multicellularity – highlighting the necessity of a unicellular perspective to understand early animal evolution.

DOI:http://dx.doi.org/10.7554/eLife.26036.001

Hundreds of millions of years ago, some single-celled organisms gained the ability to work together and form multicellular organisms. This transition was a major step in evolution and took place at separate times in several parts of the tree of life, including in animals, plants, fungi and algae.

Animals are some of the most complex organisms on Earth. Their single-celled ancestors were also quite genetically complex themselves and their genomes (the complete set of the organism’s DNA) already contained many genes that now coordinate the activity of the cells in a multicellular organism.

The genome of an animal typically has certain features: it is large, diverse and contains many segments (called introns) that are not genes. By seeing if the single-celled relatives of animals share these traits, it is possible to learn more about when specific genetic features first evolved, and whether they are linked to the origin of animals.

Now, Grau-Bové et al. have studied the genomes of several of the animal kingdom’s closest single-celled relatives using a technique called whole genome sequencing. This revealed that there was a period of rapid genetic change in the single-celled ancestors of animals during which their genes became much more diverse. Another ‘explosion’ of diversity happened after animals had evolved. Furthermore, the overall amount of the genomic content inside cells and the number of introns found in the genome rapidly increased in separate, independent events in both animals and their single-celled ancestors.

Future research is needed to investigate whether other multicellular life forms – such as plants, fungi and algae – originated in the same way as animal life. Understanding how the genetic material of animals evolved also helps us to understand the genetic structures that affect our health. For example, genes that coordinate the behavior of cells (and so are important for multicellular organisms) also play a role in cancer, where cells break free of this regulation to divide uncontrollably.

DOI:http://dx.doi.org/10.7554/eLife.26036.002

The Kinome of Pacific Oyster Crassostrea gigas, Its Expression during Development and in Response to Environmental Factors

Oysters play an important role in estuarine and coastal marine habitats, where the majority of humans live. In these ecosystems, environmental degradation is substantial, and oysters must cope with highly dynamic and stressful environmental constraints during their lives in the intertidal zone. The availability of the genome sequence of the Pacific oyster Crassostrea gigas represents a unique opportunity for a comprehensive assessment of the signal transduction pathways that the species has developed to deal with this unique habitat. We performed an in silico analysis to identify, annotate and classify protein kinases in C. gigas, according to their kinase domain taxonomy classification, and compared with kinome already described in other animal species. The C. gigas kinome consists of 371 protein kinases, making it closely related to the sea urchin kinome, which has 353 protein kinases. The absence of gene redundancy in some groups of the C. gigas kinome may simplify functional studies of protein kinases. Through data mining of transcriptomes in C. gigas, we identified part of the kinome which may be central during development and may play a role in response to various environmental factors. Overall, this work contributes to a better understanding of key sensing pathways that may be central for adaptation to a highly dynamic marine environment.

The Genomes of Three Uneven Siblings: Footprints of the Lifestyles of Three Trichoderma Species

The genus Trichoderma contains fungi with high relevance for humans, with applications in enzyme production for plant cell wall degradation and use in biocontrol. Here, we provide a broad, comprehensive overview of the genomic content of these species for "hot topic" research aspects, including CAZymes, transport, transcription factors, and development, along with a detailed analysis and annotation of less-studied topics, such as signal transduction, genome integrity, chromatin, photobiology, or lipid, sulfur, and nitrogen metabolism in T. reesei, T. atroviride, and T. virens, and we open up new perspectives to those topics discussed previously. In total, we covered more than 2,000 of the predicted 9,000 to 11,000 genes of each Trichoderma species discussed, which is >20% of the respective gene content. Additionally, we considered available transcriptome data for the annotated genes. Highlights of our analyses include overall carbohydrate cleavage preferences due to the different genomic contents and regulation of the respective genes. We found light regulation of many sulfur metabolic genes. Additionally, a new Golgi 1,2-mannosidase likely involved in N-linked glycosylation was detected, as were indications for the ability of Trichoderma spp. to generate hybrid galactose-containing N-linked glycans. The genomic inventory of effector proteins revealed numerous compounds unique to Trichoderma, and these warrant further investigation. We found interesting expansions in the Trichoderma genus in several signaling pathways, such as G-protein-coupled receptors, RAS GTPases, and casein kinases. A particularly interesting feature absolutely unique to T. atroviride is the duplication of the alternative sulfur amino acid synthesis pathway.

Evolutionary and structural analyses of heterodimeric proteins composed of subunits with same fold

Heterodimeric proteins with homologous subunits of same fold are involved in various biological processes. The objective of this study is to understand the evolution of structural and functional features of such heterodimers. Using a non-redundant dataset of 70 such heterodimers of known 3D structure and an independent dataset of 173 heterodimers from yeast, we note that the mean sequence identity between interacting homologous subunits is only 23-24% suggesting that, generally, highly diverged paralogues assemble to form such a heterodimer. We also note that the functional roles of interacting subunits/domains are generally quite different. This suggests that, though the interacting subunits/domains are homologous, the high evolutionary divergence characterize their high functional divergence which contributes to a gross function for the heterodimer considered as a whole. The inverse relationship between sequence identity and RMSD of interacting homologues in heterodimers is not followed. We also addressed the question of formation of homodimers of the subunits of heterodimers by generating models of fictitious homodimers on the basis of the 3D structures of the heterodimers. Interaction energies associated with these homodimers suggests that, in overwhelming majority of the cases, such homodimers are unlikely to be stable. Majority of the homologues of heterodimers of known structures form heterodimers (51.8%) and a small proportion (14.6%) form homodimers. Comparison of 3D structures of heterodimers with homologous homodimers suggests that interfacial nature of residues is not well conserved. In over 90% of the cases we note that the interacting subunits of heterodimers are co-localized in the cell.

Small molecule substrate phosphorylation site inhibitors of protein kinases: approaches and challenges

Protein kinases are important mediators of cellular communication and attractive drug targets for many diseases. Although success has been achieved with developing ATP-competitive kinase inhibitors, the disadvantages of ATP-competitive inhibitors have led to increased interest in targeting sites outside of the ATP binding pocket. Kinase inhibitors with substrate-competitive, ATP-noncompetitive binding modes are promising due to the possibility of increased selectivity and better agreement between biochemical and in vitro potency. However, the difficulty of identifying these types of inhibitors has resulted in significantly fewer small molecule substrate phosphorylation site inhibitors being reported compared to ATP-competitive inhibitors. This review surveys reported substrate phosphorylation site inhibitors and methods that can be applied to the discovery of such inhibitors, including a discussion of the challenges inherent to these screening methods.

Sweet taste receptor signaling network: possible implication for cognitive functioning

Sweet taste receptors are transmembrane protein network specialized in the transmission of information from special "sweet" molecules into the intracellular domain. These receptors can sense the taste of a range of molecules and transmit the information downstream to several acceptors, modulate cell specific functions and metabolism, and mediate cell-to-cell coupling through paracrine mechanism. Recent reports indicate that sweet taste receptors are widely distributed in the body and serves specific function relative to their localization. Due to their pleiotropic signaling properties and multisubstrate ligand affinity, sweet taste receptors are able to cooperatively bind multiple substances and mediate signaling by other receptors. Based on increasing evidence about the role of these receptors in the initiation and control of absorption and metabolism, and the pivotal role of metabolic (glucose) regulation in the central nervous system functioning, we propose a possible implication of sweet taste receptor signaling in modulating cognitive functioning.

Typical and atypical domain combinations in human protein kinases: functions, disease causing mutations and conservation in other primates

A twist in the evolution of human kinases resulting in kinases with hybrid and rogue properties.

The eukaryotic ancestor had a complex ubiquitin signaling system of archaeal origin

The origin of the eukaryotic cell is one of the most important transitions in the history of life. However, the emergence and early evolution of eukaryotes remains poorly understood. Recent data have shown that the last eukaryotic common ancestor (LECA) was much more complex than previously thought. The LECA already had the genetic machinery encoding the endomembrane apparatus, spliceosome, nuclear pore, and myosin and kinesin cytoskeletal motors. It is unclear, however, when the functional regulation of these cellular components evolved. Here, we address this question by analyzing the origin and evolution of the ubiquitin (Ub) signaling system, one of the most important regulatory layers in eukaryotes. We delineated the evolution of the whole Ub, Small-Ub-related MOdifier (SUMO), and Ub-fold modifier 1 (Ufm1) signaling networks by analyzing representatives from all major eukaryotic, bacterial, and archaeal lineages. We found that the Ub toolkit had a pre-eukaryotic origin and is present in three extant archaeal groups. The pre-eukaryotic Ub toolkit greatly expanded during eukaryogenesis, through massive gene innovation and diversification of protein domain architectures. This resulted in a LECA with essentially all of the Ub-related genes, including the SUMO and Ufm1 Ub-like systems. Ub and SUMO signaling further expanded during eukaryotic evolution, especially labeling and delabeling enzymes responsible for substrate selection. Additionally, we analyzed protein domain architecture evolution and found that multicellular lineages have the most complex Ub systems in terms of domain architectures. Together, we demonstrate that the Ub system predates the origin of eukaryotes and that a burst of innovation during eukaryogenesis led to a LECA with complex posttranslational regulation.

Neem Limonoids as Anticancer Agents: Modulation of Cancer Hallmarks and Oncogenic Signaling

Neem (Azadirachta indica A. Juss) is one of the most versatile medicinal plants, widely distributed in the Indian subcontinent. Neem is a rich source of limonoids that are endowed with potent medicinal properties predominantly antioxidant, anti-inflammatory, and anticancer activities. Azadirachtin, gedunin, and nimbolide are more extensively investigated relative to other neem limonoids. Accumulating evidence indicates that the anticancer effects of neem limonoids are mediated through the inhibition of hallmark capabilities of cancer such as cell proliferation, apoptosis evasion, inflammation, invasion, and angiogenesis. The neem limonoids have been demonstrated to target oncogenic signaling kinases and transcription factors chiefly, NF-κB, Wnt/β-catenin, PI3K/Akt, MAPK, and JAK/STAT signaling pathways. Neem limonoids that target multiple pathways that are aberrant in cancer are ideal candidates for cancer chemoprevention and therapy.

Tethering preferences of domain families co-occurring in multi-domain proteins

Genomic data of several organisms have revealed the presence of a vast repertoire of multi-domain proteins. The role played by individual domains in a multi-domain protein has a profound influence on the overall function of the protein. In the present analysis an attempt has been made to better understand the tethering preferences of domain families that occur in multi-domain proteins. The analysis has been carried out on an exhaustive dataset of 2 961 898 sequences of proteins from 930 organisms, where 741 274 proteins are comprised of at least two domain families. For every domain family, the number of other domain families with which it co-occurs within a protein in this dataset has been enumerated and is referred to as the tethering number of the domain family. It was found that, in the general dataset, the AAA ATPase family and the family of Ser/Thr kinases have the highest tethering numbers of 450 and 444 respectively. Further analysis reveals significant correlation between the number of members in a family and its tethering number. Positive correlation was also observed for the extent of a sequence and functional diversity within a family and the tethering numbers of domain families. Domain families that are present ubiquitously in diverse organisms tend to have large tethering numbers, while organism/kingdom-specific families have low tethering numbers. Thus, the analysis uncovers how domain families recombine and evolve to give rise to multi-domain proteins.

Design principles underpinning the regulatory diversity of protein kinases

Protein phosphorylation in eukaryotes is carried out by a large and diverse family of protein kinases, which display remarkable diversity and complexity in their modes of regulation. The complex modes of regulation have evolved as a consequence of natural selection operating on protein kinase sequences for billions of years. Here we describe how quantitative comparisons of protein kinase sequences from diverse organisms, in particular prokaryotes, have contributed to our understanding of the structural organization and evolution of allosteric regulation in the protein kinase domain. An emerging view from these studies is that regulatory diversity and complexity in the protein kinase domain evolved in a ‘modular’ fashion through elaboration of an ancient core component, which existed before the emergence of eukaryotes. The core component provided the conformational flexibility required for ATP binding and phosphoryl transfer in prokaryotic kinases, but evolved into a highly regulatable domain in eukaryotes through the addition of exaggerated structural features that facilitated tight allosteric control. Family and group-specific features are built upon the core component in eukaryotes to provide additional layers of control. We propose that ‘modularity’ and ‘conformational flexibility’ are key evolvable traits of the protein kinase domain that contributed to its extensive regulatory diversity and complexity.

Structural and mechanistic insights into the bifunctional enzyme isocitrate dehydrogenase kinase/phosphatase AceK

The switch between the Krebs cycle and the glyoxylate bypass is controlled by isocitrate dehydrogenase kinase/phosphatase (AceK). AceK, a bifunctional enzyme, phosphorylates and dephosphorylates isocitrate dehydrogenase (IDH) with its unique active site that harbours both the kinase and ATP/ADP-dependent phosphatase activities. AceK was the first example of prokaryotic phosphorylation identified, and the recent characterization of the structures of AceK and its complex with its protein substrate, IDH, now offers a new understanding of both previous and future endeavours. AceK is structurally similar to the eukaryotic protein kinase superfamily, sharing many of the familiar catalytic and regulatory motifs, demonstrating a close evolutionary relationship. Although the active site is shared by both the kinase and phosphatase functions, the catalytic residues needed for phosphatase function are readily seen when compared with the DXDX(T/V) family of phosphatases, despite the fact that the phosphatase function of AceK is strictly ATP/ADP-dependent. Structural analysis has also allowed a detailed look at regulation and its stringent requirements for interacting with IDH.

Diversification of a protein kinase cascade: IME-2 is involved in nonself recognition and programmed cell death in Neurospora crassa

Kinase cascades and the modification of proteins by phosphorylation are major mechanisms for cell signaling and communication, and evolution of these signaling pathways can contribute to new developmental or environmental response pathways. The Saccharomyces cerevisiae kinase Ime2 has been well characterized for its role in meiosis. However, recent studies have revealed alternative functions for Ime2 in both S. cerevisiae and other fungi. In the filamentous fungus Neurospora crassa, the IME2 homolog (ime-2) is not required for meiosis. Here we determine that ime-2 interacts genetically with a transcription factor vib-1 during nonself recognition and programmed cell death (PCD). Mutations in vib-1 (Δvib-1) suppress PCD due to nonself recognition events; however, a Δvib-1 Δime-2 mutant restored wild-type levels of cell death. A role for ime-2 in the post-translational processing and localization of a mitochondrial matrix protein was identified, which may implicate mitochondria in N. crassa nonself recognition and PCD. Further, Δvib-1 strains do not produce extracellular proteases, but protease secretion reverted to near wild-type levels in a Δvib-1 Δime-2 strain. Mass spectrometry analysis revealed that the VIB-1 protein is phosphorylated at several sites, including a site that matches the IME-2 consensus. The genetic and biochemical data for ime-2 and vib-1 indicate that IME-2 is a negative regulator of VIB-1 and suggest parallel negative regulation by IME-2 of a cell death pathway in N. crassa that functions in concert with the VIB-1 cell death pathway. Thus, IME2 kinase function has evolved following the divergence of S. cerevisiae and N. crassa and provides insight into the evolution of kinases and their regulatory targets.

Tracing the origin and evolutionary history of plant nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes

Plant disease resistance genes (R genes) encode proteins that function to monitor signals indicating pathogenic infection, thus playing a critical role in the plant's defense system. Although many studies have been performed to explore the functional details of these important genes, their origin and evolutionary history remain unclear. In this study, focusing on the largest group of R genes, the nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes, we conducted an extensive genome-wide survey of 38 representative model organisms and obtained insights into the evolutionary stage and timing of NBS-LRR genes. Our data show that the two major domains, NBS and LRR, existed before the split of prokaryotes and eukaryotes but their fusion was observed only in land plant lineages. The Toll/interleukin-1 receptor (TIR) class of NBS-LRR genes probably had an earlier origin than its nonTIR counterpart. The similarities of the innate immune systems of plants and animals are likely to have been shaped by convergent evolution after their independent origins. Our findings start to unravel the evolutionary history of these important genes from the perspective of comparative genomics and also highlight the important role of reorganizing pre-existing building blocks in generating evolutionary novelties.

The secret life of kinases: functions beyond catalysis

Protein phosphorylation participates in the regulation of all fundamental biological processes, and protein kinases have been intensively studied. However, while the focus was on catalytic activities, accumulating evidence suggests that non-catalytic properties of protein kinases are essential, and in some cases even sufficient for their functions. These non-catalytic functions include the scaffolding of protein complexes, the competition for protein interactions, allosteric effects on other enzymes, subcellular targeting, and DNA binding. This rich repertoire often is used to coordinate phosphorylation events and enhance the specificity of substrate phosphorylation, but also can adopt functions that do not rely on kinase activity. Here, we discuss such kinase independent functions of protein and lipid kinases focussing on kinases that play a role in the regulation of cell proliferation, differentiation, apoptosis, and motility.

Functional Diversity of the Schistosoma mansoni Tyrosine Kinases

Schistosoma mansoni, one of the causative agents of schistosomiasis, has a complex life cycle infecting over 200 million people worldwide. Such a successful and prolific parasite life cycle has been shown to be dependent on the adaptive interaction between the parasite and hosts. Tyrosine kinases (TKs) play a key role in signaling pathways as demonstrated by a large body of experimental work in eukaryotes. Furthermore, comparative genomics have allowed the identification of TK homologs and provided insights into the functional role of TKs in several biological systems. Finally, TK structural biology has provided a rational basis for obtaining selective inhibitors directed to the treatment of human diseases. This paper covers the important aspects of the phospho-tyrosine signaling network in S. mansoni, Caenorhabditis elegans, and humans, the main process of functional diversification of TKs, that is, protein-domain shuffling, and also discusses TKs as targets for the development of new anti-schistosome drugs.

Eukaryotic protein kinases (ePKs) of the helminth parasite Schistosoma mansoni

BACKGROUND

Schistosomiasis remains an important parasitic disease and a major economic problem in many countries. The Schistosoma mansoni genome and predicted proteome sequences were recently published providing the opportunity to identify new drug candidates. Eukaryotic protein kinases (ePKs) play a central role in mediating signal transduction through complex networks and are considered druggable targets from the medical and chemical viewpoints. Our work aimed at analyzing the S. mansoni predicted proteome in order to identify and classify all ePKs of this parasite through combined computational approaches. Functional annotation was performed mainly to yield insights into the parasite signaling processes relevant to its complex lifestyle and to select some ePKs as potential drug targets.

RESULTS

We have identified 252 ePKs, which corresponds to 1.9% of the S. mansoni predicted proteome, through sequence similarity searches using HMMs (Hidden Markov Models). Amino acid sequences corresponding to the conserved catalytic domain of ePKs were aligned by MAFFT and further used in distance-based phylogenetic analysis as implemented in PHYLIP. Our analysis also included the ePK homologs from six other eukaryotes. The results show that S. mansoni has proteins in all ePK groups. Most of them are clearly clustered with known ePKs in other eukaryotes according to the phylogenetic analysis. None of the ePKs are exclusively found in S. mansoni or belong to an expanded family in this parasite. Only 16 S. mansoni ePKs were experimentally studied, 12 proteins are predicted to be catalytically inactive and approximately 2% of the parasite ePKs remain unclassified. Some proteins were mentioned as good target for drug development since they have a predicted essential function for the parasite.

CONCLUSIONS

Our approach has improved the functional annotation of 40% of S. mansoni ePKs through combined similarity and phylogenetic-based approaches. As we continue this work, we will highlight the biochemical and physiological adaptations of S. mansoni in response to diverse environments during the parasite development, vector interaction, and host infection.

The kinome of Phytophthora infestans reveals oomycete-specific innovations and links to other taxonomic groups

BACKGROUND

Oomycetes are a large group of economically and ecologically important species. Its most notorious member is Phytophthora infestans, the cause of the devastating potato late blight disease. The life cycle of P. infestans involves hyphae which differentiate into spores used for dispersal and host infection. Protein phosphorylation likely plays crucial roles in these stages, and to help understand this we present here a genome-wide analysis of the protein kinases of P. infestans and several relatives. The study also provides new insight into kinase evolution since oomycetes are taxonomically distant from organisms with well-characterized kinomes.

RESULTS

Bioinformatic searches of the genomes of P. infestans, P. ramorum, and P. sojae reveal they have similar kinomes, which for P. infestans contains 354 eukaryotic protein kinases (ePKs) and 18 atypical kinases (aPKs), equaling 2% of total genes. After refining gene models, most were classifiable into families seen in other eukaryotes. Some ePK families are nevertheless unusual, especially the tyrosine kinase-like (TKL) group which includes large oomycete-specific subfamilies. Also identified were two tyrosine kinases, which are rare in non-metazoans. Several ePKs bear accessory domains not identified previously on kinases, such as cyclin-dependent kinases with integral cyclin domains. Most ePKs lack accessory domains, implying that many are regulated transcriptionally. This was confirmed by mRNA expression-profiling studies that showed that two-thirds vary significantly between hyphae, sporangia, and zoospores. Comparisons to neighboring taxa (apicomplexans, ciliates, diatoms) revealed both clade-specific and conserved features, and multiple connections to plant kinases were observed. The kinome of Hyaloperonospora arabidopsidis, an oomycete with a simpler life cycle than P. infestans, was found to be one-third smaller. Some differences may be attributable to gene clustering, which facilitates subfamily expansion (or loss) through unequal crossing-over.

CONCLUSION

The large sizes of the Phytophthora kinomes imply that phosphorylation plays major roles in their life cycles. Their kinomes also include many novel ePKs, some specific to oomycetes or shared with neighboring groups. Little experimentation to date has addressed the biological functions of oomycete kinases, but this should be stimulated by the structural, evolutionary, and expression data presented here. This may lead to targets for disease control.

Classification of protein kinases on the basis of both kinase and non-kinase regions

Background

Protein phosphorylation is a generic way to regulate signal transduction pathways in all kingdoms of life. In many organisms, it is achieved by the large family of Ser/Thr/Tyr protein kinases which are traditionally classified into groups and subfamilies on the basis of the amino acid sequence of their catalytic domains. Many protein kinases are multi-domain in nature but the diversity of the accessory domains and their organization are usually not taken into account while classifying kinases into groups or subfamilies.

Methodology

Here, we present an approach which considers amino acid sequences of complete gene products, in order to suggest refinements in sets of pre-classified sequences. The strategy is based on alignment-free similarity scores and iterative Area Under the Curve (AUC) computation. Similarity scores are computed by detecting common patterns between two sequences and scoring them using a substitution matrix, with a consistent normalization scheme. This allows us to handle full-length sequences, and implicitly takes into account domain diversity and domain shuffling. We quantitatively validate our approach on a subset of 212 human protein kinases. We then employ it on the complete repertoire of human protein kinases and suggest few qualitative refinements in the subfamily assignment stored in the KinG database, which is based on catalytic domains only. Based on our new measure, we delineate 37 cases of potential hybrid kinases: sequences for which classical classification based entirely on catalytic domains is inconsistent with the full-length similarity scores computed here, which implicitly consider multi-domain nature and regions outside the catalytic kinase domain. We also provide some examples of hybrid kinases of the protozoan parasite Entamoeba histolytica.

Conclusions

The implicit consideration of multi-domain architectures is a valuable inclusion to complement other classification schemes. The proposed algorithm may also be employed to classify other families of enzymes with multi-domain architecture.