Prediction of 492 human protein kinase substrate specificities

Background

Complex intracellular signaling networks monitor diverse environmental inputs to evoke appropriate and coordinated effector responses. Defective signal transduction underlies many pathologies, including cancer, diabetes, autoimmunity and about 400 other human diseases. Therefore, there is high impetus to define the composition and architecture of cellular communications networks in humans. The major components of intracellular signaling networks are protein kinases and protein phosphatases, which catalyze the reversible phosphorylation of proteins. Here, we have focused on identification of kinase-substrate interactions through prediction of the phosphorylation site specificity from knowledge of the primary amino acid sequence of the catalytic domain of each kinase.

Results

The presented method predicts 488 different kinase catalytic domain substrate specificity matrices in 478 typical and 4 atypical human kinases that rely on both positive and negative determinants for scoring individual phosphosites for their suitability as kinase substrates. This represents a marked advancement over existing methods such as those used in NetPhorest (179 kinases in 76 groups) and NetworKIN (123 kinases), which consider only positive determinants for kinase substrate prediction. Comparison of our predicted matrices with experimentally-derived matrices from about 9,000 known kinase-phosphosite substrate pairs revealed a high degree of concordance with the established preferences of about 150 well studied protein kinases. Furthermore for many of the better known kinases, the predicted optimal phosphosite sequences were more accurate than the consensus phosphosite sequences inferred by simple alignment of the phosphosites of known kinase substrates.

Conclusions

Application of this improved kinase substrate prediction algorithm to the primary structures of over 23, 000 proteins encoded by the human genome has permitted the identification of about 650, 000 putative phosphosites, which are posted on the open source PhosphoNET website (http://www.phosphonet.ca).

     

Using multitask classification methods to investigate the kinase-specific phosphorylation sites

Background

Identification of phosphorylation sites by computational methods is becoming increasingly important because it reduces labor-intensive and costly experiments and can improve our understanding of the common properties and underlying mechanisms of protein phosphorylation.

Methods

A multitask learning framework for learning four kinase families simultaneously, instead of studying each kinase family of phosphorylation sites separately, is presented in the study. The framework includes two multitask classification methods: the Multi-Task Least Squares Support Vector Machines (MTLS-SVMs) and the Multi-Task Feature Selection (MT-Feat3).

Results

Using the multitask learning framework, we successfully identify 18 common features shared by four kinase families of phosphorylation sites. The reliability of selected features is demonstrated by the consistent performance in two multi-task learning methods.

Conclusions

The selected features can be used to build efficient multitask classifiers with good performance, suggesting they are important to protein phosphorylation across 4 kinase families.

     

Recurrent Nicolau syndrome associated with subcutaneous glatiramer acetate injection—a case report

Background

Glatiramer acetate is worldwide used as first line treatment in relapsing remitting multiple sclerosis. Local skin reactions associated with glatiramer acetate are common, however, only isolated cases of severe local injection site reactions known as Nicolau Syndrome have been reported so far.

Case presentation

We describe the case of a recurrent Nicolau Syndrome occurred during longstanding glatiramer acetate treatment in a woman with multiple sclerosis. The haemorrhagic patch necrotized and was treated locally as a deep second degree burn with excision of dead skin tissue and was healed. Treatment with glatiramer acetate was definitely suspended.

Conclusions

GA injections can be complicated by isolated or recurrent Nicolau Syndrome, a potentially life-threatening condition of which neurologists should be aware.

     

Triple A patient cells suffering from mitotic defects fail to localize PGRMC1 to mitotic kinetochore fibers

Background

Membrane-associated progesterone receptors are restricted to the endoplasmic reticulum and are shown to regulate the activity of cytochrome P450 enzymes which are involved in steroidogenesis or drug detoxification. PGRMC1 and PGRMC2 belong to the membrane-associated progesterone receptor family and are of interest due to their suspected role during cell cycle. PGRMC1 and PGRMC2 are thought to bind to each other; thereby suppressing entry into mitosis. We could previously report that PGRMC2 interacts with the nucleoporin ALADIN which when mutated results in the autosomal recessive disorder triple A syndrome. ALADIN is a novel regulator of mitotic controller Aurora kinase A and depletion of this nucleoporin leads to microtubule instability.

Results

In the current study, we present that proliferation is decreased when ALADIN, PGRMC1 or PGRMC2 are over-expressed. Furthermore, we find that depletion of ALADIN results in mislocalization of Aurora kinase A and PGRMC1 in metaphase cells. Additionally, PGRMC2 is over-expressed in triple A patient fibroblasts.

Conclusion

Our results emphasize the possibility that loss of the regulatory association between ALADIN and PGRMC2 gives rise to a depletion of PGRMC1 at kinetochore fibers. This observation may explain part of the symptoms seen in triple A syndrome patients.

     

Global phosphoproteomic analysis identifies SRMS-regulated secondary signaling intermediates

Background

The non-receptor tyrosine kinase, SRMS (Src-related kinase lacking C-terminal regulatory tyrosine and N-terminal myristoylation sites) is a member of the BRK family kinases (BFKs) which represents an evolutionarily conserved relative of the Src family kinases (SFKs). Tyrosine kinases are known to regulate a number of cellular processes and pathways via phosphorylating substrate proteins directly and/or by partaking in signaling cross-talks leading to the indirect modulation of various signaling intermediates. In a previous study, we profiled the tyrosine-phosphoproteome of SRMS and identified multiple candidate substrates of the kinase. The broader cellular signaling intermediates of SRMS are unknown.

Methods

In order to uncover the broader SRMS-regulated phosphoproteome and identify the SRMS-regulated indirect signaling intermediates, we performed label-free global phosphoproteomics analysis on cells expressing wild-type SRMS. Using computational database searching and bioinformatics analyses we characterized the dataset.

Results

Our analyses identified 60 hyperphosphorylated (phosphoserine/phosphothreonine) proteins mapped from 140 hyperphosphorylated peptides. Bioinfomatics analyses identified a number of significantly enriched biological and cellular processes among which DNA repair pathways were found to be upregulated while apoptotic pathways were found to be downregulated. Analyses of motifs derived from the upregulated phosphosites identified Casein kinase 2 alpha (CK2α) as one of the major potential kinases contributing to the SRMS-dependent indirect regulation of signaling intermediates.

Conclusions

Overall, our phosphoproteomics analyses identified serine/threonine phosphorylation dynamics as important secondary events of the SRMS-regulated phosphoproteome with implications in the regulation of cellular and biological processes.

     

Weekly oral alendronate in mevalonate kinase deficiency

Background

Mevalonate kinase deficiency (MKD) is caused by mutations in the MVK gene, encoding the second enzyme of mevalonate pathway, which results in subsequent shortage of downstream compounds, and starts in childhood with febrile attacks, skin, joint, and gastrointestinal symptoms, sometimes induced by vaccinations.

Methods

For a history of early-onset corticosteroid-induced reduction of bone mineral density in a 14-year-old boy with MKD, who also had presented three bone fractures, we administered weekly oral alendronate, a drug widely used in the management of osteoporosis and other high bone turnover diseases, which blocks mevalonate and halts the prenylation process.

Results

All of the patient’s MKD clinical and laboratory abnormalities were resolved after starting alendronate treatment.

Conclusions

This observation appears enigmatic, since alendronate should reinforce the metabolic block characterizing MKD, but is crucial because of the ultimate improvement shown by this patient. The anti-inflammatory properties of bisphosphonates are a new question for debate among physicians across various specialties, and requires further biochemical and clinical investigation.

     

Optimization of fecal sample preparation for untargeted LC-HRMS based metabolomics

Introduction

Collecting feces is easy. It offers direct outcome to endogenous and microbial metabolites.

Objectives

In a context of lack of consensus about fecal sample preparation, especially in animal species, we developed a robust protocol allowing untargeted LC-HRMS fingerprinting.

Methods

The conditions of extraction (quantity, preparation, solvents, dilutions) were investigated in bovine feces.

Results

A rapid and simple protocol involving feces extraction with methanol (1/3, M/V) followed by centrifugation and a step filtration (10 kDa) was developed.

Conclusion

The workflow generated repeatable and informative fingerprints for robust metabolome characterization.

     

Choroidal eye metastases from (recurrent) primary peritoneal carcinoma: case report and review of the literature

Background

Choroidal metastases from gynaecological primary are extremely rare. There is no documented case in the literature of choroid metastasis in a patient with primary peritoneal carcinoma (PPC).

Methods & Results

We describe the first case of a 54-year-old woman with a history of borderline mucinous tumour who presented 17 months later with PPC and 21 months after with recurrent disease metastatic to the eye, and review pertinent literature.

Conclusion

High index of suspicion is warranted when patients with history of primary peritoneal carcinoma present with visual complaints in order to treat and/or relieve symptomatology from metastatic eye disease.

     

The repertoire of protein kinases encoded in the draft version of the human genome: atypical variations and uncommon domain combinations

Background

Phosphorylation by protein kinases is central to cellular signal transduction. Abnormal functioning of kinases has been implicated in developmental disorders and malignancies. Their activity is regulated by second messengers and by the binding of associated domains, which are also influential in translocating the catalytic component to their substrate sites, in mediating interaction with other proteins and carrying out their biological roles.

Result

Using sensitive profile-search methods and manual analysis, the human genome has been surveyed for protein kinases. A set of 448 sequences, which show significant similarity to protein kinases and contain the critical residues essential for kinase function, have been selected for an analysis of domain combinations after classifying the kinase domains into subfamilies. The unusual domain combinations in particular kinases suggest their involvement in ubiquitination pathways and alternative modes of regulation for mitogen-activated protein kinase kinases (MAPKKs) and cyclin-dependent kinase (CDK)-like kinases. Previously unexplored kinases have been implicated in osteoblast differentiation and embryonic development on the basis of homology with kinases of known functions from other organisms. Kinases potentially unique to vertebrates are involved in highly evolved processes such as apoptosis, protein translation and tyrosine kinase signaling. In addition to coevolution with the kinase domain, duplication and recruitment of non-catalytic domains is apparent in signaling domains such as the PH, DAG-PE, SH2 and SH3 domains.

Conclusions

Expansion of the functional repertoire and possible existence of alternative modes of regulation of certain kinases is suggested by their uncommon domain combinations. Experimental verification of the predicted implications of these kinases could enhance our understanding of their biological roles.

     

Retrieving hidden atrial repolarization waves from standard surface ECGs

Background

This study estimates atrial repolarization activities (Ta waves), which are typically hidden most of the time from body surface electrocardiography when diagnosing cardiovascular diseases. The morphology of Ta waves has been proven to be an important marker for the early sign of inferior injury, such as acute atrial infarction, or arrhythmia, such as atrial fibrillation. However, Ta waves are usually unseen except during conduction system malfunction, such as long QT interval or atrioventricular block. Therefore, justifying heart diseases based on atrial repolarization becomes impossible in sinus rhythm.

Methods

We obtain TMPs in the atrial part of the myocardium which reflects the correct excitation sequence starting from the atrium to the end of the apex.

Results

The resulting TMP shows the hidden atrial part of ECG waves.

Conclusions

This extraction makes many diseases, such as acute atrial infarction or arrhythmia, become easily diagnosed.

     

Iridoid and phenylethanoid/phenylpropanoid metabolite profiles of <Emphasis Type="Italic">Scrophularia</Emphasis> and <Emphasis Type="Italic">Verbascum</Emphasis> species used medicinally in North America

Introduction

Botanicals containing iridoid and phenylethanoid/phenylpropanoid glycosides are used worldwide for the treatment of inflammatory musculoskeletal conditions that are primary causes of human years lived with disability, such as arthritis and lower back pain.

Objectives

We report the analysis of candidate anti-inflammatory metabolites of several endemic Scrophularia species and Verbascum thapsus used medicinally by peoples of North America.

Methods

Leaves, stems, and roots were analyzed by ultra-performance liquid chromatography-mass spectrometry (UPLC-MS) and partial least squares-discriminant analysis (PLS-DA) was performed in MetaboAnalyst 3.0 after processing the datasets in Progenesis QI.

Results

Comparison of the datasets revealed significant and differential accumulation of iridoid and phenylethanoid/phenylpropanoid glycosides in the tissues of the endemic Scrophularia species and Verbascum thapsus.

Conclusions

Our investigation identified several species of pharmacological interest as good sources for harpagoside and other important anti-inflammatory metabolites.

     

NS-kNN: a modified <Emphasis Type="Italic">k</Emphasis>-nearest neighbors approach for imputing metabolomics data

Introduction

A common problem in metabolomics data analysis is the existence of a substantial number of missing values, which can complicate, bias, or even prevent certain downstream analyses. One of the most widely-used solutions to this problem is imputation of missing values using a k-nearest neighbors (kNN) algorithm to estimate missing metabolite abundances. kNN implicitly assumes that missing values are uniformly distributed at random in the dataset, but this is typically not true in metabolomics, where many values are missing because they are below the limit of detection of the analytical instrumentation.

Objectives

Here, we explore the impact of nonuniformly distributed missing values (missing not at random, or MNAR) on imputation performance. We present a new model for generating synthetic missing data and a new algorithm, No-Skip kNN (NS-kNN), that accounts for MNAR values to provide more accurate imputations.

Methods

We compare the imputation errors of the original kNN algorithm using two distance metrics, NS-kNN, and a recently developed algorithm KNN-TN, when applied to multiple experimental datasets with different types and levels of missing data.

Results

Our results show that NS-kNN typically outperforms kNN when at least 20–30% of missing values in a dataset are MNAR. NS-kNN also has lower imputation errors than KNN-TN on realistic datasets when at least 50% of missing values are MNAR.

Conclusion

Accounting for the nonuniform distribution of missing values in metabolomics data can significantly improve the results of imputation algorithms. The NS-kNN method imputes missing metabolomics data more accurately than existing kNN-based approaches when used on realistic datasets.

     

The radial organization of neuronal primary cilia is acutely disrupted by seizure and ischemic brain injury

Background

Neuronal primary cilia are sensory organelles that are critically involved in the proper growth, development, and function of the central nervous system (CNS). Recent work also suggests that they signal in the context of CNS injury, and that abnormal ciliary signaling may be implicated in neurological diseases.

Methods

We quantified the distribution of neuronal primary cilia alignment throughout the normal adult mouse brain by immunohistochemical staining for the primary cilia marker adenylyl cyclase III (ACIII) and measuring the angles of primary cilia with respect to global and local coordinate planes. We then introduced two different models of acute brain insult—temporal lobe seizure and cerebral ischemia, and re-examined neuronal primary cilia distribution, as well as ciliary lengths and the proportion of neurons harboring cilia.

Results

Under basal conditions, cortical cilia align themselves radially with respect to the cortical surface, while cilia in the dentate gyrus align themselves radially with respect to the granule cell layer. Cilia of neurons in the striatum and thalamus, by contrast, exhibit a wide distribution of ciliary arrangements. In both cases of acute brain insult, primary cilia alignment was significantly disrupted in a region-specific manner, with areas affected by the insult preferentially disrupted. Further, the two models promoted differential effects on ciliary lengths, while only the ischemia model decreased the proportion of ciliated cells.

Conclusions

These findings provide evidence for the regional anatomical organization of neuronal primary cilia in the adult brain and suggest that various brain insults may disrupt this organization.

     

The C-terminal domain of TPX2 is made of alpha-helical tandem repeats

Background

TPX2 (Targeting Protein for Xklp2) is essential for spindle assembly, activation of the mitotic kinase Aurora A and for triggering microtubule nucleation. Homologs of TPX2 in Chordata and plants were previously identified. Currently, proteins of the TPX2 family have little structural information and only small parts are covered by defined protein domains.

Methods

We have used computational sequence analyses and structural predictions of proteins of the TPX2 family, supported with Circular Dichroism (CD) measurements.

Results

Here, we report our finding that the C-terminal domain of TPX2, which is responsible of its microtubule nucleation capacity and is conserved in all members of the family, is actually formed by tandem repeats, covering well above 2/3 of the protein. We propose that this region forms a flexible solenoid involved in protein-protein interactions. Structural prediction and molecular modeling, combined with Circular Dichroism (CD) measurements reveal a predominant alpha-helical content. Furthermore, we identify full length homologs in fungi and shorter homologs with a different domain organization in diptera (including a paralogous expansion in Drosophila).

Conclusions

Our results, represent the first computational and biophysical analysis of the TPX2 proteins family and help understand the structure and evolution of this conserved protein family to direct future structural studies.

     

Incidence of T315I mutation in BCR/ABL-positive CML and ALL patients

Objectives

Targeted therapy of Philadelphia-positive ALL and CML patients using imatinib (IM) has caused significant changes in treatment course and has increased the survival of patients. A small group of patients show resistance to IM. Acquired mutations in tyrosine kinase domain of BCR-ABL protein are a mechanism for development of resistance. T315I is one of the most common acquired mutations in this domain, which occurs in ATP binding site and inhibits the formation of hydrogen bond with IM. The aim of this study was to evaluate the prevalence of this mutation in BCR/ABL-positive CML and ALL patients.

Methods

To conduct this study, 60 BCR-ABL-positive patients (including 50 CML and 10 ALL patients) who were subject to treatment with IM were selected. After taking the samples, presence of T315I mutation was assessed using ARMS-PCR on cDNA and its polymorphism was evaluated by sequencing.

Results

The results showed that among 60 patients, only three patients had T315I mutation, which was detected using ARMS technique. The three patients bearing mutation were afflicted with CML and no significant association was found between blood parameters with duration of treatment in presence of mutation.

Conclusions

The mutation was found in three CML patients, which indicated lower likelihood and diagnostic value of this mutation in ALL patients. Given the negative direct sequencing results in T315I patients, it can be concluded that ARMS-PCR is a more sensitive technique when the number of cancer cells is low in patients during treatment.

     

Crystal structure of <Emphasis Type="Italic">Escherichia coli</Emphasis> protein ybgI,a toroidal structure with a dinuclear metal site

Background

The protein encoded by the gene ybgI was chosen as a target for a structural genomics project emphasizing the relation of protein structure to function.

Results

The structure of the ybgI protein is a toroid composed of six polypeptide chains forming a trimer of dimers. Each polypeptide chain binds two metal ions on the inside of the toroid.

Conclusion

The toroidal structure is comparable to that of some proteins that are involved in DNA metabolism. The di-nuclear metal site could imply that the specific function of this protein is as a hydrolase-oxidase enzyme.

     

Calcium is the switch in the moonlighting dual function of the ligand-activated receptor kinase phytosulfokine receptor 1

Background

A number of receptor kinases contain guanylate cyclase (GC) catalytic centres encapsulated in the cytosolic kinase domain. A prototypical example is the phytosulfokine receptor 1 (PSKR1) that is involved in regulating growth responses in plants. PSKR1 contains both kinase and GC activities however the underlying mechanisms regulating the dual functions have remained elusive.

Findings

Here, we confirm the dual activity of the cytoplasmic domain of the PSKR1 receptor. We show that mutations within the guanylate cyclase centre modulate the GC activity while not affecting the kinase catalytic activity. Using physiologically relevant Ca²⁺ levels, we demonstrate that its GC activity is enhanced over two-fold by Ca²⁺ in a concentration-dependent manner. Conversely, increasing Ca²⁺ levels inhibits kinase activity up to 500-fold at 100 nM Ca²⁺.

Conclusions

Changes in calcium at physiological levels can regulate the kinase and GC activities of PSKR1. We therefore propose a functional model of how calcium acts as a bimodal switch between kinase and GC activity in PSKR1 that could be relevant to other members of this novel class of ligand-activated receptor kinases.

     

NMRProcFlow: a graphical and interactive tool dedicated to 1D spectra processing for NMR-based metabolomics

Introduction

Concerning NMR-based metabolomics, 1D spectra processing often requires an expert eye for disentangling the intertwined peaks.

Objectives

The objective of NMRProcFlow is to assist the expert in this task in the best way without requirement of programming skills.

Methods

NMRProcFlow was developed to be a graphical and interactive 1D NMR (¹H & ¹³C) spectra processing tool.

Results

NMRProcFlow (http://nmrprocflow.org), dedicated to metabolic fingerprinting and targeted metabolomics, covers all spectra processing steps including baseline correction, chemical shift calibration and alignment.

Conclusion

Biologists and NMR spectroscopists can easily interact and develop synergies by visualizing the NMR spectra along with their corresponding experimental-factor levels, thus setting a bridge between experimental design and subsequent statistical analyses.

     

The D Domain of LRRC4 anchors ERK1/2 in the cytoplasm and competitively inhibits MEK/ERK activation in glioma cells

Background

As a well-characterized key player in various signal transduction networks, extracellular-signal-regulated kinase (ERK1/2) has been widely implicated in the development of many malignancies. We previously found that Leucine-rich repeat containing 4 (LRRC4) was a tumor suppressor and a negative regulator of the ERK/MAPK pathway in glioma tumorigenesis. However, the precise molecular role of LRRC4 in ERK signal transmission is unclear.

Methods

The interaction between LRRC4 and ERK1/2 was assessed by co-immunoprecipitation and GST pull-down assays in vivo and in vitro. We also investigated the interaction of LRRC4 and ERK1/2 and the role of the D domain in ERK activation in glioma cells.

Results

Here, we showed that LRRC4 and ERK1/2 interact via the D domain and CD domain, respectively. Following EGF stimuli, the D domain of LRRC4 anchors ERK1/2 in the cytoplasm and abrogates ERK1/2 activation and nuclear translocation. In glioblastoma cells, ectopic LRRC4 expression competitively inhibited the interaction of endogenous mitogen-activated protein kinase (MEK) and ERK1/2. Mutation of the D domain decreased the LRRC4-mediated inhibition of MAPK signaling and its anti-proliferation and anti-invasion roles.

Conclusions

Our results demonstrated that the D domain of LRRC4 anchors ERK1/2 in the cytoplasm and competitively inhibits MEK/ERK activation in glioma cells. These findings identify a new mechanism underlying glioblastoma progression and suggest a novel therapeutic strategy by restoring the activity of LRRC4 to decrease MAPK cascade activation.