In-silico cleavage of polypeptides using the cleaver package
26 October 2021
Abstract
This vignette describes the in-silico cleavage of polypeptides using thecleaver package.
Package
cleaver 1.32.0
1 Introduction
Most proteomics experiments need protein (peptide) separation and cleavage procedures before these molecules could be analyzed or identified by mass spectrometry or other analytical tools.
cleaver allows in-silico cleavage of polypeptide sequences to e.g. create theoretical mass spectrometry data.
The cleavage rules are taken from the ExPASy PeptideCutter tool (Gasteiger et al. 2005).
2 Simple Usage
Loading the cleaver package:
library("cleaver")Getting help and list all available cleavage rules:
help("cleave")Cleaving of Gastric juice peptide 1 (P01358) using Trypsin:
## cleave it
cleave("LAAGKVEDSD", enzym="trypsin")## $LAAGKVEDSD
## [1] "LAAGK" "VEDSD"## get the cleavage ranges
cleavageRanges("LAAGKVEDSD", enzym="trypsin")## $LAAGKVEDSD
## start end
## [1,] 1 5
## [2,] 6 10## get only cleavage sites
cleavageSites("LAAGKVEDSD", enzym="trypsin")## $LAAGKVEDSD
## [1] 5Sometimes cleavage is not perfect and the enzym miss some cleavage positions:
## miss one cleavage position
cleave("LAAGKVEDSD", enzym="trypsin", missedCleavages=1)## $LAAGKVEDSD
## [1] "LAAGKVEDSD"cleavageRanges("LAAGKVEDSD", enzym="trypsin", missedCleavages=1)## $LAAGKVEDSD
## start end
## [1,] 1 10## miss zero or one cleavage positions
cleave("LAAGKVEDSD", enzym="trypsin", missedCleavages=0:1)## $LAAGKVEDSD
## [1] "LAAGK" "VEDSD" "LAAGKVEDSD"cleavageRanges("LAAGKVEDSD", enzym="trypsin", missedCleavages=0:1)## $LAAGKVEDSD
## start end
## [1,] 1 5
## [2,] 6 10
## [3,] 1 10Combine cleaver and Biostrings (Pages et al., n.d.):
## create AAStringSet object
p <- AAStringSet(c(gaju="LAAGKVEDSD", pnm="AGEPKLDAGV"))
## cleave it
cleave(p, enzym="trypsin")## AAStringSetList of length 2
## [["gaju"]] LAAGK VEDSD
## [["pnm"]] AGEPK LDAGVcleavageRanges(p, enzym="trypsin")## IRangesList object of length 2:
## $gaju
## IRanges object with 2 ranges and 0 metadata columns:
## start end width
## <integer> <integer> <integer>
## [1] 1 5 5
## [2] 6 10 5
##
## $pnm
## IRanges object with 2 ranges and 0 metadata columns:
## start end width
## <integer> <integer> <integer>
## [1] 1 5 5
## [2] 6 10 5cleavageSites(p, enzym="trypsin")## $gaju
## [1] 5
##
## $pnm
## [1] 53 Insulin & Somatostatin Example
Downloading Insulin (P01308) and Somatostatin (P61278) sequences from the UniProt (The UniProt Consortium 2012) database using UniProt.ws (Carlson, n.d.).
## load UniProt.ws library
library("UniProt.ws")
## select species Homo sapiens
UniProt.ws <- UniProt.ws(taxId=9606)
## download sequences of Insulin/Somatostatin
s <- select(UniProt.ws, keys=c("P01308", "P61278"), columns=c("SEQUENCE"))## Getting extra data for P01308, P61278## 'select()' returned 1:1 mapping between keys and columns## fetch only sequences
sequences <- setNames(s$SEQUENCE, s$UNIPROTKB)
## remove whitespaces
sequences <- gsub(pattern="[[:space:]]", replacement="", x=sequences)Cleaving using Pepsin:
cleave(sequences, enzym="pepsin")## $P01308
## [1] "MA" "L" "W" "MRLLP"
## [5] "LL" "A" "WGPDPAAA" "F"
## [9] "VNQH" "CGSH" "VEA" "Y"
## [13] "VCGERG" "FF" "YTPKTRREAED" "QVGQVE"
## [17] "GGGPGAGS" "LQP" "LA" "EGS"
## [21] "QKRGIVEQCCTSICS" "Q" "EN" "CN"
##
## $P61278
## [1] "ML" "SCRL" "QCA"
## [4] "L" "AA" "SIV"
## [7] "A" "GCVTGAPSDPRL" "RQ"
## [10] "FL" "QKS" "LAAAAGKQEL"
## [13] "AK" "Y" "AE"
## [16] "SEPNQTENDA" "LEPED" "SQAAEQDEMRL"
## [19] "EL" "QRSANSNPAMAPRERKAGCKN" "FF"
## [22] "W" "KT" "FTSC"4 Isotopic Distribution Of Tryptic Digested Insulin
A common use case of in-silico cleavage is the calculation of the isotopic distribution of peptides (which were enzymatic digested in the in-vitro experimental workflow). Here BRAIN (Claesen et al. 2012; Dittwald et al. 2013) is used to calculate the isotopic distribution of cleaver’s output. (please note: it is only a toy example, e.g. the relation of intensity values between peptides isn’t correct).
## load BRAIN library
library("BRAIN")
## cleave insulin
cleavedInsulin <- cleave(sequences[1], enzym="trypsin")[[1]]
## create empty plot area
plot(NA, xlim=c(150, 4300), ylim=c(0, 1),
xlab="mass", ylab="relative intensity",
main="tryptic digested insulin - isotopic distribution")
## loop through peptides
for (i in seq(along=cleavedInsulin)) {
## count C, H, N, O, S atoms in current peptide
atoms <- BRAIN::getAtomsFromSeq(cleavedInsulin[[i]])
## calculate isotopic distribution
d <- useBRAIN(atoms)
## draw peaks
lines(d$masses, d$isoDistr, type="h", col=2)
}
5 Session Information
## R version 4.1.1 (2021-08-10)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 20.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.14-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.14-bioc/R/lib/libRlapack.so
##
## locale:
## [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
## [3] LC_TIME=en_GB LC_COLLATE=C
## [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8
## [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
## [9] LC_ADDRESS=C LC_TELEPHONE=C
## [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
##
## attached base packages:
## [1] stats4 stats graphics grDevices utils datasets methods
## [8] base
##
## other attached packages:
## [1] BRAIN_1.40.0 lattice_0.20-45 PolynomF_2.0-3
## [4] UniProt.ws_2.34.0 RCurl_1.98-1.5 RSQLite_2.2.8
## [7] cleaver_1.32.0 Biostrings_2.62.0 GenomeInfoDb_1.30.0
## [10] XVector_0.34.0 IRanges_2.28.0 S4Vectors_0.32.0
## [13] BiocGenerics_0.40.0 BiocStyle_2.22.0
##
## loaded via a namespace (and not attached):
## [1] Rcpp_1.0.7 png_0.1-7 assertthat_0.2.1
## [4] digest_0.6.28 utf8_1.2.2 BiocFileCache_2.2.0
## [7] R6_2.5.1 evaluate_0.14 highr_0.9
## [10] httr_1.4.2 pillar_1.6.4 zlibbioc_1.40.0
## [13] rlang_0.4.12 curl_4.3.2 jquerylib_0.1.4
## [16] blob_1.2.2 magick_2.7.3 rmarkdown_2.11
## [19] stringr_1.4.0 bit_4.0.4 compiler_4.1.1
## [22] xfun_0.27 pkgconfig_2.0.3 htmltools_0.5.2
## [25] tidyselect_1.1.1 KEGGREST_1.34.0 tibble_3.1.5
## [28] GenomeInfoDbData_1.2.7 bookdown_0.24 fansi_0.5.0
## [31] crayon_1.4.1 dplyr_1.0.7 dbplyr_2.1.1
## [34] bitops_1.0-7 rappdirs_0.3.3 grid_4.1.1
## [37] jsonlite_1.7.2 lifecycle_1.0.1 DBI_1.1.1
## [40] magrittr_2.0.1 stringi_1.7.5 cachem_1.0.6
## [43] bslib_0.3.1 ellipsis_0.3.2 filelock_1.0.2
## [46] generics_0.1.1 vctrs_0.3.8 tools_4.1.1
## [49] bit64_4.0.5 Biobase_2.54.0 glue_1.4.2
## [52] purrr_0.3.4 fastmap_1.1.0 yaml_2.2.1
## [55] AnnotationDbi_1.56.0 BiocManager_1.30.16 memoise_2.0.0
## [58] knitr_1.36 sass_0.4.0References
Carlson, Marc. n.d. UniProt.ws: R Interface to Uniprot Web Services.
Claesen, Jürgen, Piotr Dittwald, Tomasz Burzykowski, and Dirk Valkenborg. 2012. “An Efficient Method to Calculate the Aggregated Isotopic Distribution and Exact Center-Masses.” Journal of the American Society for Mass Spectrometry 23 (4): 753–63.
Dittwald, Piotr, Jürgen Claesen, Tomasz Burzykowski, Dirk Valkenborg, and Anna Gambin. 2013. “BRAIN: A Universal Tool for High-Throughput Calculations of the Isotopic Distribution for Mass Spectrometry.” Analytical Chemistry 85 (4): 1991–4.
Gasteiger, Elisabeth, Christine Hoogland, Alexandre Gattiker, S’everine Duvaud, Marc R. Wilkins, Ron D. Appel, and Amos Bairoch. 2005. “Protein Identification and Analysis Tools on the Expasy Server.” In The Proteomics Protocols Handbook, edited by John M. Walker, 571–607. Humana Press. https://doi.org/10.1385/1-59259-890-0:571.
Pages, H., P. Aboyoun, R. Gentleman, and S. DebRoy. n.d. Biostrings: String Objects Representing Biological Sequences, and Matching Algorithms.
The UniProt Consortium. 2012. “Reorganizing the Protein Space at the Universal Protein Resource (Uniprot).” Nucleic Acids Research 40 (D1): D71–D75. https://doi.org/10.1093/nar/gkr981.