Antonio Ramos-Fernández

Improved Method for Differential Expression Proteomics Using Trypsin-catalyzed 18O Labeling with a Correction for Labeling Efficiency

Quantitative strategies relying on stable isotope labeling and isotope dilution mass spectrometry have proven to be a very robust alternative to the well established gel-based techniques for the study of the dynamic proteome. Postdigestion 18O labeling is becoming very popular mainly due to the simplicity of the enzyme-catalyzed exchange reaction, the peptide handling and storage procedures, and the flexibility and versatility introduced by decoupling protein digestion from peptide labeling. Despite recent progresses, peptide quantification by postdigestion 18O labeling still involves several computational problems. In this work we analyzed the behavior of large collections of peptides when they were subjected to postdigestion labeling and concluded that this process can be explained by a universal kinetic model. On the basis of this observation, we developed an advanced quantification algorithm for this kind of labeling. Our method fits the entire isotopic envelope to parameters related with the kinetic exchange model, allowing at the same time an accurate calculation of the relative proportion of peptides in the original samples and of the specific labeling efficiency of each one of the peptides. We demonstrated that the new method eliminates artifacts produced by incomplete oxygen exchange in subsets of peptides that have a relatively low labeling efficiency and that may be considered indicative of false protein ratio deviations. Finally using a rigorous statistical analysis based on the calculation of error rates associated with false expression changes, we showed the validity of the method in the practice by detecting significant expression changes, produced by the activation of a model preparation of T cells, with only 5 μg of protein in three proteins among a pool of more than 100. By allowing a full control over potential artifacts, our method may improve automation of the procedures for relative protein quantification using this labeling strategy.

Properties of Average Score Distributions of SEQUEST The Probability Ratio Method

High throughput identification of peptides in databases from tandem mass spectrometry data is a key technique in modern proteomics. Common approaches to interpret large scale peptide identification results are based on the statistical analysis of average score distributions, which are constructed from the set of best scores produced by large collections of MS/MS spectra by using searching engines such as SEQUEST. Other approaches calculate individual peptide identification probabilities on the basis of theoretical models or from single-spectrum score distributions constructed by the set of scores produced by each MS/MS spectrum. In this work, we study the mathematical properties of average SEQUEST score distributions by introducing the concept of spectrum quality and expressing these average distributions as compositions of single-spectrum distributions. We predict and demonstrate in the practice that average score distributions are dominated by the quality distribution in the spectra collection, except in the low probability region, where it is possible to predict the dependence of average probability on database size. Our analysis leads to a novel indicator, the probability ratio, which takes optimally into account the statistical information provided by the first and second best scores. The probability ratio is a non-parametric and robust indicator that makes spectra classification according to parameters such as charge state unnecessary and allows a peptide identification performance, on the basis of false discovery rates, that is better than that obtained by other empirical statistical approaches. The probability ratio also compares favorably with statistical probability indicators obtained by the construction of single-spectrum SEQUEST score distributions. These results make the robustness, conceptual simplicity, and ease of automation of the probability ratio algorithm a very attractive alternative to determine peptide identification confidences and error rates in high throughput experiments.

Generalized Method for Probability-based Peptide and Protein Identification from Tandem Mass Spectrometry Data and Sequence Database Searching

Tandem mass spectrometry-based proteomics is currently in great demand of computational methods that facilitate the elimination of likely false positives in peptide and protein identification. In the last few years, a number of new peptide identification programs have been described, but scores or other significance measures reported by these programs cannot always be directly translated into an easy to interpret error rate measurement such as the false discovery rate. In this work we used generalized lambda distributions to model frequency distributions of database search scores computed by MASCOT, X!TANDEM with k-score plug-in, OMSSA, and InsPecT. From these distributions, we could successfully estimate p values and false discovery rates with high accuracy. From the set of peptide assignments reported by any of these engines, we also defined a generic protein scoring scheme that enabled accurate estimation of protein-level p values by simulation of random score distributions that was also found to yield good estimates of protein-level false discovery rate. The performance of these methods was evaluated by searching four freely available data sets ranging from 40,000 to 285,000 MS/MS spectra.

A DERL3 -associated defect in the degradation of SLC2A1 mediates the Warburg effect

Cancer cells possess aberrant proteomes that can arise by the disruption of genes involved in physiological protein degradation. Here we demonstrate the presence of promoter CpG island hypermethylation-linked inactivation of DERL3 (Derlin-3), a key gene in the endoplasmic reticulum-associated protein degradation pathway, in human tumours. The restoration of in vitro and in vivo DERL3 activity highlights the tumour suppressor features of the gene. Using the stable isotopic labelling of amino acids in cell culture workflow for differential proteome analysis, we identify SLC2A1 (glucose transporter 1, GLUT1) as a downstream target of DERL3. Most importantly, SLC2A1 overexpression mediated by DERL3 epigenetic loss contributes to the Warburg effect in the studied cells and pinpoints a subset of human tumours with greater vulnerability to drugs targeting glycolysis.

Increased Diversity of the HLA-B40 Ligandome by the Presentation of Peptides Phosphorylated at Their Main Anchor Residue

Human leukocyte antigen (HLA) class I molecules bind peptides derived from the intracellular degradation of endogenous proteins and present them to cytotoxic T lymphocytes, allowing the immune system to detect transformed or virally infected cells. It is known that HLA class I–associated peptides may harbor posttranslational modifications. In particular, phosphorylated ligands have raised much interest as potential targets for cancer immunotherapy. By combining affinity purification with high-resolution mass spectrometry, we identified more than 2000 unique ligands bound to HLA-B40. Sequence analysis revealed two major anchor motifs: aspartic or glutamic acid at peptide position 2 (P2) and methionine, phenylalanine, or aliphatic residues at the C terminus. The use of immobilized metal ion and TiO2 affinity chromatography allowed the characterization of 85 phosphorylated ligands. We further confirmed every sequence belonging to this subset by comparing its experimental MS2 spectrum with that obtained upon fragmentation of the corresponding synthetic peptide. Remarkably, three phospholigands lacked a canonical anchor residue at P2, containing phosphoserine instead. Binding assays showed that these peptides bound to HLA-B40 with high affinity. Together, our data demonstrate that the peptidome of a given HLA allotype can be broadened by the presentation of peptides with posttranslational modifications at major anchor positions. We suggest that ligands with phosphorylated residues at P2 might be optimal targets for T-cell-based cancer immunotherapy.

A molecular basis for the presentation of phosphorylated peptides by HLA-B antigens

As aberrant protein phosphorylation is a hallmark of tumor cells, the display of tumor-specific phosphopeptides by Human Leukocyte Antigen (HLA) class I molecules can be exploited in the treatment of cancer by T-cell-based immunotherapy. Yet, the characterization and prediction of HLA-I phospholigands is challenging as the molecular determinants of the presentation of such post-translationally modified peptides are not fully understood. Here, we employed a peptidomic workflow to identify 256 unique phosphorylated ligands associated with HLA-B*40, -B*27, -B*39, or -B*07. Remarkably, these phosphopeptides showed similar molecular features. Besides the specific anchor motifs imposed by the binding groove of each allotype, the predominance of phosphorylation at peptide position 4 (P4) became strikingly evident, as was the enrichment of basic residues at P1. To determine the structural basis of this observation, we carried out a series of peptide binding assays and solved the crystal structures of HLA-B*40 in complex with a phosphorylated ligand or its nonphosphorylated counterpart. Overall, our data provide a clear explanation to the common motif found in the phosphopeptidomes associated to different HLA-B molecules. The high prevalence of phosphorylation at P4 is dictated by the presence of the conserved residue Arg62 in the heavy chain, a structural feature shared by most HLA-B alleles. In contrast, the preference for basic residues at P1 is allotype-dependent and might be linked to the structure of the A pocket. This molecular understanding of the presentation of phosphopeptides by HLA-B molecules provides a base for the improved prediction and identification of phosphorylated neo-antigens, as potentially used for cancer immunotherapy.

Inverse regulation in the metabolic genes pckA and metE revealed by proteomic analysis of the Salmonella RcsCDB regulon.

The RcsC, RcsD, and RcsB proteins compose a system used by enteric bacteria to sense envelope stress. Signal transmission occurs from the sensor RcsC to the transcriptional regulator RcsB. Accessory proteins, such as IgaA, are known to adjust the response level. In a previous transcriptomic study, we uncovered 85 genes differentially expressed in Salmonella enterica serovar Typhimurium igaA mutants. Here, we extended these observations to proteomics by performing differential isotope-coded protein labeling (ICPL) followed by liquid chromatography-electrospray ionization tandem mass spectrometry. Five-hundred five proteins were identified and quantified, with 75 of them displaying significant changes in response to alterations in the RcsCDB system. Divergent expression at the RNA and protein level was observed for the metabolic genes pckA and metE, involved in gluconeogenesis and methionine synthesis, respectively. When analyzed in diverse environmental conditions, including the intracellular niche of eukaryotic cells, inverse regulation was more evident for metE and in bacteria growing in defined minimal medium or to stationary phase. The RcsCDB system was also shown to repress the synthesis of the small RNA FnrS, previously reported to modulate metE expression. Collectively, these findings provide new insights into post-transcriptional regulatory mechanisms involving the RcsCDB system and its control over metabolic functions.

Comparative Analysis of the Endogenous Peptidomes Displayed by HLA-B*27 and Mamu-B*08: Two MHC Class I Alleles Associated with Elite Control of HIV/SIV Infection

Indian rhesus macaques are arguably the most reliable animal models in AIDS research. In this species the MHC class I allele Mamu-B*08, among others, is associated with elite control of SIV replication. A similar scenario is observed in humans where the expression of HLA-B*27 or HLA-B*57 has been linked to slow or no progression to AIDS after HIV infection. Despite having large differences in their primary structure, it has been reported that HLA-B*27 and Mamu-B*08 display peptides with sequence similarity. To fine-map the Mamu-B*08 binding motif and assess its similarities with that of HLA-B*27, we affinity purified the peptidomes bound to these MHC class I molecules and analyzed them by LC-MS, identifying several thousands of endogenous ligands. Sequence analysis of both sets of peptides revealed a degree of similarity in their binding motifs, especially at peptide position 2 (P2), where arginine was present in the vast majority of ligands of both allotypes. In addition, several differences emerged from this analysis: (i) ligands displayed by Mamu-B*08 tended to be shorter and to have lower molecular weight, (ii) Mamu-B*08 showed a higher preference for glutamine at P2 as a suboptimal binding motif, and (iii) the second major anchor position, found at PΩ, was much more restrictive in Mamu-B*08. In this regard, HLA-B*27 bound efficiently peptides with aliphatic, aromatic (including tyrosine), and basic C-terminal residues while Mamu-B*08 preferred peptides with leucine and phenylalanine in this position. Finally, in silico estimations of binding efficiency and competitive binding assays to Mamu-B*08 of several selected peptides revealed a good correlation between the characterized anchor motif and binding affinity. These results deepen our understanding of the molecular basis of the presentation of peptides by Mamu-B*08 and can contribute to the detection of novel SIV epitopes restricted by this allotype.

Identification of the Missing Protein Hyaluronan Synthase 1 in Human Mesenchymal Stem Cells Derived from Adipose Tissue or Umbilical Cord

Currently, 14% of the human proteome is made up of proteins whose existence is not confirmed by mass spectrometry. We performed a proteomic profiling of human mesenchymal stem cells derived from adipose tissue or umbilical cord (PRIDE accession number: PXD009893) and identified peptides derived from 13 of such missing proteins. Remarkably, we found compelling evidence of the expression of hyaluronan synthase 1 (NX_Q92839-1) and confirmed its identification by the fragmentation of four heavy-labeled peptides that coeluted with their endogenous light counterparts. Our data also suggest that mesenchymal stem cells constitute a promising source for the detection of missing proteins.

Corrigendum: A DERL3 -associated defect in the degradation of SLC2A1 mediates the Warburg effect

Nature Communications 5: Article number: 4608 (2014); Published: 3 April 2014; Updated: 2 November 2016 During the preparation of the PCR gel panels in Fig. 1b,c, duplicate or otherwise irrelevant lanes were excised; these excisions were not noted in the published figure as per the Nature journal policy.