GPS 5.0: An Update on the Prediction of Kinase-specific Phosphorylation Sites in Proteins

The global phosphorylation landscape of mouse oocytes during meiotic maturation

Phosphorylation is a key post-translational modification regulating protein function and biological outcomes. However, the phosphorylation dynamics orchestrating mammalian oocyte development remains poorly understood. In the present study, we apply high-resolution mass spectrometry-based phosphoproteomics to obtain the first global in vivo quantification of mouse oocyte phosphorylation. Of more than 8000 phosphosites, 75% significantly oscillate and 64% exhibit marked upregulation during meiotic maturation, indicative of the dominant regulatory role. Moreover, we identify numerous novel phosphosites on oocyte proteins and a few highly conserved phosphosites in oocytes from different species. Through functional perturbations, we demonstrate that phosphorylation status of specific sites participates in modulating critical events including metabolism, translation, and RNA processing during meiosis. Finally, we combine inhibitor screening and enzyme-substrate network prediction to discover previously unexplored kinases and phosphatases that are essential for oocyte maturation. In sum, our data define landscape of the oocyte phosphoproteome, enabling in-depth mechanistic insights into developmental control of germ cells.

Ketogenic diet reshapes cancer metabolism through lysine β-hydroxybutyrylation

Lysine β-hydroxybutyrylation (Kbhb) is a post-translational modification induced by the ketogenic diet (KD), a diet showing therapeutic effects on multiple human diseases. Little is known how cellular processes are regulated by Kbhb. Here we show that protein Kbhb is strongly affected by the KD through a multi-omics analysis of mouse livers. Using a small training dataset with known functions, we developed a bioinformatics method for the prediction of functionally important lysine modification sites (pFunK), which revealed functionally relevant Kbhb sites on various proteins, including aldolase B (ALDOB) Lys108. KD consumption or β-hydroxybutyrate supplementation in hepatocellular carcinoma cells increases ALDOB Lys108bhb and inhibits the enzymatic activity of ALDOB. A Kbhb-mimicking mutation (p.Lys108Gln) attenuates ALDOB activity and its binding to substrate fructose-1,6-bisphosphate, inhibits mammalian target of rapamycin signalling and glycolysis, and markedly suppresses cancer cell proliferation. Our study reveals a critical role of Kbhb in regulating cancer cell metabolism and provides a generally applicable algorithm for predicting functionally important lysine modification sites.

Distinctive expression and cellular localisation of zinc homeostasis-related proteins in breast and prostate cancer cells

BACKGROUND

Zinc transport proteins (ZIP and ZnT), metallothioneins (MT) and protein kinase CK2 are involved in dysregulation of zinc homeostasis in breast and prostate cancer cells. Following up our previous research, we targeted ZIP12, ZnT1, MT2A and CK2 in this study by investigating their expression levels and protein localisation.

METHODS

Quantitative reverse transcription polymerase chain reaction (qRT-PCR) and immunofluorescence confocal microscopy were employed to quantify the expression of ZIP12, ZnT1, MT2A and CK2 subunits in a panel of breast and prostate cell lines without or with extracellular zinc exposure. The cellular localisations of these target proteins were also examined by immunofluorescence confocal microscopy.

RESULTS

In response to the extracellular zinc exposure, the gene expression was elevated for SLC39A12 (ZIP12), SLC30A1 (ZnT1) and MT2A (MT2A) in normal prostate epithelial cells (RWPE-1) in contrast to their cancerous counterparts (PC3 and DU145), whilst the gene expression was higher for SLC39A12 (ZIP12) and SLC30A1 (ZnT1) in both normal (MCF10A) and basal breast cancer cells (MDA-MB-231) compared to luminal breast cancer cells (MCF-7). At the protein level, the expression for both ZIP12 and ZnT1 was trending lower in the time course for the breast cancer cells whilst their expression was remained constant in the normal breast epithelial cells. The expression of ZIP12 in prostate cancer cells was higher than the normal prostate cells. The protein expression for CK2 α/αꞌ and CK2β was markedly higher in prostate cancer cells than the normal prostate cells. Upon extracellular zinc exposure, ZIP12 was, for the first time, conspicuously localised in the plasma membrane of breast cancer cells but not in normal breast epithelial cells and prostate cells. ZnT1 is only localised in the plasma membrane of breast cancer cells. MT2A is distinctively seen close to the plasma membrane in breast cancer cells. CK2 is also for the first time shown to be localised in proximity to the plasma membrane of breast cancer cells.

CONCLUSION

The findings, particularly the localisation of ZIP12 and CK2, are novel and significant for our understanding of zinc homeostasis in breast and prostate cancer cells.

Insights into the role of FoxL2 in tebuconazole-induced male- biased sex differentiation of zebrafish

Tebuconazole (TEB) is a commonly used fungicide that inhibits the aromatase Cyp19A and downregulates the transcription factor forkhead box L2 (FoxL2), leading to male-biased sex differentiation in zebrafish larvae. However, the specific mechanism by which FoxL2 functions following TEB exposure remains unclear. In this study, the phosphorylation sites and kinase-specific residues in zebrafish FoxL2 protein (zFoxL2) were predicted. Subsequently, recombinant zFoxL2 was prepared via prokaryotic expression, and a polyclonal rabbit-anti-zFoxL2 antibody was generated. Zebrafish fibroblast (ZF4) cells were exposed to 100-μM TEB alone for 8 h, after which changes in the expression of genes involved in the foxl2 regulatory pathway (akt1, pi3k, cyp19a1b, c/ebpb and sox9a) were detected. When co-exposed to 1-μM estradiol and 100-μM TEB, the expression of these key genes tended to be restored. Interestingly, TEB did not affect the expression of the foxl2 gene or protein but it significantly suppressed the phosphorylation of FoxL2 (pFoxL2) at serine 238 (decreased by 43.64 %, p = 0.009). Co-immunoprecipitation assays showed that, following exposure to 100-μM TEB, the total precipitated proteins in ZF4 cells decreased by 17.02 % (p = 0.029) and 31.39 % (p = 0.027) in the anti-zFoxL2 antibody group and anti-pFoxL2 (ser238) antibody group, respectively, indicating that TEB suppressed the capacity of the FoxL2 protein to bind to other proteins via repression of its own phosphorylation. The pull-down assay confirmed this conclusion. This study preliminarily elucidated that the foxl2 gene functions via post-translational regulation through hypophosphorylation of its encoded protein during TEB-induced male-biased sex differentiation.

Computational screening of pathogenic missense nsSNPs in heme oxygenase 1 (HMOX1) gene and their structural and functional consequences

Heme Oxygenase 1 (HMOX1) is a cytoprotective enzyme, exhibiting the highest activity in the spleen, catalyzing the heme ring breakdown into products of biological significance- biliverdin, CO, and Fe2+. In vascular cells, HMOX1 possesses strong anti-apoptotic, antioxidant, anti-proliferative, anti-inflammatory, and immunomodulatory actions. The majority of these activities are crucial for the prevention of atherogenesis. Single amino acid substitutions in proteins generated by missense non-synonymous single nucleotide polymorphism (nsSNPs) in the protein-encoding regions of genes are potent enough to cause significant medical challenges due to the alteration of protein structure and function. The current study aimed at characterizing and analyzing high-risk nsSNPs associated with the human HMOX1 gene. Preliminary screening of the total available 288 missense SNPs was performed through the lens of deleteriousness and stability prediction tools. Finally, a total of seven nsSNPs (Y58D, A131T, Y134H, F166S, F167S, R183S and M186V) were found to be most deleterious by all tools that are present at highly conserved positions. Molecular dynamics simulations (MDS) analysis explained the mutational effects on the dynamic action of the wild-type and mutant proteins. In a nutshell, R183S (rs749644285) was identified as a highly detrimental mutation that could significantly render the enzymatic activity of HMOX1. The finding of this computational analysis might help subject the experimental confirmatory analysis to characterize the role of nsSNPs in HMOX1.Communicated by Ramaswamy H. Sarma.

In silico transcriptome screens identify epidermal growth factor receptor inhibitors as therapeutics for noise-induced hearing loss

Noise-induced hearing loss (NIHL) is a common sensorineural hearing impairment that lacks U.S. Food and Drug Administration-approved drugs. To fill the gap in effective screening models, we used an in silico transcriptome-based drug screening approach, identifying 22 biological pathways and 64 potential small molecule treatments for NIHL. Two of these, afatinib and zorifertinib [epidermal growth factor receptor (EGFR) inhibitors], showed efficacy in zebrafish and mouse models. Further tests with EGFR knockout mice and EGF-morpholino zebrafish confirmed their protective role against NIHL. Molecular studies in mice highlighted EGFR's crucial involvement in NIHL and the protective effect of zorifertinib. When given orally, zorifertinib was found in the perilymph with favorable pharmacokinetics. In addition, zorifertinib combined with AZD5438 (a cyclin-dependent kinase 2 inhibitor) synergistically prevented NIHL in zebrafish. Our results underscore the potential for in silico transcriptome-based drug screening in diseases lacking efficient models and suggest EGFR inhibitors as potential treatments for NIHL, meriting clinical trials.

Atypical Protein Kinase C Promotes its own Asymmetric Localisation by Phosphorylating Cdc42 in the C. elegans zygote

Atypical protein kinase C (aPKC) is a major regulator of cell polarity. Acting in conjunction with Par6, Par3 and the small GTPase Cdc42, aPKC becomes asymmetrically localised and drives the polarisation of cells. aPKC activity is crucial for its own asymmetric localisation, suggesting a hitherto unknown feedback mechanism contributing to polarisation. Here we show in the C. elegans zygote that the feedback relies on aPKC phosphorylation of Cdc42 at serine 71. The turnover of CDC-42 phosphorylation ensures optimal aPKC asymmetry and activity throughout polarisation by tuning Par6/aPKC association with Par3 and Cdc42. Moreover, turnover of Cdc42 phosphorylation regulates actomyosin cortex dynamics that are known to drive aPKC asymmetry. Given the widespread role of aPKC and Cdc42 in cell polarity, this form of self-regulation of aPKC may be vital for the robust control of polarisation in many cell types.

Lysophosphatidic Acid Receptor 3 (LPA3): Signaling and Phosphorylation Sites

LPA3 receptors were expressed in TREx HEK 293 cells, and their signaling and phosphorylation were studied. The agonist, lysophosphatidic acid (LPA), increased intracellular calcium and ERK phosphorylation through pertussis toxin-insensitive processes. Phorbol myristate acetate, but not LPA, desensitizes LPA3-mediated calcium signaling, the agonists, and the phorbol ester-induced LPA3 internalization. Pitstop 2 (clathrin heavy chain inhibitor) markedly reduced LPA-induced receptor internalization; in contrast, phorbol ester-induced internalization was only delayed. LPA induced rapid β-arrestin-LPA3 receptor association. The agonist and the phorbol ester-induced marked LPA3 receptor phosphorylation, and phosphorylation sites were detected using mass spectrometry. Phosphorylated residues were detected in the intracellular loop 3 (S221, T224, S225, and S229) and in the carboxyl terminus (S321, S325, S331, T333, S335, Y337, and S343). Interestingly, phosphorylation sites are within sequences predicted to constitute β-arrestin binding sites. These data provide insight into LPA3 receptor signaling and regulation.

Understanding Genetic Risks: Computational Exploration of Human β-Synuclein nsSNPs and their Potential Impact on Structural Alteration

Synucleins are pivotal in neurodegenerative conditions. Beta-synuclein (β-synuclein) is part of the synuclein protein family alongside alpha-synuclein (α-synuclein) and gamma-synuclein (γ-synuclein). These proteins, found mainly in brain tissue and cancers, are soluble and unstructured. β-synuclein shares significant similarity with α-synuclein, especially in their N-terminus, with a 90% match. However, their aggregation tendencies differ significantly. While α-synuclein aggregation is believed to be counteracted by β-synuclein, which occurs in conditions like Parkinson's disease, β-synuclein may counteract α-synuclein's toxic effects on the nervous system, offering potential treatment for neurodegenerative diseases. Under normal circumstances, β-synuclein may guard against disease by interacting with α-synuclein. Yet, in pathological environments with heightened levels or toxic substances, it might contribute to disease. Our research aims to explore potential harmful mutations in the β-synuclein using computational tools to predict their destabilizing impact on protein structure. Consensus analysis revealed rs1207608813 (A63P), rs1340051870 (S72F), and rs1581178262 (G36C) as deleterious. These findings highlight the intricate relationship between nsSNPs and protein function, shedding light on their potential implications in disease pathways. Understanding the structural consequences of nsSNPs is crucial for elucidating their role in pathogenesis and developing targeted therapeutic interventions. Our results offer a robust computational framework for identifying neurodegenerative disorder-related mutations from SNP datasets, potentially reducing the costs associated with experimental characterization.

S-PLM: Structure-aware Protein Language Model via Contrastive Learning between Sequence and Structure

Proteins play an essential role in various biological and engineering processes. Large protein language models (PLMs) present excellent potential to reshape protein research by accelerating the determination of protein function and the design of proteins with the desired functions. The prediction and design capacity of PLMs relies on the representation gained from the protein sequences. However, the lack of crucial 3D structure information in most PLMs restricts the prediction capacity of PLMs in various applications, especially those heavily dependent on 3D structures. To address this issue, we introduce S-PLM, a 3D structure-aware PLM that utilizes multi-view contrastive learning to align the sequence and 3D structure of a protein in a coordinated latent space. S-PLM applies Swin-Transformer on AlphaFold-predicted protein structures to embed the structural information and fuses it into sequence-based embedding from ESM2. Additionally, we provide a library of lightweight tuning tools to adapt S-PLM for diverse protein property prediction tasks. Our results demonstrate S-PLM's superior performance over sequence-only PLMs on all protein clustering and classification tasks, achieving competitiveness comparable to state-of-the-art methods requiring both sequence and structure inputs. S-PLM and its lightweight tuning tools are available at https://github.com/duolinwang/S-PLM/ .

Alternative splicing of clock transcript mediates the response of circadian clocks to temperature changes

Circadian clocks respond to temperature changes over the calendar year, allowing organisms to adjust their daily biological rhythms to optimize health and fitness. In Drosophila, seasonal adaptations and temperature compensation are regulated by temperature-sensitive alternative splicing (AS) of period (per) and timeless (tim) genes that encode key transcriptional repressors of clock gene expression. Although clock (clk) gene encodes the critical activator of clock gene expression, AS of its transcripts and its potential role in temperature regulation of clock function have not been explored. We therefore sought to investigate whether clk exhibits AS in response to temperature and the functional changes of the differentially spliced transcripts. We observed that clk transcripts indeed undergo temperature-sensitive AS. Specifically, cold temperature leads to the production of an alternative clk transcript, hereinafter termed clk-cold, which encodes a CLK isoform with an in-frame deletion of four amino acids proximal to the DNA binding domain. Notably, serine 13 (S13), which we found to be a CK1α-dependent phosphorylation site, is among the four amino acids deleted in CLK-cold protein. Using a combination of transgenic fly, tissue culture, and in vitro experiments, we demonstrated that upon phosphorylation at CLK(S13), CLK-DNA interaction is reduced, thus decreasing CLK occupancy at clock gene promoters. This is in agreement with our findings that CLK occupancy at clock genes and transcriptional output are elevated at cold temperature, which can be explained by the higher amounts of CLK-cold isoforms that lack S13 residue. This study provides new insights into the complex collaboration between AS and phospho-regulation in shaping temperature responses of the circadian clock.

KinPred-RNA-kinase activity inference and cancer type classification using machine learning on RNA-seq data

Kinases as important enzymes can transfer phosphate groups from high-energy and phosphate-donating molecules to specific substrates and play essential roles in various cellular processes. Existing algorithms for kinase activity from phosphorylated proteomics data are often costly, requiring valuable samples. Moreover, methods to extract kinase activities from bulk RNA sequencing data remain undeveloped. In this study, we propose a computational framework KinPred-RNA to derive kinase activities from bulk RNA-sequencing data in cancer samples. KinPred-RNA framework, using the extreme gradient boosting (XGBoost) regression model, outperforms random forest regression, multiple linear regression, and support vector machine regression models in predicting kinase activities from cancer-related RNA sequencing data. Efficient gene signatures from the LINCS-L1000 dataset were used as inputs for KinPred-RNA. The results highlight its potential to be related to biological function. In conclusion, KinPred RNA constitutes a significant advance in cancer research by potentially facilitating the identification of cancer.

Unveiling orphan receptor-like kinases in plants: novel client discovery using high-confidence library predictions in the Kinase-Client (KiC) assay

Plants are remarkable in their ability to adapt to changing environments, with receptor-like kinases (RLKs) playing a pivotal role in perceiving and transmitting environmental cues into cellular responses. Despite extensive research on RLKs from the plant kingdom, the function and activity of many kinases, i.e., their substrates or "clients", remain uncharted. To validate a novel client prediction workflow and learn more about an important RLK, this study focuses on P2K1 (DORN1), which acts as a receptor for extracellular ATP (eATP), playing a crucial role in plant stress resistance and immunity. We designed a Kinase-Client (KiC) assay library of 225 synthetic peptides, incorporating previously identified P2K phosphorylated peptides and novel predictions from a deep-learning phosphorylation site prediction model (MUsite) and a trained hidden Markov model (HMM) based tool, HMMER. Screening the library against purified P2K1 cytosolic domain (CD), we identified 46 putative substrates, including 34 novel clients, 27 of which may be novel peptides, not previously identified experimentally. Gene Ontology (GO) analysis among phosphopeptide candidates revealed proteins associated with important biological processes in metabolism, structure development, and response to stress, as well as molecular functions of kinase activity, catalytic activity, and transferase activity. We offer selection criteria for efficient further in vivo experiments to confirm these discoveries. This approach not only expands our knowledge of P2K1's substrates and functions but also highlights effective prediction algorithms for identifying additional potential substrates. Overall, the results support use of the KiC assay as a valuable tool in unraveling the complexities of plant phosphorylation and provide a foundation for predicting the phosphorylation landscape of plant species based on peptide library results.

Control of cell proliferation by memories of mitosis

Mitotic duration is tightly constrained, and extended mitosis is characteristic of problematic cells prone to chromosome missegregation and genomic instability. We show here that mitotic extension leads to the formation of p53-binding protein 1 (53BP1)-ubiquitin-specific protease 28 (USP28)-p53 protein complexes that are transmitted to, and stably retained by, daughter cells. Complexes assembled through a Polo-like kinase 1-dependent mechanism during extended mitosis and elicited a p53 response in G1 that prevented the proliferation of the progeny of cells that experienced an approximately threefold extended mitosis or successive less extended mitoses. The ability to monitor mitotic extension was lost in p53-mutant cancers and some p53-wild-type (p53-WT) cancers, consistent with classification of TP53BP1 and USP28 as tumor suppressors. Cancers retaining the ability to monitor mitotic extension exhibited sensitivity to antimitotic agents.

Precise prediction of phase-separation key residues by machine learning

Understanding intracellular phase separation is crucial for deciphering transcriptional control, cell fate transitions, and disease mechanisms. However, the key residues, which impact phase separation the most for protein phase separation function have remained elusive. We develop PSPHunter, which can precisely predict these key residues based on machine learning scheme. In vivo and in vitro validations demonstrate that truncating just 6 key residues in GATA3 disrupts phase separation, enhancing tumor cell migration and inhibiting growth. Glycine and its motifs are enriched in spacer and key residues, as revealed by our comprehensive analysis. PSPHunter identifies nearly 80% of disease-associated phase-separating proteins, with frequent mutated pathological residues like glycine and proline often residing in these key residues. PSPHunter thus emerges as a crucial tool to uncover key residues, facilitating insights into phase separation mechanisms governing transcriptional control, cell fate transitions, and disease development.

DF-Phos: Prediction of Protein Phosphorylation Sites by Deep Forest

Phosphorylation is the most important and studied post-translational modification (PTM), which plays a crucial role in protein function studies and experimental design. Many significant studies have been performed to predict phosphorylation sites using various machine-learning methods. Recently, several studies have claimed that deep learning-based methods are the best way to predict the phosphorylation sites because deep learning as an advanced machine learning method can automatically detect complex representations of phosphorylation patterns from raw sequences and thus offers a powerful tool to improve phosphorylation site prediction. In this study, we report DF-Phos, a new phosphosite predictor based on the Deep Forest to predict phosphorylation sites. In DF-Phos, the feature vector taken from the CkSAApair method is as input for a Deep Forest framework for predicting phosphorylation sites. The results of 10-fold cross-validation show that the Deep Forest method has the highest performance among other available methods. We implemented a Python program of DF-Phos, which is freely available for non-commercial use at https://github.com/zahiriz/DF-Phos Moreover, users can use it for various PTM predictions.

Plasma p-tau212 antemortem diagnostic performance and prediction of autopsy verification of Alzheimer's disease neuropathology

Blood phosphorylated tau (p-tau) biomarkers, including p-tau217, show high associations with Alzheimer's disease (AD) neuropathologic change and clinical stage. Certain plasma p-tau217 assays recognize tau forms phosphorylated additionally at threonine-212, but the contribution of p-tau212 alone to AD is unknown. We developed a blood-based immunoassay that is specific to p-tau212 without cross-reactivity to p-tau217. Here, we examined the diagnostic utility of plasma p-tau212. In five cohorts (n = 388 participants), plasma p-tau212 showed high performances for AD diagnosis and for the detection of both amyloid and tau pathology, including at autopsy as well as in memory clinic populations. The diagnostic accuracy and fold changes of plasma p-tau212 were similar to those for p-tau217 but higher than p-tau181 and p-tau231. Immunofluorescent staining of brain tissue slices showed prominent p-tau212 reactivity in neurofibrillary tangles that co-localized with p-tau217 and p-tau202/205. These findings support plasma p-tau212 as a peripherally accessible biomarker of AD pathophysiology.

Illuminating the Dark Cancer Phosphoproteome Through a Machine-Learned Co-Regulation Map of 26,280 Phosphosites

Mass spectrometry-based phosphoproteomics offers a comprehensive view of protein phosphorylation, but limited knowledge about the regulation and function of most phosphosites restricts our ability to extract meaningful biological insights from phosphoproteomics data. To address this, we combine machine learning and phosphoproteomic data from 1,195 tumor specimens spanning 11 cancer types to construct CoPheeMap, a network mapping the co-regulation of 26,280 phosphosites. Integrating network features from CoPheeMap into a machine learning model, CoPheeKSA, we achieve superior performance in predicting kinase-substrate associations. CoPheeKSA reveals 24,015 associations between 9,399 phosphosites and 104 serine/threonine kinases, including many unannotated phosphosites and under-studied kinases. We validate the accuracy of these predictions using experimentally determined kinase-substrate specificities. By applying CoPheeMap and CoPheeKSA to phosphosites with high computationally predicted functional significance and cancer-associated phosphosites, we demonstrate the effectiveness of these tools in systematically illuminating phosphosites of interest, revealing dysregulated signaling processes in human cancer, and identifying under-studied kinases as putative therapeutic targets.

Endothelial β-catenin upregulation and Y142 phosphorylation drive diabetic angiogenesis via upregulating KDR/HDAC9

BACKGROUND

Diabetic angiogenesis is closely associated with disabilities and death caused by diabetic microvascular complications. Advanced glycation end products (AGEs) are abnormally accumulated in diabetic patients and are a key pathogenic factor for diabetic angiogenesis. The present study focuses on understanding the mechanisms underlying diabetic angiogenesis and identifying therapeutic targets based on these mechanisms.

METHODS

In this study, AGE-induced angiogenesis serves as a model to investigate the mechanisms underlying diabetic angiogensis. Mouse aortic rings, matrigel plugs, and HUVECs or 293T cells were employed as research objects to explore this pathological process by using transcriptomics, gene promoter reporter assays, virtual screening and so on.

RESULTS

Here, we found that AGEs activated Wnt/β-catenin signaling pathway and enhanced the β-catenin protein level by affecting the expression of β-catenin degradation-related genes, such as FZDs (Frizzled receptors), LRPs (LDL Receptor Related Proteins), and AXIN1. AGEs could also mediate β-catenin Y142 phosphorylation through VEGFR1 isoform5. These dual effects of AGEs elevated the nuclear translocation of β-catenin and sequentially induced the expression of KDR (Kinase Insert Domain Receptor) and HDAC9 (Histone Deacetylase 9) by POU5F1 and NANOG, respectively, thus mediating angiogenesis. Finally, through virtual screening, Bioymifi, an inhibitor that blocks VEGFR1 isoform5-β-catenin complex interaction and alleviates AGE-induced angiogenesis, was identified.

CONCLUSION

Collectively, this study offers insight into the pathophysiological functions of β-catenin in diabetic angiogenesis.

Trioxidized cysteine in the aging proteome mimics the structural dynamics and interactome of phosphorylated serine

Aging is the primary risk factor for the development of numerous human chronic diseases. On a molecular level, it significantly impacts the regulation of protein modifications, leading to the accumulation of degenerative protein modifications (DPMs) such as aberrant serine phosphorylation (p-Ser) and trioxidized cysteine (t-Cys) within the proteome. The altered p-Ser is linked to abnormal cell signaling, while the accumulation of t-Cys is associated with chronic diseases induced by oxidative stress. Despite this, the potential cross-effects and functional interplay between these two critical molecular factors of aging remain undisclosed. This study analyzes the aging proteome of wild-type C57BL/6NTac mice over 2 years using advanced proteomics and bioinformatics. Our objective is to provide a comprehensive analysis of how t-Cys affects cell signaling and protein structure in the aging process. The results obtained indicate that t-Cys residues accumulate in the aging proteome, interact with p-Ser interacting enzymes, as validated in vitro, and alter their structures similarly to p-Ser. These findings have significant implications for understanding the interplay of oxidative stress and phosphorylation in the aging process. Additionally, they open new venues for further research on the role(s) of these protein modifications in various human chronic diseases and aging, wherein exacerbated oxidation and aberrant phosphorylation are implicated.

The role of RNF138 in DNA end resection is regulated by ubiquitylation and CDK phosphorylation

Double-strand breaks (DSBs) are DNA lesions that pose a significant threat to genomic stability. The repair of DSBs by the homologous recombination (HR) pathway is preceded by DNA end resection, the 5' to 3' nucleolytic degradation of DNA away from the DSB. We and others previously identified a role for RNF138, a really interesting new gene finger E3 ubiquitin ligase, in stimulating DNA end resection and HR. Yet, little is known about how RNF138's function is regulated in the context of DSB repair. Here, we show that RNF138 is phosphorylated at residue T27 by cyclin-dependent kinase (CDK) activity during the S and G2 phases of the cell cycle. We also observe that RNF138 is ubiquitylated constitutively, with ubiquitylation occurring in part on residue K158 and rising during the S/G2 phases. Interestingly, RNF138 ubiquitylation decreases upon genotoxic stress. By mutating RNF138 at residues T27, K158, and the previously identified S124 ataxia telangiectasia mutated phosphorylation site (Han et al., 2016, ref. 22), we find that post-translational modifications at all three positions mediate DSB repair. Cells expressing the T27A, K158R, and S124A variants of RNF138 are impaired in DNA end resection, HR activity, and are more sensitive to ionizing radiation compared to those expressing wildtype RNF138. Our findings shed more light on how RNF138 activity is controlled by the cell during HR.

The Quantitative Biotinylproteomics Studies Reveal a WInd-Related Kinase 1 (Raf-Like Kinase 36) Functioning as an Early Signaling Component in Wind-Induced Thigmomorphogenesis and Gravitropism

Wind is one of the most prevalent environmental forces entraining plants to develop various mechano-responses, collectively called thigmomorphogenesis. Largely unknown is how plants transduce these versatile wind force signals downstream to nuclear events and to the development of thigmomorphogenic phenotype or anemotropic response. To identify molecular components at the early steps of the wind force signaling, two mechanical signaling-related phosphoproteins, identified from our previous phosphoproteomic study of Arabidopsis touch response, mitogen-activated protein kinase kinase 1 (MKK1) and 2 (MKK2), were selected for performing in planta TurboID (ID)-based quantitative proximity-labeling (PL) proteomics. This quantitative biotinylproteomics was separately performed on MKK1-ID and MKK2-ID transgenic plants, respectively, using the genetically engineered TurboID biotin ligase expression transgenics as a universal control. This unique PTM proteomics successfully identified 11 and 71 MKK1 and MKK2 putative interactors, respectively. Biotin occupancy ratio (BOR) was found to be an alternative parameter to measure the extent of proximity and specificity between the proximal target proteins and the bait fusion protein. Bioinformatics analysis of these biotinylprotein data also found that TurboID biotin ligase favorably labels the loop region of target proteins. A WInd-Related Kinase 1 (WIRK1), previously known as rapidly accelerated fibrosarcoma (Raf)-like kinase 36 (RAF36), was found to be a putative common interactor for both MKK1 and MKK2 and preferentially interacts with MKK2. Further molecular biology studies of the Arabidopsis RAF36 kinase found that it plays a role in wind regulation of the touch-responsive TCH3 and CML38 gene expression and the phosphorylation of a touch-regulated PATL3 phosphoprotein. Measurement of leaf morphology and shoot gravitropic response of wirk1 (raf36) mutant revealed that the WIRK1 gene is involved in both wind-triggered rosette thigmomorphogenesis and gravitropism of Arabidopsis stems, suggesting that the WIRK1 (RAF36) protein probably functioning upstream of both MKK1 and MKK2 and that it may serve as the crosstalk point among multiple mechano-signal transduction pathways mediating both wind mechano-response and gravitropism.

Hypomagnesemia with Secondary Hypoparathyroidism and Hypocalcemia due to Novel Variants in the Transient Receptor Potential Cation Channel Subfamily M Member 6 ( TRPM6 ) Gene

HOMG1 (hypomagnesemia 1, intestinal) or hypomagnesemia with secondary hypocalcemia is a rare autosomal recessive disorder of magnesium metabolism, characterized by impaired magnesium absorption. This disorder may mimic other conditions presenting with neonatal seizures. Here, we report an infant diagnosed to have hypomagnesemia with secondary hypocalcemia due to novel variants in TRPM6 gene.

Association between missense variants of uncertain significance in the CHEK2 gene and hereditary breast cancer: a cosegregation and bioinformatics analysis

Inherited mutations in the CHEK2 gene have been associated with an increased lifetime risk of developing breast cancer (BC). We aim to identify in the study population the prevalence of mutations in the CHEK2 gene in diagnosed BC patients, evaluate the phenotypic characteristics of the tumor and family history, and predict the deleteriousness of the variants of uncertain significance (VUS). A genetic study was performed, from May 2016 to April 2020, in 396 patients diagnosed with BC at the University Hospital Lozano Blesa of Zaragoza, Spain. Patients with a genetic variant in the CHEK2 gene were selected for the study. We performed a descriptive analysis of the clinical variables, a bibliographic review of the variants, and a cosegregation study when possible. Moreover, an in-depth bioinformatics analysis of CHEK2 VUS was carried out. We identified nine genetic variants in the CHEK2 gene in 10 patients (two pathogenic variants and seven VUS). This supposes a prevalence of 0.75% and 1.77%, respectively. In all cases, there was a family history of BC in first- and/or second-degree relatives. We carried out a cosegregation study in two families, being positive in one of them. The bioinformatics analyses predicted the pathogenicity of six of the VUS. In conclusion, CHEK2 mutations have been associated with an increased risk for BC. This risk is well-established for foundation variants. However, the risk assessment for other variants is unclear. The incorporation of bioinformatics analysis provided supporting evidence of the pathogenicity of VUS.

PP2A Affects Angiogenesis via Its Interaction with a Novel Phosphorylation Site of TSP1

Alterations in angiogenic properties play a pivotal role in the manifestation and onset of various pathologies, including vascular diseases and cancer. Thrombospondin-1 (TSP1) protein is one of the master regulators of angiogenesis. This study unveils a novel aspect of TSP1 regulation through reversible phosphorylation. The silencing of the B55α regulatory subunit of protein phosphatase 2A (PP2A) in endothelial cells led to a significant decrease in TSP1 expression. Direct interaction between TSP1 and PP2A-B55α was confirmed via various methods. Truncated TSP1 constructs were employed to identify the phosphorylation site and the responsible kinase, ultimately pinpointing PKC as the enzyme phosphorylating TSP1 on Ser93. The biological effects of B55α-TSP1 interaction were also analyzed. B55α silencing not only counteracted the increase in TSP1 expression during wound closure but also prolonged wound closure time. Although B55α silenced cells initiated tube-like structures earlier than control cells, their spheroid formation was disrupted, leading to disintegration. Cells transfected with phosphomimic TSP1 S93D exhibited smaller spheroids and reduced effectiveness in tube formation, revealing insights into the effects of TSP1 phosphorylation on angiogenic properties. In this paper, we introduce a new regulatory mechanism of angiogenesis by reversible phosphorylation on TSP1 S93 by PKC and PP2A B55α.

Deep Learning Prediction Boosts Phosphoproteomics-Based Discoveries Through Improved Phosphopeptide Identification

Shotgun phosphoproteomics enables high-throughput analysis of phosphopeptides in biological samples. One of the primary challenges associated with this technology is the relatively low rate of phosphopeptide identification during data analysis. This limitation hampers the full realization of the potential offered by shotgun phosphoproteomics. Here we present DeepRescore2, a computational workflow that leverages deep learning-based retention time and fragment ion intensity predictions to improve phosphopeptide identification and phosphosite localization. Using a state-of-the-art computational workflow as a benchmark, DeepRescore2 increases the number of correctly identified peptide-spectrum matches by 17% in a synthetic dataset and identifies 19% to 46% more phosphopeptides in biological datasets. In a liver cancer dataset, 30% of the significantly altered phosphosites between tumor and normal tissues and 60% of the prognosis-associated phosphosites identified from DeepRescore2-processed data could not be identified based on the state-of-the-art workflow. Notably, DeepRescore2-processed data uniquely identifies EGFR hyperactivation as a new target in poor-prognosis liver cancer, which is validated experimentally. Integration of deep learning prediction in DeepRescore2 improves phosphopeptide identification and facilitates biological discoveries.

Molecular profiling of sponge deflation reveals an ancient relaxant-inflammatory response

A hallmark of animals is the coordination of whole-body movement. Neurons and muscles are central to this, yet coordinated movements also exist in sponges that lack these cell types. Sponges are sessile animals with a complex canal system for filter-feeding. They undergo whole-body movements resembling "contractions" that lead to canal closure and water expulsion. Here, we combine live 3D optical coherence microscopy, pharmacology, and functional proteomics to elucidate the sequence and detail of shape changes, the tissues and molecular physiology involved, and the control of these movements. Morphometric analysis and targeted perturbation suggest that the movement is driven by the relaxation of actomyosin stress fibers in epithelial canal cells, which leads to whole-body deflation via collapse of the incurrent and expansion of the excurrent canal system. Thermal proteome profiling and quantitative phosphoproteomics confirm the control of cellular relaxation by an Akt/NO/PKG/PKA pathway. Agitation-induced deflation leads to differential phosphorylation of proteins forming epithelial cell junctions, implying their mechanosensitive role. Unexpectedly, untargeted metabolomics detect a concomitant decrease in antioxidant molecules during deflation, reflecting an increase in reactive oxygen species. Together with the secretion of proteinases, cytokines, and granulin, this indicates an inflammation-like state of the deflating sponge reminiscent of vascular endothelial cells experiencing oscillatory shear stress. These results suggest the conservation of an ancient relaxant-inflammatory response of perturbed fluid-carrying systems in animals and offer a possible mechanism for whole-body coordination through diffusible paracrine signals and mechanotransduction.

A computational study on structural and functional consequences of nsSNPs in human dopa decarboxylase

The Dopa Decarboxylase (DDC) gene plays an important role in the synthesis of biogenic amines such as dopamine, serotonin, and histamine. Non-synonymous single nucleotide polymorphisms (nsSNPs) in the DDC gene have been linked with various neurodegenerative disorders. In this study, a comprehensive in silico analysis of nsSNPs in the DDC gene was conducted to assess their potential functional consequences and associations with disease outcomes. Using publicly available databases, a complete list of nsSNPs in the DDC gene was obtained. 29 computational tools and algorithms were used to characterise the effects of these nsSNPs on protein structure, function, and stability. In addition, the population-based association studies were performed to investigate possible associations between specific nsSNPs and arthritis. Our research identified four novel DDC gene nsSNPs that have a major impact on the structure and function of proteins. Through molecular dynamics simulations (MDS), we observed changes in the stability of the DDC protein induced by specific nsSNPs. Furthermore, population-based association studies have revealed potential associations between certain DDC nsSNPs and various neurological disorders, including Parkinson's disease and dementia. The in silico approach used in this study offers insightful information about the functional effects of nsSNPs in the DDC gene. These discoveries provide insight into the cellular processes that underlie cognitive disorders. Furthermore, the detection of disease-associated nsSNPs in the DDC gene may facilitate the development of tailored and targeted therapy approaches.Communicated by Ramaswamy H. Sarma.

Troponin I Tyrosine Phosphorylation Beneficially Accelerates Diastolic Function

BACKGROUND

A healthy heart is able to modify its function and increase relaxation through post-translational modifications of myofilament proteins. While there are known examples of serine/threonine kinases directly phosphorylating myofilament proteins to modify heart function, the roles of tyrosine (Y) phosphorylation to directly modify heart function have not been demonstrated. The myofilament protein TnI (troponin I) is the inhibitory subunit of the troponin complex and is a key regulator of cardiac contraction and relaxation. We previously demonstrated that TnI-Y26 phosphorylation decreases calcium-sensitive force development and accelerates calcium dissociation, suggesting a novel role for tyrosine kinase-mediated TnI-Y26 phosphorylation to regulate cardiac relaxation. Therefore, we hypothesize that increasing TnI-Y26 phosphorylation will increase cardiac relaxation in vivo and be beneficial during pathological diastolic dysfunction.

METHODS

The signaling pathway involved in TnI-Y26 phosphorylation was predicted in silico and validated by tyrosine kinase activation and inhibition in primary adult murine cardiomyocytes. To investigate how TnI-Y26 phosphorylation affects cardiac muscle, structure, and function in vivo, we developed a novel TnI-Y26 phosphorylation-mimetic mouse that was subjected to echocardiography, pressure-volume loop hemodynamics, and myofibril mechanical studies. TnI-Y26 phosphorylation-mimetic mice were further subjected to the nephrectomy/DOCA (deoxycorticosterone acetate) model of diastolic dysfunction to investigate the effects of increased TnI-Y26 phosphorylation in disease.

RESULTS

Src tyrosine kinase is sufficient to phosphorylate TnI-Y26 in cardiomyocytes. TnI-Y26 phosphorylation accelerates in vivo relaxation without detrimental structural or systolic impairment. In a mouse model of diastolic dysfunction, TnI-Y26 phosphorylation is beneficial and protects against the development of disease.

CONCLUSIONS

We have demonstrated that tyrosine kinase phosphorylation of TnI is a novel mechanism to directly and beneficially accelerate myocardial relaxation in vivo.

Proteomics-based mass spectrometry profiling of SARS-CoV-2 infection from human nasopharyngeal samples

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the on-going global pandemic of coronavirus disease 2019 (COVID-19) that continues to pose a significant threat to public health worldwide. SARS-CoV-2 encodes four structural proteins namely membrane, nucleocapsid, spike, and envelope proteins that play essential roles in viral entry, fusion, and attachment to the host cell. Extensively glycosylated spike protein efficiently binds to the host angiotensin-converting enzyme 2 initiating viral entry and pathogenesis. Reverse transcriptase polymerase chain reaction on nasopharyngeal swab is the preferred method of sample collection and viral detection because it is a rapid, specific, and high-throughput technique. Alternate strategies such as proteomics and glycoproteomics-based mass spectrometry enable a more detailed and holistic view of the viral proteins and host-pathogen interactions and help in detection of potential disease markers. In this review, we highlight the use of mass spectrometry methods to profile the SARS-CoV-2 proteome from clinical nasopharyngeal swab samples. We also highlight the necessity for a comprehensive glycoproteomics mapping of SARS-CoV-2 from biological complex matrices to identify potential COVID-19 markers.

HEI10 is subject to phase separation and mediates RPA1a degradation during meiotic interference-sensitive crossover formation

Reciprocal exchanges of DNA between homologous chromosomes during meiosis, or crossovers (COs), shuffle genetic information in gametes and progeny. In many eukaryotes, the majority of COs (class I COs) are sensitive to a phenomenon called interference, which influences the occurrence of closely spaced double COs. Class I COs depend on a group of factors called ZMM (Zip, Msh, Mer) proteins including HEI10 (Human Enhancer of Invasion-10). However, how these proteins are recruited to class I CO sites is unclear. Here, we show that HEI10 forms foci on chromatin via a liquid-liquid phase separation (LLPS) mechanism that relies on residue Ser70. A HEI10S70F allele results in LLPS failure and a defect in class I CO formation. We further used immunoprecipitation-mass spectrometry to identify RPA1a (Replication Protein A 1) as a HEI10 interacting protein. Surprisingly, we find that RPA1a also undergoes phase separation and its ubiquitination and degradation are directly regulated by HEI10. We also show that HEI10 is required for the condensation of other class I CO factors. Thus, our results provide mechanistic insight into how meiotic class I CO formation is controlled by HEI10 coupling LLPS and ubiquitination.

Plasma p-tau212: antemortem diagnostic performance and prediction of autopsy verification of Alzheimer's disease neuropathology

Blood phosphorylated tau (p-tau) biomarkers, including p-tau217, show high associations with Alzheimer's disease (AD) neuropathologic change and clinical stage. Certain plasma p-tau217 assays recognize tau forms phosphorylated additionally at threonine-212, but the contribution of p-tau212 alone to AD is unknown. We developed a blood-based immunoassay that is specific to p-tau212 without cross-reactivity to p-tau217. Thereafter, we examined the diagnostic utility of plasma p-tau212. In five cohorts (n=388 participants), plasma p-tau212 showed high performances for AD diagnosis and for the detection of both amyloid and tau pathology, including at autopsy as well as in memory clinic populations. The diagnostic accuracy and fold changes of plasma p-tau212 were similar to those for p-tau217 but higher than p-tau181 and p-tau231. Immunofluorescent staining of brain tissue slices showed prominent p-tau212 reactivity in neurofibrillary tangles that co-localized with p-tau217 and p-tau202/205. These findings support plasma p-tau212 as a novel peripherally accessible biomarker of AD pathophysiology.

EGFR promotes ALKBH5 nuclear retention to attenuate N6-methyladenosine and protect against ferroptosis in glioblastoma

Growth factor receptors rank among the most important oncogenic pathways, but pharmacologic inhibitors often demonstrate limited benefit as monotherapy. Here, we show that epidermal growth factor receptor (EGFR) signaling repressed N6-methyladenosine (m6A) levels in glioblastoma stem cells (GSCs), whereas genetic or pharmacologic EGFR targeting elevated m6A levels. Activated EGFR induced non-receptor tyrosine kinase SRC to phosphorylate the m6A demethylase, AlkB homolog 5 (ALKBH5), thereby inhibiting chromosomal maintenance 1 (CRM1)-mediated nuclear export of ALKBH5 to permit sustained mRNA m6A demethylation in the nucleus. ALKBH5 critically regulated ferroptosis through m6A modulation and YTH N6-methyladenosine RNA binding protein (YTHDF2)-mediated decay of the glutamate-cysteine ligase modifier subunit (GCLM). Pharmacologic targeting of ALKBH5 augmented the anti-tumor efficacy of EGFR and GCLM inhibitors, supporting an EGFR-ALKBH5-GCLM oncogenic axis. Collectively, EGFR reprograms the epitranscriptomic landscape through nuclear retention of the ALKBH5 demethylase to protect against ferroptosis, offering therapeutic paradigms for the treatment of lethal cancers.

SnapKin: a snapshot deep learning ensemble for kinase-substrate prediction from phosphoproteomics data

A major challenge in mass spectrometry-based phosphoproteomics lies in identifying the substrates of kinases, as currently only a small fraction of substrates identified can be confidently linked with a known kinase. Machine learning techniques are promising approaches for leveraging large-scale phosphoproteomics data to computationally predict substrates of kinases. However, the small number of experimentally validated kinase substrates (true positive) and the high data noise in many phosphoproteomics datasets together limit their applicability and utility. Here, we aim to develop advanced kinase-substrate prediction methods to address these challenges. Using a collection of seven large phosphoproteomics datasets, and both traditional and deep learning models, we first demonstrate that a 'pseudo-positive' learning strategy for alleviating small sample size is effective at improving model predictive performance. We next show that a data resampling-based ensemble learning strategy is useful for improving model stability while further enhancing prediction. Lastly, we introduce an ensemble deep learning model ('SnapKin') by incorporating the above two learning strategies into a 'snapshot' ensemble learning algorithm. We propose SnapKin, an ensemble deep learning method, for predicting substrates of kinases from large-scale phosphoproteomics data. We demonstrate that SnapKin consistently outperforms existing methods in kinase-substrate prediction. SnapKin is freely available at https://github.com/PYangLab/SnapKin.

Comprehensive analysis of non-synonymous missense SNPs of human galactose mutarotase (GALM) gene: an integrated computational approach

Missense Non-synonymous single nucleotide polymorphisms (nsSNPs) of Galactose Mutarotase (GALM) are associated with the Novel type of Galactosemia (Galactosemia type 4) together with symptoms such as high blood galactose levels and eye cataracts. The objective of the present study was to identify deleterious nsSNPs of GALM recorded on the dbSNP database through comprehensive insilico analysis. Among the 319 missense nsSNPs reported, various insilco tools predicted R78S, R82G, A163E, P210S, Y281C, E307G and F339C as the most deleterious mutations. Structural analysis, PTM analysis and molecular dynamics simulations (MDS) were carried out to understand the effect of these mutations on the structural and physicochemical properties of the GALM protein. The residues R82G and E307G were found to be part of the binding site that resulted in decreased surface accessibility. Replacing the charged wild-type residue with a neutral mutant type affected its substrate binding. All 7 mutations were found to increase the rigidity of the protein structure, which is unfavorable during ligand binding. The mutation F339E made the protein structure more rigid than all the other mutations. Y281 is a phosphorylated site, and therefore, less significant structural changes were observed when compared to other mutations; however, it may have significant differences in the usual functioning of the protein. In summary, the structural and functional analysis of missense SNPs of GALM is important to reduce the number of potential mutations to be evaluated in vitro to understand the association with some genetic diseases.Communicated by Ramaswamy H. Sarma.

Scribble promotes fibrosis-dependent mechanisms of hepatocarcinogenesis by p53/PUMA-mediated glycolysis

BACKGROUNDS AND AIMS

Liver cancer is the sixth most common type of cancer and the fifth leading cause of cancer mortality worldwide. Scribble has been shown to function as a neoplastic tumor suppressor gene in most tumors. Our previous studies reported that down-regulation or mislocalization of Scribble was sufficient to initiate mammary tumorigenesis and NSCLC. Recently, it was reported that Scribble was highly expressed in hepatocellular carcinoma (HCC). We aim to study how it was up-regulated and the contradictory role of Scribble in HCC.

METHODS AND RESULTS

Using a mouse model of carbon tetrachloride (CCl4)-induced liver fibrosis system, we showed that Scribble was over-expressed and which may protect the mice against hepatic fibrosis. Unexpectedly, we found out the potential for Scribble to act as a tumor driver at the advanced stage of N-nitrosodiethylamine (DEN) plus CCl4 induced HCC mice model in vivo. In addition, we observed even higher expression of Scribble in HCC tumors harboring elevated levels of wild-type p53. Most importantly, nuclear translocated Scribble could interact with p53, which lead to enhanced stability and transcriptional activity of p53. Mechanistically, our data suggested that Scribble might drive HCC progression by promoting metabolic regulation of p53 through p53-upregulated modulator of apoptosis (PUMA)-mediated Warburg effect.

CONCLUSIONS

Our data identified the molecular basis of hepatic fibrosis-specific gene expression of polarity gene, such as Scribble. Interestingly, with the progression from fibrosis to cirrhosis to HCC, its nuclear translocation promoted a wild-type p53-mediated cancer metabolic switch and tumor progression in HCC. Taken together, we demonstrated that Scribble was up-regulated and served a protective role in liver fibrosis, while also apparently acting as a tumor driver in fibrosis-dependent hepatocarcinogenesis.

CDK13 promotes lipid deposition and prostate cancer progression by stimulating NSUN5-mediated m5C modification of ACC1 mRNA

Cyclin-dependent kinases (CDKs) regulate cell cycle progression and the transcription of a number of genes, including lipid metabolism-related genes, and aberrant lipid metabolism is involved in prostate carcinogenesis. Previous studies have shown that CDK13 expression is upregulated and fatty acid synthesis is increased in prostate cancer (PCa). However, the molecular mechanisms linking CDK13 upregulation and aberrant lipid metabolism in PCa cells remain largely unknown. Here, we showed that upregulation of CDK13 in PCa cells increases the fatty acyl chains and lipid classes, leading to lipid deposition in the cells, which is positively correlated with the expression of acetyl-CoA carboxylase (ACC1), the first rate-limiting enzyme in fatty acid synthesis. Gain- and loss-of-function studies showed that ACC1 mediates CDK13-induced lipid accumulation and PCa progression by enhancing lipid synthesis. Mechanistically, CDK13 interacts with RNA-methyltransferase NSUN5 to promote its phosphorylation at Ser327. In turn, phosphorylated NSUN5 catalyzes the m5C modification of ACC1 mRNA, and then the m5C-modified ACC1 mRNA binds to ALYREF to enhance its stability and nuclear export, thereby contributing to an increase in ACC1 expression and lipid deposition in PCa cells. Overall, our results disclose a novel function of CDK13 in regulating the ACC1 expression and identify a previously unrecognized CDK13/NSUN5/ACC1 pathway that mediates fatty acid synthesis and lipid accumulation in PCa cells, and targeting this newly identified pathway may be a novel therapeutic option for the treatment of PCa.

A Review of Machine Learning and Algorithmic Methods for Protein Phosphorylation Site Prediction

Post-translational modifications (PTMs) have key roles in extending the functional diversity of proteins and, as a result, regulating diverse cellular processes in prokaryotic and eukaryotic organisms. Phosphorylation modification is a vital PTM that occurs in most proteins and plays a significant role in many biological processes. Disorders in the phosphorylation process lead to multiple diseases, including neurological disorders and cancers. The purpose of this review is to organize this body of knowledge associated with phosphorylation site (p-site) prediction to facilitate future research in this field. At first, we comprehensively review all related databases and introduce all steps regarding dataset creation, data preprocessing, and method evaluation in p-site prediction. Next, we investigate p-site prediction methods, which are divided into two computational groups: algorithmic and machine learning (ML). Additionally, it is shown that there are basically two main approaches for p-site prediction by ML: conventional and end-to-end deep learning methods, both of which are given an overview. Moreover, this review introduces the most important feature extraction techniques, which have mostly been used in p-site prediction. Finally, we create three test sets from new proteins related to the released version of the database of protein post-translational modifications (dbPTM) in 2022 based on general and human species. Evaluating online p-site prediction tools on newly added proteins introduced in the dbPTM 2022 release, distinct from those in the dbPTM 2019 release, reveals their limitations. In other words, the actual performance of these online p-site prediction tools on unseen proteins is notably lower than the results reported in their respective research papers.

DeepMPSF: A Deep Learning Network for Predicting General Protein Phosphorylation Sites Based on Multiple Protein Sequence Features

Phosphorylation, as one of the most important post-translational modifications, plays a key role in various cellular physiological processes and disease occurrences. In recent years, computer technology has been gradually applied to the prediction of protein phosphorylation sites. However, most existing methods rely on simple protein sequence features that provide limited contextual information. To overcome this limitation, we propose DeepMPSF, a phosphorylation site prediction model based on multiple protein sequence features. There are two types of features: sequence semantic features, which comprise protein residue type information and relative position information within protein sequence, and protein background biophysical features, which include global semantic information containing more comprehensive protein background information obtained from pretrained models. To extract these features, DeepMPSF employs two separate subnetworks: the S71SFE module and the BBFE module, which automatically extract high-level semantic features. Our model incorporates a learning strategy for handling imbalanced datasets through ensemble learning during training and prediction. DeepMPSF is trained and evaluated on a well-established dataset of human proteins. Comparing the analysis with other benchmark methods reveals that DeepMPSF outperforms in predicting both S/T residues and Y residues. In particular, DeepMPSF showed excellent generalization performance in cross-species blind test performance, with an average improvement of 5.63%/5.72%, 22.28%/25.94%, 20.11%/17.49%, and 26.40%/28.33% for Mus musculus/Rattus norvegicus test sets in area under curves (AUCs) of ROC curve, AUC of the PR curve, F1-score, and MCC metrics, respectively. Furthermore, it also shows excellent performance in the latest updated case of natural proteins with functional phosphorylation sites. Through an ablation study and visual analysis, we uncover that the design of different feature modules significantly contributes to the accurate classification of DeepMPSF, which provides valuable insights for predicting phosphorylation sites and offers effective support for future downstream research.

Schnurri-3 drives tumor growth and invasion in cancer cells expressing interleukin-13 receptor alpha 2

Interleukin 13 receptor alpha 2 (IL13Rα2) is a relevant therapeutic target in glioblastoma (GBM) and other tumors associated with tumor growth and invasion. In a previous study, we demonstrated that protein tyrosine phosphatase 1B (PTP1B) is a key mediator of the IL-13/IL13Rα2 signaling pathway. PTP1B regulates cancer cell invasion through Src activation. However, PTP1B/Src downstream signaling mechanisms that modulate the invasion process remain unclear. In the present research, we have characterized the PTP1B interactome and the PTP1B-associated phosphoproteome after IL-13 treatment, in different cellular contexts, using proteomic strategies. PTP1B was associated with proteins involved in signal transduction, vesicle transport, and with multiple proteins from the NF-κB signaling pathway, including Tenascin-C (TNC). PTP1B participated with NF-κB in TNC-mediated proliferation and invasion. Analysis of the phosphorylation patterns obtained after PTP1B activation with IL-13 showed increased phosphorylation of the transcription factor Schnurri-3 (SHN3), a reported competitor of NF-κB. SHN3 silencing caused a potent inhibition in cell invasion and proliferation, associated with a down-regulation of the Wnt/β-catenin pathway, an extensive decline of MMP9 expression and the subsequent inhibition of tumor growth and metastasis in mouse models. Regarding clinical value, high expression of SHN3 was associated with poor survival in GBM, showing a significant correlation with the classical and mesenchymal subtypes. In CRC, SHN3 expression showed a preferential association with the mesenchymal subtypes CMS4 and CRIS-B. Moreover, SHN3 expression strongly correlated with IL13Rα2 and MMP9-associated poor prognosis in different cancers. In conclusion, we have uncovered the participation of SNH3 in the IL-13/IL13Rα2/PTP1B pathway to promote tumor growth and invasion. These findings support a potential therapeutic value for SHN3.

Emerging trends in post-translational modification: Shedding light on Glioblastoma multiforme

Recent multi-omics studies, including proteomics, transcriptomics, genomics, and metabolomics have revealed the critical role of post-translational modifications (PTMs) in the progression and pathogenesis of Glioblastoma multiforme (GBM). Further, PTMs alter the oncogenic signaling events and offer a novel avenue in GBM therapeutics research through PTM enzymes as potential biomarkers for drug targeting. In addition, PTMs are critical regulators of chromatin architecture, gene expression, and tumor microenvironment (TME), that play a crucial function in tumorigenesis. Moreover, the implementation of artificial intelligence and machine learning algorithms enhances GBM therapeutics research through the identification of novel PTM enzymes and residues. Herein, we briefly explain the mechanism of protein modifications in GBM etiology, and in altering the biologics of GBM cells through chromatin remodeling, modulation of the TME, and signaling pathways. In addition, we highlighted the importance of PTM enzymes as therapeutic biomarkers and the role of artificial intelligence and machine learning in protein PTM prediction.

SARS-CoV-2 variants evolve convergent strategies to remodel the host response

SARS-CoV-2 variants of concern (VOCs) emerged during the COVID-19 pandemic. Here, we used unbiased systems approaches to study the host-selective forces driving VOC evolution. We discovered that VOCs evolved convergent strategies to remodel the host by modulating viral RNA and protein levels, altering viral and host protein phosphorylation, and rewiring virus-host protein-protein interactions. Integrative computational analyses revealed that although Alpha, Beta, Gamma, and Delta ultimately converged to suppress interferon-stimulated genes (ISGs), Omicron BA.1 did not. ISG suppression correlated with the expression of viral innate immune antagonist proteins, including Orf6, N, and Orf9b, which we mapped to specific mutations. Later Omicron subvariants BA.4 and BA.5 more potently suppressed innate immunity than early subvariant BA.1, which correlated with Orf6 levels, although muted in BA.4 by a mutation that disrupts the Orf6-nuclear pore interaction. Our findings suggest that SARS-CoV-2 convergent evolution overcame human adaptive and innate immune barriers, laying the groundwork to tackle future pandemics.

PIM1 controls GBP1 activity to limit self-damage and to guard against pathogen infection

Disruption of cellular activities by pathogen virulence factors can trigger innate immune responses. Interferon-γ (IFN-γ)-inducible antimicrobial factors, such as the guanylate binding proteins (GBPs), promote cell-intrinsic defense by attacking intracellular pathogens and by inducing programmed cell death. Working in human macrophages, we discovered that GBP1 expression in the absence of IFN-γ killed the cells and induced Golgi fragmentation. IFN-γ exposure improved macrophage survival through the activity of the kinase PIM1. PIM1 phosphorylated GBP1, leading to its sequestration by 14-3-3σ, which thereby prevented GBP1 membrane association. During Toxoplasma gondii infection, the virulence protein TgIST interfered with IFN-γ signaling and depleted PIM1, thereby increasing GBP1 activity. Although infected cells can restrain pathogens in a GBP1-dependent manner, this mechanism can protect uninfected bystander cells. Thus, PIM1 can provide a bait for pathogen virulence factors, guarding the integrity of IFN-γ signaling.

Survey of Acetylation for Thermoanaerobacter tengcongensis

Non-histone protein acetylation is involved in key cellular processes both in eukaryotes and prokaryotes. Acetylation in bacteria is used to modify proteins involved in metabolism and allow the bacteria to adapt to their environment. TTE (Thermoanaerobacter tengcongensis) is an anaerobic, thermophilic saccharolytic bacterium that grows at extreme temperature range between 50 and 80 ℃. The annotated TTE proteome contains less than 3000 proteins. We analyzed the proteome and acetylome of TTE using 2DLC-MS/MS (2-dimensional liquid chromatography mass spectrum). We evaluated the ability of mass spectrometry technology to cover a relatively small proteome as much as possible. And we also observed wide spread of acetylation in TTE, which changed under different temperatures. A total of 2082 proteins were identified, which accounts for about 82% of the database. A total of 2050 (~ 98%) proteins were quantified in at least one culture condition and 1818 proteins were quantified in all 4 conditions. The result also consisted 3457 acetylation sites corresponding to 827 distinct proteins, which covered 40% of the proteins identified. Bioinformatics analysis reported that proteins related to replication, recombination, repair, and extracellular structure cell wall biogenesis had more than half members acetylated, while energy production, carbohydrate transport, and metabolism related proteins were least acetylated. Our result suggested that acetylation affects the ATP-related energy metabolism and energy-dependent biosynthesis process. Comparing the enzymes related with lysine acetylation and acetyl-CoA (acetyl-coenzyme A) metabolism, we suggested that the acetylation of TTE took a non-enzymatic mechanism and affected by abundance of acetyl-CoA.

Phosphorylation of Schizosaccharomyces pombe Dss1 mediates direct binding to the ubiquitin-ligase Dma1 in vitro

Intrinsically disordered proteins (IDPs) are often multifunctional and frequently posttranslationally modified. Deleted in split hand/split foot 1 (Dss1-Sem1 in budding yeast) is a highly multifunctional IDP associated with a range of protein complexes. However, it remains unknown if the different functions relate to different modified states. In this work, we show that Schizosaccharomyces pombe Dss1 is a substrate for casein kinase 2 in vitro, and we identify three phosphorylated threonines in its linker region separating two known disordered ubiquitin-binding motifs. Phosphorylations of the threonines had no effect on ubiquitin-binding but caused a slight destabilization of the C-terminal α-helix and mediated a direct interaction with the forkhead-associated (FHA) domain of the RING-FHA E3-ubiquitin ligase defective in mitosis 1 (Dma1). The phosphorylation sites are not conserved and are absent in human Dss1. Sequence analyses revealed that the Txx(E/D) motif, which is important for phosphorylation and Dma1 binding, is not linked to certain branches of the evolutionary tree. Instead, we find that the motif appears randomly, supporting the mechanism of ex nihilo evolution of novel motifs. In support of this, other threonine-based motifs, although frequent, are nonconserved in the linker, pointing to additional functions connected to this region. We suggest that Dss1 acts as an adaptor protein that docks to Dma1 via the phosphorylated FHA-binding motifs, while the C-terminal α-helix is free to bind mitotic septins, thereby stabilizing the complex. The presence of Txx(D/E) motifs in the disordered regions of certain septin subunits may be of further relevance to the formation and stabilization of these complexes.

Phosphorylation of KRT8 (keratin 8) by excessive mechanical load-activated PKN (protein kinase N) impairs autophagosome initiation and contributes to disc degeneration

Excessive mechanical load (overloading) is a well-documented pathogenetic factor for many mechano stress-induced pathologies, i.e. intervertebral disc degeneration (IDD). Under overloading, the balance between anabolism and catabolism within nucleus pulposus (NP) cells are badly thrown off, and NP cells undergo apoptosis. However, little is known about how the overloading is transduced to the NP cells and contributes to disc degeneration. The current study shows that conditional knockout of Krt8 (keratin 8) within NP aggravates load-induced IDD in vivo, and overexpression of Krt8 endows NP cells greater resistance to overloading-induced apoptosis and degeneration in vitro. Discovery-driven experiments shows that phosphorylation of KRT8 on Ser43 by overloading activated RHOA-PKN (protein kinase N) impedes trafficking of Golgi resident small GTPase RAB33B, suppresses the autophagosome initiation and contributes to IDD. Overexpression of Krt8 and knockdown of Pkn1 and Pkn2, at an early stage of IDD, ameliorates disc degeneration; yet only knockdown of Pkn1 and Pkn2, when treated at late stage of IDD, shows a therapeutic effect. This study validates a protective role of Krt8 during overloading-induced IDD and demonstrates that targeting overloading activation of PKNs could be a novel and effective approach to mechano stress-induced pathologies with a wider window of therapeutic opportunity.Abbreviations: AAV: adeno-associated virus; AF: anulus fibrosus; ANOVA: analysis of variance; ATG: autophagy related; BSA: bovine serum albumin; cDNA: complementary deoxyribonucleic acid; CEP: cartilaginous endplates; CHX: cycloheximide; cKO: conditional knockout; Cor: coronal plane; CT: computed tomography; Cy: coccygeal vertebra; D: aspartic acid; DEG: differentially expressed gene; DHI: disc height index; DIBA: dot immunobinding assay; dUTP: 2'-deoxyuridine 5'-triphosphate; ECM: extracellular matrix; EDTA: ethylene diamine tetraacetic acid; ER: endoplasmic reticulum; FBS: fetal bovine serum; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; GPS: group-based prediction system; GSEA: gene set enrichment analysis; GTP: guanosine triphosphate; HE: hematoxylin-eosin; HRP: horseradish peroxidase; IDD: intervertebral disc degeneration; IF: immunofluorescence staining; IL1: interleukin 1; IVD: intervertebral disc; KEGG: Kyoto encyclopedia of genes and genomes; KRT8: keratin 8; KD: knockdown; KO: knockout; L: lumbar vertebra; LBP: low back pain; LC/MS: liquid chromatograph mass spectrometer; LSI: mouse lumbar instability model; MAP1LC3/LC3: microtubule associated protein 1 light chain 3; MMP3: matrix metallopeptidase 3; MRI: nuclear magnetic resonance imaging; NC: negative control; NP: nucleus pulposus; PBS: phosphate-buffered saline; PE: p-phycoerythrin; PFA: paraformaldehyde; PI: propidium iodide; PKN: protein kinase N; OE: overexpression; PTM: post translational modification; PVDF: polyvinylidene fluoride; qPCR: quantitative reverse-transcriptase polymerase chain reaction; RHOA: ras homolog family member A; RIPA: radio immunoprecipitation assay; RNA: ribonucleic acid; ROS: reactive oxygen species; RT: room temperature; TCM: rat tail compression-induced IDD model; TCS: mouse tail suturing compressive model; S: serine; Sag: sagittal plane; SD rats: Sprague-Dawley rats; shRNA: short hairpin RNA; siRNA: small interfering RNA; SOFG: safranin O-fast green; SQSTM1: sequestosome 1; TUNEL: terminal deoxynucleotidyl transferase dUTP nick end labeling; VG/ml: viral genomes per milliliter; WCL: whole cell lysate.

Mitochondrial polymorphism m.3017C>T of SHLP6 relates to heterothermy

Heterothermic thermoregulation requires intricate regulation of metabolic rate and activation of pro-survival factors. Eliciting these responses and coordinating the necessary energy shifts likely involves retrograde signalling by mitochondrial-derived peptides (MDPs). Members of the group were suggested before to play a role in heterothermic physiology, a key component of hibernation and daily torpor. Here we studied the mitochondrial single-nucleotide polymorphism (SNP) m.3017C>T that resides in the evolutionarily conserved gene MT-SHLP6. The substitution occurring in several mammalian orders causes truncation of SHLP6 peptide size from twenty to nine amino acids. Public mass spectrometric (MS) data of human SHLP6 indicated a canonical size of 20 amino acids, but not the use of alternative translation initiation codons that would expand the peptide. The shorter isoform of SHLP6 was found in heterothermic rodents at higher frequency compared to homeothermic rodents (p < 0.001). In heterothermic mammals it was associated with lower minimal body temperature (T b, p < 0.001). In the thirteen-lined ground squirrel, brown adipose tissue-a key organ required for hibernation, showed dynamic changes of the steady-state transcript level of mt-Shlp6. The level was significantly higher before hibernation and during interbout arousal and lower during torpor and after hibernation. Our finding argues to further explore the mode of action of SHLP6 size isoforms with respect to mammalian thermoregulation and possibly mitochondrial retrograde signalling.

Loss of function of the ALS-associated NEK1 kinase disrupts microtubule homeostasis and nuclear import

Loss-of-function variants in NIMA-related kinase 1 (NEK1) constitute a major genetic cause of amyotrophic lateral sclerosis (ALS), accounting for 2 to 3% of all cases. However, how NEK1 mutations cause motor neuron (MN) dysfunction is unknown. Using mass spectrometry analyses for NEK1 interactors and NEK1-dependent expression changes, we find functional enrichment for proteins involved in the microtubule cytoskeleton and nucleocytoplasmic transport. We show that α-tubulin and importin-β1, two key proteins involved in these processes, are phosphorylated by NEK1 in vitro. NEK1 is essential for motor control and survival in Drosophila models in vivo, while using several induced pluripotent stem cell (iPSC)-MN models, including NEK1 knockdown, kinase inhibition, and a patient mutation, we find evidence for disruptions in microtubule homeostasis and nuclear import. Notably, stabilizing microtubules with two distinct classes of drugs restored NEK1-dependent deficits in both pathways. The capacity of NEK1 to modulate these processes that are critically involved in ALS pathophysiology renders this kinase a formidable therapeutic candidate.

Molecular profiling of sponge deflation reveals an ancient relaxant-inflammatory response

A hallmark of animals is the coordination of whole-body movement. Neurons and muscles are central to this, yet coordinated movements also exist in sponges that lack these cell types. Sponges are sessile animals with a complex canal system for filter-feeding. They undergo whole-body movements resembling "contractions" that lead to canal closure and water expulsion. Here, we combine 3D optical coherence microscopy, pharmacology, and functional proteomics to elucidate anatomy, molecular physiology, and control of these movements. We find them driven by the relaxation of actomyosin stress fibers in epithelial canal cells, which leads to whole-body deflation via collapse of the incurrent and expansion of the excurrent system, controlled by an Akt/NO/PKG/A pathway. A concomitant increase in reactive oxygen species and secretion of proteinases and cytokines indicate an inflammation-like state reminiscent of vascular endothelial cells experiencing oscillatory shear stress. This suggests an ancient relaxant-inflammatory response of perturbed fluid-carrying systems in animals.

Phosphorylation of AHR by PLK1 promotes metastasis of LUAD via DIO2-TH signaling

Metastasis of Lung adenocarcinoma (LUAD) is a major cause of death in patients. Aryl hydrocarbon receptor (AHR) is an important transcription factor involved in the initiation and progression of lung cancer. Polo-like kinase 1 (PLK1), a serine/threonine kinase, is an oncogene that promotes the malignancy of multiple cancer types. Nonetheless, the interaction between these two factors and significance in lung cancer remains to be determined. Here, we demonstrate that PLK1 phosphorylates AHR at S489 in LUAD, which leads to epithelial-mesenchymal transition (EMT) and metastatic events. RNA-seq analyses show that type 2 deiodinase (DIO2) is responsible for EMT and enhanced metastatic potential. DIO2 converts tetraiodothyronine (T4) to triiodothyronine (T3), which then activates thyroid hormone signaling. In vitro and in vivo experiments demonstrate that treatment with T3 or T4 promotes the metastasis of LUAD, whereas depletion of DIO2 or deiodinase inhibitor disrupts this property. Taken together, our results identify the phosphorylation of AHR by PLK1 as a mechanism leading to the progression of LUAD and provide possible therapeutic interventions for this event.