Search

Abstract Fructosamine-3-kinases (FN3Ks) are a conserved family of repair enzymes that phosphorylate reactive sugars attached to lysine residues in peptides and proteins. Although FN3Ks are present across the Tree of Life and share detectable sequence similarity to eukaryotic protein kinases, the biological processes regulated by these kinases are largely unknown. To address this knowledge gap, we leveraged the FN3K CRISPR Knock-Out (KO) HepG2 cell line alongside an integrative multi-omics study combining transcriptomics, metabolomics, and interactomics to place these enzymes in a pathway context. The integrative analyses revealed the enrichment of pathways related to oxidative stress response, lipid biosynthesis (cholesterol and fatty acids), and carbon and co-factor metabolism. Moreover, enrichment of nicotinamide adenine dinucleotide (NAD) binding proteins and localization of human FN3K (HsFN3K) to mitochondria suggests potential links between FN3K and NAD-mediated energy metabolism and redox balance. We report specific binding of HsFN3K to NAD compounds in a metal and concentration-dependent manner and provide insight into their binding mode using modeling and experimental site-directed mutagenesis. Our studies provide a framework for targeting these understudied kinases in diabetic complications and metabolic disorders where redox balance and NAD-dependent metabolic processes are altered.

Abstract The necroptosis pathway is a lytic, pro-inflammatory mode of cell death that is widely implicated in human disease, including renal, pulmonary, gut and skin inflammatory pathologies. The precise mechanism of the terminal steps in the pathway, where the RIPK3 kinase phosphorylates and triggers a conformation change and oligomerization of the terminal pathway effector, MLKL, are only emerging. Here, we structurally identify RIPK3-mediated phosphorylation of the human MLKL activation loop as a cue for MLKL pseudokinase domain dimerization. MLKL pseudokinase domain dimerization subsequently drives formation of elongated homotetramers. Negative stain electron microscopy and modelling support nucleation of the MLKL tetramer assembly by a central coiled coil formed by the extended, ~80 Å brace helix that connects the pseudokinase and executioner four-helix bundle domains. Mutational data assert MLKL tetramerization as an essential prerequisite step to enable the release and reorganization of four-helix bundle domains for membrane permeabilization and cell death.

Abstract Autophosphorylation controls the transition between discrete functional and conformational states in protein kinases, yet the structural and molecular determinants underlying this fundamental process remain unclear. Here we show that c-terminal Tyr 530 is a de facto c-Src autophosphorylation site with slow time-resolution kinetics and a strong intermolecular component. On the contrary, activation-loop Tyr 419 undergoes faster kinetics and a cis-to-trans phosphorylation switch that controls c-terminal Tyr 530 autophosphorylation, enzyme specificity, and strikingly, c-Src non-catalytic function as a substrate. In line with this, we visualize by X-ray crystallography a snapshot of Tyr 530 intermolecular autophosphorylation. In an asymmetric arrangement of both catalytic domains, a c-terminal palindromic phospho-motif flanking Tyr 530 on the substrate molecule engages the G-loop of the active kinase adopting a position ready for entry into the catalytic cleft. Perturbation of the phospho-motif accounts for c-Src dysfunction as indicated by viral and colorectal cancer (CRC)-associated c-terminal deleted variants. We show that c-terminal residues 531 to 536 are required for c-Src Tyr 530 autophosphorylation, and such a detrimental effect is caused by the substrate molecule inhibiting allosterically the active kinase. Our work reveals a crosstalk between the activation and c-terminal segments that control the allosteric interplay between substrate- and enzyme-acting kinases during autophosphorylation.

The Protein Kinase Ontology (ProKinO) is an integrated knowledge graph that conceptualizes the complex relationships among protein kinase sequence, structure, function, and disease in a human and machine-readable format. In this study, we have significantly expanded ProKinO by incorporating additional data on expression patterns and drug interactions. Furthermore, we have developed a completely new browser from the ground up to render the knowledge graph visible and interactive on the web. We have enriched ProKinO with new classes and relationships that capture information on kinase ligand binding sites, expression patterns, and functional features. These additions extend ProKinO’s capabilities as a discovery tool, enabling it to uncover novel insights about understudied members of the protein kinase family. We next demonstrate the application of ProKinO. Specifically, through graph mining and aggregate SPARQL queries, we identify the p21-activated protein kinase 5 (PAK5) as one of the most frequently mutated dark kinases in human cancers with abnormal expression in multiple cancers, including a previously unappreciated role in acute myeloid leukemia. We have identified recurrent oncogenic mutations in the PAK5 activation loop predicted to alter substrate binding and phosphorylation. Additionally, we have identified common ligand/drug binding residues in PAK family kinases, underscoring ProKinO’s potential application in drug discovery. The updated ontology browser and the addition of a web component, ProtVista, which enables interactive mining of kinase sequence annotations in 3D structures and Alphafold models, provide a valuable resource for the signaling community. The updated ProKinO database is accessible at https://prokino.uga.edu.

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The doublecortin-like kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory ‘tail’ segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations, and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to ‘supercharge’ the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other calcium calmodulin kinases (CAMKs), and a ‘Swiss Army’ assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for autoregulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically divergent DCLK1 modulators, stabilizers, or degraders.

The 534 protein kinases encoded in the human genome constitute a large druggable class of proteins that include both well-studied and understudied “dark” members. Accurate prediction of dark kinase functions is a major bioinformatics challenge. Here, we employ a graph mining approach that uses the evolutionary and functional context encoded in knowledge graphs (KGs) to predict protein and pathway associations for understudied kinases. We propose a new scalable graph embedding approach, RegPattern2Vec, which employs regular pattern constrained random walks to sample diverse aspects of node context within a KG flexibly. RegPattern2Vec learns functional representations of kinases, interacting partners, post-translational modifications, pathways, cellular localization, and chemical interactions from a kinase-centric KG that integrates and conceptualizes data from curated heterogeneous data resources. By contextualizing information relevant to prediction, RegPattern2Vec improves accuracy and efficiency in comparison to other random walk-based graph embedding approaches. We show that the predictions produced by our model overlap with pathway enrichment data produced using experimentally validated Protein-Protein Interaction (PPI) data from both publicly available databases and experimental datasets not used in training. Our model also has the advantage of using the collected random walks as biological context to interpret the predicted protein-pathway associations. We provide high-confidence pathway predictions for 34 dark kinases and present three case studies in which analysis of meta-paths associated with the prediction enables biological interpretation. Overall, RegPattern2Vec efficiently samples multiple node types for link prediction on biological knowledge graphs and the predicted associations between understudied kinases, pseudokinases, and known pathways serve as a conceptual starting point for hypothesis generation and testing.

ABSTRACT Pneumococcal pneumonia remains a WHO high-priority disease despite multivalent conjugate vaccines administered in clinical practice worldwide. A protein-based, serotype-independent vaccine has long-promised comprehensive coverage of most clinical isolates of the pneumococcus. Along with numerous pneumococcal surface protein immunogens, the pneumococcal serine-rich repeat protein (PsrP) has been investigated as a potential vaccine target due to its surface exposure and functions toward bacterial virulence and lung infection. Three critical criteria for its vaccine potential — the clinical prevalence, serotype distribution, and sequence homology of PsrP — have yet to be well characterized. Here, we used genomes of 13,454 clinically isolated pneumococci from the Global Pneumococcal Sequencing project to investigate PsrP presence among isolates, distribution among serotypes, and interrogate its homology as a protein across species. These isolates represent all age groups, countries worldwide, and types of pneumococcal infection. We found PsrP present in at least 50% of all isolates across all determined serotypes and nontypeable (NT) clinical isolates. Using a combination of peptide matching and HMM profiles built on full-length and individual PsrP domains, we identified novel variants that expand PsrP diversity and prevalence. We also observed sequence variability in its basic region (BR) between isolates and serotypes. PsrP has a strong vaccine potential due to its breadth of coverage, especially in nonvaccine serotypes (NVTs) when exploiting its regions of conservation in vaccine design. IMPORTANCE An updated outlook on PsrP prevalence and serotype distribution sheds new light on the comprehensiveness of a PsrP-based protein vaccine. The protein is present in all vaccine serotypes and highly present in the next wave of potentially disease-causing serotypes not included in the current multivalent conjugate vaccines. Furthermore, PsrP is strongly correlated with clinical isolates harboring pneumococcal disease as opposed to pneumococcal carriage. PsrP is also highly present in strains and serotypes from Africa, where the need for a protein-based vaccine is the greatest, giving new reasoning to pursue PsrP as a protein vaccine.

Abstract Background Protein kinases are among the largest druggable family of signaling proteins, involved in various human diseases, including cancers and neurodegenerative disorders. Despite their clinical relevance, nearly 30% of the 545 human protein kinases remain highly understudied. Comparative genomics is a powerful approach for predicting and investigating the functions of understudied kinases. However, an incomplete knowledge of kinase orthologs across fully sequenced kinomes severely limits the application of comparative genomics approaches for illuminating understudied kinases. Here, we introduce KinOrtho, a query- and graph-based orthology inference method that combines full-length and domain-based approaches to map one-to-one kinase orthologs across 17 thousand species. Results Using multiple metrics, we show that KinOrtho performed better than existing methods in identifying kinase orthologs across evolutionarily divergent species and eliminated potential false positives by flagging sequences without a proper kinase domain for further evaluation. We demonstrate the advantage of using domain-based approaches for identifying domain fusion events, highlighting a case between an understudied serine/threonine kinase TAOK1 and a metabolic kinase PIK3C2A with high co-expression in human cells. We also identify evolutionary fission events involving the understudied OBSCN kinase domains, further highlighting the value of domain-based orthology inference approaches. Using KinOrtho-defined orthologs, Gene Ontology annotations, and machine learning, we propose putative biological functions of several understudied kinases, including the role of TP53RK in cell cycle checkpoint(s), the involvement of TSSK3 and TSSK6 in acrosomal vesicle localization, and potential functions for the ULK4 pseudokinase in neuronal development. Conclusions In sum, KinOrtho presents a novel query-based tool to identify one-to-one orthologous relationships across thousands of proteomes that can be applied to any protein family of interest. We exploit KinOrtho here to identify kinase orthologs and show that its well-curated kinome ortholog set can serve as a valuable resource for illuminating understudied kinases, and the KinOrtho framework can be extended to any protein-family of interest.

Glycosyltransferases (GT) are proteins that display extensive sequence and functional variation on a subset of 3D folds. Here, the authors use interpretable deep learning to predict 3D folds from sequence without the need for sequence alignment, which also enables the prediction of GTs with new folds.

Mixed Lineage Kinase Domain-Like (MLKL) pseudokinase is phosphorylated by RIPK3 kinase prior to cell death by necroptosis. Here, the authors use monobodies that bind to the MLKL pseudokinase domain as tools, which allowed them to determine the crystal structures of the MLKL pseudokinase domain in two distinct conformations. By combining their structural data with cell signalling assays and MD simulations they provide evidence that endogenous MLKL preassociates with its upstream regulator RIPK3, and that MLKL disengages from RIPK3 following the induction of necroptosis.

The viral Protein Kinase-1 (PK-1) phosphorylates the regulatory protein p6.9, which facilitates baculoviral genome release. Here, the authors combine X-ray crystallography with biophysical and biochemical analyses as well as molecular dynamics simulations to characterize Cydia pomenella granulovirus PK-1, which forms a dimer with a parallel side-to-side arrangement of the kinase domains and furthermore, they provide insights into its catalytic mechanism and evolutionary relationships with other kinases.

Aim: The present study was planned to clinically determine the effectiveness of colchicine as a monotherapy in oral submucous fibrosis (OSMF) patients. Materials and Methods: Thirty OSMF patients were enrolled in the study and randomly divided into two groups by permuted-block randomization. Group A: Patients were treated with 0.5 mg colchicine (Goutnil) tablets twice daily for 3 months. Whereas, in Group B, intralesional injection of dexamethasone 2 mL, hyaluronidase 1,500 IU with 2 mL lignocaine HCl biweekly for 4–6 weeks were given. For all the individuals, baseline parameters like mouth opening, buccal mucosal flexibility, burning sensation, and tongue protrusion were recorded and reassessment was done at 1, 3, and 6 months followed by statistical analysis. Results: A significant decrease in the Visual Analogue Scale (VAS) score for severity of burning sensation was observed in Group A whereas mouth opening and tongue protrusion were higher in Group B patients. Statistically significant differences in the buccal mucosal flexibility were appreciable in the intragroup comparisons of both the groups from baseline to 6-month follow-up. Conclusion : Colchicine can be suggested as an adjuvant drug for reducing the burning sensation for the OSMF patients in whom steroids are contraindicated.

Abstract Background Protein kinases are a large family of druggable proteins that are genomically and proteomically altered in many human cancers. Kinase-targeted drugs are emerging as promising avenues for personalized medicine because of the differential response shown by altered kinases to drug treatment in patients and cell-based assays. However, an incomplete understanding of the relationships connecting genome, proteome and drug sensitivity profiles present a major bottleneck in targeting kinases for personalized medicine. Results In this study, we propose a multi-component Quantitative Structure–Mutation–Activity Relationship Tests (QSMART) model and neural networks framework for providing explainable models of protein kinase inhibition and drug response ( $$\hbox {IC}_{50}$$ IC 50 ) profiles in cell lines. Using non-small cell lung cancer as a case study, we show that interaction terms that capture associations between drugs, pathways, and mutant kinases quantitatively contribute to the response of two EGFR inhibitors (afatinib and lapatinib). In particular, protein–protein interactions associated with the JNK apoptotic pathway, associations between lung development and axon extension, and interaction terms connecting drug substructures and the volume/charge of mutant residues at specific structural locations contribute significantly to the observed $$\hbox {IC}_{50}$$ IC 50 values in cell-based assays. Conclusions By integrating multi-omics data in the QSMART model, we not only predict drug responses in cancer cell lines with high accuracy but also identify features and explainable interaction terms contributing to the accuracy. Although we have tested our multi-component explainable framework on protein kinase inhibitors, it can be extended across the proteome to investigate the complex relationships connecting genotypes and drug sensitivity profiles.

Protein prenylation by farnesyltransferase (FTase) is often described as the targeting of a cysteine-containing motif (CaaX) that is enriched for aliphatic amino acids at the a1 and a2 positions, while quite flexible at the X position. Prenylation prediction methods often rely on these features despite emerging evidence that FTase has broader target specificity than previously considered. Using a machine learning approach and training sets based on canonical (prenylated, proteolyzed, and carboxymethylated) and recently identified shunted motifs (prenylation only), this study aims to improve prenylation predictions with the goal of determining the full scope of prenylation potential among the 8000 possible Cxxx sequence combinations. Further, this study aims to subdivide the prenylated sequences as either shunted (i.e., uncleaved) or cleaved (i.e., canonical). Predictions were determined for Saccharomyces cerevisiae FTase and compared to results derived using currently available prenylation prediction methods. In silico predictions were further evaluated using in vivo methods coupled to two yeast reporters, the yeast mating pheromone a-factor and Hsp40 Ydj1p, that represent proteins with canonical and shunted CaaX motifs, respectively. Our machine learning-based approach expands the repertoire of predicted FTase targets and provides a framework for functional classification.

The necroptotic cell death pathway involves signaling through pseudokinases. Here the authors define the structural determinants of species specificity in necroptosis signaling mediated by the essential necroptotic effector pseudokinase, Mixed Lineage Kinase Domain-Like (MLKL).

Glycosyltransferases (GTs) are prevalent across the tree of life and regulate nearly all aspects of cellular functions. The evolutionary basis for their complex and diverse modes of catalytic functions remain enigmatic. Here, based on deep mining of over half million GT-A fold sequences, we define a minimal core component shared among functionally diverse enzymes. We find that variations in the common core and emergence of hypervariable loops extending from the core contributed to GT-A diversity. We provide a phylogenetic framework relating diverse GT-A fold families for the first time and show that inverting and retaining mechanisms emerged multiple times independently during evolution. Using evolutionary information encoded in primary sequences, we trained a machine learning classifier to predict donor specificity with nearly 90% accuracy and deployed it for the annotation of understudied GTs. Our studies provide an evolutionary framework for investigating complex relationships connecting GT-A fold sequence, structure, function and regulation.

The length of cilia is controlled by a poorly understood mechanism that involves members of the conserved RCK kinase group, and among them, the LF4/MOK kinases. The multiciliated protist model, Tetrahymena, carries two types of cilia (oral and locomotory) and the length of the locomotory cilia is dependent on their position with the cell. In Tetrahymena, loss of an LF4/MOK ortholog, LF4A, lengthened the locomotory cilia, but also reduced their number. Without LF4A, cilia assembled faster and showed signs of increased intraflagellar transport (IFT). Consistently, overproduced LF4A shortened cilia and downregulated IFT. GFP-tagged LF4A, expressed in the native locus and imaged by total internal reflection microscopy, was enriched at the basal bodies and distributed along the shafts of cilia. Within cilia, most LF4A-GFP particles were immobile and a few either diffused or moved by IFT. We suggest that the distribution of LF4/MOK along the cilium delivers a uniform dose of inhibition to IFT trains that travel from the base to the tip. In a longer cilium, the IFT machinery may experience a higher cumulative dose of inhibition by LF4/MOK. Thus, LF4/MOK activity could be a readout of cilium length that helps to balance the rate of IFT-driven assembly with the rate of disassembly at steady state. We used a forward genetic screen to identify a CDK-related kinase, CDKR1, whose loss-of-function suppressed the shortening of cilia caused by overexpression of LF4A, by reducing its kinase activity. Loss of CDKR1 alone lengthened both the locomotory and oral cilia. CDKR1 resembles other known ciliary CDK-related kinases: LF2 of Chlamydomonas, mammalian CCRK and DYF-18 of C. elegans, in lacking the cyclin-binding motif and acting upstream of RCKs. The new genetic tools we developed here for Tetrahymena have potential for further dissection of the principles of cilia length regulation in multiciliated cells.

ABSTRACT As the protozoan parasite Toxoplasma gondii disseminates through its host, it responds to environmental changes by altering its gene expression, metabolism, and other processes. Oxygen is one variable environmental factor, and properly adapting to changes in oxygen levels is critical to prevent the accumulation of reactive oxygen species and other cytotoxic factors. Thus, oxygen-sensing proteins are important, and among these, 2-oxoglutarate-dependent prolyl hydroxylases are highly conserved throughout evolution. Toxoplasma expresses two such enzymes, TgPHYa, which regulates the SCF-ubiquitin ligase complex, and TgPHYb. To characterize TgPHYb, we created a Toxoplasma strain that conditionally expresses TgPHYb and report that TgPHYb is required for optimal parasite growth under normal growth conditions. However, exposing TgPHYb-depleted parasites to extracellular stress leads to severe decreases in parasite invasion, which is likely due to decreased abundance of parasite adhesins. Adhesin protein abundance is reduced in TgPHYb-depleted parasites as a result of inactivation of the protein synthesis elongation factor eEF2 that is accompanied by decreased rates of translational elongation. In contrast to most other oxygen-sensing proteins that mediate cellular responses to low O2, TgPHYb is specifically required for parasite growth and protein synthesis at high, but not low, O2 tensions as well as resistance to reactive oxygen species. In vivo, reduced TgPHYb expression leads to lower parasite burdens in oxygen-rich tissues. Taken together, these data identify TgPHYb as a sensor of high O2 levels, in contrast to TgPHYa, which supports the parasite at low O2. IMPORTANCE Because oxygen plays a key role in the growth of many organisms, cells must know how much oxygen is available. O2-sensing proteins are therefore critical cellular factors, and prolyl hydroxylases are the best-studied type of O2-sensing proteins. In general, prolyl hydroxylases trigger cellular responses to decreased oxygen availability. But, how does a cell react to high levels of oxygen? Using the protozoan parasite Toxoplasma gondii, we discovered a prolyl hydroxylase that allows the parasite to grow at elevated oxygen levels and does so by regulating protein synthesis. Loss of this enzyme also reduces parasite burden in oxygen-rich tissues, indicating that sensing both high and low levels of oxygen impacts the growth and physiology of Toxoplasma.

[This corrects the article DOI: 10.1371/journal.pgen.1005885.].

Protein tyrosine kinases (PTKs) are a group of closely related enzymes that have evolutionarily diverged from serine/threonine kinases (STKs) to regulate pathways associated with multi-cellularity. Evolutionary divergence of PTKs from STKs has occurred through accumulation of mutations in the active site as well as in the commonly conserved hydrophobic core. While the functional significance of active site variations is well understood, relatively little is known about how hydrophobic core variations contribute to PTK evolutionary divergence. Here, using a combination of statistical sequence comparisons, molecular dynamics simulations, mutational analysis and in vitro thermostability and kinase assays, we investigate the structural and functional significance of key PTK-specific variations in the kinase core. We find that the nature of residues and interactions in the hydrophobic core of PTKs is strikingly different from other protein kinases, and PTK-specific variations in the core contribute to functional divergence by altering the stability and dynamics of the kinase domain. In particular, a functionally critical STK-conserved histidine that stabilizes the regulatory spine in STKs is selectively mutated to an alanine, serine or glutamate in PTKs, and this loss-of-function mutation is accommodated, in part, through compensatory PTK-specific interactions in the core. In particular, a PTK-conserved phenylalanine in the I-helix appears to structurally and functionally compensate for the loss of STK-histidine by interacting with the regulatory spine, which has far-reaching effects on enzyme activity, inhibitor sensing, and stability. We propose that hydrophobic core variations provide a selective advantage during PTK evolution by increasing the conformational flexibility, and therefore the allosteric potential of the kinase domain. Our studies also suggest that Tyrosine Kinase Like kinases such as RAF are intermediates in PTK evolutionary divergence inasmuch as they share features of both PTKs and STKs in the core. Finally, our studies provide an evolutionary framework for identifying and characterizing disease and drug resistance mutations in the kinase core.

A novel form of copy number control (CNC) helps maintain a low number of Ty1 retrovirus-like transposons in the Saccharomyces genome. Ty1 produces an alternative transcript that encodes p22, a trans-dominant negative inhibitor of Ty1 retrotransposition whose sequence is identical to the C-terminal half of Gag. The level of p22 increases with copy number and inhibits normal Ty1 virus-like particle (VLP) assembly and maturation through interactions with full length Gag. A forward genetic screen for CNC-resistant (CNCR) mutations in Ty1 identified missense mutations in GAG that restore retrotransposition in the presence of p22. Some of these mutations map within a predicted UBN2 domain found throughout the Ty1/copia family of long terminal repeat retrotransposons, and others cluster within a central region of Gag that is referred to as the CNCR domain. We generated multiple alignments of yeast Ty1-like Gag proteins and found that some Gag proteins, including those of the related Ty2 elements, contain non-Ty1 residues at multiple CNCR sites. Interestingly, the Ty2-917 element is resistant to p22 and does not undergo a Ty1-like form of CNC. Substitutions conferring CNCR map within predicted helices in Ty1 Gag that overlap with conserved sequence in Ty1/copia, suggesting that p22 disturbs a central function of the capsid during VLP assembly. When hydrophobic residues within predicted helices in Gag are mutated, Gag level remains unaffected in most cases yet VLP assembly and maturation is abnormal. Gag CNCR mutations do not alter binding to p22 as determined by co-immunoprecipitation analyses, but instead, exclude p22 from Ty1 VLPs. These findings suggest that the CNCR alleles enhance retrotransposition in the presence of p22 by allowing productive Gag-Gag interactions during VLP assembly. Our work also expands the strategies used by retroviruses for developing resistance to Gag-like restriction factors to now include retrotransposons.

The epidermal growth factor receptor (EGFR) is overexpressed in multiple carcinomas and is the focus of a variety of targeted therapies. Here we report the design of peptide-based compounds that mimic the EGFR dimerization arm and inhibit allosteric activation of EGFR. These peptides are modified to contain a triazolyl bridge between the peptide strands to constrain the EGFR dimerization arm β-loop. In this study, we demonstrate that these peptides have significantly improved proteolytic stability over the non-modified peptide sequence, and their inhibitory effects are dependent on the number of the methylene units and orientation of the introduced triazolyl bridge. We identified a peptide, EDA2, which downregulates receptor phosphorylation and dimerization and reduces cell viability. This is the first example of a biologically active triazolyl-bridged peptide targeting the EGFR dimerization interface that effectively downregulates EGFR activation.

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.

Cancer is a genetic disease that develops through a series of somatic mutations, a subset of which drive cancer progression. Although cancer genome sequencing studies are beginning to reveal the mutational patterns of genes in various cancers, identifying the small subset of "causative" mutations from the large subset of "non-causative" mutations, which accumulate as a consequence of the disease, is a challenge. In this article, we present an effective machine learning approach for identifying cancer-associated mutations in human protein kinases, a class of signaling proteins known to be frequently mutated in human cancers. We evaluate the performance of 11 well known supervised learners and show that a multiple-classifier approach, which combines the performances of individual learners, significantly improves the classification of known cancer-associated mutations. We introduce several novel features related specifically to structural and functional characteristics of protein kinases and find that the level of conservation of the mutated residue at specific evolutionary depths is an important predictor of oncogenic effect. We consolidate the novel features and the multiple-classifier approach to prioritize and experimentally test a set of rare unconfirmed mutations in the epidermal growth factor receptor tyrosine kinase (EGFR). Our studies identify T725M and L861R as rare cancer-associated mutations inasmuch as these mutations increase EGFR activity in the absence of the activating EGF ligand in cell-based assays.

Eukaryotic protein kinases (EPKs) regulate numerous signaling processes by phosphorylating targeted substrates through the highly conserved catalytic domain. Our previous computational studies proposed a model stating that a properly assembled nonlinear motif termed the Regulatory (R) spine is essential for catalytic activity of EPKs. Here we define the required intramolecular interactions and biochemical properties of the R-spine and the newly identified "Shell" that surrounds the R-spine using site-directed mutagenesis and various in vitro phosphoryl transfer assays using cyclic AMP-dependent protein kinase as a representative of the entire kinome. Analysis of the 172 available Apo EPK structures in the protein data bank (PDB) revealed four unique structural conformations of the R-spine that correspond with catalytic inactivation of various EPKs. Elucidating the molecular entities required for the catalytic activation of EPKs and the identification of these inactive conformations opens new avenues for the design of efficient therapeutic EPK inhibitors.

BACKGROUND: Protein kinases are a large and diverse family of enzymes that are genomically altered in many human cancers. Targeted cancer genome sequencing efforts have unveiled the mutational profiles of protein kinase genes from many different cancer types. While mutational data on protein kinases is currently catalogued in various databases, integration of mutation data with other forms of data on protein kinases such as sequence, structure, function and pathway is necessary to identify and characterize key cancer causing mutations. Integrative analysis of protein kinase data, however, is a challenge because of the disparate nature of protein kinase data sources and data formats. RESULTS: Here, we describe ProKinO, a protein kinase-specific ontology, which provides a controlled vocabulary of terms, their hierarchy, and relationships unifying sequence, structure, function, mutation and pathway information on protein kinases. The conceptual representation of such diverse forms of information in one place not only allows rapid discovery of significant information related to a specific protein kinase, but also enables large-scale integrative analysis of protein kinase data in ways not possible through other kinase-specific resources. We have performed several integrative analyses of ProKinO data and, as an example, found that a large number of somatic mutations (∼288 distinct mutations) associated with the haematopoietic neoplasm cancer type map to only 8 kinases in the human kinome. This is in contrast to glioma, where the mutations are spread over 82 distinct kinases. We also provide examples of how ontology-based data analysis can be used to generate testable hypotheses regarding cancer mutations. CONCLUSION: We present an integrated framework for large-scale integrative analysis of protein kinase data. Navigation and analysis of ontology data can be performed using the ontology browser available at: http://vulcan.cs.uga.edu/prokino.

BACKGROUND: The epidermal growth factor receptor kinases, or ErbB kinases, belong to a large sub-group of receptor tyrosine kinases (RTKs), which share a conserved catalytic core. The catalytic core of ErbB kinases have functionally diverged from other RTKs in that they are activated by a unique allosteric mechanism that involves specific interactions between the kinase core and the flanking Juxtamembrane (JM) and COOH-terminal tail (C-terminal tail). Although extensive studies on ErbB and related tyrosine kinases have provided important insights into the structural basis for ErbB kinase functional divergence, the sequence features that contribute to the unique regulation of ErbB kinases have not been systematically explored. METHODOLOGY/PRINCIPAL FINDINGS: In this study, we use a Bayesian approach to identify the selective sequence constraints that most distinguish ErbB kinases from other receptor tyrosine kinases. We find that strong ErbB kinase-specific constraints are imposed on residues that tether the JM and C-terminal tail to key functional regions of the kinase core. A conserved RIxKExE motif in the JM-kinase linker region and a glutamine in the inter-lobe linker are identified as two of the most distinguishing features of the ErbB family. While the RIxKExE motif tethers the C-terminal tail to the N-lobe of the kinase domain, the glutamine tethers the C-terminal tail to hinge regions critical for inter-lobe movement. Comparison of the active and inactive crystal structures of ErbB kinases indicates that the identified residues are conformationally malleable and can potentially contribute to the cis regulation of the kinase core by the JM and C-terminal tail. ErbB3, and EGFR orthologs in sponges and parasitic worms, diverge from some of the canonical ErbB features, providing insights into sub-family and lineage-specific functional specialization. CONCLUSION/SIGNIFICANCE: Our analysis pinpoints key residues for mutational analysis, and provides new clues to cancer mutations that alter the canonical modes of ErbB kinase regulation.

The eukaryotic protein kinase (ePK) domain mediates the majority of signaling and coordination of complex events in eukaryotes. By contrast, most bacterial signaling is thought to occur through structurally unrelated histidine kinases, though some ePK-like kinases (ELKs) and small molecule kinases are known in bacteria. Our analysis of the Global Ocean Sampling (GOS) dataset reveals that ELKs are as prevalent as histidine kinases and may play an equally important role in prokaryotic behavior. By combining GOS and public databases, we show that the ePK is just one subset of a diverse superfamily of enzymes built on a common protein kinase-like (PKL) fold. We explored this huge phylogenetic and functional space to cast light on the ancient evolution of this superfamily, its mechanistic core, and the structural basis for its observed diversity. We cataloged 27,677 ePKs and 18,699 ELKs, and classified them into 20 highly distinct families whose known members suggest regulatory functions. GOS data more than tripled the count of ELK sequences and enabled the discovery of novel families and classification and analysis of all ELKs. Comparison between and within families revealed ten key residues that are highly conserved across families. However, all but one of the ten residues has been eliminated in one family or another, indicating great functional plasticity. We show that loss of a catalytic lysine in two families is compensated by distinct mechanisms both involving other key motifs. This diverse superfamily serves as a model for further structural and functional analysis of enzyme evolution.

There have been significant developments in the use of siRNA in the silencing of cancer-allied target proteins with substantial apoptotic effects. Nevertheless, the challenges regarding siRNA delivery for cancer therapy remain a major concern for taking these therapies successfully from laboratory to in vivo studies. Biomaterials are defined traditionally as any material that is used for either medical or dental applications that contact the host cells in any form, such as a drug carrier a device, or as a prosthesis towards the replacement of damaged tissues. The main issues to be fulfilled by the biomaterial for clinical applications are biocompatibility, bioactivity, ability to carry drug to target site, inflammatory responses and other factors based on its application. The present study focuses on the in vitro inflammatory response to the CS nanoparticles using RAW 264.7 and bone marrow derived macrophage cells. Additionally, the in vitro release kinetics of siRNA with varying concentrations and pH, transfection efficacy and biocompatibility were also investigated. The results of siRNA cumulative release increased at pH 5 and 3, which may be corresponding to the protonation, and a delayed release was seen at 7, which was ascribed to unprotonated amine groups inside the CS. The results of release kinetics confirmed a sustained release of siRNA from CS NPs. Considering that CS is a biocompatible polymer, it typically has little impact/damage on cells, as numerous researchers have observed during in vitro experiments. Inflammatory studies were carried out in vitro with RAW 264.7 and BMC cells derived from mice. The gene and protein expression studies showed that the materials might cause some slight inflammation on exposure with both RAW 264.7 and BMC cells in vitro, which is completely negligible. However, putting together the overall data it can be concluded that CS NPs can be a promising material for in vivo applications, which is in agreement with the results of other researchers, but the only concern being its ability to carry siRNA and protect it from nuclease and other enzymatic attacks. [ABSTRACT FROM AUTHOR]

Chaperones TAPBPR and tapasin associate with class I major histocompatibility complexes (MHC-I) to promote optimization (editing) of peptide cargo. Here, we use solution NMR to investigate the mechanism of peptide exchange. We identify TAPBPR-induced conformational changes on conserved MHC-I molecular surfaces, consistent with our independently determined X-ray structure of the complex. Dynamics present in the empty MHC-I are stabilized by TAPBPR and become progressively dampened with increasing peptide occupancy. Incoming peptides are recognized according to the global stability of the final pMHC-I product and anneal in a native-like conformation to be edited by TAPBPR. Our results demonstrate an inverse relationship between MHC-I peptide occupancy and TAPBPR binding affinity, wherein the lifetime and structural features of transiently bound peptides control the regulation of a conformational switch located near the TAPBPR binding site, which triggers TAPBPR release. These results suggest a similar mechanism for the function of tapasin in the peptide-loading complex.

The molecular mechanism through which the interaction of a clonotypic αβ T-cell receptor (TCR) with a peptide-loaded major histocompatibility complex (p/MHC) leads to T-cell activation is not yet fully understood. Here we exploit a high-affinity TCR (B4.2.3) to examine the structural changes that accompany binding to its p/MHC ligand (P18-I10/H2-Dd). In addition to conformational changes in complementarity-determining regions (CDRs) of the TCR seen in comparison of unliganded and bound X-ray structures, NMR characterization of the TCR β-chain dynamics reveals significant chemical shift effects in sites removed from the MHC-binding site. Remodelling of electrostatic interactions near the Cβ H3 helix at the membrane-proximal face of the TCR, a region implicated in interactions with the CD3 co-receptor, suggests a possible role for an allosteric mechanism in TCR signalling. The contribution of these TCR residues to signal transduction is supported by mutagenesis and T-cell functional assays.