Toward computer-based cleavage site prediction of cysteine endopeptidases

Prediction and design of protease enzyme specificity using a structure-aware graph convolutional network

Site-specific proteolysis by the enzymatic cleavage of small linear sequence motifs is a key posttranslational modification involved in physiology and disease. The ability to robustly and rapidly predict protease-substrate specificity would also enable targeted proteolytic cleavage by designed proteases. Current methods for predicting protease specificity are limited to sequence pattern recognition in experimentally derived cleavage data obtained for libraries of potential substrates and generated separately for each protease variant. We reasoned that a more semantically rich and robust model of protease specificity could be developed by incorporating the energetics of molecular interactions between protease and substrates into machine learning workflows. We present Protein Graph Convolutional Network (PGCN), which develops a physically grounded, structure-based molecular interaction graph representation that describes molecular topology and interaction energetics to predict enzyme specificity. We show that PGCN accurately predicts the specificity landscapes of several variants of two model proteases. Node and edge ablation tests identified key graph elements for specificity prediction, some of which are consistent with known biochemical constraints for protease:substrate recognition. We used a pretrained PGCN model to guide the design of protease libraries for cleaving two noncanonical substrates, and found good agreement with experimental cleavage results. Importantly, the model can accurately assess designs featuring diversity at positions not present in the training data. The described methodology should enable the structure-based prediction of specificity landscapes of a wide variety of proteases and the construction of tailor-made protease editors for site-selectively and irreversibly modifying chosen target proteins.

Prediction and Design of Protease Enzyme Specificity Using a Structure-Aware Graph Convolutional Network

Site-specific proteolysis by the enzymatic cleavage of small linear sequence motifs is a key post-translational modification involved in physiology and disease. The ability to robustly and rapidly predict protease substrate specificity would also enable targeted proteolytic cleavage - editing - of a target protein by designed proteases. Current methods for predicting protease specificity are limited to sequence pattern recognition in experimentally-derived cleavage data obtained for libraries of potential substrates and generated separately for each protease variant. We reasoned that a more semantically rich and robust model of protease specificity could be developed by incorporating the three-dimensional structure and energetics of molecular interactions between protease and substrates into machine learning workflows. We present Protein Graph Convolutional Network (PGCN), which develops a physically-grounded, structure-based molecular interaction graph representation that describes molecular topology and interaction energetics to predict enzyme specificity. We show that PGCN accurately predicts the specificity landscapes of several variants of two model proteases: the NS3/4 protease from the Hepatitis C virus (HCV) and the Tobacco Etch Virus (TEV) proteases. Node and edge ablation tests identified key graph elements for specificity prediction, some of which are consistent with known biochemical constraints for protease:substrate recognition. We used a pre-trained PGCN model to guide the design of TEV protease libraries for cleaving two non-canonical substrates, and found good agreement with experimental cleavage results. Importantly, the model can accurately assess designs featuring diversity at positions not present in the training data. The described methodology should enable the structure-based prediction of specificity landscapes of a wide variety of proteases and the construction of tailor-made protease editors for site-selectively and irreversibly modifying chosen target proteins.

Cell Death Related Proteins Beyond Apoptosis in the CNS

Cell death related (CDR) proteins are a diverse group of proteins whose original function was ascribed to apoptotic cell death signaling. Recently, descriptions of non-apoptotic functions for CDR proteins have increased. In this minireview, we comment on recent studies of CDR proteins outside the field of apoptosis in the CNS, encompassing areas such as the inflammasome and non-apoptotic cell death, cytoskeleton reorganization, synaptic plasticity, mitophagy, neurodegeneration and calcium signaling among others. Furthermore, we discuss the evolution of proteomic techniques used to predict caspase substrates that could potentially explain their non-apoptotic roles. Finally, we address new concepts in the field of non-apoptotic functions of CDR proteins that require further research such the effect of sexual dimorphism on non-apoptotic CDR protein function and the emergence of zymogen-specific caspase functions.

An automated protocol for modelling peptide substrates to proteases

BACKGROUND

Proteases are key drivers in many biological processes, in part due to their specificity towards their substrates. However, depending on the family and molecular function, they can also display substrate promiscuity which can also be essential. Databases compiling specificity matrices derived from experimental assays have provided valuable insights into protease substrate recognition. Despite this, there are still gaps in our knowledge of the structural determinants. Here, we compile a set of protease crystal structures with bound peptide-like ligands to create a protocol for modelling substrates bound to protease structures, and for studying observables associated to the binding recognition.

RESULTS

As an application, we modelled a subset of protease-peptide complexes for which experimental cleavage data are available to compare with informational entropies obtained from protease-specificity matrices. The modelled complexes were subjected to conformational sampling using the Backrub method in Rosetta, and multiple observables from the simulations were calculated and compared per peptide position. We found that some of the calculated structural observables, such as the relative accessible surface area and the interaction energy, can help characterize a protease's substrate recognition, giving insights for the potential prediction of novel substrates by combining additional approaches.

CONCLUSION

Overall, our approach provides a repository of protease structures with annotated data, and an open source computational protocol to reproduce the modelling and dynamic analysis of the protease-peptide complexes.

Twenty years of bioinformatics research for protease-specific substrate and cleavage site prediction: a comprehensive revisit and benchmarking of existing methods

The roles of proteolytic cleavage have been intensively investigated and discussed during the past two decades. This irreversible chemical process has been frequently reported to influence a number of crucial biological processes (BPs), such as cell cycle, protein regulation and inflammation. A number of advanced studies have been published aiming at deciphering the mechanisms of proteolytic cleavage. Given its significance and the large number of functionally enriched substrates targeted by specific proteases, many computational approaches have been established for accurate prediction of protease-specific substrates and their cleavage sites. Consequently, there is an urgent need to systematically assess the state-of-the-art computational approaches for protease-specific cleavage site prediction to further advance the existing methodologies and to improve the prediction performance. With this goal in mind, in this article, we carefully evaluated a total of 19 computational methods (including 8 scoring function-based methods and 11 machine learning-based methods) in terms of their underlying algorithm, calculated features, performance evaluation and software usability. Then, extensive independent tests were performed to assess the robustness and scalability of the reviewed methods using our carefully prepared independent test data sets with 3641 cleavage sites (specific to 10 proteases). The comparative experimental results demonstrate that PROSPERous is the most accurate generic method for predicting eight protease-specific cleavage sites, while GPS-CCD and LabCaS outperformed other predictors for calpain-specific cleavage sites. Based on our review, we then outlined some potential ways to improve the prediction performance and ease the computational burden by applying ensemble learning, deep learning, positive unlabeled learning and parallel and distributed computing techniques. We anticipate that our study will serve as a practical and useful guide for interested readers to further advance next-generation bioinformatics tools for protease-specific cleavage site prediction.

Toward more accurate prediction of caspase cleavage sites: a comprehensive review of current methods, tools and features

As one of the few irreversible protein posttranslational modifications, proteolytic cleavage is involved in nearly all aspects of cellular activities, ranging from gene regulation to cell life-cycle regulation. Among the various protease-specific types of proteolytic cleavage, cleavages by casapses/granzyme B are considered as essential in the initiation and execution of programmed cell death and inflammation processes. Although a number of substrates for both types of proteolytic cleavage have been experimentally identified, the complete repertoire of caspases and granzyme B substrates remains to be fully characterized. To tackle this issue and complement experimental efforts for substrate identification, systematic bioinformatics studies of known cleavage sites provide important insights into caspase/granzyme B substrate specificity, and facilitate the discovery of novel substrates. In this article, we review and benchmark 12 state-of-the-art sequence-based bioinformatics approaches and tools for caspases/granzyme B cleavage prediction. We evaluate and compare these methods in terms of their input/output, algorithms used, prediction performance, validation methods and software availability and utility. In addition, we construct independent data sets consisting of caspases/granzyme B substrates from different species and accordingly assess the predictive power of these different predictors for the identification of cleavage sites. We find that the prediction results are highly variable among different predictors. Furthermore, we experimentally validate the predictions of a case study by performing caspase cleavage assay. We anticipate that this comprehensive review and survey analysis will provide an insightful resource for biologists and bioinformaticians who are interested in using and/or developing tools for caspase/granzyme B cleavage prediction.

Proteolysis mediated by cysteine cathepsins and legumain-recent advances and cell biological challenges

Proteases play essential roles in protein degradation, protein processing, and extracellular matrix remodeling in all cell types and tissues. They are also involved in protein turnover for maintenance of homeostasis and protein activation or inactivation for cell signaling. Proteases range in function and specificity, with some performing distinct substrate cleavages, while others accomplish proteolysis of a wide range of substrates. As such, different cell types use specialized molecular mechanisms to regulate the localization of proteases and their function within the compartments to which they are destined. Here, we focus on the cysteine family of cathepsin proteases and legumain, which act predominately within the endo-lysosomal pathway. In particular, recent knowledge on cysteine cathepsins and their primary regulator legumain is scrutinized in terms of their trafficking to endo-lysosomal compartments and other less recognized cellular locations. We further explore the mechanisms that regulate these processes and point to pathological cases which arise from detours taken by these proteases. Moreover, the emerging biological roles of specific forms and variants of cysteine cathepsins and legumain are discussed. These may be decisive, pathogenic, or even deadly when localizing to unusual cellular compartments in their enzymatically active form, because they may exert unexpected effects by alternative substrate cleavage. Hence, we propose future perspectives for addressing the actions of cysteine cathepsins and legumain as well as their specific forms and variants. The increasing knowledge in non-canonical aspects of cysteine cathepsin- and legumain-mediated proteolysis may prove valuable for developing new strategies to utilize these versatile proteases in therapeutic approaches.

Structure of the kidney slit diaphragm adapter protein CD2-associated protein as determined with electron microscopy

CD2-associated protein (CD2AP) is a multidomain scaffolding protein that has a critical role in renal function. CD2AP is expressed in glomerular podocytes at the slit diaphragm, a modified adherens junction that comprises the protein filtration barrier of the kidney, and interacts with a number of protein ligands involved in cytoskeletal remodeling, membrane trafficking, cell motility, and cell survival. The structure of CD2AP is unknown. We used electron microscopy and single particle image analysis to determine the three-dimensional structure of recombinant full-length CD2AP and found that the protein is a tetramer in solution. Image reconstruction of negatively stained protein particles generated a structure at 21 Å resolution. The protein assumed a roughly spherical, very loosely packed structure. Analysis of the electron density map revealed that CD2AP consists of a central coiled-coil domain, which forms the tetramer interface, surrounded by four symmetry-related motifs, each containing three globular domains corresponding to the three SH3 domains. The spatial organization exposes the binding sites of all 12 SH3 domains in the tetramer, allowing simultaneous binding to multiple targets. Determination of the structure of CD2AP provides novel insights into the biology of this slit diaphragm protein and lays the groundwork for characterizing the interactions between key molecules of the slit diaphragm that control glomerular filtration.

Macrobrachium rosenbergii cathepsin L: molecular characterization and gene expression in response to viral and bacterial infections

Cathepsin L (MrCathL) was identified from a constructed cDNA library of freshwater prawn Macrobrachium rosenbergii. MrCathL full-length cDNA is 1161 base pairs (bp) with an ORF of 1026bp which encodes a polypeptide of 342 amino acid (aa) long. The eukaryotic cysteine proteases, histidine and asparagine active site residues were identified in the aa sequence of MrCathL at 143-154, 286-296 and 304-323, respectively. The pair wise clustalW analysis of MrCathL showed the highest similarity (97%) with the homologous cathepsin L from Macrobrachium nipponense and the lowest similarity (70%) from human. Phylogenetic analysis revealed two distinct clusters of the invertebrates and vertebrates cathepsin L in the phylogenetic tree. MrCathL and cathepsin L from M. nipponense were clustered together, formed a sister group to cathepsin L of Penaeus monodon, and finally clustered to Lepeophtheirus salmonis. High level of (P<0.05) MrCathL gene expression was noticed in haemocyte and lowest in eyestalk. Furthermore, the MrCathL gene expression in M. rosenbergii was up-regulated in haemocyte by virus [M. rosenbergii nodovirus (MrNV) and white spot syndrome baculovirus (WSBV)] and bacteria (Vibrio harveyi and Aeromonas hydrophila). The recombinant MrCathL exhibited a wide range of activity in various pH between 3 and 10 and highest at pH 7.5. Cysteine proteinase (stefin A, stefin B and antipain) showed significant influence (100%) on recombinant MrCathL enzyme activity. The relative activity and residual activity of recombinant MrCathL against various metal ions or salts and detergent tested at different concentrations. These results indicated that the metal ions, salts and detergent had an influence on the proteinase activity of recombinant MrCathL. Conclusively, the results of this study imply that MrCathL has high pH stability and is fascinating object for further research on the function of cathepsin L in prawn innate immune system.

Cathepsin B is secreted apically from Xenopus 2F3 cells and cleaves the epithelial sodium channel (ENaC) to increase its activity

The epithelial sodium channel (ENaC) plays an important role in regulating sodium balance, extracellular volume, and blood pressure. Evidence suggests the α and γ subunits of ENaC are cleaved during assembly before they are inserted into the apical membranes of epithelial cells, and maximal activity of ENaC depends on cleavage of the extracellular loops of α and γ subunits. Here, we report that Xenopus 2F3 cells apically express the cysteine protease cathepsin B, as indicated by two-dimensional gel electrophoresis and mass spectrometry analysis. Recombinant GST ENaC α, β, and γ subunit fusion proteins were expressed in Escherichia coli and then purified and recovered from bacterial inclusion bodies. In vitro cleavage studies revealed the full-length ENaC α subunit fusion protein was cleaved by active cathepsin B but not the full-length β or γ subunit fusion proteins. Both single channel patch clamp studies and short circuit current experiments show ENaC activity decreases with the application of a cathepsin B inhibitor directly onto the apical side of 2F3 cells. We suggest a role for the proteolytic cleavage of ENaC by cathepsin B, and we suggest two possible mechanisms by which cathepsin B could regulate ENaC. Cathepsin B may cleave ENaC extracellularly after being secreted or intracellularly, while ENaC is present in the Golgi or in recycling endosomes.

Activation of the Nipah virus fusion protein in MDCK cells is mediated by cathepsin B within the endosome-recycling compartment

Proteolytic activation of the fusion protein of the highly pathogenic Nipah virus (NiV F) is a prerequisite for the production of infectious particles and for virus spread via cell-to-cell fusion. Unlike other paramyxoviral fusion proteins, functional NiV F activation requires endocytosis and pH-dependent cleavage at a monobasic cleavage site by endosomal proteases. Using prototype Vero cells, cathepsin L was previously identified to be a cleavage enzyme. Compared to Vero cells, MDCK cells showed substantially higher F cleavage rates in both NiV-infected and NiV F-transfected cells. Surprisingly, this could not be explained either by an increased F endocytosis rate or by elevated cathepsin L activities. On the contrary, MDCK cells did not display any detectable cathepsin L activity. Though we could confirm cathepsin L to be responsible for F activation in Vero cells, inhibitor studies revealed that in MDCK cells, cathepsin B was required for F-protein cleavage and productive replication of pathogenic NiV. Supporting the idea of an efficient F cleavage in early and recycling endosomes of MDCK cells, endocytosed F proteins and cathepsin B colocalized markedly with the endosomal marker proteins early endosomal antigen 1 (EEA-1), Rab4, and Rab11, while NiV F trafficking through late endosomal compartments was not needed for F activation. In summary, this study shows for the first time that endosomal cathepsin B can play a functional role in the activation of highly pathogenic NiV.

Profiling protease activities by dynamic proteomics workflows

Proteases play prominent roles in many physiological processes and the pathogenesis of various diseases, which makes them interesting drug targets. To fully understand the functional role of proteases in these processes, it is necessary to characterize the target specificity of the enzymes, identify endogenous substrates and cleavage products as well as protease activators and inhibitors. The complexity of these proteolytic networks presents a considerable analytic challenge. To comprehensively characterize these systems, quantitative methods that capture the spatial and temporal distributions of the network members are needed. Recently, activity-based workflows have come to the forefront to tackle the dynamic aspects of proteolytic processing networks in vitro, ex vivo and in vivo. In this review, we will discuss how mass spectrometry-based approaches can be used to gain new insights into protease biology by determining substrate specificities, profiling the activity-states of proteases, monitoring proteolysis in vivo, measuring reaction kinetics and defining in vitro and in vivo proteolytic events. In addition, examples of future aspects of protease research that go beyond mass spectrometry-based applications are given.

Developing a powerful in silico tool for the discovery of novel caspase-3 substrates: a preliminary screening of the human proteome

Background

Caspases are a family of cysteinyl proteases that regulate apoptosis and other biological processes. Caspase-3 is considered the central executioner member of this family with a wide range of substrates. Identification of caspase-3 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 200 caspase-3 substrates have been identified experimentally. However, many are still awaiting discovery.

Results

Here, we describe a powerful bioinformatics tool that can predict the presence of caspase-3 cleavage sites in a given protein sequence using a Position-Specific Scoring Matrix (PSSM) approach. The present tool, which we call CAT3, was built using 227 confirmed caspase-3 substrates that were carefully extracted from the literature. Assessing prediction accuracy using 10 fold cross validation, our method shows AUC (area under the ROC curve) of 0.94, sensitivity of 88.83%, and specificity of 89.50%. The ability of CAT3 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, CAT3 showed a significant decrease in the false positive rate. This cost effective and powerful feature makes CAT3 an ideal tool for high-throughput screening to identify novel caspase-3 substrates.

The developed tool, CAT3, was used to screen 13,066 human proteins with assigned gene ontology terms. The analyses revealed the presence of many potential caspase-3 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.

Conclusions

CAT3 is a powerful tool that is a clear improvement over existing similar tools, especially in reducing the false positive rate. Human proteome screening, using CAT3, indicate the presence of a large number of possible caspase-3 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.

A review of statistical methods for prediction of proteolytic cleavage

A fundamental component of systems biology, proteolytic cleavage is involved in nearly all aspects of cellular activities: from gene regulation to cell lifecycle regulation. Current sequencing technologies have made it possible to compile large amount of cleavage data and brought greater understanding of the underlying protein interactions. However, the practical impossibility to exhaustively retrieve substrate sequences through experimentation alone has long highlighted the need for efficient computational prediction methods. Such methods must be able to quickly mark substrate candidates and putative cleavage sites for further analysis. Available methods and expected reliability depend heavily on the type and complexity of proteolytic action, as well as the availability of well-labelled experimental data sets: factors varying greatly across enzyme families. For this review, we chose to give a quick overview of the general issues and challenges in cleavage prediction methods followed by a more in-depth presentation of major techniques and implementations, with a focus on two particular families of cysteine proteases: caspases and calpains. Through their respective differences in proteolytic specificity (high for caspases, broader for calpains) and data availability (much lower for calpains), we aimed to illustrate the strengths and limitations of techniques ranging from position-based matrices and decision trees to more flexible machine-learning methods such as hidden Markov models and Support Vector Machines. In addition to a technical overview for each family of algorithms, we tried to provide elements of evaluation and performance comparison across methods.

ApoptoProteomics, an integrated database for analysis of proteomics data obtained from apoptotic cells

Apoptosis is the most commonly described form of programmed cell death, and dysfunction is implicated in a large number of human diseases. Many quantitative proteome analyses of apoptosis have been performed to gain insight in proteins involved in the process. This resulted in large and complex data sets that are difficult to evaluate. Therefore, we developed the ApoptoProteomics database for storage, browsing, and analysis of the outcome of large scale proteome analyses of apoptosis derived from human, mouse, and rat. The proteomics data of 52 publications were integrated and unified with protein annotations from UniProt-KB, the caspase substrate database homepage (CASBAH), and gene ontology. Currently, more than 2300 records of more than 1500 unique proteins were included, covering a large proportion of the core signaling pathways of apoptosis. Analysis of the data set revealed a high level of agreement between the reported changes in directionality reported in proteomics studies and expected apoptosis-related function and may disclose proteins without a current recognized involvement in apoptosis based on gene ontology. Comparison between induction of apoptosis by the intrinsic and the extrinsic apoptotic signaling pathway revealed slight differences. Furthermore, proteomics has significantly contributed to the field of apoptosis in identifying hundreds of caspase substrates. The database is available at http://apoptoproteomics.uio.no.

Proteomic identification of protease cleavage sites characterizes prime and non-prime specificity of cysteine cathepsins B, L, and S

Cysteine cathepsins mediate proteome homeostasis and have pivotal functions in diseases such as cancer. To better understand substrate recognition by cathepsins B, L, and S, we applied proteomic identification of protease cleavage sites (PICS) for simultaneous profiling of prime and non-prime specificity. PICS profiling of cathepsin B endopeptidase specificity highlights strong selectivity for glycine in P3' due to an occluding loop blocking access to the primed subsites. In P1', cathepsin B has a partial preference for phenylalanine, which is not found for cathepsins L and S. Occurrence of P1' phenylalanine often coincides with aromatic residues in P2. For cathepsin L, PICS identifies 845 cleavage sites, representing the most comprehensive PICS profile to date. Cathepsin L specificity is dominated by the canonical preference for aromatic residues in P2 with limited contribution of prime-site selectivity determinants. Profiling of cathepsins B and L with a shorter incubation time (4 h instead of 16 h) did not reveal time-dependency of individual specificity determinants. Cathepsin S specificity was profiled at pH 6.0 and 7.5. The PICS profiles at both pH values display a high degree of similarity. Cathepsin S specificity is primarily guided by aliphatic residues in P2 with limited importance of prime-site residues.

CD2AP in mouse and human podocytes controls a proteolytic program that regulates cytoskeletal structure and cellular survival

Kidney podocytes are highly differentiated epithelial cells that form interdigitating foot processes with bridging slit diaphragms (SDs) that regulate renal ultrafiltration. Podocyte injury results in proteinuric kidney disease, and genetic deletion of SD-associated CD2-associated protein (CD2AP) leads to progressive renal failure in mice and humans. Here, we have shown that CD2AP regulates the TGF-β1-dependent translocation of dendrin from the SD to the nucleus. Nuclear dendrin acted as a transcription factor to promote expression of cytosolic cathepsin L (CatL). CatL proteolyzed the regulatory GTPase dynamin and the actin-associated adapter synaptopodin, leading to a reorganization of the podocyte microfilament system and consequent proteinuria. CD2AP itself was proteolyzed by CatL, promoting sustained expression of the protease during podocyte injury, and in turn increasing the apoptotic susceptibility of podocytes to TGF-β1. Our study identifies CD2AP as the gatekeeper of the podocyte TGF-β response through its regulation of CatL expression and defines a molecular mechanism underlying proteinuric kidney disease.

Specialized roles for cysteine cathepsins in health and disease

Cathepsins were originally identified as proteases that act in the lysosome. Recent work has uncovered nontraditional roles for cathepsins in the extracellular space as well as in the cytosol and nucleus. There is strong evidence that subspecialized and compartmentalized cathepsins participate in many physiologic and pathophysiologic cellular processes, in which they can act as both digestive and regulatory proteases. In this review, we discuss the transcriptional and translational control of cathepsin expression, the regulation of intracellular sorting of cathepsins, and the structural basis of cathepsin activation and inhibition. In particular, we highlight the emerging roles of various cathepsin forms in disease, particularly those of the cardiac and renal systems.

Pripper: prediction of caspase cleavage sites from whole proteomes

BACKGROUND

Caspases are a family of proteases that have central functions in programmed cell death (apoptosis) and inflammation. Caspases mediate their effects through aspartate-specific cleavage of their target proteins, and at present almost 400 caspase substrates are known. There are several methods developed to predict caspase cleavage sites from individual proteins, but currently none of them can be used to predict caspase cleavage sites from multiple proteins or entire proteomes, or to use several classifiers in combination. The possibility to create a database from predicted caspase cleavage products for the whole genome could significantly aid in identifying novel caspase targets from tandem mass spectrometry based proteomic experiments.

RESULTS

Three different pattern recognition classifiers were developed for predicting caspase cleavage sites from protein sequences. Evaluation of the classifiers with quality measures indicated that all of the three classifiers performed well in predicting caspase cleavage sites, and when combining different classifiers the accuracy increased further. A new tool, Pripper, was developed to utilize the classifiers and predict the caspase cut sites from an arbitrary number of input sequences. A database was constructed with the developed tool, and it was used to identify caspase target proteins from tandem mass spectrometry data from two different proteomic experiments. Both known caspase cleavage products as well as novel cleavage products were identified using the database demonstrating the usefulness of the tool. Pripper is not restricted to predicting only caspase cut sites, but it gives the possibility to scan protein sequences for any given motif(s) and predict cut sites once a suitable cut site prediction model for any other protease has been developed. Pripper is freely available and can be downloaded from http://users.utu.fi/mijopi/Pripper.

CONCLUSIONS

We have developed Pripper, a tool for reading an arbitrary number of proteins in FASTA format, predicting their caspase cleavage sites and outputting the cleaved sequences to a new FASTA format sequence file. We show that Pripper is a valuable tool in identifying novel caspase target proteins from modern proteomics experiments.

Enzymatic disease of the podocyte

Proteinuria is an early sign of kidney disease and has gained increasing attention over the past decade because of its close association with cardio-vascular and renal morbidity and mortality. Podocytes have emerged as the cell type that is critical in maintaining proper functioning of the kidney filter. A few genes have been identified that explain genetic glomerular failure and recent insights shed light on the pathogenesis of acquired proteinuric diseases. This review highlights the unique role of the cysteine protease cathepsin L as a regulatory rather than a digestive protease and its action on podocyte structure and function. We provide arguments why many glomerular diseases can be regarded as podocyte enzymatic disorders.

Prediction of protease substrates using sequence and structure features

MOTIVATION

Granzyme B (GrB) and caspases cleave specific protein substrates to induce apoptosis in virally infected and neoplastic cells. While substrates for both types of proteases have been determined experimentally, there are many more yet to be discovered in humans and other metazoans. Here, we present a bioinformatics method based on support vector machine (SVM) learning that identifies sequence and structural features important for protease recognition of substrate peptides and then uses these features to predict novel substrates. Our approach can act as a convenient hypothesis generator, guiding future experiments by high-confidence identification of peptide-protein partners.

RESULTS

The method is benchmarked on the known substrates of both protease types, including our literature-curated GrB substrate set (GrBah). On these benchmark sets, the method outperforms a number of other methods that consider sequence only, predicting at a 0.87 true positive rate (TPR) and a 0.13 false positive rate (FPR) for caspase substrates, and a 0.79 TPR and a 0.21 FPR for GrB substrates. The method is then applied to approximately 25 000 proteins in the human proteome to generate a ranked list of predicted substrates of each protease type. Two of these predictions, AIF-1 and SMN1, were selected for further experimental analysis, and each was validated as a GrB substrate.

AVAILABILITY

All predictions for both protease types are publically available at http://salilab.org/peptide. A web server is at the same site that allows a user to train new SVM models to make predictions for any protein that recognizes specific oligopeptide ligands.

Toward the development of podocyte-specific drugs

Most kidney diseases that ultimately lead to end-stage renal failure originate within the glomerulus and are associated with proteinuria. Treatment options are unspecific and offer partial cures at best because available therapies do not primarily treat glomerular cells but rather act systemically and thus cause many side effects. Most glomerulopathies directly stem from injury to podocytes, cells that have a key role in the maintenance of the glomerular filter. Thus, these cells constitute an obvious and promising target for the development of novel kidney-protective drugs. During the last decade, enormous advances have been made in the understanding of podocyte structure and function. A number of pathways that are altered during glomerular diseases may be targeted by novel small- and large-molecule drugs as well as biologicals that have been identified in nephrology and other areas of drug development. Cultured podocytes provide a valuable model for high-throughput drug screening assays. Furthermore, podocytes have been shown to possess many features that make them particularly good target cells for renal protection. This mini-review discusses some of the most recent promising data related to potential drug therapy for proteinuria and kidney disease through direct podocyte targeting.

Cascleave: towards more accurate prediction of caspase substrate cleavage sites

MOTIVATION

The caspase family of cysteine proteases play essential roles in key biological processes such as programmed cell death, differentiation, proliferation, necrosis and inflammation. The complete repertoire of caspase substrates remains to be fully characterized. Accordingly, systematic computational screening studies of caspase substrate cleavage sites may provide insight into the substrate specificity of caspases and further facilitating the discovery of putative novel substrates.

RESULTS

In this article we develop an approach (termed Cascleave) to predict both classical (i.e. following a P(1) Asp) and non-typical caspase cleavage sites. When using local sequence-derived profiles, Cascleave successfully predicted 82.2% of the known substrate cleavage sites, with a Matthews correlation coefficient (MCC) of 0.667. We found that prediction performance could be further improved by incorporating information such as predicted solvent accessibility and whether a cleavage sequence lies in a region that is most likely natively unstructured. Novel bi-profile Bayesian signatures were found to significantly improve the prediction performance and yielded the best performance with an overall accuracy of 87.6% and a MCC of 0.747, which is higher accuracy than published methods that essentially rely on amino acid sequence alone. It is anticipated that Cascleave will be a powerful tool for predicting novel substrate cleavage sites of caspases and shedding new insights on the unknown caspase-substrate interactivity relationship.

AVAILABILITY

http://sunflower.kuicr.kyoto-u.ac.jp/ approximately sjn/Cascleave/

CONTACT

jiangning.song@med.monash.edu.au; takutsu@kuicr.kyoto-u.ac.jp; james; whisstock@med.monash.edu.au

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Caspase substrates: easily caught in deep waters?

Caspases are key players in various cellular processes, such as apoptosis, proliferation and differentiation, and in pathological conditions including cancer and inflammation. Although caspases preferentially cleave C-terminal of aspartic acid residues, their action is restricted generally to one or a few sites per protein substrate. Caspase-specific substrate recognition appears to be determined by the substrate sequences adjacent to the scissile bond. Knowledge of these substrates and the generated fragments is crucial for a thorough understanding of the functional implications of caspase-mediated proteolysis. In addition, insight into the cleavage specificity might assist in designing inhibitors that target disease-related caspase activities. Here, we critically review recently published procedures used to generate a proteome-wide view of caspase substrates.

Comparative docking study of anibamine as the first natural product CCR5 antagonist in CCR5 homology models

Anibamine, a novel pyridine quaternary alkaloid recently isolated from Aniba sp., has been found to effectively bind to the chemokine receptor CCR5 with an IC(50) at 1 microM in competition with (125)I-gp120, an HIV viral envelope protein binding to CCR5 with high affinity. Since CCR5, a G-protein-coupled receptor, is an essential coreceptor for the human immunodeficiency virus type I (HIV-1) entry to host cells, a CCR5 antagonist that inhibits the cellular entry of HIV-1 provides a new therapy choice for the treatment of HIV. Anibamine provides a novel structural skeleton that is remarkably different from all lead compounds previously identified as CCR5 antagonists. Here, we report comparative docking studies of anibamine with several other known CCR5 antagonists in two CCR5 homology models built based on the crystal structures of bovine rhodopsin and human beta(2)-adrenergic receptor. The binding pocket of anibamine has some common features shared with other high affinity CCR5 antagonists, suggesting that they may bind in similar binding sites and/or modes. At the same time, several unique binding features of anibamine were identified, and it will likely prove beneficial in future molecular design of novel CCR5 antagonists based on the anibamine scaffold.

The actin cytoskeleton of kidney podocytes is a direct target of the antiproteinuric effect of cyclosporine A

The immunosuppressive action of the calcineurin inhibitor cyclosporine A (CsA) stems from the inhibition of nuclear factor of activated T cells (NFAT) signaling in T cells. CsA is also used for the treatment of proteinuric kidney diseases. As it stands, the antiproteinuric effect of CsA is attributed to its immunosuppressive action. Here we show that the beneficial effect of CsA on proteinuria is not dependent on NFAT inhibition in T cells, but rather results from the stabilization of the actin cytoskeleton in kidney podocytes. CsA blocks the calcineurin-mediated dephosphorylation of synaptopodin, a regulator of Rho GTPases in podocytes, thereby preserving the phosphorylation-dependent synaptopodin-14-3-3 beta interaction. Preservation of this interaction, in turn, protects synaptopodin from cathepsin L-mediated degradation. These results represent a new view of calcineurin signaling and shed further light on the treatment of proteinuric kidney diseases. Novel calcineurin substrates such as synaptopodin may provide promising starting points for antiproteinuric drugs that avoid the serious side effects of long-term CsA treatment.

Basement membrane derived fibulin-1 and fibulin-5 function as angiogenesis inhibitors and suppress tumor growth

The fibulins are a family of secreted glycoproteins that are characterized by repeated epidermal-growth-factor-like domains and a unique C-terminus structure. Fibulins modulate cell morphology, growth, adhesion, and motility. Our initial basement membrane degradome screen using Cathepsin D, a tumor microenvironment-associated protease, contained fragments of fibulin-1 and full length fibulin-5. In this report, we evaluate the antiangiogenic activity of fibulin-1 and fibulin-5. Tumor studies demonstrate that both fibulin-1 and fibulin-5 suppress HT1080 tumor growth. CD31 labeling and TUNEL assay further reveal that fibulin-1 suppression of HT1080 tumor growth is associated with diminished angiogenesis and also enhanced apoptosis of endothelial cells and tumor cells. In contrast, fibulin-5 inhibits tumor angiogenesis with a minimal anti-apoptotic affect. Cathepsin D digestion of fibulin-1 produces a fragment with nearly the same molecular weight as fibulin-5, and this fragment (named Neostatin) inhibits endothelial cell proliferation. Additionally, degradation of basement membrane by cathepsin D liberates both fibulin-1 fragments and fibulin-5, which function to inhibit angiogenesis.

Proteolytic processing of dynamin by cytoplasmic cathepsin L is a mechanism for proteinuric kidney disease

Kidney podocytes and their foot processes maintain the ultrafiltration barrier and prevent urinary protein loss (proteinuria). Here we show that the GTPase dynamin is essential for podocyte function. During proteinuric kidney disease, induction of cytoplasmic cathepsin L leads to cleavage of dynamin at an evolutionary conserved site, resulting in reorganization of the podocyte actin cytoskeleton and proteinuria. Dynamin mutants that lack the cathepsin L site, or render the cathepsin L site inaccessible through dynamin self-assembly, are resistant to cathepsin L cleavage. When delivered into mice, these mutants restored podocyte function and resolve proteinuria. Our study identifies dynamin as a critical regulator of renal permselectivity that is specifically targeted by proteolysis under pathological conditions.

Proteolytic processing of dynamin by cytoplasmic cathepsin L is a mechanism for proteinuric kidney disease

Kidney podocytes and their foot processes maintain the ultrafiltration barrier and prevent urinary protein loss (proteinuria). Here we show that the GTPase dynamin is essential for podocyte function. During proteinuric kidney disease, induction of cytoplasmic cathepsin L leads to cleavage of dynamin at an evolutionary conserved site, resulting in reorganization of the podocyte actin cytoskeleton and proteinuria. Dynamin mutants that lack the cathepsin L site, or render the cathepsin L site inaccessible through dynamin self-assembly, are resistant to cathepsin L cleavage. When delivered into mice, these mutants restored podocyte function and resolve proteinuria. Our study identifies dynamin as a critical regulator of renal permselectivity that is specifically targeted by proteolysis under pathological conditions.

Cysteine cathepsins: cellular roadmap to different functions

Cysteine cathepsins belong to the papain-like family C1 of clan CA cysteine peptidases. These enzymes are ubiquitously expressed and exert their proteolytic activity mainly, but not exclusively within the compartments along the endocytic pathway. Moreover, cysteine cathepsins are active in pericellular environments as soluble enzymes or bound to cell surface receptors at the plasma membrane, and possibly even within secretory vesicles, the cytosol, mitochondria, and within the nuclei of eukaryotic cells. Proteolytic actions performed by cysteine cathepsins are essential in the maintenance of homeostasis and depend heavily upon their correct sorting and trafficking within cells. As a consequence, the numerous and diverse approaches to identification, qualitative and quantitative determination, and visualization of cysteine cathepsin functions in vitro, in situ, and in vivo cover the entire spectrum of biochemistry, molecular and cell biology. This review focuses upon the transport pathways directing cysteine cathepsins to their points of action and thus emphasizes the broader role and functionality of cysteine cathepsins in a number of specific cellular locales. Such understanding will provide a foundation for future research investigating the involvement of these peptidases with their substrates, inhibitors, and the intertwined proteolytic networks at the hubs of complex biological systems.

CASVM: web server for SVM-based prediction of caspase substrates cleavage sites

Caspases belong to a unique class of cysteine proteases which function as critical effectors of apoptosis, inflammation and other important cellular processes. Caspases cleave substrates at specific tetrapeptide sites after a highly conserved aspartic acid residue. Prediction of such cleavage sites will complement structural and functional studies on substrates cleavage as well as discovery of new substrates. We have recently developed a support vector machines (SVM) method to address this issue. Our algorithm achieved an accuracy ranging from 81.25 to 97.92%, making it one of the best methods currently available. CASVM is the web server implementation of our SVM algorithms, written in Perl and hosted on a Linux platform. The server can be used for predicting non-canonical caspase substrate cleavage sites. We have also included a relational database containing experimentally verified caspase substrates retrievable using accession IDs, keywords or sequence similarity.

AVAILABILITY

http://www.casbase.org/casvm/index.html

Cathepsin B mediates cleavage of herpes simplex virus type 1 origin binding protein (OBP) to yield OBPC-1, and cleavage is dependent upon viral DNA replication

Although the seven viral proteins required for herpes simplex virus type 1 (HSV-1) DNA replication have been identified, the mechanism by which viral DNA synthesis is regulated is unclear. HSV-1 DNA replication is thought to occur in two stages: origin-dependent DNA replication (stage I) mediated by the origin binding protein (OBP), followed by origin- and OBP-independent DNA replication (stage II). The mechanism that facilitates the switch from stage I to stage II is unknown; however, it must involve the loss of OBP function or OBP itself from the replication initiation complex. Previous studies from this laboratory identified a transcript (UL8.5) and protein (OBPC) that are in frame with and comprise the C terminus of the gene specifying OBP. Because of its DNA binding ability, OBPC has been hypothesized to mediate the switch from stage I to stage II. Here, we identify a second protein (OBPC-2) that is also in frame with the C terminus of OBP but comprises a smaller portion of the protein. We demonstrate that the protein originally identified (OBPC-1) is a cathepsin B-mediated cleavage product of OBP, while OBPC-2 may be the product of the UL8.5 transcript. We further demonstrate that the cleavage of OBP to yield OBPC-1 is dependent upon viral DNA replication. These results suggest that cleavage may be a mechanism by which OBP levels and/or activity are regulated during infection.

Processing sites are different in the generation of HLA-A2.1-restricted, T cell reactive tumor antigen epitopes and viral epitopes

In order to improve the processing efficiency of T cell tumor antigen epitopes, this bioinformatic study compares proteolytic sites in the generation of 47 experimentally identified HLA-A2.1-restricted immunodominant tumor antigen epitopes to those of 52 documented HLA-A2.1-restricted immunodominant viral antigen epitopes. Our results show that the amino acid frequencies in the C-terminal cleavage sites of the tumor antigen epitopes, as well as several positions within the 10 amino acid (aa) flanking regions, are significantly different from those of the viral antigen epitopes. In the 9 amino acid epitope region, frequencies differed somewhat in the secondary-anchored amino acid residues on E3 (the third aa of the epitope), E4, E6, E7 and E8; however, frequencies in the primary-anchored positions, on E2 and E9, for binding in the HLA-A2.1 groove, remained almost identical. The most frequently occurring amino acid pairs in both N-terminal and C-terminal cleavage sites in the generation of tumor antigen epitopes were different from those of the viral antigen epitopes. Our findings demonstrate for the first time that these two groups of epitopes may be cleaved by distinct sets of proteasomes and peptidases or similar enzymes with lower efficiencies for tumor epitopes. In the future, in order to more effectively generate tumor antigen epitopes, targeted activation of the immunoproteasomes and peptidases that mediate the cleavage of viral epitopes could be achieved, thus enhancing our potential for antigen-specific tumor immunotherapy.

SVM-based prediction of caspase substrate cleavage sites

BACKGROUND

Caspases belong to a class of cysteine proteases which function as critical effectors in apoptosis and inflammation by cleaving substrates immediately after unique sites. Prediction of such cleavage sites will complement structural and functional studies on substrates cleavage as well as discovery of new substrates. Recently, different computational methods have been developed to predict the cleavage sites of caspase substrates with varying degrees of success. As the support vector machines (SVM) algorithm has been shown to be useful in several biological classification problems, we have implemented an SVM-based method to investigate its applicability to this domain.

RESULTS

A set of unique caspase substrates cleavage sites were obtained from literature and used for evaluating the SVM method. Datasets containing (i) the tetrapeptide cleavage sites, (ii) the tetrapeptide cleavage sites, augmented by two adjacent residues, P1' and P2' amino acids and (iii) the tetrapeptide cleavage sites with ten additional upstream and downstream flanking sequences (where available) were tested. The SVM method achieved an accuracy ranging from 81.25% to 97.92% on independent test sets. The SVM method successfully predicted the cleavage of a novel caspase substrate and its mutants.

CONCLUSION

This study presents an SVM approach for predicting caspase substrate cleavage sites based on the cleavage sites and the downstream and upstream flanking sequences. The method shows an improvement over existing methods and may be useful for predicting hitherto undiscovered cleavage sites.

Caspase substrates

The relatively common occurrence of sequences within proteins that match the consensus substrate specificity of caspases in intracellular proteins suggests a multitude of substrates in vivo - somewhere in the order of several hundred in humans alone. Indeed, the list of proteins that are reported to be cleaved by caspases in vitro proliferates rapidly. However, only a few of these proteins have been rigorously established as biologically or pathologically relevant, bona fide substrates in vivo. Many of them probably simply represent 'innocent bystanders' or erroneous assignments. In this review we discuss concepts of caspase substrate recognition and specificity, give resources for the discovery and annotation of caspase substrates, and highlight some specific human or mouse proteins where there is strong evidence for biologic or pathologic relevance.

PoPS: a computational tool for modeling and predicting protease specificity

Proteases play a fundamental role in the control of intra- and extracellular processes by binding and cleaving specific amino acid sequences. Identifying these targets is extremely challenging. Current computational attempts to predict cleavage sites are limited, representing these amino acid sequences as patterns or frequency matrices. Here we present PoPS, a publicly accessible bioinformatics tool (http://pops.csse.monash.edu.au/) which provides a novel method for building computational models of protease specificity that, while still being based on these amino acid sequences, can be built from any experimental data or expert knowledge available to the user. PoPS specificity models can be used to predict and rank likely cleavages within a single substrate, and within entire proteomes. Other factors, such as the secondary or tertiary structure of the substrate, can be used to screen unlikely sites. Furthermore, the tool also provides facilities to infer, compare and test models, and to store them in a publicly accessible database.

GraBCas: a bioinformatics tool for score-based prediction of Caspase- and Granzyme B-cleavage sites in protein sequences

Caspases and granzyme B are proteases that share the primary specificity to cleave at the carboxyl terminal of aspartate residues in their substrates. Both, caspases and granzyme B are enzymes that are involved in fundamental cellular processes and play a central role in apoptotic cell death. Although various targets are described, many substrates still await identification and many cleavage sites of known substrates are not identified or experimentally verified. A more comprehensive knowledge of caspase and granzyme B substrates is essential to understand the biological roles of these enzymes in more detail. The relatively high variability in cleavage site recognition sequence often complicates the identification of cleavage sites. As of yet there is no software available that allows identification of caspase and/or granzyme with cleavage sites differing from the consensus sequence. Here, we present a bioinformatics tool ‘GraBCas’ that provides score-based prediction of potential cleavage sites for the caspases 1–9 and granzyme B including an estimation of the fragment size. We tested GraBCas on already known substrates and showed its usefulness for protein sequence analysis. GraBCas is available at .