Block-based characterization of protease specificity from substrate sequence profile
© The Author(s). 2017
Received: 15 June 2017
Accepted: 26 September 2017
Published: 3 October 2017
The mechanism of action of proteases has been widely studied based on substrate specificity. Prior research has been focused on the amino acids at a single amino acid site, but rarely on combinations of amino acids around the cleavage bond.
We propose a novel block-based approach to reveal the potential combinations of amino acids which may regulate the action of proteases. Using the entropies of eight blocks centered at a cleavage bond, we created a distance matrix for 61 proteases to compare their specificities. After quantitative analysis, we discovered a number of prominent blocks, each of which consists of successive amino acids near a cleavage bond, intuitively characterizing the site cooperation of the substrate sequences.
This approach will help in the discovery of specific substrate sequences which may bridge between proteases and cleavage substrate as more substrate information becomes available.
Proteases are a category of enzymes capable of hydrolyzing peptide bonds and irreversibly modifying functions of substrate proteins. These hydrolyzations and modifications are essential for cell growth and differentiation [1, 2]. Recognition of the target substrate of a protease depends partly on the complementation between the protease active site and the sequence surrounding the scissile bond in the substrate. Proteases have pockets that accommodate substrate residues. Substrate sequences that bind the pockets are indexed by P4, P3, P2, P1, P1’, P2’, P3’, P4’ in order from N-terminal to C-terminal following the convention of Schechter and Berger .
Some proteases show strict specificities on the cleavage sequences of the substrates. For example, trypsin 1 requires Lys and Arg at the P1 site , and granzyme B shows strict specificity for Asp at the P1 site . The specificity of protease has been widely used not only in identifying the biologically relevant substrates, but also in applying protease to site-specific proteolysis [6, 7]. Proteases participate in various disease processes, exhibiting a potentially huge future application in the design of new drug targets for enzyme [8, 9] and protease inhibitors . Although all the proteases function in hydrolyzing peptide bonds, almost all are linked to a particular cleavage pattern .
The MEROPS database is a manually curated information resource for peptidases . According to MEROPS, more than 10,000 known substrates are profiled for some proteases , so it is necessary to develop an approach to map the abundant substrate-sequence information to specificities of proteases to highlight the enzymatic preferences, especially for specific catalytic types . Integrating features of substrate sequences characteristics, PoPS  and PROSPER  are proposed to predict protease substrate cleavage sites. A well-designed approach of identifying the specificity of the protease will contribute to a better method of predicting the substrate cleavage site.
Previous analyses of protease cleavage data, such as visualized sequence logos , iceLogo , heat maps  and several techniques [20–22], have been focused on qualitative interpretation. Using LC-MS/MS sequencing , a simple and rapid multiplex substrate profiling method was presented to demonstrate the substrate specificity. Further measures include using fluorogenic substrates , specific labeling techniques of N-terminal [24, 25], and proteome-derived peptide libraries [26–30]. Fuchs  developed a method to quantify protease specificity and rank proteases with the cleavage entropy of a single position. Several quantitative measures were developed [32–35], in which the specificities of proteases were shown by the occurrences of amino acids at the binding sites. As mentioned by Schilling and Overall  in profiling the specificity of the MMP2, the preferred amino acid residues at different site may cooperate in the hydrolysis process, therefore, it is critically important to elucidate the hydrolyzation process by the closely cooperative relationship of successive positions on the substrate sequences.
In this study, we designed a novel approach to present the protease specificity based on blocks which are composed by successive amino acids from the substrate sequence. The essential difference between our approach and previous ones lies in that we characterize the specificity of proteases based on successive amino acids rather than a single binding site. This new approach could more reliably identify protease specificity by considering cooperation among the successive sites of the substrate peptides during the hydrolyzation process.
The dataset is composed of 61 proteases for analysis as described by Fuchs . The cleavage information from all experimental sources is obtained from the MEROPS database  and is updated according to MEROPS 10.0. This study focuses on the protease specificity directly on the active sites, ignoring differences in allosteric sites and exosite interactions. Among the data, signal peptidase complex (XS26.001) has been deleted from the dataset since the complex contains two peptidases, and it is not possible to assign a particular cleavage to one activity due to not a single component.
Greedy algorithm for filtering the data
First, the substrate sequence with less than two amino acids is filtered out. Then all of substrate sequences left primarily are aligned pairwise. Starting from sequences with the maximum number of similar amino acids, we remove redundant sequences by greedy algorithm  to make sure that there is no pair of sequences whose similarity is greater than or equal to 0.875. Therefore, there are at least two different amino acid residues between any two remaining substrate sequences.
Construction of blocks
Calculation of entropy
Calculation of distance matrix
Where E i (P) and E i ’(P) are the entropies of blocks B i and B i ′ of protease P respectively. This yields a symmetric distance matrix. The elements on the diagonal are 0, which is the distance of identical proteases.
Principal components analysis
All the eight blocks for each of 61 proteases are used for principal components analysis (PCA). Kaiser-Meyer-Olkin (KMO) Measure of Sampling Adequacy is computed as 0.733 which indicates that the sample size is sufficient for the application of PCA. The PCA is performed in SPSS 19.0 (SPSS Inc., Chicago, IL, USA) with the correlation method and Varimax with Kaiser Normalization as the rotation method.
Fisher’s exact test
2 × 2 contingency table for Fisher’s exact test
a + b
c + d
a + c
b + d
a + b + c + d
Creation of sequence profile
To depict the substrate preferences at different sites, the data of substrate sequences after removing the redundancy is submitted to Weblogo [17, 40] to generate sequence profiles of substrate cleavage site.
Distance character of 61 proteases
The entropies of eight blocks B4, B3, B2, B1, B1’, B2’, B3’, B4’ are calculated and denoted as E4, E3, E2, E1, E1’, E2’, E3’, E4’ correspondingly (Additional file 2: Table S1). There are three blocks with entropy 0, including, caspase 6 with E1 = 0 implying the unique amino acid Asp at site P1; peptidyl-Lys metallopeptidase with E1’ = 0 implying the unique amino acid Lys at site P1’; lysyl peptidase with E1 = 0 implying the unique amino acid Lys at site P1.
Principal components analysis
The entropies of eight blocks reveal the complexity of different combination types. In order to mine the blocks which play the crucial role in the specificity recognition of substrate sequence, we used principal components analysis.
The distribution of eight different eigenvalues is shown in a scree plot (Additional file 2: Figure S1.). Three principal components (PC1: the first principal component; PC2: the second principal component; PC3: the third principal component) are obtained according to the principle of eigenvalues more than 1. Among the three principal components, PC1, PC2 and PC3 contribute 57.938%, 23.284% and 15.960% to the total variance respectively, and the cumulative contribution of three principal components is 97.183% (Additional file 2: Table S3). Thus, the three principal components may represent the main features in the recognition of substrate specificities of different proteases.
Block-based sequence profile
There are a few prominent blocks from prime side. For instance, except for strict specificity for Lys at the P1’ site, peptidyl-Lys metallopeptidase has block B2’ with LysGlu = 179 from 1869 substrates, and signal peptidase 1 has block B3’ with AlaGluAla = 19 from 297 substrates (The number behind the equal sign represents the amount of combination of amino acids in the corresponding block).
Meanwhile, a few blocks from non-prime side show the specificity. For example, kexin has block B2 LysArg = 147 from 171 substrates. With caspase 3 having 571 substrates, besides the prominent block B3 GluValAsp = 43, we still find the prominent block B4 AspGluValAsp = 19.
The top prominent B2 blocks of proteases listed in Fig. 6
(a) The top prominent B2 blocks of proteases listed in Fig. 6a
(b) The top prominent B2 blocks of proteases listed in Fig. 6b
Val, Glu, Ile
Ala, Ser, Gly
Ala, Gly, Asn
Ala, Gly, Pro
Gly, Ala, Pro
Some blocks B2’ show the similar combination property as in blocks B2. For example, the top three amino acids of HIV-1 retropepsin at sites P1’ and P2’ are Leu, Val, Phe and Glu, Val, Ala, respectively, yet the prominent block B2’ with the highest number of combination are LeuAla = 33. For LAST_MAM peptidase, amino acids on the top at sites P1’ and P2’ are Asp, Ala, Glu and Pro, Ala, Glu respectively, yet the top one Asp at the site P1’ shows no preference of the top amino acid Pro at the site P2’, and the prominent block B2’ with the highest number of combination is AlaPro = 44 from 429 substrates. For the proteases which cleavage sites possess two or more preferred residues, the prominent combinations in the blocks reflect the cooperation of the residues in one position with other positions, characterizing the specificity of proteases detailedly.
Some specificities of certain proteases have been determined, such as trypsin 1 , caspase 3 , kexin , furin  and so on. However, by focusing on single positions and not taking into consideration the interaction of adjacent amino acids, the study of substrate specificity is too limited.
Taking the cooperation of amino acids into account, we propose a quantitative method to characterize substrates specificity of different proteases. By calculating entropies of different blocks, some distinctions of substrates between different proteases can be conceived (Fig. 2a). The principal component analysis gives evidence on the existence of blocks which play the crucial role in the specificity recognition of substrate sequence, and most of them are block B2, B1 and B1’. This is confirmed by the statistical analysis showing the ratios at B2, B1, and B1’ are higher than those in other blocks (Fig. 5). With Fisher’s exact test, a number of prominent blocks of different proteases have been discovered. For example, blocks B2 in kexin and furin are consistent with the previous discovery that both of the proteases cleave after dibasic residues . Other block B2, e.g. GluLeu in HIV-1 retropepsin, AlaAla in MMP2 and ProGly in MMP 9, are more likely to reflect the preferences and the cooperation of the successive amino acids in the substrate sequences which could not be found previously.
Cathepsin B is an endopeptidase and as an exopeptidase acts as a peptidyl-dipeptidase, releasing a dipeptide from the C-terminus of a protein or peptide. As no distinction is made in MEROPS between cleavages resulting from either activity, a view of the endopeptidase activity would be clear if the substrates of the exopeptidase activity were filtered out.
From the specificity matrix in MEROPS and the heat map , the preference of the protease is shown by the amino acids at one single binding site. However, it will not show the combinations of amino acids if proteases show multiple preferences at each binding site. Our method indicates interactions of different compositions of successive amino acids which can’t be obtained previously. For example, MMP9 has preferences for Ala, Gly and Pro at the site P2, Gly, Ala and Pro at the site P1 from the specificity matrix, yet the combination is clear using our method, such as ProGly, AlaAla in block B2. Whether a prominent combination exists in a block is obviously presented in the heat map of prominent combinations in each block (Fig. 4). These findings of specific blocks will shed light on future experiments and further investigation of proteolytic specificity.
Although in this study we only focused on the specificity of selected proteases, the method would be applicable to other proteases for mining the specificity pattern of substrates. In conclusion, we can obtain more substrate specificity patterns by site cooperation as more and more substrates data becomes available. Further investigations of the substrate specificity will be important to reveal the hydrolyzation mechanism of proteases.
Generally, the design of experiments and the description of the specificity of the protease are based on the assumption that the process of binding amino acid residue to the corresponding subsite is independent. However, it is not exactly true and the binding of amino acid residues at one site can more or less influence the binding at other subsites. It is essential to take the site cooperation into consideration for understanding fully the active site.
Our approach provides a new framework for dealing with the specificity pattern of substrates of the proteases. The combinations of site cooperation in the substrates offer a new sight in mining the specificity of the protease. We successfully find the significant blocks B2 in kexin and furin which are consistent with the previous discovery that both of the proteases cleave after dibasic residues. Other significant combinations found by the new approach could be more reliable to capture the activity of the active site. In principle, this method is useful for the further research relying on the substrate dataset, such as the identification of the novel substrate and the design of the inhibitor for the protease.
The authors thank the editorial staff for their help in editing this manuscript and thank the anonymous reviewers for their suggestions and comments to improve the manuscript.
This work was supported by National Science Foundation (NSF Grant No.1553680), and National Science Foundation of China (NSFC Grant No. 61432010, 61,272,016 and 31,571,354).
Availability of data and materials
The datasets used in this research are available at http://www.merops.ac.uk.
E.Q and G.L conceived and designed the approach. E.Q, B.G and Y.L implemented the software. D.W and Y.L performed the data analysis. E.Q wrote the manuscript. G.L contributed to revise the manuscript. All author approved the final version of this manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Turk B, Turk D, Turk V. Protease signalling: the cutting edge. EMBO J. 2012;31(7):1630–43.View ArticlePubMedPubMed CentralGoogle Scholar
- López-Otín C, Bond JS. Proteases: multifunctional enzymes in life and disease. J Biol Chem. 2008;283(45):30433–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Schechter I, Berger A. On the size of the active site in proteases. I Papain Biochem Bioph Res Co. 1967;27(2):157–62.View ArticleGoogle Scholar
- Harris JL, Backes BJ, Leonetti F, Mahrus S, Ellman JA, Craik CS. Rapid and general profiling of protease specificity by using combinatorial fluorogenic substrate libraries. Proc Natl Acad Sci U S A. 2000;97(14):7754–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Waugh SM, Harris JL, Fletterick R, Craik CS. The Structure of the Pro-Apoptotic Protease Granzyme B Reveals the Molecular Determinants of its Specificity. Nat Struct Biol. 2000;7(9):762–5.View ArticlePubMedGoogle Scholar
- Denning DW, Anderson MJ, Turner G, Latgé JP, Bennett JW. Sequencing the Aspergillus fumigatus genome. Lancet Infect Dis. 2002;2(4):251–3.View ArticlePubMedGoogle Scholar
- López-Otín C, Overall CM. Protease degradomics: a new challenge for proteomics. Nat Rev Mol Cell Bio. 2002;3(7):509–19.View ArticleGoogle Scholar
- Turk B. Targeting proteases: successes, failures and future prospects. Nat Rev Drug Discov. 2006;5(9):785–99.View ArticlePubMedGoogle Scholar
- Lopez-Otin C, Matrisian LM. Emerging roles of proteases in tumour suppression. Nat Rev Cancer. 2007;7(10):800–8.View ArticlePubMedGoogle Scholar
- Liu H, Shi X, Guo D, Zhao Z, Yimin. Feature Selection Combined with Neural Network Structure Optimization for HIV-1 Protease Cleavage Site Prediction. Biomed Res Int. 2015;2015:263586.PubMedPubMed CentralGoogle Scholar
- Hedstrom L.Introduction: proteases. 2002;102(12):4429.Google Scholar
- Rawlings ND, Barrett AJ, Finn R. Twenty years of the MEROPS database of proteolytic enzymes, their substrates and inhibitors. Nucleic Acids Res Nucleic Acids Res. 2016;44(D1):D343–50.View ArticlePubMedGoogle Scholar
- Rawlings ND. Peptidase specificity from the substrate cleavage collection in the MEROPS database and a tool to measure cleavage site conservation. Biochimie. 2016;122:5–30.View ArticlePubMedPubMed CentralGoogle Scholar
- Song J, Tan H, Boyd SE, Shen H, Mahmood K, Webb GI, Akutsu T, Whisstock JC, Pike RN. Bioinformatic approaches for predicting substrates of proteases. J Bioinforma Comput Biol. 2011;9(1):149–78.View ArticleGoogle Scholar
- Boyd SE, Pike RN, Rudy GB, Whisstock JC, Garcia de la Banda M. PoPS: a computational tool for modeling and predicting protease specificity. J Bioinforma Comput Biol. 2005;3(3):551–85.View ArticleGoogle Scholar
- Song J, Tan H, Perry AJ, Akutsu T, Webb GI, Whisstock JC, Pike RNPROSPER. an integrated feature-based tool for predicting protease substrate cleavage sites. PLoS One. 2012;7(11):e50300.View ArticlePubMedPubMed CentralGoogle Scholar
- Schneider TD, Stephens RM. Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 1990;18(20):6097–100.View ArticlePubMedPubMed CentralGoogle Scholar
- Colaert N, Helsens K, Martens L, Vandekerckhove J, Gevaert K. Improved visualization of protein consensus sequences by iceLogo. Nat Methods. 2009;6(11):786–7.View ArticlePubMedGoogle Scholar
- Schilling O, Overall CM. database-searchable peptide libraries for identifying protease cleavage sites. Nat Biotechnol. 2008;26(6):685–94.View ArticlePubMedGoogle Scholar
- Poreba M, Drag M. Current strategies for probing substrate specificity of proteases. Curr Med Chem. 2010;17(33):3968–95.View ArticlePubMedGoogle Scholar
- Huesgen PF, Overall CM. N- and C-terminal degradomics: new approaches to reveal biological roles for plant proteases from substrate identification. Physiol Plant. 2012;145(1):5–17.View ArticlePubMedGoogle Scholar
- Van Damme P, Staes A, Bronsoms S, Helsens K, Colaert N, Timmerman E, Aviles FX, Vandekerckhove J, Gevaert K. Complementary positional proteomics for screening substrates of endo- and exoproteases. Nat Methods. 2010;7(7):512–5.View ArticlePubMedGoogle Scholar
- O'Donoghue AJ, Eroy-Reveles AA, Knudsen GM, Ingram J, Zhou M, Statnekov JB, Greninger AL, Hostetter DR, Qu G, Maltby DA, Anderson MO, Derisi JL, McKerrow JH, Burlingame AL, Craik CS. Nat Methods 2012;9(11):1095–100.Google Scholar
- Mahrus S, Trinidad JC, Barkan DT, Sali A, Burlingame AL, Wells JA. Global sequencing of proteolytic cleavage sites in apoptosis by specific labeling of protein N termini. Cell. 2008;134(5):866–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Kleifeld O, Doucet A, Prudova A, Auf dem Keller U, Gioia M, Kizhakkedathu JN, Overall CM. Identifying and quantifying proteolytic events and the natural N terminome by terminal amine isotopic labeling of substrates. Nat Protoc. 2011;6(10):1578–611.View ArticlePubMedGoogle Scholar
- Boulware KT, Daugherty PS. Protease specificity determination by using cellular libraries of peptide substrates (CLiPS). Proc Natl Acad Sci U S A. 2006;103(20):7583–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Turk BE, Huang LL, Piro ET, Cantley LC. Determination of protease cleavage site motifs using mixture-based oriented peptide libraries. Nat Biotechnol. 2001;19(7):661–7.View ArticlePubMedGoogle Scholar
- Schilling O, Huesgen PF, Barré O, Auf dem Keller U, Overall CM. Characterization of the prime and non-prime active site specificities of proteases by proteome-derived peptide libraries and tandem mass spectrometry. Nat Protoc. 2011;6(1):111–20.View ArticlePubMedGoogle Scholar
- Wang C, Ye M, Bian Y, Liu F, Cheng K, Dong M, Dong J, Zou H. Determination of CK2 specificity and substrates by proteome-derived peptide libraries. J Proteome Res. 2013;12(8):3813–21.View ArticlePubMedGoogle Scholar
- Tucher J, Linke D, Koudelka T, Cassidy L, Tredup C, Wichert R, Pietrzik C, Becker-Pauly C, Tholey A. LC-MS based cleavage site profiling of the proteases ADAM10 and ADAM17 using proteome-derived peptide libraries. J Proteome Res. 2014;13(4):2205–14.View ArticlePubMedGoogle Scholar
- Fuchs JE, von Grafenstein S, Huber RG, Margreiter MA, Spitzer GM, Wallnoefer HG, Liedl KR. Cleavage entropy as quantitative measure of protease specificity. PLoS Comput Biol. 2013;9(4):e1003007.View ArticlePubMedPubMed CentralGoogle Scholar
- Julien O, Zhuang M, Wiita AP, O'Donoghue AJ, Knudsen GM, Craik CS, Wells JA. Quantitative MS-based enzymology of caspases reveals distinct protein substrate specificities, hierarchies, and cellular roles. Proc Natl Acad Sci U S A. 2016;113(14):E2001–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Schauperl M, Fuchs JE, Waldner BJ, Huber RG, Kramer C, Liedl KR. Characterizing protease specificity: how many substrates do we need? PLoS One. 2015;10(11):e0142658.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu J, Duan X, Sun J, Yin Y, Li G, Wang L, Liu B. Bi-factor analysis based on noise-reduction (BIFANR): a new algorithm for detecting coevolving amino acid sites in proteins. PLoS One. 2013;8(11):e79764.View ArticlePubMedPubMed CentralGoogle Scholar
- Fuchs JE, von Grafenstein S, Huber RG, Kramer C, Liedl KR. Substrate-driven mapping of the degradome by comparison of sequence logos. PLoS Comput Biol. 2013;9(11):e1003353.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang Z, Schwartz S, Wagner L, Miller WA. greedy algorithm for aligning DNA sequences. J Comput Biol. 2000;7(1–2):203–14.View ArticlePubMedGoogle Scholar
- Shannon CEA. mathematical theory of communication. Bell Syst Tech J. 1948;27(3):379–423.View ArticleGoogle Scholar
- Fisher RA. On the Interpretation of χ2 from Contingency Tables, and the Calculation of P. J R Stat Soc. 1922;85(1):87–94.View ArticleGoogle Scholar
- Miller RG. Simultaneous statistical inference. 2nd ed. New York: Springer; 1981.View ArticleGoogle Scholar
- Crooks GE, Hon G, Chandonia JM, Brenner SE. WebLogo: a sequence logo generator. Genome Res. 2004;14(6):1188–90.View ArticlePubMedPubMed CentralGoogle Scholar
- Oliveira V, Campos M, Melo RL, Ferro ES, Camargo AC, Juliano MA, Juliano L. Substrate specificity characterization of recombinant metallo oligo-peptidases thimet oligopeptidase and neurolysin. Biochemistry. 2001;40(14):4417–25.View ArticlePubMedGoogle Scholar
- Demon D, Van Damme P, Vanden Berghe T, Deceuninck A, Van Durme J, Verspurten J, Helsens K, Impens F, Wejda M, Schymkowitz J, Rousseau F, Madder A, Vandekerckhove J, Declercq W, Gevaert K, Vandenabeele P. Proteome-wide substrate analysis indicates substrate exclusion as a mechanism to generate caspase-7 versus caspase-3 specificity. Mol Cell Proteomics. 2009;8(12):2700–14.View ArticlePubMedPubMed CentralGoogle Scholar
- Bader O, Krauke Y, Hube B. Processing of predicted substrates of fungal Kex2 proteinases from Candida albicans, C. glabrata, Saccharomyces cerevisiae and Pichia pastoris. BMC Microbiol. 2008;8:116.View ArticlePubMedPubMed CentralGoogle Scholar
- Remacle AG, Shiryaev SA, ES O, Cieplak P, Srinivasan A, Wei G, Liddington RC, Ratnikov BI, Parent A, Desjardins R, Day R, Smith JW, Lebl M, Strongin AY. Substrate cleavage analysis of furin and related proprotein convertases, A comparative study. J Biol Chem. 2008;283(30):20897–906.View ArticlePubMedPubMed CentralGoogle Scholar
- Page MJ, Di Cera E. Serine peptidases: classification, structure and function. Cell Mol Life Sci. 2008;65(7–8):1220–36.View ArticlePubMedGoogle Scholar