BACKGROUND
caspases are a family of intracellular cysteinyl aspartate-specific proteases that are highly conserved in multicellular organisms and are key regulators of apoptosis initiation and execution.
at least  <dig> members of the caspase family have been identified in mammals and they are grouped into two major sub-families, namely inflammatory caspases and apoptotic caspases.
apoptosis associated caspases can be further classified into two groups: initiator caspases, including caspase- <dig>  - <dig>  and - <dig>  which are present upstream of apoptosis signalling pathways; and executioner  caspases - <dig>  - <dig>  and - <dig>  <cit> .
initiator capsases- <dig> and - <dig> are activated through an auto-cleavage process that is mediated by large adaptor-caspase complexes known respectively as the death-inducing signalling complex  and apoptosome.
these complexes are usually formed in response to an intrinsic or extrinsic cell death stimulus  <cit> .
the main targets of the activated initiator caspases are the executioner procaspases.
it is interesting to notice that substrates of initiator caspases are limited to their own precursors, executioner procaspases- <dig>  - <dig>  and - <dig>  and few more proteins  <cit> .
on the other hand, executioner caspases target a large number of cellular proteins to control the dismantling process of the cell  <cit> .
in addition to their essential role in apoptosis, recent accumulated evidence demonstrates various non-apoptotic functions of executioner caspases including: regulation of the immune response, cell proliferation, differentiation and motility  <cit> .
caspases are characterized by high substrate selectivity.
they recognize a specific sequence signal in their target proteins.
resolving substrate specificity for caspases was initially investigated using a combinatorial approach with positional scanning of synthetic tetrapeptidyl-aminomethyl coumarin derivatives.
the results of this approach determined the absolute requirements for aspartic acid at position p <dig>  <cit> .
in addition, p <dig> to p <dig> positions demonstrate high preference for certain amino acids.
based on positional scanning of synthetic tetrapeptides, the preferred recognition sequences for caspases - <dig>  - <dig>  and - <dig> were determined to be ehd, whereas caspases - <dig>  and - <dig> recognize the sequence devd, while caspases-  <dig>  - <dig>  - <dig>  and - <dig> recognize the sequence ed.
it is important to emphasize that the in vitro caspase substrate specificity, determined by the synthetic tetrapeptide method, is not absolutely representative of the cleavage conditions in vivo.
the cleavage specificities of caspases in vivo are influenced by sequence-dependent conformational features, flanking the cleavage tetrapeptide motif, which can control the molecular electrostatic potential and the steric accessibility of the enzyme to its target protein.
for example, and in spite of their identical preference for the devd tetrapeptide cleavage motif, caspase- <dig> and  <dig> show a clear differential preference for various natural substrates  <cit> .
demon et al  <cit>  demonstrated that in addition to the tetrapeptide cleavage core, devd , several amino acid positions located outside this core such as p <dig>  p <dig>  p2' and p3' are critical in the discrimination of caspase- <dig> and caspase- <dig> for their specific substrates.
shen et al  <cit>  reported another interesting example which proves that the relatively similar tetrapeptide cleavage motif of caspase- <dig> and caspase- <dig>  which are functionally distinct, does not imply a similar recognition preference for their natural substrates.
via a thorough statistical analysis of a window size of p10-p10' for a collection of caspase- <dig> and caspase- <dig> natural substrates, shen et al have determined the significance of various amino acids and/or certain physiochemical properties at certain positions outside the canonical tetrapeptide motif  <cit> .
among executioner caspases, caspase- <dig> is considered the major enzyme with a wide array of cellular substrates.
while immunodepletion of caspase- <dig> abolishes the majority of proteolytic events observed during apoptosis, immunodepletion of other executioner caspases shows a minimal impact on apoptosis markers and its proteolytic cleavage outcomes  <cit> .
in the last decade, extensive research on caspases led to the identification of more than  <dig> caspase- <dig> substrates and the list is still growing.
with the increasing number of proteins that have been discovered, thanks to the sequencing of the human genome and the genomes of many other organisms, there is a need for efficient methods that can help in discovering new caspase- <dig> substrates.
the identification of new cellular substrates for caspase- <dig> would lead to further insights into the cellular mechanisms that regulate apoptosis, proliferation, and other biological processes.
bioinformatics tools would allow high-throughput analyses of proteomic data in order to screen for putative caspase- <dig> substrates.
in addition, such tools can provide researchers with an accurate map of the potential cleavage site for a given sequence of interest.
in the last few years, several computer-based tools were developed with the aim of predicting caspase substrates.
prediction of endopeptidase substrates   <cit>  was among the initial tools and it was developed in order to predict putative caspase- <dig>  cathepsin b and cathepsin l cleavage sites using cleavage site scoring matrices .
peptidecutter  <cit>  is another tool that was designed with the objective of predicting cleavage sites for a wide range of proteases including various caspases.
grabcas  <cit>  is a tool that uses a position specific scoring matrix for caspases 1- <dig> and granzyme b, based on substrate specificities that were determined by positional scanning of synthetic peptides.
caspredictor  <cit>  was developed based on the assumption that sequences rich in the amino acids ser , thr , pro , glu or asp   are favoured caspase cleavage sites.
caspredictor was built based on  <dig> experimentally verified natural substrates for caspase- <dig>  - <dig>  - <dig>  - <dig>  - <dig>  - <dig>  - <dig>  and - <dig>
in addition to the previously aforementioned scoring-matrix based approaches, several groups recently reported the development of tools that were mainly built up using the support vector machine  technique.
wee et al.
<cit>  described an svm based approach  using  <dig> substrates for different caspases from various organisms.
cascleave  <cit>  is an interesting tool that was recently developed utilizing primary sequence, as well as secondary structure features, of the cleaved sites based on svm approach.
cascleave was built using a dataset of  <dig> substrates of the different caspases.
piippo et al  <cit>  described another tool  using three different pattern recognition classifiers, namely: svm, a decision tree based method known as j <dig>  and the random forest classifier.
the three classifiers were trained on  <dig> different caspase cleavage sites.
li et al  <cit>  proposed a hybrid svm-pssm method based on an extended dataset.
unfortunately, some of these tools are not available for testing and comparison purposes.
despite the substantial efforts to develop in silico systems to predict sites of caspase cleavage, the accuracy of such tools is still a challenging issue.
the major drawback of the early tools is the use of training datasets that represent synthetic peptides or limited natural substrates for various proteases including caspases.
on the other hand, the recently developed svm-based tools were built using a mixture of heterogeneous data that represent cleavage sites of the different caspases including non-canonical sites as well as some unverified caspase substrates.
in general the svm-based tools such as cascleave achieve good levels of sensitivity, yet they suffer from high rates of false positive results.
it is generally expected that training a prediction tool on data representing distinctive patterns can lead to overgeneralization and hence a high rate of false positive results.
it is important to recall that although different caspases share the primary sequence-requirement to cleave at the carboxyl terminal of aspartate residues in their protein targets, each one of these proteases recognizes a unique context surrounding the cleavage position.
even the caspases that appear to have identical tetrapeptide cleavage specificities such as caspase- <dig> and - <dig> are actually distinct in terms of the amino acids preferences outside the tetrapeptide core sequence  <cit> .
based on this assumption, we decided to develop a prediction tool focusing on data that represent substrates of a single caspase.
caspase- <dig> was selected for this objective as it represents the major executioner caspase with a considerable number of substrates.
in this work, we present a novel tool designated caspase analysis tool  <dig> , which was developed based on an extensive and highly curated dataset of caspase- <dig> substrates.
cat <dig> showed an obvious improvement in the overall prediction accuracy as well as a marked reduction of the false positive rate.
using cat <dig>  a high-throughput screening was performed on a large set of human proteins with assigned gene ontology  annotations.
the screening results reveal the existence of a large number of potential caspase- <dig> substrates.
implementation
methods
caspase- <dig> substrates
the pubmed literature database  <cit>  was used to search for papers that describe human, mouse, and rat caspase- <dig> substrates.
each paper was critically analyzed to determine the experimentally demonstrated cleavage position in the relevant protein.
the confirmed caspase- <dig> substrates were  <dig> proteins with a total of  <dig> cleavage sites.
the amino acid sequences of the proteins were then obtained from the universal protein resource knowledgebase   <cit> .
of the  <dig> cleavage sites,  <dig> sequences were sorted randomly to be used later to compare the performance of cat <dig> versus existing similar tools; the remaining  <dig> sequences, which we refer to as the positive  peptide data, were used for training and validation of the cat <dig>
definition of study controls
the following datasets were established as controls in this study:
the negative uncleaved peptides
this data set consists of all the peptides that contain aspartic acid residues and are presumed to be uncleaved.
this control group was established based on the assumption that any d residue, in a caspase- <dig> substrate, apart from the mapped cleavage site is most likely uncleaved.
after excluding the positive peptides that exist in training data, the remaining  <dig> d residues were used to create this negative  control.
amino acids natural frequency
this control represents the frequency of each one of the  <dig> amino acids in a group of  <dig>  human proteins that were available as reviewed proteins in the uniprotkb  <cit>  as of march  <dig>
amino acid r-groups frequency
this control represents the frequencies of the different r-groups of amino acids; acidic , basic , polar , and non-polar .
the frequencies were calculated based on the above mentioned  <dig>  reviewed proteins.
physiochemical characteristics flanking the cleavage site
the positive peptide sequences were aligned in reference to the cleaved aspartic acid residues.
the resulting multiple sequence alignment was divided into three regions: the central tetrapeptide cleavage motif , the n-terminal region preceding the motif and the c-terminal region following the motif that was designated "before-motif" and "after-motif", respectively.
the analyses for the regions flanking the motif were made serially:  <dig>   <dig>   <dig>   <dig>  and  <dig> amino acids before and after the motif .
the analysis included: the frequencies of amino acids represented by their r-groups , the frequencies of hydrophobic and hydrophilic amino acids and finally the frequencies for each single amino acid.
in the case of tetrapeptide motif analysis, the different frequencies were calculated separately for each position: p <dig>  p <dig>  p <dig> and p <dig>  however, the frequencies within  <dig>   <dig>   <dig>   <dig>  and  <dig> amino acids, before and after the motif, were calculated collectively for each region.
establishment of scoring matrices
the peptides that fulfil the length criterion, p9-p5', which means having  <dig> amino acids before and  <dig> amino acids after the aspartate residue of interest, were used to build the scoring matrices.
both the positive and the negative peptide data sets were used to build the scoring matrices.
the first step was to generate position specific frequency matrices from the multiple sequence alignments of the relevant set of peptides.
each matrix consisted of  <dig> rows, representing positions p9p8p7p6p <dig> p <dig> p <dig> p <dig> p1p1' p2' p3'p4'p5', where a d amino acid is at the position p <dig>  the  <dig> columns of the matrix represent the frequencies of each amino acid.
from the multiple sequence alignment of the positive peptides, we noticed the presence of two possibly different patterns; the first pattern has a d at p <dig>  and the second has any amino acid except d at p <dig> .
to represent this subtle difference, we decided to construct amino acid frequency matrices to represent each sub-pattern.
two weighting systems were used in order to correct the probability of overrepresented and underrepresented amino acids in the frequency matrices so as to establish the scoring matrices:
i) calculating log odd ratio: this weighting system involves calculating log odd ratio for each element in the frequency matrix by dividing the observed frequency of a given amino acid over its corresponding natural frequency .
ii) subtraction of negative control background: instead of relying only on the common log odd weighting system and in order to minimize scoring bias, we decided to add a second normalization approach.
the method relies on comparing the positive peptides with the negative peptides to further remove the noise signals around the cleaved aspartate residues.
four scoring matrices are involved in the overall calculation of the final score of cat <dig> tool.
we propose the following notation to define each scoring matrix and the overall score.
first, let fm1+ denote the frequency matrix that was constructed from all the positive  peptides.
the corresponding scoring matrix 'a' is defined as:
a=log2fm1+Ω
where Ω is the natural frequencies of the amino acids.
in addition to the above scoring matrix we define fm1- as the frequency matrix generated from the negative  peptides.
a new frequency matrix 'b' is defined as:
b=log2fm1+-fm1-Ω
let fmc denote a frequency matrix calculated from a subset of peptides that fulfil the constraint 'c'.
here,  is either + or - as explained before.
therefore, we define the following scoring matrices:
c1=log2fm1p4=d+-fm1p4=d-Ω
and
c2=log2fm1p4≠d+-fm1p4≠d-Ω
cat <dig> implementation and scoring
cat <dig> tool was built using perl language.
the input protein can be entered either as a fasta format sequence or as a text file.
once a p <dig> peptide with a d residue at p <dig> is identified, it is analyzed to calculate the final score 's' as follows:
s=a+b+c <dig>
where 'a' and 'b' are scores generated from the scoring matrices 'a' and 'b' in equation  <dig> and equation  <dig>  respectively.
the 'c' score is generated either from the scoring matrix 'c1'or 'c2' as follows:
c=c1ifp4=dc2ifp4≠d
we refer to the scoring matrix 'c1' if the peptide contains the amino acid d at p <dig> or the scoring matrix 'c2' if the amino acid at p <dig> is not d. the three scores  are normalized to a 100% score by dividing each score by the maximum score that could be obtained from each formula.
cat <dig> validation
to examine the prediction power of cat <dig> a k =  <dig> fold cross validation was performed.
the positive data were the actual cleavage sites, whereas the negative data were obtained from the uncleaved dataset.
in each fold four pssm matrices were created from 9/10th of the positive substrates.
then, the remaining 1/10th positive and negative substrates are used for testing.
since the number of the negative peptides was much larger than the positive peptides, an equal number of the negative peptides were randomly obtained.
the whole  <dig> fold cross validation experiment was repeated  <dig> times to ensure a good coverage of the negative dataset.
the sensitivity , specificity , positive predictive value , negative predictive value , accuracy  and the matthew's correlation coefficient  were calculated as in  <cit> .
the areas under the receiver operating characteristic  curves were calculated by plotting the sensitivity against the corresponding 1-specificity.
the optimal cut-off point was defined as that measurement that corresponded to the point on the roc curve closest to the top left corner, i.e., closest to having sensitivity =  <dig> and specificity =  <dig>
local window size
the most appropriate local window size of amino acid sequence encompassing the cleaved aspartate and other critical exosite residues was determined based on comparing the prediction performance of the following peptides: p3-p <dig>  p4-p <dig>  p5-p1', p6-p2', p7-p3', p8-p4', p9-p5', p11-p7', p14-p10', p17-p13', p20-p16', p23-p19'.
the performance of each window size was evaluated using the same cross validation approach as mentioned above.
the obtained area under curve  and matthew's correlation coefficient  for the different experiments were plotted for comparison purposes.
performance comparison
a performance comparison was carried out for cat <dig> versus two recently published prediction tools, namely casvm and cascleave  <cit> .
the aim of the test was to assess how accurate the three tools were in predicting caspase- <dig> cleavage sites.
the comparison was made on  <dig> caspase- <dig> substrates that were randomly excluded from the training dataset.
since cat <dig> is a prediction tool specific for caspase- <dig> cleavage sites, whereas casvm and cascleave were developed to predict cleavage sites of different caspases, there is a possibility to misjudge true positive sites of other caspases by assigning them to the false positive category of casvm and cascleave.
to avoid such unfair comparison, the  <dig> substrates were carefully inspected to find all caspase cleavage sites.
the search was performed using the pubmed database, google searching engine, the caspase substrate database homepage   <cit> , and merpos - the peptidases database  <cit> .
the protein sequences of the  <dig> substrates were analyzed individually and the prediction results for each tool were counted according to software default parameters.
the true positives are the positively predicted caspase- <dig> cleavage sites, whereas the false positives are the positively predicted aspartates that are actually not recognized by any caspase.
high-throughput screening
the uniprotkb  <cit>  was used to retrieve human proteins with known biological processes.
two filters of the advance search option were used: the first was organism: homo sapiens  and the second was gene ontology go: biological process .
after excluding the experimentally verified  <dig> human caspase- <dig> substrates, a total of  <dig> reviewed human proteins with defined gene ontology  were obtained.
the protein sequences were analyzed by cat <dig> to screen for potential novel caspase- <dig> substrates.
only results of scores ≥  <dig> were considered for further analyses.
proteins that were predicted as potential caspase- <dig> substrates were further analyzed using toppgene suite tool  <cit>  to retrieve the most significant gene ontology  terms.
RESULTS
caspase  <dig> substrates
our search in the pubmed literature database for caspase- <dig> substrates revealed the presence of  <dig> proteins:  <dig> of human origin,  <dig> of mouse origin, and  <dig> of rat origin.
all the substrates were experimentally verified as natural substrates and their cleavage sites were mapped.
of the  <dig> substrates, the cleavage sites of  <dig> proteins were mapped by site-directed mutagenesis technique, while the remaining  <dig> were mapped by different high-throughput proteomic screening approaches.
the full list and description of the obtained substrates are available in the additional materials .
the obtained caspase- <dig> substrates as well as other caspase substrates will be available in the caspase substrates comprehensive database  that has been developed at our institute .
tetrapeptide cleavage motif analysis
the tetrapeptide cleavage motifs  of the training group were analyzed to determine physiochemical properties and frequencies of amino acids at each position.
the examination of amino acid frequencies within the tetrapeptide motif revealed a unique distribution pattern of hydrophobic and hydrophilic amino acids .
hydrophilic amino acids are  <dig>  times more frequent in p <dig> than hydrophobic amino acids.
interestingly, p <dig> has an opposite pattern to p <dig>  in p <dig> hydrophilic amino acids are nearly two times more frequent than hydrophobic amino acids, whereas in p <dig> the converse is true.
features surrounding the cleavage site
the amino acid sequences surrounding the tetrapeptide cleavage motifs were thoroughly analyzed to identify necessary feature for caspase- <dig> recognition.
the analyses include: secondary structure, amino acids physicochemical properties and amino acid composition.
the secondary structure prediction method gor <dig>  <cit>  was used to investigate the cleavage motif and its flanking regions for any common secondary structure.
the analysis of gor <dig> results showed that the majority  of the cleaved sites are located within unstructured context, while 18% are located within alpha helical regions, and only 2% are located in beta sheets.
we then analyzed the biochemical properties of amino acids that flank the tetrapeptide cleavage motif to determine amino acids preferences for caspase- <dig> substrate recognition.
no significant differences in the frequencies of acidic, basic, polar and non-polar amino acids between the tested region and the corresponding control group were found when examining  <dig>   <dig>   <dig>  and  <dig> amino acids before and after the tetrapeptide cleavage motif.
when testing the region of  <dig> amino acids before the cleavage motif a slightly higher percentage of acidic amino acids were noticed, while basic and polar amino acids were strongly unfavoured.
in the region of  <dig> amino acids after the cleavage motif, lower percentages of acidic and basic amino acids were noticed .
to further explore the characteristic biochemical properties, we examined individual amino acid frequencies for the entire cleavage vicinity: the 5-amino acids before and after the tetrapeptide motif.
as shown in figure  <dig>  the frequencies of glycine, alanine, serine and proline have altered distributions in regions before and after the tetrapeptide motif that may indicate size and charge requirements for caspase- <dig> recognition and binding to the substrates.
position specific scoring matrices
in order to determine the most appropriate window size to construct efficient scoring matrices for cat <dig>  a series of gradually increasing window sizes ranging from p3-p <dig> to p23-p19' were evaluated .
as can be seen, it is obvious that window sizes equal to or shorter than the tetrapeptide cleavage motif  are not adequate to develop a reliable prediction tool.
despite of a marginally higher mcc at the window size of p6-p2', the overall prediction efficiency of the window sizes ranging from p5-p1' to p9-p5' are apparently quite similar.
however, the efficiency seems to decrease gradually when extending the window size beyond p9-p5'.
we actually preferred the window size p9-p5' over other seemingly comparable shorter alternatives for several reasons.
first, the critical analysis of amino acids over-/under-representation scores demonstrated the significance of all the positions in this extended window p9-p5' .
second, a careful analysis of natural caspase substrates, available in merops database, with cleavage positions near to n- or c- terminals, indicates that minimal adequate n- and c- terminal spacers comparable to the length of p <dig> and p5', respectively, are required for efficient recognition.
therefore, our scoring matrices were developed by calculating the weight of each amino acid in the 14-mer peptide sequence from p <dig> to p5'.
to evaluate the contribution of the different amino acids at the positions surrounding the cleaved aspartate, the scoring matrix 'a'  was drawn as a heat-map .
analyzing the heat-map shows that apart from the tetrapeptide cleavage motif, all positions have either overrepresented or underrepresented amino acids.
however, the positions including p <dig>  p <dig>  p1', p3', p4' and p5' have a remarkable rejection or preference for certain amino acids.
prediction power of cat3
a  <dig> fold cross validation was used to evaluate the predictive power of cat <dig>  figure  <dig> shows the roc curve that represents the average of  <dig> different experiments of the 10-fold cross validation.
the optimal cut-off score was found to be  <dig>  at this cut-off point the prediction statistical measures are shown in table  <dig>
the values of the statistical measures represent the average of  <dig> different experiments of the  <dig> fold cross validation test.
the optimal cut-off score was found to be  <dig>  auc: area under roc curve; spe: specificity; sen: sensitivity; ppv: positive predictive value; npv: negative predictive value; acc: accuracy; mcc: matthew's correlation coefficient.
to demonstrate specificity of cat <dig> for caspase- <dig> cleavage sites, a group of  <dig> non-caspase- <dig> substrates were examined by cat <dig>  the substrates included  <dig>   <dig> and  <dig> cleavage sites of caspase- <dig>  caspase- <dig> and caspase- <dig>  respectively.
we avoided using any cleavage site that was known to be a shared target with caspase- <dig>  interestingly,  <dig> of the  <dig> cleavage sites  showed cat <dig> scores below  <dig>  which is the minimum cut-off for predicting a caspase- <dig> cleavage site.
this result provides a clear evidence to substantiate the very high specificity of cat <dig> for predicting caspase- <dig> substrates.
the detailed results of this experiment are available in the additional materials .
evaluating the performance of different binary classifiers is frequently made by comparing their reported statistical measures such as specificity, sensitivity etc.
which are usually calculated under different conditions.
we have avoided the use of such a comparison as it can lead to a biased conclusion.
instead, we compared the performance of cat <dig> versus two recently reported tools, namely casvm and cascleave, on a group of  <dig> caspase- <dig> substrates that were initially excluded from our training data.
it is worth mentioning that some of these substrates could have been used in the training of the other tools, which could offer unfair advantage to the other two tools versus cat <dig>  a thorough examination using different databases revealed that the  <dig> substrates contain a total of  <dig> aspartate residues, of which  <dig> are caspase- <dig> cleavage sites,  <dig> are cleavage sites of other caspases and  <dig> aspartate residues that are evidently not cleaved by any caspase.
out of the  <dig> actual caspase- <dig> cleavage sites, the predicted true positive results for the three tools were as follows: cat <dig> 14/ <dig> , casvm 8/ <dig> , and cascleave 16/ <dig> .
however, cat <dig> was the best tool in reducing the false positive results ; out of  <dig> actual uncleaved aspartic acid containing peptide, cat <dig> predicted only  <dig> false positives , while casvm predicted  <dig> false positive  and cascleave had  <dig>  false positives.
figure  <dig> shows the result of comparing cat <dig> versus casvm and cascleave.
it is noteworthy that both casvm and cascleave correctly predicted two of the  <dig> non-caspase- <dig> cleavage sites.
the detailed results of the comparison are available in the additional materials .
high-throughput screening for novel caspase- <dig> substrates
screening of  <dig> reviewed human proteins with ascribed gene ontology  using our cat <dig> tool showed that  <dig> proteins are predicted to be caspase- <dig> substrates with a total of  <dig> potential caspase- <dig> cleavage sites .
to further investigate the function of these potential substrates we used toppfun: an annotations based gene list functional enrichment analysis tool  <cit> .
out of the  <dig> genes only  <dig> had annotations in toppfun.
the analysis revealed a group of  <dig> biological processes that showed significant enrichment in predicted proteins versus the whole human genome as a control .
a careful analysis of these biological processes was performed to shortlist the most significant biological processes by excluding general roots  and detailed leaves  of go terms.
the most significant biological processes are shown in table  <dig>
the go terms of the listed biological processes were manually filtered, to reduce redundancy, by removing general roots  and detailed leaves  of the enriched go terms that were obtained by toppfun tool.
the p-value of each go term in the predicted caspase- <dig> proteins was derived by random sampling of the whole genome.
toppfun was also used to examine the enriched cellular component go terms.
interestingly, the majority of the predicted proteins are located in different nuclear components, cytoskeleton, cell projection, membrane fraction, cell junction, and extracellular matrix, where most typical apoptotic morphological and biochemical changes are observed.
the detailed list of the enriched cellular component go terms is available in the additional materials .
discussion
in addition to its well known key function in apoptosis, caspase- <dig> has been shown to play a crucial role in the regulation of various biological processes such as cell differentiation, adhesion, neurodevelopment and neuronal signalling  <cit> .
recognition of caspase- <dig> substrates is becoming a vital need to understand molecular mechanisms behind many disorders including cancer, autoimmune and neurodegenerative diseases.
currently, most of the known caspase- <dig> substrates have been identified using in vitro proteolytic cleavage assays, coupled with site-directed mutagenesis to determine the exact cleavage position.
in recognition of the physiological importance of caspase- <dig>  many labs began to perform high-throughput proteomic screening to identify novel substrates of this major caspase  <cit> .
such techniques are relatively expensive and cumbersome.
in addition to the high cost, the number of identified substrates is usually limited to the proteins that are relatively abundant in the examined cell type.
recently, several computer-based prediction tools such as casvm, cascleave, and pripper were developed in order to help discover novel caspase substrates  <cit> .
these tools were trained on data that represent substrates of different caspases and in some cases non-caspase endopeptidase.
although different caspases share a primary sequence-requirement, to cleave at the carboxyl terminal of aspartate residues in their protein targets, each one of these proteases needs a special context surrounding the cleavage position.
even the caspases that appear to have identical tetrapeptide cleavage specificities such as caspase- <dig> and - <dig> are actually distinct in terms of the amino acids preferences outside the tetrapeptide core sequence  <cit> .
therefore, we believe that building a single algorithm for predicting the cleavage of multi-caspases would likely have low prediction specificity.
based on this hypothesis, we decided to develop an algorithm focusing only on substrates of caspase- <dig>  which is the major executioner caspase with a considerable number of targets.
our caspase-3-specific approach  has indeed outperformed other multi-caspases prediction tools on an independent comparison-dataset .
cat <dig> has three distinctive features.
firstly, it was developed using pssm instead of other relatively complex approaches.
pssm is known to be practical, require low computation power and is able to represent the statistical weights of amino acids at each position.
in addition, it can be easily combined with other machine learning tools to generate hybrid approaches that might enhance the prediction performance.
secondly, instead of using data for different caspases, which are actually a mixture of heterogeneous patterns, we used an extended set of highly-curated caspase- <dig> natural substrates.
we believe that inclusion of data that represents other proteases or cleavage sites of caspases with very few substrates and/or cleavage positions representing non-canonical patterns can lead to overgeneralization.
in this situation, the classification model is required to loosen its decision boundary to increase sensitivity, but at the cost of having more false positive results.
it is therefore generally accepted that improvement in prediction accuracy is more likely to be associated with the good quality of the used data rather than the complexity of the classification method.
cat <dig>  indeed, showed a very low rate of false positive results in comparison to existing state-of-the-art tools, namely, casvm and cascleave.
thirdly, cat <dig> is a straightforward sequence-based scoring system that offers an easy to use reference scale to determine the potential cleavage site instead of offering a yes-no answer or providing many suggested cleavage sites without any score to rank them.
in contrast to other tools that can execute a single sequence per query, cat <dig> is a fast system that can process both single and multiple sequence inputs: a feature that would assist biologists to perform large scale in silico screening to identify novel caspase- <dig> substrates.
our secondary structure analyses of caspase- <dig> substrates showed that regardless of cleavage patterns, aspartic acid residues are predominantly located in unstructured regions and to a lesser extent within alpha-helices.
in addition, we found the amino acids d, e, a, g and s appear more frequently in natively unstructured regions no matter whether they lie within or outside cleavage motifs.
these findings are in agreement with various reports that used statistical analysis to determine the natural distribution of these amino acids and their influence on secondary structure  <cit> .
this interesting result raises a question about the benefit of using local secondary structure properties of the cleavage sites as additional features to enhance the discrimination between cleaved and noncleaved patterns  <cit> .
careful evaluation of amino acid preference, at the positions surrounding the tetrapeptide cleavage motif, points to a general trend where the unfavourable amino acids have greater weight than the favoured ones, especially at p <dig>  p <dig>  p1', p3', p4', and p5' .
nevertheless, p1' has a remarkable preference for specific amino acids, namely glycine and serine.
in addition to their role in determining the molecular electrostatic potential and the steric accessibility of the enzyme, the post-translational modification potential of these two amino acids is vital for determining the timing and functional consequences of cleavage.
tözsér et al.
<cit>  demonstrated that phosphorylation of serine residues in close proximity to the tetrapeptide cleavage core can determine caspase- <dig> cleavability.
on the other hand, the high preference for glycine at p1' can be crucial to the acquisition of a myristic acid at this residue.
myristoylation is a co-translational reaction that occurs after the removal of the initiator methionine residue.
it can also occur as a post-translational modification when internal glycine residues become exposed after caspase cleavage.
the addition of a myristate moiety can alter subcellular localization of the cleaved proteins by facilitating their attachment to membranes and other proteins  <cit> .
by using cat <dig>  we carried out a large scale proteomic screen to identify novel potential caspase- <dig> substrates.
the initial screening showed that  <dig> human proteins can be potential caspase- <dig> substrates.
even after normalizing this result by excluding the noise coming from the presumed false positive rate, the percentage of potential caspase- <dig> substrates in the human proteome would be roughly ~14%.
this means that only a small fraction  of caspase- <dig> substrates has so far been discovered.
the results of go term enrichment analysis using toppfun showed that the majority of the predicted caspase- <dig> substrates are involved in cell adhesion, signal transduction, cell cycle, cytoskeleton organization, chromosome organization, neurogenesis, embryo development, cell morphogenesis, dna metabolism .
it is interesting to note the direct association of some of these processes to the biochemical events that lead to characteristic morphological changes in an apoptotic cell.
these changes include the breakup of the nuclear envelope and actin filaments in the cytoskeleton, blebbing of the plasma membrane, cell shrinkage, nuclear fragmentation, chromatin condensation, and chromosomal dna fragmentation  <cit> .
it is interesting to notice the remarkable presentation of some biological processes that are not related to apoptosis such as cell development and neurogenesis.
the careful analysis of the enriched biological processes go terms demonstrates a possible significant role of caspases- <dig> in the development and differentiation of nerve cells.
in fact, several reports have shown a strong expression of non-apoptotic active caspase- <dig> in various proliferating and differentiating neuronal cells  <cit> .
further investigation focusing on the role of caspases in nerve tissue may reveal new pathways that are necessary for the development and differentiation of nerve cells.
the results of enriched cellular component go terms showed that most of the predicted substrates are distributed to nuclear components , cytoplasmic components , and plasma membrane part .
this distribution is correlated with the normal subcellular localization of caspase- <dig>  although the procaspase- <dig> is localized in the cytoplasm, active caspase- <dig> plays essential roles both in the cytoplasm and nucleus  <cit> .
feng et al.
<cit>  have shown that the activated caspase- <dig> is first observed close to the inside surface of the cellular membrane, then transferred to the cytoplasm, and finally translocated to the nucleus.
an interesting fraction of the predicted caspase- <dig> substrates are proteins of the extracellular matrix.
the cleavage of such proteins can be achieved through their cytoplasmic embedded domains.
further investigations are needed to shed light on the biological importance of extracellular matrix proteins and their association with apoptotic and non-apoptotic roles of caspases.
CONCLUSIONS
in this work, we introduce a significant improvement to the in silico prediction approach of caspase substrates.
based on our results and in order to increase prediction specificity, we suggest the caspase-specific approach instead of that based upon considering the different caspases' substrates as having one common pattern.
cat <dig> can be considered a prototype system that would be easily utilized in developing prediction tools for other caspases and endopeptidases.
the predicted cellular targets of cat <dig> might be used to explore new pathways to gain further insight into the cellular mechanisms that regulate apoptosis, proliferation, and other biological processes.
in addition, the discovery of such targets might have significant implications for the development of drugs for various diseases including cancer, autoimmune disorders and neurodegenerative pathologies.
availability and requirements
project name: cat3
operating system: windows
programming language: perl
other requirements: none
any restrictions to use by non-academics: none
note: cat- <dig> v  <dig>  software that can process both single and multiple sequences is provided in the additional materials .
list of abbreviations
pssm: position specific scoring matrix; auc: area under curve; roc: receiver operating characteristic; svm: support vector machine; uniprotkb: universal protein resource knowledgebase; go: gene ontology; sen: sensitivity; spe: specificity; ppv: positive predictive value; npv: negative predictive value; acc: accuracy; mcc: matthew's correlation coefficient
authors' contributions
ma and ya jointly performed data collection and verification, ma performed the sequence analysis, wrote the cat <dig> perl code, and helped to draft the methodology section of the manuscript.
ht carried out the cross validation experiments, wrote the code to run the high-throughput screening, and helped to draft the validation and implantation sections of the manuscript.
ya designed and supervised the study, carried out the cleavage pattern analysis, performed the analysis of the high-throughput screening, and wrote the manuscript.
all authors read and approved the final manuscript.
supplementary material
additional file 1
full list of caspase- <dig> substrates.
this table shows the description of the obtained  <dig> caspase- <dig> substrates.
the cleavage evidence refers to the experimental method through which cleavage was identified.
sdm stands for site directed mutagenesis, while proteomics refers to experiments of high-throughput proteomic screening.
click here for file
additional file 2
heat-map representing caspase- <dig> cleavage pattern.
this heat-map represents the scores of the  <dig> amino acids in the scoring matrix 'a' .
the colour intensities reflect the magnitude of amino acid scores.
the blue scale denotes the positive scores, while the yellow to red scale denotes the negative values.
click here for file
additional file 3
prediction results of non-caspase- <dig> substrates.
the three excel sheets show the prediction results for  <dig> non-caspase- <dig> substrates .
the pub med id reference and the cat <dig> score for each cleavage site are shown.
click here for file
additional file 4
detailed results of the performance comparison.
this table shows the prediction results of cat <dig> versus casvm and cascleave for  <dig> caspase substrates  that were not originally included in the training of cat <dig>  the cleavage position in all these substrates were confirmed by site directed mutagenesis.
click here for file
additional file 5
list of predicted potential human caspase- <dig> substrates.
this table shows the list of  <dig> human proteins that were predicted to be caspase- <dig> substrates.
several proteins have more than one predicted site leading to a total of  <dig> potential cleavage sites.
only cleavage sites of scores ≥  <dig> are shown.
click here for file
additional file 6
list of the significantly enriched biological process go terms.
this table shows the full list of the significantly enriched terms of go: biological process among the  <dig> predicted caspase- <dig> substrates.
click here for file
additional file 7
chart representing biological processes among predicted substrates.
this figure shows the graphical representation of the enriched go terms: biological process among the  <dig> predicted caspase- <dig> substrates.
the red bars indicate the gene count  per each go term: biological process .
the blue dotted line shows the p-value for each go term that was derived by random sampling from the whole genome analysis.
click here for file
additional file 8
list of the significantly enriched cellular component go terms.
this table shows the full list of the significantly enriched terms of go: cellular component among the  <dig> predicted caspase- <dig> substrates.
click here for file
additional file 9
installation files for cat- <dig> tool.
this compressed folder contains a total of  <dig> files: final_cat <dig> exe, process.exe, algmodule.pm, and submodules.pm that are needed for cat <dig> software to work on windows operating system.
the two files: single_sequence-test.txt and multiple_sequences-test.txt are provided for testing purposes.
to install cat- <dig> v <dig>  on your pc, copy the  <dig> files to one folder and thereafter you can simply start working by clicking on the cat <dig> exe file.
click here for file
acknowledgements
the authors wish to thank mohamad amro, amjad alkhatib and hasan altaradah for their technical help and for developing the gui of cat <dig>  dr. fawzi alrazem for the careful reading of the manuscript.
we also wish to express our sincere thanks to dr. robin abu ghazaleh for his helpful comments and english proofreading.

 caspases are a family of intracellular cysteinyl aspartate-specific proteases that are highly conserved in multicellular organisms and are key regulators of apoptosis initiation and execution. at least members of the caspase family have been identified in mammals and they are grouped into two major sub-families, namely inflammatory caspases and apoptotic caspases. apoptosis associated caspases can be further classified into two groups: initiator caspases, including caspase- - and - which are present upstream of apoptosis signalling pathways; and executioner caspases - - and -. initiator capsases- and - are activated through an auto-cleavage process that is mediated by large adaptor-caspase complexes known respectively as the death-inducing signalling complex and apoptosome. these complexes are usually formed in response to an intrinsic or extrinsic cell death stimulus. the main targets of the activated initiator caspases are the executioner procaspases. it is interesting to notice that substrates of initiator caspases are limited to their own precursors, executioner procaspases- - and - and few more proteins. on the other hand, executioner caspases target a large number of cellular proteins to control the dismantling process of the cell. in addition to their essential role in apoptosis, recent accumulated evidence demonstrates various non-apoptotic functions of executioner caspases including: regulation of the immune response, cell proliferation, differentiation and motility. caspases are characterized by high substrate selectivity. they recognize a specific sequence signal in their target proteins. resolving substrate specificity for caspases was initially investigated using a combinatorial approach with positional scanning of synthetic tetrapeptidyl-aminomethyl coumarin derivatives. the results of this approach determined the absolute requirements for aspartic acid at position p. in addition, p to p positions demonstrate high preference for certain amino acids. based on positional scanning of synthetic tetrapeptides, the preferred recognition sequences for caspases - - and - were determined to be ehd, whereas caspases - and - recognize the sequence devd, while caspases- - - and - recognize the sequence ed. it is important to emphasize that the in vitro caspase substrate specificity, determined by the synthetic tetrapeptide method, is not absolutely representative of the cleavage conditions in vivo. the cleavage specificities of caspases in vivo are influenced by sequence-dependent conformational features, flanking the cleavage tetrapeptide motif, which can control the molecular electrostatic potential and the steric accessibility of the enzyme to its target protein. for example, and in spite of their identical preference for the devd tetrapeptide cleavage motif, caspase- and show a clear differential preference for various natural substrates. demon et al demonstrated that in addition to the tetrapeptide cleavage core, devd, several amino acid positions located outside this core such as p p p2' and p3' are critical in the discrimination of caspase- and caspase- for their specific substrates. shen et al reported another interesting example which proves that the relatively similar tetrapeptide cleavage motif of caspase- and caspase- which are functionally distinct, does not imply a similar recognition preference for their natural substrates. via a thorough statistical analysis of a window size of p10-p10' for a collection of caspase- and caspase- natural substrates, shen et al have determined the significance of various amino acids and/or certain physiochemical properties at certain positions outside the canonical tetrapeptide motif. among executioner caspases, caspase- is considered the major enzyme with a wide array of cellular substrates. while immunodepletion of caspase- abolishes the majority of proteolytic events observed during apoptosis, immunodepletion of other executioner caspases shows a minimal impact on apoptosis markers and its proteolytic cleavage outcomes. in the last decade, extensive research on caspases led to the identification of more than caspase- substrates and the list is still growing. with the increasing number of proteins that have been discovered, thanks to the sequencing of the human genome and the genomes of many other organisms, there is a need for efficient methods that can help in discovering new caspase- substrates. the identification of new cellular substrates for caspase- would lead to further insights into the cellular mechanisms that regulate apoptosis, proliferation, and other biological processes. bioinformatics tools would allow high-throughput analyses of proteomic data in order to screen for putative caspase- substrates. in addition, such tools can provide researchers with an accurate map of the potential cleavage site for a given sequence of interest. in the last few years, several computer-based tools were developed with the aim of predicting caspase substrates. prediction of endopeptidase substrates was among the initial tools and it was developed in order to predict putative caspase- cathepsin b and cathepsin l cleavage sites using cleavage site scoring matrices. peptidecutter is another tool that was designed with the objective of predicting cleavage sites for a wide range of proteases including various caspases. grabcas is a tool that uses a position specific scoring matrix for caspases 1- and granzyme b, based on substrate specificities that were determined by positional scanning of synthetic peptides. caspredictor was developed based on the assumption that sequences rich in the amino acids ser, thr, pro, glu or asp are favoured caspase cleavage sites. caspredictor was built based on experimentally verified natural substrates for caspase- - - - - - - and - in addition to the previously aforementioned scoring-matrix based approaches, several groups recently reported the development of tools that were mainly built up using the support vector machine technique. wee et al. described an svm based approach using substrates for different caspases from various organisms. cascleave is an interesting tool that was recently developed utilizing primary sequence, as well as secondary structure features, of the cleaved sites based on svm approach. cascleave was built using a dataset of substrates of the different caspases. piippo et al described another tool using three different pattern recognition classifiers, namely: svm, a decision tree based method known as j and the random forest classifier. the three classifiers were trained on different caspase cleavage sites. li et al proposed a hybrid svm-pssm method based on an extended dataset. unfortunately, some of these tools are not available for testing and comparison purposes. despite the substantial efforts to develop in silico systems to predict sites of caspase cleavage, the accuracy of such tools is still a challenging issue. the major drawback of the early tools is the use of training datasets that represent synthetic peptides or limited natural substrates for various proteases including caspases. on the other hand, the recently developed svm-based tools were built using a mixture of heterogeneous data that represent cleavage sites of the different caspases including non-canonical sites as well as some unverified caspase substrates. in general the svm-based tools such as cascleave achieve good levels of sensitivity, yet they suffer from high rates of false positive results. it is generally expected that training a prediction tool on data representing distinctive patterns can lead to overgeneralization and hence a high rate of false positive results. it is important to recall that although different caspases share the primary sequence-requirement to cleave at the carboxyl terminal of aspartate residues in their protein targets, each one of these proteases recognizes a unique context surrounding the cleavage position. even the caspases that appear to have identical tetrapeptide cleavage specificities such as caspase- and - are actually distinct in terms of the amino acids preferences outside the tetrapeptide core sequence. based on this assumption, we decided to develop a prediction tool focusing on data that represent substrates of a single caspase. caspase- was selected for this objective as it represents the major executioner caspase with a considerable number of substrates. in this work, we present a novel tool designated caspase analysis tool, which was developed based on an extensive and highly curated dataset of caspase- substrates. cat showed an obvious improvement in the overall prediction accuracy as well as a marked reduction of the false positive rate. using cat a high-throughput screening was performed on a large set of human proteins with assigned gene ontology annotations. the screening results reveal the existence of a large number of potential caspase- substrates. implementation methods caspase- substrates the pubmed literature database was used to search for papers that describe human, mouse, and rat caspase- substrates. each paper was critically analyzed to determine the experimentally demonstrated cleavage position in the relevant protein. the confirmed caspase- substrates were proteins with a total of cleavage sites. the amino acid sequences of the proteins were then obtained from the universal protein resource knowledgebase. of the cleavage sites, sequences were sorted randomly to be used later to compare the performance of cat versus existing similar tools; the remaining sequences, which we refer to as the positive peptide data, were used for training and validation of the cat definition of study controls the following datasets were established as controls in this study: the negative uncleaved peptides this data set consists of all the peptides that contain aspartic acid residues and are presumed to be uncleaved. this control group was established based on the assumption that any d residue, in a caspase- substrate, apart from the mapped cleavage site is most likely uncleaved. after excluding the positive peptides that exist in training data, the remaining d residues were used to create this negative control. amino acids natural frequency this control represents the frequency of each one of the amino acids in a group of human proteins that were available as reviewed proteins in the uniprotkb as of march amino acid r-groups frequency this control represents the frequencies of the different r-groups of amino acids; acidic, basic, polar, and non-polar. the frequencies were calculated based on the above mentioned reviewed proteins. physiochemical characteristics flanking the cleavage site the positive peptide sequences were aligned in reference to the cleaved aspartic acid residues. the resulting multiple sequence alignment was divided into three regions: the central tetrapeptide cleavage motif, the n-terminal region preceding the motif and the c-terminal region following the motif that was designated "before-motif" and "after-motif", respectively. the analyses for the regions flanking the motif were made serially: and amino acids before and after the motif. the analysis included: the frequencies of amino acids represented by their r-groups, the frequencies of hydrophobic and hydrophilic amino acids and finally the frequencies for each single amino acid. in the case of tetrapeptide motif analysis, the different frequencies were calculated separately for each position: p p p and p however, the frequencies within and amino acids, before and after the motif, were calculated collectively for each region. establishment of scoring matrices the peptides that fulfil the length criterion, p9-p5', which means having amino acids before and amino acids after the aspartate residue of interest, were used to build the scoring matrices. both the positive and the negative peptide data sets were used to build the scoring matrices. the first step was to generate position specific frequency matrices from the multiple sequence alignments of the relevant set of peptides. each matrix consisted of rows, representing positions p9p8p7p6p p p p p1p1' p2' p3'p4'p5', where a d amino acid is at the position p the columns of the matrix represent the frequencies of each amino acid. from the multiple sequence alignment of the positive peptides, we noticed the presence of two possibly different patterns; the first pattern has a d at p and the second has any amino acid except d at p. to represent this subtle difference, we decided to construct amino acid frequency matrices to represent each sub-pattern. two weighting systems were used in order to correct the probability of overrepresented and underrepresented amino acids in the frequency matrices so as to establish the scoring matrices: i) calculating log odd ratio: this weighting system involves calculating log odd ratio for each element in the frequency matrix by dividing the observed frequency of a given amino acid over its corresponding natural frequency. ii) subtraction of negative control background: instead of relying only on the common log odd weighting system and in order to minimize scoring bias, we decided to add a second normalization approach. the method relies on comparing the positive peptides with the negative peptides to further remove the noise signals around the cleaved aspartate residues. four scoring matrices are involved in the overall calculation of the final score of cat tool. we propose the following notation to define each scoring matrix and the overall score. first, let fm1+ denote the frequency matrix that was constructed from all the positive peptides. the corresponding scoring matrix'a' is defined as: a=log2fm1+Ω where Ω is the natural frequencies of the amino acids. in addition to the above scoring matrix we define fm1- as the frequency matrix generated from the negative peptides. a new frequency matrix'b' is defined as: b=log2fm1+-fm1-Ω let fmc denote a frequency matrix calculated from a subset of peptides that fulfil the constraint'c'. here, is either + or - as explained before. therefore, we define the following scoring matrices: c1=log2fm1p4=d+-fm1p4=d-Ω and c2=log2fm1p4≠d+-fm1p4≠d-Ω cat implementation and scoring cat tool was built using perl language. the input protein can be entered either as a fasta format sequence or as a text file. once a p peptide with a d residue at p is identified, it is analyzed to calculate the final score's' as follows: s=a+b+c where'a' and'b' are scores generated from the scoring matrices'a' and'b' in equation and equation respectively. the'c' score is generated either from the scoring matrix'c1'or'c2' as follows: c=c1ifp4=dc2ifp4≠d we refer to the scoring matrix'c1' if the peptide contains the amino acid d at p or the scoring matrix'c2' if the amino acid at p is not d. the three scores are normalized to a 100% score by dividing each score by the maximum score that could be obtained from each formula. cat validation to examine the prediction power of cat a k = fold cross validation was performed. the positive data were the actual cleavage sites, whereas the negative data were obtained from the uncleaved dataset. in each fold four pssm matrices were created from 9/10th of the positive substrates. then, the remaining 1/10th positive and negative substrates are used for testing. since the number of the negative peptides was much larger than the positive peptides, an equal number of the negative peptides were randomly obtained. the whole fold cross validation experiment was repeated times to ensure a good coverage of the negative dataset. the sensitivity, specificity, positive predictive value, negative predictive value, accuracy and the matthew's correlation coefficient were calculated as in. the areas under the receiver operating characteristic curves were calculated by plotting the sensitivity against the corresponding 1-specificity. the optimal cut-off point was defined as that measurement that corresponded to the point on the roc curve closest to the top left corner, i.e., closest to having sensitivity = and specificity = local window size the most appropriate local window size of amino acid sequence encompassing the cleaved aspartate and other critical exosite residues was determined based on comparing the prediction performance of the following peptides: p3-p p4-p p5-p1', p6-p2', p7-p3', p8-p4', p9-p5', p11-p7', p14-p10', p17-p13', p20-p16', p23-p19'. the performance of each window size was evaluated using the same cross validation approach as mentioned above. the obtained area under curve and matthew's correlation coefficient for the different experiments were plotted for comparison purposes. performance comparison a performance comparison was carried out for cat versus two recently published prediction tools, namely casvm and cascleave. the aim of the test was to assess how accurate the three tools were in predicting caspase- cleavage sites. the comparison was made on caspase- substrates that were randomly excluded from the training dataset. since cat is a prediction tool specific for caspase- cleavage sites, whereas casvm and cascleave were developed to predict cleavage sites of different caspases, there is a possibility to misjudge true positive sites of other caspases by assigning them to the false positive category of casvm and cascleave. to avoid such unfair comparison, the substrates were carefully inspected to find all caspase cleavage sites. the search was performed using the pubmed database, google searching engine, the caspase substrate database homepage, and merpos - the peptidases database. the protein sequences of the substrates were analyzed individually and the prediction results for each tool were counted according to software default parameters. the true positives are the positively predicted caspase- cleavage sites, whereas the false positives are the positively predicted aspartates that are actually not recognized by any caspase. high-throughput screening the uniprotkb was used to retrieve human proteins with known biological processes. two filters of the advance search option were used: the first was organism: homo sapiens and the second was gene ontology go: biological process. after excluding the experimentally verified human caspase- substrates, a total of reviewed human proteins with defined gene ontology were obtained. the protein sequences were analyzed by cat to screen for potential novel caspase- substrates. only results of scores ≥ were considered for further analyses. proteins that were predicted as potential caspase- substrates were further analyzed using toppgene suite tool to retrieve the most significant gene ontology terms.  caspase substrates our search in the pubmed literature database for caspase- substrates revealed the presence of proteins: of human origin, of mouse origin, and of rat origin. all the substrates were experimentally verified as natural substrates and their cleavage sites were mapped. of the substrates, the cleavage sites of proteins were mapped by site-directed mutagenesis technique, while the remaining were mapped by different high-throughput proteomic screening approaches. the full list and description of the obtained substrates are available in the additional materials. the obtained caspase- substrates as well as other caspase substrates will be available in the caspase substrates comprehensive database that has been developed at our institute. tetrapeptide cleavage motif analysis the tetrapeptide cleavage motifs of the training group were analyzed to determine physiochemical properties and frequencies of amino acids at each position. the examination of amino acid frequencies within the tetrapeptide motif revealed a unique distribution pattern of hydrophobic and hydrophilic amino acids. hydrophilic amino acids are times more frequent in p than hydrophobic amino acids. interestingly, p has an opposite pattern to p in p hydrophilic amino acids are nearly two times more frequent than hydrophobic amino acids, whereas in p the converse is true. features surrounding the cleavage site the amino acid sequences surrounding the tetrapeptide cleavage motifs were thoroughly analyzed to identify necessary feature for caspase- recognition. the analyses include: secondary structure, amino acids physicochemical properties and amino acid composition. the secondary structure prediction method gor was used to investigate the cleavage motif and its flanking regions for any common secondary structure. the analysis of gor results showed that the majority of the cleaved sites are located within unstructured context, while 18% are located within alpha helical regions, and only 2% are located in beta sheets. we then analyzed the biochemical properties of amino acids that flank the tetrapeptide cleavage motif to determine amino acids preferences for caspase- substrate recognition. no significant differences in the frequencies of acidic, basic, polar and non-polar amino acids between the tested region and the corresponding control group were found when examining and amino acids before and after the tetrapeptide cleavage motif. when testing the region of amino acids before the cleavage motif a slightly higher percentage of acidic amino acids were noticed, while basic and polar amino acids were strongly unfavoured. in the region of amino acids after the cleavage motif, lower percentages of acidic and basic amino acids were noticed. to further explore the characteristic biochemical properties, we examined individual amino acid frequencies for the entire cleavage vicinity: the 5-amino acids before and after the tetrapeptide motif. as shown in figure the frequencies of glycine, alanine, serine and proline have altered distributions in regions before and after the tetrapeptide motif that may indicate size and charge requirements for caspase- recognition and binding to the substrates. position specific scoring matrices in order to determine the most appropriate window size to construct efficient scoring matrices for cat a series of gradually increasing window sizes ranging from p3-p to p23-p19' were evaluated. as can be seen, it is obvious that window sizes equal to or shorter than the tetrapeptide cleavage motif are not adequate to develop a reliable prediction tool. despite of a marginally higher mcc at the window size of p6-p2', the overall prediction efficiency of the window sizes ranging from p5-p1' to p9-p5' are apparently quite similar. however, the efficiency seems to decrease gradually when extending the window size beyond p9-p5'. we actually preferred the window size p9-p5' over other seemingly comparable shorter alternatives for several reasons. first, the critical analysis of amino acids over-/under-representation scores demonstrated the significance of all the positions in this extended window p9-p5'. second, a careful analysis of natural caspase substrates, available in merops database, with cleavage positions near to n- or c- terminals, indicates that minimal adequate n- and c- terminal spacers comparable to the length of p and p5', respectively, are required for efficient recognition. therefore, our scoring matrices were developed by calculating the weight of each amino acid in the 14-mer peptide sequence from p to p5'. to evaluate the contribution of the different amino acids at the positions surrounding the cleaved aspartate, the scoring matrix'a' was drawn as a heat-map. analyzing the heat-map shows that apart from the tetrapeptide cleavage motif, all positions have either overrepresented or underrepresented amino acids. however, the positions including p p p1', p3', p4' and p5' have a remarkable rejection or preference for certain amino acids. prediction power of cat3 a fold cross validation was used to evaluate the predictive power of cat figure shows the roc curve that represents the average of different experiments of the 10-fold cross validation. the optimal cut-off score was found to be at this cut-off point the prediction statistical measures are shown in table the values of the statistical measures represent the average of different experiments of the fold cross validation test. the optimal cut-off score was found to be auc: area under roc curve; spe: specificity; sen: sensitivity; ppv: positive predictive value; npv: negative predictive value; acc: accuracy; mcc: matthew's correlation coefficient. to demonstrate specificity of cat for caspase- cleavage sites, a group of non-caspase- substrates were examined by cat the substrates included and cleavage sites of caspase- caspase- and caspase- respectively. we avoided using any cleavage site that was known to be a shared target with caspase- interestingly, of the cleavage sites showed cat scores below which is the minimum cut-off for predicting a caspase- cleavage site. this result provides a clear evidence to substantiate the very high specificity of cat for predicting caspase- substrates. the detailed results of this experiment are available in the additional materials. evaluating the performance of different binary classifiers is frequently made by comparing their reported statistical measures such as specificity, sensitivity etc. which are usually calculated under different conditions. we have avoided the use of such a comparison as it can lead to a biased conclusion. instead, we compared the performance of cat versus two recently reported tools, namely casvm and cascleave, on a group of caspase- substrates that were initially excluded from our training data. it is worth mentioning that some of these substrates could have been used in the training of the other tools, which could offer unfair advantage to the other two tools versus cat a thorough examination using different databases revealed that the substrates contain a total of aspartate residues, of which are caspase- cleavage sites, are cleavage sites of other caspases and aspartate residues that are evidently not cleaved by any caspase. out of the actual caspase- cleavage sites, the predicted true positive results for the three tools were as follows: cat 14/, casvm 8/, and cascleave 16/. however, cat was the best tool in reducing the false positive results; out of actual uncleaved aspartic acid containing peptide, cat predicted only false positives, while casvm predicted false positive and cascleave had false positives. figure shows the result of comparing cat versus casvm and cascleave. it is noteworthy that both casvm and cascleave correctly predicted two of the non-caspase- cleavage sites. the detailed results of the comparison are available in the additional materials. high-throughput screening for novel caspase- substrates screening of reviewed human proteins with ascribed gene ontology using our cat tool showed that proteins are predicted to be caspase- substrates with a total of potential caspase- cleavage sites. to further investigate the function of these potential substrates we used toppfun: an annotations based gene list functional enrichment analysis tool. out of the genes only had annotations in toppfun. the analysis revealed a group of biological processes that showed significant enrichment in predicted proteins versus the whole human genome as a control. a careful analysis of these biological processes was performed to shortlist the most significant biological processes by excluding general roots and detailed leaves of go terms. the most significant biological processes are shown in table the go terms of the listed biological processes were manually filtered, to reduce redundancy, by removing general roots and detailed leaves of the enriched go terms that were obtained by toppfun tool. the p-value of each go term in the predicted caspase- proteins was derived by random sampling of the whole genome. toppfun was also used to examine the enriched cellular component go terms. interestingly, the majority of the predicted proteins are located in different nuclear components, cytoskeleton, cell projection, membrane fraction, cell junction, and extracellular matrix, where most typical apoptotic morphological and biochemical changes are observed. the detailed list of the enriched cellular component go terms is available in the additional materials. discussion in addition to its well known key function in apoptosis, caspase- has been shown to play a crucial role in the regulation of various biological processes such as cell differentiation, adhesion, neurodevelopment and neuronal signalling. recognition of caspase- substrates is becoming a vital need to understand molecular mechanisms behind many disorders including cancer, autoimmune and neurodegenerative diseases. currently, most of the known caspase- substrates have been identified using in vitro proteolytic cleavage assays, coupled with site-directed mutagenesis to determine the exact cleavage position. in recognition of the physiological importance of caspase- many labs began to perform high-throughput proteomic screening to identify novel substrates of this major caspase. such techniques are relatively expensive and cumbersome. in addition to the high cost, the number of identified substrates is usually limited to the proteins that are relatively abundant in the examined cell type. recently, several computer-based prediction tools such as casvm, cascleave, and pripper were developed in order to help discover novel caspase substrates. these tools were trained on data that represent substrates of different caspases and in some cases non-caspase endopeptidase. although different caspases share a primary sequence-requirement, to cleave at the carboxyl terminal of aspartate residues in their protein targets, each one of these proteases needs a special context surrounding the cleavage position. even the caspases that appear to have identical tetrapeptide cleavage specificities such as caspase- and - are actually distinct in terms of the amino acids preferences outside the tetrapeptide core sequence. therefore, we believe that building a single algorithm for predicting the cleavage of multi-caspases would likely have low prediction specificity. based on this hypothesis, we decided to develop an algorithm focusing only on substrates of caspase- which is the major executioner caspase with a considerable number of targets. our caspase-3-specific approach has indeed outperformed other multi-caspases prediction tools on an independent comparison-dataset. cat has three distinctive features. firstly, it was developed using pssm instead of other relatively complex approaches. pssm is known to be practical, require low computation power and is able to represent the statistical weights of amino acids at each position. in addition, it can be easily combined with other machine learning tools to generate hybrid approaches that might enhance the prediction performance. secondly, instead of using data for different caspases, which are actually a mixture of heterogeneous patterns, we used an extended set of highly-curated caspase- natural substrates. we believe that inclusion of data that represents other proteases or cleavage sites of caspases with very few substrates and/or cleavage positions representing non-canonical patterns can lead to overgeneralization. in this situation, the classification model is required to loosen its decision boundary to increase sensitivity, but at the cost of having more false positive results. it is therefore generally accepted that improvement in prediction accuracy is more likely to be associated with the good quality of the used data rather than the complexity of the classification method. cat indeed, showed a very low rate of false positive results in comparison to existing state-of-the-art tools, namely, casvm and cascleave. thirdly, cat is a straightforward sequence-based scoring system that offers an easy to use reference scale to determine the potential cleavage site instead of offering a yes-no answer or providing many suggested cleavage sites without any score to rank them. in contrast to other tools that can execute a single sequence per query, cat is a fast system that can process both single and multiple sequence inputs: a feature that would assist biologists to perform large scale in silico screening to identify novel caspase- substrates. our secondary structure analyses of caspase- substrates showed that regardless of cleavage patterns, aspartic acid residues are predominantly located in unstructured regions and to a lesser extent within alpha-helices. in addition, we found the amino acids d, e, a, g and s appear more frequently in natively unstructured regions no matter whether they lie within or outside cleavage motifs. these findings are in agreement with various reports that used statistical analysis to determine the natural distribution of these amino acids and their influence on secondary structure. this interesting result raises a question about the benefit of using local secondary structure properties of the cleavage sites as additional features to enhance the discrimination between cleaved and noncleaved patterns. careful evaluation of amino acid preference, at the positions surrounding the tetrapeptide cleavage motif, points to a general trend where the unfavourable amino acids have greater weight than the favoured ones, especially at p p p1', p3', p4', and p5'. nevertheless, p1' has a remarkable preference for specific amino acids, namely glycine and serine. in addition to their role in determining the molecular electrostatic potential and the steric accessibility of the enzyme, the post-translational modification potential of these two amino acids is vital for determining the timing and functional consequences of cleavage. tözsér et al. demonstrated that phosphorylation of serine residues in close proximity to the tetrapeptide cleavage core can determine caspase- cleavability. on the other hand, the high preference for glycine at p1' can be crucial to the acquisition of a myristic acid at this residue. myristoylation is a co-translational reaction that occurs after the removal of the initiator methionine residue. it can also occur as a post-translational modification when internal glycine residues become exposed after caspase cleavage. the addition of a myristate moiety can alter subcellular localization of the cleaved proteins by facilitating their attachment to membranes and other proteins. by using cat we carried out a large scale proteomic screen to identify novel potential caspase- substrates. the initial screening showed that human proteins can be potential caspase- substrates. even after normalizing this result by excluding the noise coming from the presumed false positive rate, the percentage of potential caspase- substrates in the human proteome would be roughly ~14%. this means that only a small fraction of caspase- substrates has so far been discovered. the results of go term enrichment analysis using toppfun showed that the majority of the predicted caspase- substrates are involved in cell adhesion, signal transduction, cell cycle, cytoskeleton organization, chromosome organization, neurogenesis, embryo development, cell morphogenesis, dna metabolism. it is interesting to note the direct association of some of these processes to the biochemical events that lead to characteristic morphological changes in an apoptotic cell. these changes include the breakup of the nuclear envelope and actin filaments in the cytoskeleton, blebbing of the plasma membrane, cell shrinkage, nuclear fragmentation, chromatin condensation, and chromosomal dna fragmentation. it is interesting to notice the remarkable presentation of some biological processes that are not related to apoptosis such as cell development and neurogenesis. the careful analysis of the enriched biological processes go terms demonstrates a possible significant role of caspases- in the development and differentiation of nerve cells. in fact, several reports have shown a strong expression of non-apoptotic active caspase- in various proliferating and differentiating neuronal cells. further investigation focusing on the role of caspases in nerve tissue may reveal new pathways that are necessary for the development and differentiation of nerve cells. the results of enriched cellular component go terms showed that most of the predicted substrates are distributed to nuclear components, cytoplasmic components, and plasma membrane part. this distribution is correlated with the normal subcellular localization of caspase- although the procaspase- is localized in the cytoplasm, active caspase- plays essential roles both in the cytoplasm and nucleus. feng et al. have shown that the activated caspase- is first observed close to the inside surface of the cellular membrane, then transferred to the cytoplasm, and finally translocated to the nucleus. an interesting fraction of the predicted caspase- substrates are proteins of the extracellular matrix. the cleavage of such proteins can be achieved through their cytoplasmic embedded domains. further investigations are needed to shed light on the biological importance of extracellular matrix proteins and their association with apoptotic and non-apoptotic roles of caspases.  in this work, we introduce a significant improvement to the in silico prediction approach of caspase substrates. based on our results and in order to increase prediction specificity, we suggest the caspase-specific approach instead of that based upon considering the different caspases' substrates as having one common pattern. cat can be considered a prototype system that would be easily utilized in developing prediction tools for other caspases and endopeptidases. the predicted cellular targets of cat might be used to explore new pathways to gain further insight into the cellular mechanisms that regulate apoptosis, proliferation, and other biological processes. in addition, the discovery of such targets might have significant implications for the development of drugs for various diseases including cancer, autoimmune disorders and neurodegenerative pathologies. availability and requirements project name: cat3 operating system: windows programming language: perl other requirements: none any restrictions to use by non-academics: none note: cat- v software that can process both single and multiple sequences is provided in the additional materials. list of abbreviations pssm: position specific scoring matrix; auc: area under curve; roc: receiver operating characteristic; svm: support vector machine; uniprotkb: universal protein resource knowledgebase; go: gene ontology; sen: sensitivity; spe: specificity; ppv: positive predictive value; npv: negative predictive value; acc: accuracy; mcc: matthew's correlation coefficient authors' contributions ma and ya jointly performed data collection and verification, ma performed the sequence analysis, wrote the cat perl code, and helped to draft the methodology section of the manuscript. ht carried out the cross validation experiments, wrote the code to run the high-throughput screening, and helped to draft the validation and implantation sections of the manuscript. ya designed and supervised the study, carried out the cleavage pattern analysis, performed the analysis of the high-throughput screening, and wrote the manuscript. all authors read and approved the final manuscript. supplementary material additional file 1 full list of caspase- substrates. this table shows the description of the obtained caspase- substrates. the cleavage evidence refers to the experimental method through which cleavage was identified. sdm stands for site directed mutagenesis, while proteomics refers to experiments of high-throughput proteomic screening. click here for file additional file 2 heat-map representing caspase- cleavage pattern. this heat-map represents the scores of the amino acids in the scoring matrix'a'. the colour intensities reflect the magnitude of amino acid scores. the blue scale denotes the positive scores, while the yellow to red scale denotes the negative values. click here for file additional file 3 prediction results of non-caspase- substrates. the three excel sheets show the prediction results for non-caspase- substrates. the pub med id reference and the cat score for each cleavage site are shown. click here for file additional file 4 detailed results of the performance comparison. this table shows the prediction results of cat versus casvm and cascleave for caspase substrates that were not originally included in the training of cat the cleavage position in all these substrates were confirmed by site directed mutagenesis. click here for file additional file 5 list of predicted potential human caspase- substrates. this table shows the list of human proteins that were predicted to be caspase- substrates. several proteins have more than one predicted site leading to a total of potential cleavage sites. only cleavage sites of scores ≥ are shown. click here for file additional file 6 list of the significantly enriched biological process go terms. this table shows the full list of the significantly enriched terms of go: biological process among the predicted caspase- substrates. click here for file additional file 7 chart representing biological processes among predicted substrates. this figure shows the graphical representation of the enriched go terms: biological process among the predicted caspase- substrates. the red bars indicate the gene count per each go term: biological process. the blue dotted line shows the p-value for each go term that was derived by random sampling from the whole genome analysis. click here for file additional file 8 list of the significantly enriched cellular component go terms. this table shows the full list of the significantly enriched terms of go: cellular component among the predicted caspase- substrates. click here for file additional file 9 installation files for cat- tool. this compressed folder contains a total of files: final_cat exe, process.exe, algmodule.pm, and submodules.pm that are needed for cat software to work on windows operating system. the two files: single_sequence-test.txt and multiple_sequences-test.txt are provided for testing purposes. to install cat- v on your pc, copy the files to one folder and thereafter you can simply start working by clicking on the cat exe file. click here for file acknowledgements the authors wish to thank mohamad amro, amjad alkhatib and hasan altaradah for their technical help and for developing the gui of cat dr. fawzi alrazem for the careful reading of the manuscript. we also wish to express our sincere thanks to dr. robin abu ghazaleh for his helpful comments and english proofreading.
caspases are a family of cysteinyl proteases that regulate apoptosis and other biological processes. the majority of these proteins are involved in signal transduction, regulation of cell adhesion, cytoskeleton organization, integrity of the nucleus, and development of nerve cells. identification of caspase- cellular targets is crucial to gain further insights into the cellular mechanisms that have been implicated in various diseases including: cancer, neurodegenerative, and immunodeficiency diseases. human proteome screening, using cat indicate the presence of a large number of possible caspase- substrates that exceed the anticipated figure. apoptosiscaspase-3caspase substratescleavage site predictionposition-specific scoring matrix human proteomebioinformatic toolpattern recognition in contrast to other tools which were trained on cleavage sites of various caspases as well as other similar proteases, cat showed a significant decrease in the false positive rate. caspase- is considered the central executioner member of this family with a wide range of substrates. assessing prediction accuracy using fold cross validation, our method shows auc of, sensitivity of %, and specificity of %. cat is a powerful tool that is a clear improvement over existing similar tools, especially in reducing the false positive rate. the ability of cat in predicting the precise cleavage site was demonstrated in comparison to existing state-of-the-art tools. here, we describe a powerful bioinformatics tool that can predict the presence of caspase- cleavage sites in a given protein sequence using a position-specific scoring matrix approach. the developed tool, cat was used to screen human proteins with assigned gene ontology terms. in addition to their involvement in various expected functions such as cytoskeleton organization, nuclear integrity and adhesion, a large number of the predicted substrates are remarkably associated with the development of nerve tissues.
 caspases are a family of cysteinyl proteases that regulate apoptosis and other biological processes. caspase- is considered the central executioner member of this family with a wide range of substrates. identification of caspase- cellular targets is crucial to gain further insights into the cellular mechanisms that have been implicated in various diseases including: cancer, neurodegenerative, and immunodeficiency diseases. to date, over caspase- substrates have been identified experimentally. however, many are still awaiting discovery.  here, we describe a powerful bioinformatics tool that can predict the presence of caspase- cleavage sites in a given protein sequence using a position-specific scoring matrix approach. the present tool, which we call cat was built using confirmed caspase- substrates that were carefully extracted from the literature. assessing prediction accuracy using fold cross validation, our method shows auc of, sensitivity of %, and specificity of %. the ability of cat in predicting the precise cleavage site was demonstrated in comparison to existing state-of-the-art tools. in contrast to other tools which were trained on cleavage sites of various caspases as well as other similar proteases, cat showed a significant decrease in the false positive rate. this cost effective and powerful feature makes cat an ideal tool for high-throughput screening to identify novel caspase- substrates. the developed tool, cat was used to screen human proteins with assigned gene ontology terms. the analyses revealed the presence of many potential caspase- substrates that are not yet described. the majority of these proteins are involved in signal transduction, regulation of cell adhesion, cytoskeleton organization, integrity of the nucleus, and development of nerve cells.  cat is a powerful tool that is a clear improvement over existing similar tools, especially in reducing the false positive rate. human proteome screening, using cat indicate the presence of a large number of possible caspase- substrates that exceed the anticipated figure. in addition to their involvement in various expected functions such as cytoskeleton organization, nuclear integrity and adhesion, a large number of the predicted substrates are remarkably associated with the development of nerve tissues. apoptosiscaspase-3caspase substratescleavage site predictionposition-specific scoring matrix human proteomebioinformatic toolpattern recognition 