Digging for Missing Proteins Using Low-Molecular-Weight Protein Enrichment and a “Mirror Protease” Strategy

The utility of proteases in proteomics, from sequence profiling to structure and function analysis

In mass spectrometry (MS)-based bottom-up proteomics, protease digestion plays an essential role in profiling both proteome sequences and post-translational modifications (PTMs). Trypsin is the gold standard in digesting intact proteins into small-size peptides, which are more suitable for high-performance liquid chromatography (HPLC) separation and tandem MS (MS/MS) characterization. However, protein sequences lacking Lys and Arg cannot be cleaved by trypsin and may be missed in conventional proteomic analysis. Proteases with cleavage sites complementary to trypsin are widely applied in proteomic analysis to greatly improve the coverage of proteome sequences and PTM sites. In this review, we survey the common and newly emerging proteases used in proteomics analysis mainly in the last 5 years, focusing on their unique cleavage features and specific proteomics applications such as missing protein characterization, new PTM discovery, and de novo sequencing. In addition, we summarize the applications of proteases in structural proteomics and protein function analysis in recent years. Finally, we discuss the future development directions of new proteases and applications in proteomics.

Discovery of Missing Proteins from an Aneuploidy Cell Line Using a Proteogenomic Approach

With the steadfast development of proteomic technology, the number of missing proteins (MPs) has been continuously shrinking, with approximately 1470 MPs that have not been explored yet. Due to this phenomenon, the discovery of MPs has been increasingly more difficult and elusive. In order to face this challenge, we have hypothesized that a stable aneuploid cell line with increased chromosomes serves as a useful material for assisting MP exploration. Ker-CT cell line with trisomy at chromosome 5 and 20 was selected for this purpose. With a combination strategy of RNA-Seq and LC-MS/MS, a total of 22 178 transcripts and 8846 proteins were identified in Ker-CT. Although the transcripts corresponding to 15 and 15 MP genes located at chromosome 5 and 20 were detected, none of the MPs were found in Ker-CT. Surprisingly, 3 MPs containing at least two unique non-nest peptides of length ≥9 amino acids were identified in Ker-CT, whose genes are located on chromosome 3 and 10, respectively. Furthermore, the 3 MPs were verified using the method of parallel reaction monitoring (PRM). These results suggest that the abnormal status of chromosomes may not only impact the expression of the corresponding genes in trisomy chromosomes, but also influence that of other chromosomes, which benefits MP discovery. The data obtained in this study are available via ProteomeXchange (PXD028647) and PeptideAtlas (PASS01700), respectively.

Proteogenomics Integrating Novel Junction Peptide Identification Strategy Discovers Three Novel Protein Isoforms of Human NHSL1 and EEF1B2

In eukaryotes, alternative pre-mRNA splicing allows a single gene to encode different protein isoforms that function in many biological processes, and they are used as biomarkers or therapeutic targets for diseases. Although protein isoforms in the human genome are well annotated, we speculate that some low-abundance protein isoforms may still be under-annotated because most genes have a primary coding product and alternative protein isoforms tend to be under-expressed. A peptide coencoded by a novel exon and an annotated exon separated by an intron is known as a novel junction peptide. In the absence of known transcripts and homologous proteins, traditional whole-genome six-frame translation-based proteogenomics cannot identify novel junction peptides, and it cannot capture novel alternative splice sites. In this article, we first propose a strategy and tool for identifying novel junction peptides, called CJunction, which we then integrate into a proteogenomics process specifically designed for novel protein isoform discovery and apply to the analysis of a deep-coverage HeLa mass spectrometry data set with identifier PXD004452 in ProteomeXchange. We succeeded in identifying and validating three novel protein isoforms of two functionally important genes, NHSL1 (causative gene of Nance-Horan syndrome) and EEF1B2 (translation elongation factor), which validate our hypothesis. These novel protein isoforms have significant sequence differences from the annotated gene-coding products introduced by the novel N-terminal, suggesting that they may play importantly different functions.

Proteogenomics Study of Blastobotrys adeninivorans TMCC 70007-A Dominant Yeast in the Fermentation Process of Pu-erh Tea

Blastobotrys adeninivorans plays an essential role in pile-fermenting of Pu-erh tea. Its ability to assimilate various carbon and nitrogen sources makes it available for application in a wide range of industry sectors. The genome of B. adeninivorans TMCC 70007 isolated from pile-fermented Pu-erh tea was sequenced and assembled. Proteomics analysis indicated that 4900 proteins in TMCC 70007 were expressed under various culture conditions. Proteogenomics mapping revealed 48 previously unknown genes and corrected 118 gene models predicted by GeneMark-ES. Ortho-proteogenomics analysis identified 17 previously unidentified genes in B. adeninivorans LS3, the first strain with a sequenced genome among the genus Blastobotrys as well. More importantly, five species specific genes were identified from TMCC 70007, which could serve as a barcode for strain typing and were applicable for fermentation process protection of this industrial species. The datasets generated from tea aqueous extract culture not only increased the proteome coverage and accuracy but also contributed to the identification of proteins related to polyphenols and caffeine, which were considered to change greatly during the microbial fermentation of Pu-erh tea. This study provides a proteome perspective on TMCC 70007, which was considered to be an important strain in the production of Pu-erh tea. The systematic proteogenomics analysis not only made a better annotation on the genome of B. adeninivorans TMCC 70007 as previous proteogenomics study but also provided solution for fermentation process protection on valuable industrial species with species specific genes uniquely identified from proteogenomics study.

Open-pFind Verified Four Missing Proteins from Multi-Tissues

The Chromosome-Centric Human Proteome Project (C-HPP) was launched in 2012 to perfect the annotation of human protein existence by identifying stronger evidence of the expression of missing proteins (MPs) at the protein level. After an 8 year effort all over the world, the number of MPs in the neXtProt database significantly decreased from 5511 (2012-02-24) to 1899 (2020-01-17). It is now more difficult to provide confident evidence of the remaining MPs because of their specific characteristics, including low abundance, low molecular weight, unexpected modifications, transmembrane structure, tissue-expression specificity, and so on. A higher resolution mass spectrometry (MS) interpretation engine might provide an opportunity to identify these buried MPs in complex samples by the combination with multi-tissue large-scale proteomics. In this study, open-pFind was used to dig MPs from 20 pairs of healthy human tissues by Wang et al. ( Mol. Syst. Biol. 2019, 15 (2), e8503) combined with our large-scale testis data set digested by three enzymes (Glu-C, Lys-C, and trypsin) with specificity for different amino acid residues ( J. Proteme Res. 2019, 18 (12), 4189-4196). A total of 1 535 536 peptides with 17 283 477 peptide-spectrum matches (PSMs) were mapped to 14 279 protein entries at a false discovery rate of <1% at the PSM, peptide, and protein levels. A total of 103 MP candidates were identified, among which 86 candidates had more unique peptide numbers compared with our single testis tissue. After rigorous screening, manual checks, peptide synthesis, and matching with documented peptides from PeptideAtlas, we validated four MPs, P0C7T8 (duodenum and small intestine), Q8WWZ4 (stomach and rectum), Q8IV35 (fallopian tube), and O14921 (tonsil), at the protein level. All MS raw files have been deposited to the ProteomeXchange with identifier PXD021391.

Trastuzumab-induced upregulation of a protein set in extracellular vesicles emitted by ErbB2-positive breast cancer cells correlates with their trastuzumab sensitivity

Background

ErbB2/HER2 oncoprotein often drives breast cancers (BCs) which are treated with the anti-ErbB2 antibody trastuzumab. The efficacy of trastuzumab-based metastatic BC therapies is routinely assessed by imaging studies. Trastuzumab typically becomes ineffective in the case of this disease and is then replaced by other drugs. Biomarkers of BC trastuzumab response could allow imaging studies and the switch to other drugs to occur earlier than is now possible. Moreover, bone-only BC metastases can be hard to measure, and biomarkers of their trastuzumab response could facilitate further treatment decisions. Such biomarkers are presently unavailable. In this study, we searched for proteins whose levels in BC cell-emitted extracellular vesicles (EVs) potentially correlate with BC trastuzumab sensitivity.

Methods

We isolated EVs from cultured trastuzumab-sensitive and trastuzumab-resistant human BC cells before and after trastuzumab treatment and characterized these EVs by nanoparticle tracking analysis and electron microscopy. We found previously that ErbB2 drives BC by downregulating a pro-apoptotic protein PERP. We now tested whether trastuzumab-induced PERP upregulation in EVs emitted by cultured human BC cells correlates with their trastuzumab sensitivity. We also used mass spectrometry to search for additional proteins whose levels in such EVs reflect BC cell trastuzumab sensitivity. Once we identified proteins whose EV levels correlate with this sensitivity in culture, we explored the feasibility of testing whether their levels in the blood EVs of trastuzumab-treated metastatic BC patients correlate with patients’ response to trastuzumab-based treatments.

Results

We found that neither trastuzumab nor acquisition of trastuzumab resistance by BC cells affects the size or morphology of EVs emitted by cultured BC cells. We established that EV levels of proteins PERP, GNAS2, GNA13, ITB1, and RAB10 correlate with BC cell trastuzumab response. Moreover, these proteins were upregulated during trastuzumab-based therapies in the blood EVs of a pilot cohort of metastatic BC patients that benefited from these therapies but not in those derived from patients that failed such treatments.

Conclusions

Upregulation of a protein set in EVs derived from cultured breast tumor cells correlates with tumor cell trastuzumab sensitivity. It is feasible to further evaluate these proteins as biomarkers of metastatic BC trastuzumab response.

Alternative Strategy To Explore Missing Proteins with Low Molecular Weight

Identifying more missing proteins (MPs) is an important mission of C-HPP. With the number of identified MPs being attenuated year by year (2,949 to 2,129 MPs from 2016 to 2019), we have realized that the difficulty of exploring the remaining MPs is a challenge in technique. Herein, we propose a comprehensive strategy to effectively enrich, separate, and identify proteins with low molecular weights, aiming at the discovery of MPs. Basically, a protein extract from human placenta was passed through a C18 SPE column, and the bound proteins that were eluted were further separated with an SDS-PAGE gel or a 50 kDa cutoff filter. The separated proteins were subjected to trypsin digestion, and the MS/MS signals were searched against data sets with two different digestion modes (full-trypsin and semitrypsin). The strategy was adopted, resulting in the identification of 4 MPs with 8 unique peptides (≥2 non-nested unique peptides with ≥9 amino acids). Importantly, the identification of 6 out of 8 of the unique peptides derived from the MPs was further supported by parallel reaction monitoring, which confirmed the identification of 3 MPs from human placenta tissues (Q6NT89: TMF-regulated nuclear protein 1; A0A183: late cornified envelope protein 6A; and Q6UWQ7: insulin growth factor-like family member 2, mapped to chromosomes 1, 1, and 19, respectively). The three proteins ranged in length from 80 aa to 227 aa. The study not only establishes a feasible strategy for analyzing proteins with low molecular weights but also fills a small part of a large gap in the list of MPs. The data obtained in this study are available via ProteomeXchange (PXD014083) and PeptideAtlas (PASS01389).

Utilization of the Proteome Data Deposited in SRMAtlas for Validating the Existence of the Human Missing Proteins in GPM

The Human Proteome Project (HPP) has made great efforts to clarify the existing evidence of human proteins since 2012. However, according to the recent release of neXtProt (2019-1), approximately 10% of all human genes still have inadequate or no experimental evidence of their translation at the protein level. They were categorized as missing proteins (PE2-PE4). To further the goal of HPP, we developed a two-step bioinformatic strategy addressing the utilization of the SRMAtlas synthetic peptides corresponding to the missing proteins as an exclusive reference in order to explore their natural counterparts within GPM. In the first step, we searched the GPM for the non-nested SRMAtlas peptides corresponding to the missing proteins, taking under consideration only those detected via ≥2 non-nested unitypic/proteotypic peptides "Stranded peptides" with length ≥9 amino acids in the same proteomic study. As a result, 51 missing proteins were newly detected in 35 different proteomic studies. In the second step, we validated these newly detected missing proteins based on matching the spectra of their synthetic and natural peptides in SRMAtlas and GPM, respectively. The results showed that 23 of the missing proteins with ≥2 non-nested peptides were validated by careful spectral matching.