Methods for mapping protease specificity

The study of protease specificity provides information on active-site structure and function, protein-protein interaction, regulation of intracellular and extracellular pathways, and evolution of protease and substrate genes. Peptide libraries that include fluorogenic and binding tags are often generated by solid-phase synthesis. Even larger explorations of cleavage site preferences employ positional scanning libraries or phage display. Microarrays enable presentation of individual peptides to proteases, DNA sequences for capture of peptide nucleic acid encoded peptides, or nanodroplets containing soluble peptide sequences. These new methods continue to inform the design of chemical inhibitors and the identification of substrates of proteases.     

Quantitative protease cleavage site profiling using tandem-mass-tag labeling and LC-MALDI-TOF/TOF MS/MS analysis

Knowledge of cleavage site specificity and activity are major prerequisites for understanding protease function. On the basis of a recently presented approach for proteomic identification of cleavage sites (PICS) in proteome-derived peptide libraries, we developed an isobaric labeling quantitative LC-MALDI-TOF/TOF MS/MS approach (Q-PICS) for simultaneous determination of cleavage site specificity and robust relative quantification of proteolytic events. For GluC-protease, 737 cleavage sites were identified in a yeast proteome-derived peptide library; 94.0% showed the typical GluC specificity for peptide bonds at glutamyl and aspartyl residues. The six-plex tandem mass tagging strategy allowed for three simultaneous replicates in a single run, guaranteeing high confidence and robust statistics for quantitative measurements. Using the quantitative capacity of Q-PICS, we performed a comparison of cleavage site specificity of GluC in two different buffer systems. The results support earlier findings describing that apparent difference between the buffer systems are probably caused by the inhibitory effect of bicarbonate on the overall GluC activity and that the preference for Glu-X bonds compared to Asp-X bonds is independent of the buffer system used.     

Proteome-derived, database-searchable peptide libraries for identifying protease cleavage sites

We introduce human proteome-derived, database-searchable peptide libraries for characterizing sequence-specific protein interactions. To identify endoprotease cleavage sites, we used peptides in such libraries with protected primary amines to simultaneously determine sequence preferences on the N-terminal (nonprime P) and C-terminal (prime P') sides of the scissile bond. Prime-side cleavage products were tagged with biotin, isolated and identified by tandem mass spectrometry, and the corresponding nonprime-side sequences were derived from human proteome databases using bioinformatics. Identification of hundreds to over 1,000 individual cleaved peptides allows the consensus protease cleavage site and subsite cooperativity to be readily determined from P6 to P6'. For the highly specific GluC protease, >95% of the 558 cleavage sites identified displayed the canonical selectivity. For the broad-specificity matrix metalloproteinase 2, >1,200 peptidic cleavage sites were identified. Profiling of HIV protease 1, caspase 3, caspase 7, cathepsins K and G, elastase and thrombin showed that this approach is broadly applicable to all mechanistic classes of endoproteases.     

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.     

Cleavage Specificity Analysis of Six Type II Transmembrane Serine Proteases (TTSPs) Using PICS with Proteome-Derived Peptide Libraries

Background

Type II transmembrane serine proteases (TTSPs) are a family of cell membrane tethered serine proteases with unclear roles as their cleavage site specificities and substrate degradomes have not been fully elucidated. Indeed just 52 cleavage sites are annotated in MEROPS, the database of proteases, their substrates and inhibitors.

Methodology/Principal Finding

To profile the active site specificities of the TTSPs, we applied Proteomic Identification of protease Cleavage Sites (PICS). Human proteome-derived database searchable peptide libraries were assayed with six human TTSPs (matriptase, matriptase-2, matriptase-3, HAT, DESC and hepsin) to simultaneously determine sequence preferences on the N-terminal non-prime (P) and C-terminal prime (P’) sides of the scissile bond. Prime-side cleavage products were isolated following biotinylation and identified by tandem mass spectrometry. The corresponding non-prime side sequences were derived from human proteome databases using bioinformatics. Sequencing of 2,405 individual cleaved peptides allowed for the development of the family consensus protease cleavage site specificity revealing a strong specificity for arginine in the P1 position and surprisingly a lysine in P1′ position. TTSP cleavage between R↓K was confirmed using synthetic peptides. By parsing through known substrates and known structures of TTSP catalytic domains, and by modeling the remainder, structural explanations for this strong specificity were derived.

Conclusions

Degradomics analysis of 2,405 cleavage sites revealed a similar and characteristic TTSP family specificity at the P1 and P1′ positions for arginine and lysine in unfolded peptides. The prime side is important for cleavage specificity, thus making these proteases unusual within the tryptic-enzyme class that generally has overriding non-prime side specificity.     

Assessment of substrate specificity of hepatitis G virus NS3 protease by a genetic method

The RNA genome of hepatitis G virus (HGV) encodes a large polyprotein that is processed to mature proteins by viral-encoded proteases. The HGV NS3 protease is responsible for the cleavage of the HGV polyprotein at four different locations. No conserved sequence motif has been identified for the cleavage sites of the NS3 protease. To determine the substrate specificity of the NS3 protease, amino acid sequences cleaved by the NS3 protease were obtained from randomized sequence libraries by using a screening method referred to as GASP (Genetic Assay for Site-specific Proteolysis). Based on statistical analyses of the obtained cleavable sequences, a consensus substrate sequence was deduced: Gln-Glu-Thr-Leu-Val downward arrow Ser, with the scissile bond located between Val and Ser. The relevance of this peptide as a cleavable substrate was further supported by molecular modeling of the NS3 protease. Our result would provide an insight on the molecular activity of the NS3 protease and may be useful for the design of substrate-based inhibitors.     

Altered Substrate Specificity of Drug-Resistant Human Immunodeficiency Virus Type 1 Protease

Resistance to human immunodeficiency virus type 1 protease (HIV PR) inhibitors results primarily from the selection of multiple mutations in the protease region. Because many of these mutations are selected for the ability to decrease inhibitor binding in the active site, they also affect substrate binding and potentially substrate specificity. This work investigates the substrate specificity of a panel of clinically derived protease inhibitor-resistant HIV PR variants. To compare protease specificity, we have used positional-scanning, synthetic combinatorial peptide libraries as well as a select number of individual substrates. The subsite preferences of wild-type HIV PR determined by using the substrate libraries are consistent with prior reports, validating the use of these libraries to compare specificity among a panel of HIV PR variants. Five out of seven protease variants demonstrated subtle differences in specificity that may have significant impacts on their abilities to function in viral maturation. Of these, four variants demonstrated up to fourfold changes in the preference for valine relative to alanine at position P2 when tested on individual peptide substrates. This change correlated with a common mutation in the viral NC/p1 cleavage site. These mutations may represent a mechanism by which severely compromised, drug-resistant viral strains can increase fitness levels. Understanding the altered substrate specificity of drug-resistant HIV PR should be valuable in the design of future generations of protease inhibitors as well as in elucidating the molecular basis of regulation of proteolysis in HIV.     

Factor Xa subsite mapping by proteome-derived peptide libraries improved using WebPICS, a resource for proteomic identification of cleavage sites

Proteomic identification of protease cleavage site specificity (PICS) is a recent proteomic approach for the easy mapping of protease subsite preferences that determines both the prime- and non-prime side specificity concurrently. Here we greatly facilitate user access by providing an automated and simple web-based data-analysis resource termed WebPics (http://clipserve.clip.ubc.ca/pics/). We demonstrate the utility of WebPics analysis of PICS data by determining the substrate specificity of factor Xa from P6-P6', an important blood coagulation protease that proteolytically generates thrombin from prothrombin. PICS confirms existing data on non-prime site specificity and refines our knowledge of factor Xa prime-site selectivity.     

In vivo selection of protease cleavage sites by using chimeric Sindbis virus libraries

Identifying protease cleavage sites contributes to our understanding of their specificity and biochemical properties and can help in designing specific inhibitors. One route to this end is the generation and screening of random libraries of cleavage sites. Both synthetic and phage-displayed libraries have been extensively used in vitro. We describe a novel system based on recombinant Sindbis virus which can be used to identify cleavage sites in vivo, thus eliminating the need for a purified enzyme and overcoming the problem of choosing the correct in vitro conditions. As a model we used the serine protease of the hepatitis C virus (HCV). We engineered the gene coding for this enzyme and two specific cleavage sites in the Sindbis virus structural gene and constructed libraries of viral genomes with a random sequence at either of the cleavage sites. The system was designed so that only viral genomes coding for sequences cleaved by the protease would produce viable viruses. With this system we selected viruses containing sequences mirroring those of the natural HCV protease substrates which were cleaved with comparable efficiencies.     

Active site specificity profiling of the matrix metalloproteinase family: Proteomic identification of 4300 cleavage sites by nine MMPs explored with structural and synthetic peptide cleavage analyses

Secreted and membrane tethered matrix metalloproteinases (MMPs) are key homeostatic proteases regulating the extracellular signaling and structural matrix environment of cells and tissues. For drug targeting of proteases, selectivity for individual molecules is highly desired and can be met by high yield active site specificity profiling. Using the high throughput Proteomic Identification of protease Cleavage Sites (PICS) method to simultaneously profile both the prime and non-prime sides of the cleavage sites of nine human MMPs, we identified more than 4300 cleavages from P6 to P6′ in biologically diverse human peptide libraries. MMP specificity and kinetic efficiency were mainly guided by aliphatic and aromatic residues in P1′ (with a ~ 32–93% preference for leucine depending on the MMP), and basic and small residues in P2′ and P3′, respectively. A wide differential preference for the hallmark P3 proline was found between MMPs ranging from 15 to 46%, yet when combined in the same peptide with the universally preferred P1′ leucine, an unexpected negative cooperativity emerged. This was not observed in previous studies, probably due to the paucity of approaches that profile both the prime and non-prime sides together, and the masking of subsite cooperativity effects by global heat maps and iceLogos. These caveats make it critical to check for these biologically highly important effects by fixing all 20 amino acids one-by-one in the respective subsites and thorough assessing of the inferred specificity logo changes. Indeed an analysis of bona fide MEROPS physiological substrate cleavage data revealed that of the 37 natural substrates with either a P3-Pro or a P1′-Leu only 5 shared both features, confirming the PICS data. Upon probing with several new quenched-fluorescent peptides, rationally designed on our specificity data, the negative cooperativity was explained by reduced non-prime side flexibility constraining accommodation of the rigidifying P3 proline with leucine locked in S1′. Similar negative cooperativity between P3 proline and the novel preference for asparagine in P1 cements our conclusion that non-prime side flexibility greatly impacts MMP binding affinity and cleavage efficiency. Thus, unexpected sequence cooperativity consequences were revealed by PICS that uniquely encompasses both the non-prime and prime sides flanking the proteomic-pinpointed scissile bond.     

Variable context Markov chains for HIV protease cleavage site prediction

Deciphering the knowledge of HIV protease specificity and developing computational tools for detecting its cleavage sites in protein polypeptide chain are very desirable for designing efficient and specific chemical inhibitors to prevent acquired immunodeficiency syndrome. In this study, we developed a generative model based on a generalization of variable order Markov chains (VOMC) for peptide sequences and adapted the model for prediction of their cleavability by certain proteases. The new method, called variable context Markov chains (VCMC), attempts to identify the context equivalence based on the evolutionary similarities between individual amino acids. It was applied for HIV-1 protease cleavage site prediction problem and shown to outperform existing methods in terms of prediction accuracy on a common dataset. In general, the method is a promising tool for prediction of cleavage sites of all proteases and encouraged to be used for any kind of peptide classification problem as well.     

One-bead, one-compound peptide library sequencing via high-pressure ammonia cleavage coupled to nanomanipulation/nanoelectrospray ionization mass spectrometry

Biological screening of one-bead, one-compound (OBOC) combinatorial peptide libraries is routinely carried out with the peptide remaining bound to the resin bead during screening. After a hit is identified, the bead is isolated, the peptide is cleaved from the bead, and its sequence is determined. We have developed a new technique for cleavage of peptides from resin beads whereby exposure of a 4-hydroxymethyl benzoic acid (HMBA)-linked peptide to high-pressure ammonia gas led to efficient cleavage in as little as 5 min. Here we also report a new method of extracting peptide from individual library beads for its introduction into a mass spectrometer that uses nanomanipulation combined with nanoelectrospray ionization mass spectrometry (NSI MS). Single beads analyzed by nanomanipulation/NSI MS were found to give identical MS results to those of bulk samples. Detection of 18 unique cleaved peptides 1 to 8 amino acids in length, and sequencing of 14 different peptide sequences 4 to 8 amino acids in length, was demonstrated on a combination of bulk samples and ones from individual beads of an OBOC library. The method was highly reproducible, with 100% of attempts to extract peptide resulting in high-quality MS data. This new collection of techniques allows rapid, reliable, environmentally responsible sequencing of hit beads from combinatorial peptide libraries.     

Proteomic analyses reveal an acidic prime side specificity for the astacin metalloprotease family reflected by physiological substrates

Astacins are secreted and membrane-bound metalloproteases with clear associations to many important pathological and physiological processes. Yet with only a few substrates described their biological roles are enigmatic. Moreover, the lack of knowledge of astacin cleavage site specificities hampers assay and drug development. Using PICS (proteomic identification of protease cleavage site specificity) and TAILS (terminal amine isotopic labeling of substrates) degradomics approaches >3000 cleavage sites were proteomically identified for five different astacins. Such broad coverage enables family-wide determination of specificities N- and C-terminal to the scissile peptide bond. Remarkably, meprin α, meprin β, and LAST_MAM proteases exhibit a strong preference for aspartate in the peptide (P)1' position because of a conserved positively charged residue in the active cleft subsite (S)1'. This unparalleled specificity has not been found for other families of extracellular proteases. Interestingly, cleavage specificity is also strongly influenced by proline in P2' or P3' leading to a rare example of subsite cooperativity. This specificity characterizes the astacins as unique contributors to extracellular proteolysis that is corroborated by known cleavage sites in procollagen I+III, VEGF (vascular endothelial growth factor)-A, IL (interleukin)-1β, and pro-kallikrein 7. Indeed, cleavage sites in VEGF-A and pro-kallikrein 7 identified by terminal amine isotopic labeling of substrates matched those reported by Edman degradation. Moreover, the novel substrate FGF-19 was validated biochemically and shown to exhibit altered biological activity after meprin processing.     

A cleavage enzyme-cytometric bead array provides biochemical profiling of resistance mutations in HIV-1 Gag and protease

Most protease-substrate assays rely on short, synthetic peptide substrates consisting of native or modified cleavage sequences. These assays are inadequate for interrogating the contribution of native substrate structure distal to a cleavage site that influences enzymatic cleavage or for inhibitor screening of native substrates. Recent evidence from HIV-1 isolates obtained from individuals resistant to protease inhibitors has demonstrated that mutations distal to or surrounding the protease cleavage sites in the Gag substrate contribute to inhibitor resistance. We have developed a protease-substrate cleavage assay, termed the cleavage enzyme- cytometric bead array (CE-CBA), which relies on native domains of the Gag substrate containing embedded cleavage sites. The Gag substrate is expressed as a fluorescent reporter fusion protein, and substrate cleavage can be followed through the loss of fluorescence utilizing flow cytometry. The CE-CBA allows precise determination of alterations in protease catalytic efficiency (k(cat)/K(M)) imparted by protease inhibitor resistance mutations in protease and/or gag in cleavage or noncleavage site locations in the Gag substrate. We show that the CE-CBA platform can identify HIV-1 protease present in cellular extractions and facilitates the identification of small molecule inhibitors of protease or its substrate Gag. Moreover, the CE-CBA can be readily adapted to any enzyme-substrate pair and can be utilized to rapidly provide assessment of catalytic efficiency as well as systematically screen for inhibitors of enzymatic processing of substrate.     

PACMANS: A bioinformatically informed algorithm to predict,design, and disrupt protease‐on‐protease hydrolysis

Multiple proteases in a system hydrolyze target substrates, but recent evidence indicates that some proteases will degrade other proteases as well. Cathepsin S hydrolysis of cathepsin K is one such example. These interactions may be uni‐ or bi‐directional and change the expected kinetics. To explore potential protease‐on‐protease interactions in silico, a program was developed for users to input two proteases: (1) the protease‐ase that hydrolyzes (2) the substrate, protease. This program identifies putative sites on the substrate protease highly susceptible to cleavage by the protease‐ase, using a sliding‐window approach that scores amino acid sequences by their preference in the protease‐ase active site, culled from MEROPS database. We call this PACMANS, Protease‐Ase Cleavage from MEROPS ANalyzed Specificities, and test and validate this algorithm with cathepsins S and K. PACMANS cumulative likelihood scoring identified L253 and V171 as sites on cathepsin K subject to cathepsin S hydrolysis. Mutations made at these locations were tested to block hydrolysis and validate PACMANS predictions. L253A and L253V cathepsin K mutants significantly reduced cathepsin S hydrolysis, validating PACMANS unbiased identification of these sites. Interfamilial protease interactions between cathepsin S and MMP‐2 or MMP‐9 were tested after predictions by PACMANS, confirming its utility for these systems as well. PACMANS is unique compared to other putative site cleavage programs by allowing users to define the proteases of interest and target, and can also be employed for non‐protease substrate proteins, as well as short peptide sequences.     

Substrate screening identifies a novel target sequence for the proteasomal activity regulated by ionizing radiation

The screening for treatment‐induced enzyme activities offers the opportunity to discover important regulatory mechanisms and the identification of potential targets for anti‐cancer therapies. A novel screening technique was applied to screen substrate peptide sequences for proteolytic activities up‐ or down‐regulated by ionizing radiation in tumor cells. One specific substrate sequence was cleaved in control cell extracts but to a smaller extent in irradiated cell extracts and investigated in detail. Based on protease‐class‐specific inhibitory studies and cleavage site analysis a potent warhead‐inhibitor was synthesized and used to identify the proteasome as the protease of interest. The investigated sequence shows high homology to a regulatory site of nucleoporin 50, an element of the nuclear pore complex, and site specific cleavage of nucleoporin 50 was determined in vitro suggesting a novel link between the ionizing radiation‐regulated proteasome and nuclear protein shuttling.     

Identification of inhibitors using a cell-based assay for monitoring Golgi-resident protease activity

Noninvasive real-time quantification of cellular protease activity allows monitoring of enzymatic activity and identification of activity modulators within the protease's natural milieu. We developed a protease activity assay based on differential localization of a recombinant reporter consisting of a Golgi retention signal and a protease cleavage sequence fused to alkaline phosphatase (AP). When expressed in mammalian cells, this protein localizes to Golgi bodies and, on protease-mediated cleavage, AP translocates to the extracellular medium where its activity is measured. We used this system to monitor the Golgi-associated protease furin, a pluripotent enzyme with a key role in tumorigenesis, viral propagation of avian influenza, ebola, and HIV as well as in activation of anthrax, pseudomonas, and diphtheria toxins. This technology was adapted for high-throughput screening of 39,000-compound small molecule libraries, leading to identification of furin inhibitors. Furthermore, this strategy was used to identify inhibitors of another Golgi protease, the beta-site amyloid precursor protein (APP)-cleaving enzyme (BACE). BACE cleavage of the APP leads to formation of the Abeta peptide, a key event that leads to Alzheimer's disease. In conclusion, we describe a customizable noninvasive technology for real-time assessment of Golgi protease activity used to identify inhibitors of furin and BACE.     

Selection of a CXCR4 antagonist from a human heavy chain CDR3-derived phage library

Phage display technology is a powerful selection approach to identify strong and specific binders to a large variety of targets. In this study, we compared the efficacy of a phage library displaying human heavy chain complementarity determining region 3 (HCDR3) repertoires with a set of conventional random peptide libraries for the identification of CXCR4 antagonists using a peptide corresponding to the second extracellular loop of the receptor CXCR4 as target. A total of 11 selection campaigns on this target did not result in any specific ligand from the random peptide libraries. In contrast, a single selection campaign with an HCDR3 library derived from the IgM repertoire of a nonimmunized donor resulted in nine specific peptides with lengths ranging from 10 to 19 residues. Four of these HCDR3 sequences interacted with native receptor and the most frequently isolated peptide displayed an affinity of 5.6 μm and acted as a CXCR4 antagonist (IC(50) = 23 μm). To comprehend the basis of the highly efficient HCDR3 library selection, its biochemical properties were investigated. The HCDR3 length varied from 3 to 21 residues and displayed a biased amino acid content with a predominant proportion of Tyr, Gly, Ser and Asp. Repetitive and conserved motifs were observed in the majority of the HCDR3 sequences. The strength and efficacy of the HCDR3 libraries reside in the combination of multiple size peptides and a naturally biased sequence variation. Therefore, HCDR3 libraries represent a powerful and versatile alternative to fully randomized peptide libraries, in particular for difficult targets.     

Robust production of a peptide library using methodological synchronization

Peptide libraries have proven to be useful in applications such as substrate profiling, drug candidate screening and identifying protein–protein interaction partners. However, issues of fidelity, peptide length, and purity have been encountered when peptide libraries are chemically synthesized. Biochemically produced libraries, on the other hand, circumvent many of these issues due to the fidelity of the protein synthesis machinery. Using thioredoxin as an expression partner, a stably folded peptide scaffold (avian pancreatic polypeptide) and a compatible cleavage site for human rhinovirus 3C protease, we report a method that allows robust expression of a genetically encoded peptide library, which yields peptides of high purity. In addition, we report the use of methodological synchronization, an experimental design created for the production of a library, from initial cloning to peptide characterization, within a 5-week period of time. Total peptide yields ranged from 0.8% to 16%, which corresponds to 2–70 mg of pure peptide. Additionally, no correlation was observed between the ability to be expressed or overall yield of peptide-fusions and the intrinsic chemical characteristics of the peptides, indicating that this system can be used for a wide variety of peptide sequences with a range of chemical characteristics.