In silico Proteome Cleavage Reveals Iterative Digestion Strategy for High Sequence Coverage

In the post-genome era, biologists have sought to measure the complete complement of proteins, termed proteomics. Currently, the most effective method to measure the proteome is with shotgun, or bottom-up, proteomics, in which the proteome is digested into peptides that are identified followed by protein inference. Despite continuous improvements to all steps of the shotgun proteomics workflow, observed proteome coverage is often low; some proteins are identified by a single peptide sequence. Complete proteome sequence coverage would allow comprehensive characterization of RNA splicing variants and all post translational modifications, which would drastically improve the accuracy of biological models. There are many reasons for the sequence coverage deficit, but ultimately peptide length determines sequence observability. Peptides that are too short are lost because they match many protein sequences and their true origin is ambiguous. The maximum observable peptide length is determined by several analytical challenges. This paper explores computationally how peptide lengths produced from several common proteome digestion methods limit observable proteome coverage. Iterative proteome cleavage strategies are also explored. These simulations reveal that maximized proteome coverage can be achieved by use of an iterative digestion protocol involving multiple proteases and chemical cleavages that theoretically allow 91.1% proteome coverage.

💡 Research Summary

The manuscript addresses a central limitation of contemporary shotgun (bottom‑up) proteomics: the inability to achieve comprehensive sequence coverage of the proteome. While advances in sample preparation, chromatography, and mass‑spectrometry have improved peptide identification rates, most proteins are still represented by a single peptide, leaving large portions of their primary sequences unobserved. The authors argue that peptide length is the primary determinant of observability. Peptides that are too short lack uniqueness in a database search and therefore cannot be confidently assigned, whereas peptides that are too long fall outside the optimal m/z range of most mass spectrometers and are often missed during data‑dependent acquisition.

To quantify how different digestion strategies affect the distribution of peptide lengths, the authors performed in‑silico digestions of the entire human proteome (UniProt/Swiss‑Prot) using six conventional proteases (trypsin, Lys‑C, Glu‑C, Asp‑N, Arg‑C, chymotrypsin) and two chemical cleavage reagents (cyanogen bromide, which cleaves after methionine, and N‑terminal cleavage with NTCB). For each enzyme, they calculated the resulting peptide length distribution and identified the fraction of peptides that fall within the “observable window” (5–30 amino acids), a range that balances uniqueness and mass‑spectrometer detectability.

The results show that single‑enzyme digestions achieve only 60–70 % theoretical coverage of the proteome within the observable window. Trypsin, the workhorse of most proteomics labs, produces many peptides of 7–15 residues, but still leaves a substantial proportion of the proteome either as overly short fragments (≤5 residues) or as long stretches (>30 residues) that would be difficult to detect. Other proteases shift the distribution in predictable ways: Lys‑C generates slightly longer peptides, Glu‑C produces a broader spread, and chemical reagents target residues that are rarely cleaved by proteases, thereby providing complementary coverage.

The core innovation of the paper is the concept of “Iterative Digestion,” a sequential combination of proteolytic and chemical cleavages designed to progressively reduce the size of remaining large peptides while avoiding the generation of excessively short fragments. In the simulated workflow, the authors first digest with trypsin, then subject the undigested or partially digested material to Lys‑C, followed by Glu‑C, Asp‑N, and finally the chemical reagents CNBr and NTCB. At each step, peptides already within the observable length range are filtered out, preventing unnecessary over‑digestion and reducing sample complexity for downstream LC‑MS/MS analysis.

When this multi‑step protocol is applied in silico, the theoretical proteome coverage rises dramatically to 91.1 %, a substantial improvement over any single‑enzyme approach. The authors emphasize that this gain is primarily due to the complementary cleavage specificities that together generate a dense set of uniquely mappable peptides across the entire protein sequence, including regions that are typically missed (e.g., methionine‑rich segments, N‑terminal regions, and internal stretches lacking Lys/Arg). The approach also promises better detection of splice variants and post‑translational modifications, because more of the protein backbone is sampled.

The paper does not claim that the 91.1 % figure will be directly realized in the laboratory without further optimization. It acknowledges several practical challenges: (1) enzymatic efficiency and missed cleavages, (2) incomplete or non‑specific reactions of chemical reagents, (3) sample loss and increased handling time associated with multiple digestion steps, and (4) the persisting issue of short peptides that still map to multiple proteins. The authors suggest that automation of the digestion workflow, use of high‑resolution, high‑mass‑accuracy instruments, and refined database search strategies (e.g., peptide‑level false‑discovery‑rate control) will be essential to translate the computational predictions into experimental reality.

In conclusion, the study provides a rigorous computational framework for evaluating how digestion chemistry shapes peptide length distributions and, consequently, proteome coverage. By demonstrating that a carefully designed iterative digestion scheme can theoretically achieve >90 % coverage, the work offers a roadmap for future experimental protocols aiming for near‑complete proteome characterization. If implemented successfully, such strategies could dramatically enhance our ability to map splice isoforms, quantify low‑abundance proteins, and comprehensively profile post‑translational modifications, thereby strengthening the quantitative foundation of systems biology and precision medicine.