Panoramic Perspective on Human Phosphosites.

Protein phosphorylation is the most common reversible post-translational modification of proteins and is key in the regulation of many cellular processes. Due to this importance, phosphorylation is extensively studied, resulting in the availability of a large amount of mass spectrometry-based phospho-proteomics data. Here, we leverage the information in these large-scale phospho-proteomics data sets, as contained in Scop3P, to analyze and characterize proteome-wide protein phosphorylation sites (P-sites). First, we set out to differentiate correctly observed P-sites from false-positive sites using five complementary site properties. We then describe the context of these P-sites in terms of the protein structure, solvent accessibility, structural transitions and disorder, and biophysical properties. We also investigate the relative prevalence of disease-linked mutations on and around P-sites. Moreover, we assess the structural dynamics of P-sites in their phosphorylated and unphosphorylated states. As a result, we show how large-scale reprocessing of available proteomics experiments can enable a more reliable view on proteome-wide P-sites. Furthermore, adding the structural context of proteins around P-sites helps uncover possible conformational switches upon phosphorylation. Moreover, by placing sites in different biophysical contexts, we show the differential preference in protein dynamics at phosphorylated sites when compared to the nonphosphorylated counterparts.

Sensitive and Specific Spectral Library Searching with CompOmics Spectral Library Searching Tool and Percolator.

Maintaining high sensitivity while limiting false positives is a key challenge in peptide identification from mass spectrometry data. Here, we investigate the effects of integrating the machine learning-based postprocessor Percolator into our spectral library searching tool COSS (CompOmics Spectral library Searching tool). To evaluate the effects of this postprocessing, we have used 40 data sets from 2 different projects and have searched these against the NIST and MassIVE spectral libraries. The searching is carried out using 2 spectral library search tools, COSS and MSPepSearch with and without Percolator postprocessing, and using sequence database search engine MS-GF+ as a baseline comparator. The addition of the Percolator rescoring step to COSS is effective and results in a substantial improvement in sensitivity and specificity of the identifications. COSS is freely available as open source under the permissive Apache2 license, and binaries and source code are found at https://github.com/compomics/COSS.

Bcl-xL acts as an inhibitor of IP3R channels, thereby antagonizing Ca2+-driven apoptosis.

Anti-apoptotic Bcl-2-family members not only act at mitochondria but also at the endoplasmic reticulum, where they impact Ca2+ dynamics by controlling IP3 receptor (IP3R) function. Current models propose distinct roles for Bcl-2 vs. Bcl-xL, with Bcl-2 inhibiting IP3Rs and preventing pro-apoptotic Ca2+ release and Bcl-xL sensitizing IP3Rs to low [IP3] and promoting pro-survival Ca2+ oscillations. We here demonstrate that Bcl-xL too inhibits IP3R-mediated Ca2+ release by interacting with the same IP3R regions as Bcl-2. Via in silico superposition, we previously found that the residue K87 of Bcl-xL spatially resembled K17 of Bcl-2, a residue critical for Bcl-2's IP3R-inhibitory properties. Mutagenesis of K87 in Bcl-xL impaired its binding to IP3R and abrogated Bcl-xL's inhibitory effect on IP3Rs. Single-channel recordings demonstrate that purified Bcl-xL, but not Bcl-xLK87D, suppressed IP3R single-channel openings stimulated by sub-maximal and threshold [IP3]. Moreover, we demonstrate that Bcl-xL-mediated inhibition of IP3Rs contributes to its anti-apoptotic properties against Ca2+-driven apoptosis. Staurosporine (STS) elicits long-lasting Ca2+ elevations in wild-type but not in IP3R-knockout HeLa cells, sensitizing the former to STS treatment. Overexpression of Bcl-xL in wild-type HeLa cells suppressed STS-induced Ca2+ signals and cell death, while Bcl-xLK87D was much less effective in doing so. In the absence of IP3Rs, Bcl-xL and Bcl-xLK87D were equally effective in suppressing STS-induced cell death. Finally, we demonstrate that endogenous Bcl-xL also suppress IP3R activity in MDA-MB-231 breast cancer cells, whereby Bcl-xL knockdown augmented IP3R-mediated Ca2+ release and increased the sensitivity towards STS, without altering the ER Ca2+ content. Hence, this study challenges the current paradigm of divergent functions for Bcl-2 and Bcl-xL in Ca2+-signaling modulation and reveals that, similarly to Bcl-2, Bcl-xL inhibits IP3R-mediated Ca2+ release and IP3R-driven cell death. Our work further underpins that IP3R inhibition is an integral part of Bcl-xL's anti-apoptotic function.

Massively parallel interrogation of protein fragment secretability using SECRiFY reveals features influencing secretory system transit.

While transcriptome- and proteome-wide technologies to assess processes in protein biogenesis are now widely available, we still lack global approaches to assay post-ribosomal biogenesis events, in particular those occurring in the eukaryotic secretory system. We here develop a method, SECRiFY, to simultaneously assess the secretability of >105 protein fragments by two yeast species, S. cerevisiae and P. pastoris, using custom fragment libraries, surface display and a sequencing-based readout. Screening human proteome fragments with a median size of 50-100 amino acids, we generate datasets that enable datamining into protein features underlying secretability, revealing a striking role for intrinsic disorder and chain flexibility. The SECRiFY methodology generates sufficient amounts of annotated data for advanced machine learning methods to deduce secretability patterns. The finding that secretability is indeed a learnable feature of protein sequences provides a solid base for application-focused studies.

Scop3P: A Comprehensive Resource of Human Phosphosites within Their Full Context.

Protein phosphorylation is a key post-translational modification in many biological processes and is associated to human diseases such as cancer and metabolic disorders. The accurate identification, annotation, and functional analysis of phosphosites are therefore crucial to understand their various roles. Phosphosites are mainly analyzed through phosphoproteomics, which has led to increasing amounts of publicly available phosphoproteomics data. Several resources have been built around the resulting phosphosite information, but these are usually restricted to the protein sequence and basic site metadata. What is often missing from these resources, however, is context, including protein structure mapping, experimental provenance information, and biophysical predictions. We therefore developed Scop3P: a comprehensive database of human phosphosites within their full context. Scop3P integrates sequences (UniProtKB/Swiss-Prot), structures (PDB), and uniformly reprocessed phosphoproteomics data (PRIDE) to annotate all known human phosphosites. Furthermore, these sites are put into biophysical context by annotating each phosphoprotein with per-residue structural propensity, solvent accessibility, disordered probability, and early folding information. Scop3P, available at https://iomics.ugent.be/scop3p, presents a unique resource for visualization and analysis of phosphosites and for understanding of phosphosite structure-function relationships.

The CEP5 Peptide Promotes Abiotic Stress Tolerance, As Revealed by Quantitative Proteomics, and Attenuates the AUX/IAA Equilibrium in Arabidopsis.

Peptides derived from non-functional precursors play important roles in various developmental processes, but also in (a)biotic stress signaling. Our (phospho)proteome-wide analyses of C-TERMINALLY ENCODED PEPTIDE 5 (CEP5)-mediated changes revealed an impact on abiotic stress-related processes. Drought has a dramatic impact on plant growth, development and reproduction, and the plant hormone auxin plays a role in drought responses. Our genetic, physiological, biochemical, and pharmacological results demonstrated that CEP5-mediated signaling is relevant for osmotic and drought stress tolerance in Arabidopsis, and that CEP5 specifically counteracts auxin effects. Specifically, we found that CEP5 signaling stabilizes AUX/IAA transcriptional repressors, suggesting the existence of a novel peptide-dependent control mechanism that tunes auxin signaling. These observations align with the recently described role of AUX/IAAs in stress tolerance and provide a novel role for CEP5 in osmotic and drought stress tolerance.

COSS: A Fast and User-Friendly Tool for Spectral Library Searching.

Spectral similarity searching to identify peptide-derived MS/MS spectra is a promising technique, and different spectrum similarity search tools have therefore been developed. Each of these tools, however, comes with some limitations, mainly because of low processing speed and issues with handling large databases. Furthermore, the number of spectral data formats supported is typically limited, which also creates a threshold to adoption. We have therefore developed COSS (CompOmics Spectral Searching), a new and user-friendly spectral library search tool supporting two scoring functions. COSS also includes decoy spectra generation for result validation. We have benchmarked COSS on three different spectral libraries and compared the results with established spectral searching tools and a sequence database search tool. Our comparison showed that COSS more reliably identifies spectra, is capable of handling large data sets and libraries, and is an easy to use tool that can run on low computer specifications. COSS binaries and source code can be freely downloaded from https://github.com/compomics/COSS.

Publisher Correction: Differences in antigenic sites and other functional regions between genotype A and G mumps virus surface proteins.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Bcl-2 and IP3 compete for the ligand-binding domain of IP3Rs modulating Ca2+ signaling output.

Bcl-2 proteins have emerged as critical regulators of intracellular Ca2+ dynamics by directly targeting and inhibiting the IP3 receptor (IP3R), a major intracellular Ca2+-release channel. Here, we demonstrate that such inhibition occurs under conditions of basal, but not high IP3R activity, since overexpressed and purified Bcl-2 (or its BH4 domain) can inhibit IP3R function provoked by low concentration of agonist or IP3, while fails to attenuate against high concentration of agonist or IP3. Surprisingly, Bcl-2 remained capable of inhibiting IP3R1 channels lacking the residues encompassing the previously identified Bcl-2-binding site (a.a. 1380-1408) located in the ARM2 domain, part of the modulatory region. Using a plethora of computational, biochemical and biophysical methods, we demonstrate that Bcl-2 and more particularly its BH4 domain bind to the ligand-binding domain (LBD) of IP3R1. In line with this finding, the interaction between the LBD and Bcl-2 (or its BH4 domain) was sensitive to IP3 and adenophostin A, ligands of the IP3R. Vice versa, the BH4 domain of Bcl-2 counteracted the binding of IP3 to the LBD. Collectively, our work reveals a novel mechanism by which Bcl-2 influences IP3R activity at the level of the LBD. This allows for exquisite modulation of Bcl-2's inhibitory properties on IP3Rs that is tunable to the level of IP3 signaling in cells.

Scop3D: Online Visualization of Mutation Rates on Protein Structure.

Scop3D is a tool that automatically annotates protein structure with sequence conservation starting from a set of protein sequence variants. We present a complete upgrade and rewrite of Scop3D. We have included a DNA module that allows the analysis of single nucleotide polymorphisms in relation to the structural context of the protein. Scop3D therefore forms a bridge between genomics and protein structure. Moreover, Scop3D is now also available through an intuitive web-interface that makes the tool highly user-friendly.

The study of degradation mechanisms of glyco-engineered plant produced anti-rabies monoclonal antibodies E559 and 62-71-3.

Rabies is an ancient and neglected zoonotic disease caused by the rabies virus, a neurotropic RNA virus that belongs to the Rhabdoviridae family, genus Lyssavirus. It remains an important public health problem as there are cost and health concerns imposed by the current human post exposure prophylaxis therapy. The use of monoclonal antibodies (mAbs) is therefore an attractive alternative. Rabies mostly affects people that reside in resource-limited areas where there are occasional failures in the cold-chain. These environmental changes may upset the stability of the mAbs. This study focused on mAbs 62-71-3 and E559; their structures, responses to freeze/thaw (F/T) and exposure to reactive oxygen species were therefore studied with the aid of a wide range of biophysical and in silico techniques in order to elucidate their stability and identify aggregation prone regions. E559 was found to be less stable than 62-71-3. The complementarity determining regions (CDR) contributed the most to its instability, more specifically: peptides 99EIWD102 and 92ATSPYT97 found in CDR3, Trp33 found in CDR1 and the oxidised Met34. The constant region "158SWNSGALTGHTFPAVL175" was also flagged by the special aggregation propensity (SAP) tool and F/T experiments to be highly prone to aggregation. The E559 peptides "4LQESGSVL11 from the heavy chain and 4LTQSPSSL11 from the light chain, were also highly affected by F/T. These residues may serve as good candidates for mutation, in the aim to bring forward more stable therapeutic antibodies, thus paving a way to a more safe and efficacious antibody-based cocktail treatment against rabies.

Differences in antigenic sites and other functional regions between genotype A and G mumps virus surface proteins.

The surface proteins of the mumps virus, the fusion protein (F) and haemagglutinin-neuraminidase (HN), are key factors in mumps pathogenesis and are important targets for the immune response during mumps virus infection. We compared the predicted amino acid sequences of the F and HN genes from Dutch mumps virus samples from the pre-vaccine era (1957-1982) with mumps virus genotype G strains (from 2004 onwards). Genotype G is the most frequently detected mumps genotype in recent outbreaks in vaccinated communities, especially in Western Europe, the USA and Japan. Amino acid differences between the Jeryl Lynn vaccine strains (genotype A) and genotype G strains were predominantly located in known B-cell epitopes and in N-linked glycosylation sites on the HN protein. There were eight variable amino acid positions specific to genotype A or genotype G sequences in five known B-cell epitopes of the HN protein. These differences may account for the reported antigenic differences between Jeryl Lynn and genotype G strains. We also found amino acid differences in and near sites on the HN protein that have been reported to play a role in mumps virus pathogenesis. These differences may contribute to the occurrence of genotype G outbreaks in vaccinated communities.

The mutational landscape of MYCN, Lin28b and ALKF1174L driven murine neuroblastoma mimics human disease.

Genetically engineered mouse models have proven to be essential tools for unraveling fundamental aspects of cancer biology and for testing novel therapeutic strategies. To optimally serve these goals, it is essential that the mouse model faithfully recapitulates the human disease. Recently, novel mouse models for neuroblastoma have been developed. Here, we report on the further genomic characterization through exome sequencing and DNA copy number analysis of four of the currently available murine neuroblastoma model systems (ALK, Th-MYCN, Dbh-MYCN and Lin28b). The murine tumors revealed a low number of genomic alterations - in keeping with human neuroblastoma - and a positive correlation of the number of genetic lesions with the time to onset of tumor formation was observed. Gene copy number alterations are the hallmark of both murine and human disease and frequently affect syntenic genomic regions. Despite low mutational load, the genes mutated in murine disease were found to be enriched for genes mutated in human disease. Taken together, our study further supports the validity of the tested mouse models for mechanistic and preclinical studies of human neuroblastoma.

Cross-linked peptide identification: A computational forest of algorithms.

Chemical cross-linking analyzed by mass spectrometry (XL-MS) has become an important tool in unravelling protein structure, dynamics, and complex formation. Because the analysis of cross-linked proteins with mass spectrometry results in specific computational challenges, many computational tools have been developed to identify cross-linked peptides from mass spectra and subsequently interpret the identified cross-links within their structural context. In this review, we will provide an overview of the different tools that are currently available to tackle the computational part of an XL-MS experiment. First, we give an introduction on the computational challenges encountered when processing data from a cross-linking experiment. We then discuss available tools to identify peptides that are linked by intact or MS-cleavable cross-linkers, and we provide an overview of tools to interpret cross-linked peptides in the context of protein structure. Finally, we give an outlook on data management and dissemination challenges and opportunities for cross-linking experiments.

A double point mutation at residues Ile14 and Val15 of Bcl-2 uncovers a role for the BH4 domain in both protein stability and function.

B-cell lymphoma 2 (Bcl-2) protein is the archetype apoptosis suppressor protein. The N-terminal Bcl-2-homology 4 (BH4) domain of Bcl-2 is required for the antiapoptotic function of this protein at the mitochondria and endoplasmic reticulum (ER). The involvement of the BH4 domain in Bcl-2's antiapoptotic functions has been proposed based on Gly-based substitutions of the Ile14/Val15 amino acids, two hydrophobic residues located in the center of Bcl-2's BH4 domain. Following this strategy, we recently showed that a BH4-domain-derived peptide in which Ile14 and Val15 have been replaced by Gly residues, was unable to dampen proapoptotic Ca2+ -release events from the ER. Here, we investigated the impact of these mutations on the overall structure, stability, and function of full-length Bcl-2 as a regulator of Ca2+ signaling and cell death. Our results indicate that full-length Bcl-2 Ile14Gly/Val15Gly, in contrast to wild-type Bcl-2, (a) displayed severely reduced structural stability and a shortened protein half-life; (b) failed to interact with Bcl-2-associated X protein (BAX), to inhibit the inositol 1,4,5-trisphosphate receptor (IP3 R) and to protect against Ca2+ -mediated apoptosis. We conclude that the hydrophobic face of Bcl-2's BH4 domain (Ile14, Val15) is an important structural regulatory element by affecting protein stability and turnover, thereby likely reducing Bcl-2's ability to modulate the function of its targets, like IP3 R and BAX. Therefore, Bcl-2 structure/function studies require pre-emptive and reliable determination of protein stability upon introduction of point mutations at the level of the BH4 domain.

Methods to Calculate Spectrum Similarity.

Scoring functions that assess spectrum similarity play a crucial role in many computational mass spectrometry algorithms. These functions are used to compare an experimentally acquired fragmentation (MS/MS) spectrum against two different types of target MS/MS spectra: either against a theoretical MS/MS spectrum derived from a peptide from a sequence database, or against another, previously acquired MS/MS spectrum. The former is typically encountered in database searching, while the latter is used in spectrum clustering and spectral library searching. The comparison between acquired versus theoretical MS/MS spectra is most commonly performed using cross-correlations or probability derived scoring functions, while the comparison of two acquired MS/MS spectra typically makes use of a normalized dot product, especially in spectrum library search algorithms. In addition to these scoring functions, Pearson's or Spearman's correlation coefficients, mean squared error, or median absolute deviation scores can also be used for the same purpose. Here, we describe and evaluate these scoring functions with regards to their ability to assess spectrum similarity for theoretical versus acquired, and acquired versus acquired spectra.

Xilmass: A New Approach toward the Identification of Cross-Linked Peptides.

Chemical cross-linking coupled with mass spectrometry plays an important role in unravelling protein interactions, especially weak and transient ones. Moreover, cross-linking complements several structural determination approaches such as cryo-EM. Although several computational approaches are available for the annotation of spectra obtained from cross-linked peptides, there remains room for improvement. Here, we present Xilmass, a novel algorithm to identify cross-linked peptides that introduces two new concepts: (i) the cross-linked peptides are represented in the search database such that the cross-linking sites are explicitly encoded, and (ii) the scoring function derived from the Andromeda algorithm was adapted to score against a theoretical tandem mass spectrometry (MS/MS) spectrum that contains the peaks from all possible fragment ions of a cross-linked peptide pair. The performance of Xilmass was evaluated against the recently published Kojak and the popular pLink algorithms on a calmodulin-plectin complex data set, as well as three additional, published data sets. The results show that Xilmass typically had the highest number of identified distinct cross-linked sites and also the highest number of predicted cross-linked sites.

The Lyssavirus glycoprotein: A key to cross-immunity.

Rabies is an acute viral encephalomyelitis in warm-blooded vertebrates, caused by viruses belonging to Rhabdovirus family and genus Lyssavirus. Although rabies is categorised as a neglected disease, the rabies virus (RABV) is the most studied amongst Lyssaviruses which show nearly identical infection patterns. In efforts to improving post-exposure prophylaxis, several anti-rabies monoclonal antibodies (mAbs) targeting the glycoprotein (G protein) sites I, II, III and G5 have been characterized. To explore cross-neutralization capacity of available mAbs and discover new possible B-cell epitopes, we have analyzed all available glycoprotein sequences from Lyssaviruses with a focus on sequence variation and conservation. This information was mapped on the structure of a representative G protein. We proposed several possible cross-neutralizing B-cell epitopes (GUVTTTF, WLRTV, REECLD and EHLVVEEL) in complement to the already well-characterized antigenic sites. The research could facilitate development of novel cross-reactive mAbs against RABV and even more broad, against possibly all Lyssavirus members.