Molecular Interactions Required for Activation of Complement Component C2 Include Exosites Located on the Serine Protease Domain of C1s and Mannose-Binding Lectin Associated Protease-2.

The activation of the CP/LP C3 proconvertase complex is a key event in complement activation and involves cleavage of C4 and C2 by the C1s protease (classical pathway) or the mannose-binding lectin-associated serine protease (MASP)-2 (lectin pathway). Efficient cleavage of C4 by C1s and MASP-2 involves exosites on the complement control protein and serine protease (SP) domains of the proteases. The complement control protein domain exosite is not involved in cleavage of C2 by the proteases, but the role of an anion-binding exosite (ABE) on the SP domains of the proteases has (to our knowledge) never been investigated. In this study, we have shown that the ABE on the SP of both C1s and MASP-2 is crucial for efficient cleavage of C2, with mutant forms of the proteases greatly impaired in their rate of cleavage of C2. We have additionally shown that the site of binding for the ABE of the proteases is very likely to be located on the von Willebrand factor domain of C2, with the precise area differing between the enzymes: whereas C1s requires two anionic clusters on the von Willebrand factor domain to enact efficient cleavage of C2, MASP-2 apparently only requires one. These data provide (to our knowledge) new information about the molecular determinants for efficient activation of C2 by C1s and MASP-2. The enhanced view of the molecular events underlying the early stages of complement activation provides further possible intervention points for control of this activation that is involved in a number of inflammatory diseases.

Effect of a protease-activated receptor-2 antagonist (GB88) on inflammation-related loss of alveolar bone in periodontal disease.

BACKGROUND AND OBJECTIVE: Protease-activated receptor-2 (PAR2 ), a pro-inflammatory G-protein coupled receptor, has been associated with pathogenesis of periodontitis and the resulting bone loss caused by oral pathogens, including the keystone pathogen Porphyromonas gingivalis (P. gingivalis). We hypothesised that administration of a PAR2 antagonist, GB88, might prevent inflammation and subsequent alveolar bone resorption in a mouse model of periodontal disease. METHODS: Periodontitis was induced in mice by oral inoculations with P. gingivalis for a total of eight times over 24 days. The infected mice were treated with either GB88 or vehicle for the duration of the trial. Following euthanasia on day 56, serum was collected and used for the detection of mast cell tryptase. The right maxillae were defleshed and stained with methylene blue to measure the exposed cementum in molar teeth. The left maxillae were prepared for cryosections followed by staining for tartrate-resistant acid phosphatase to identify osteoclasts or with toluidine blue to identify mast cells. Reverse transcription quantitative PCR (RT-qPCR) was used to quantify the expression of inflammatory cytokines in the gingival tissue. Supernatants of T-lymphocyte cultures isolated from the regional lymph nodes were assayed using a cytometric bead array to measure the Th1/Th2/Th17 cytokine levels. RESULTS: Measurement of the exposed cementum showed that GB88 reduced P. gingivalis-induced alveolar bone loss by up to 69%. GB88 also prevented the increase in osteoclast numbers observed in the infected mice. Serum tryptase levels were significantly elevated in both the infected groups, and not altered by treatment. RT-qPCR showed that GB88 prevented the upregulation of Il1b, Il6, Ifng and Cd11b. In T-lymphocyte supernatants, only IFNγ and IL-17A levels were increased in response to infection, but this was prevented by GB88 treatment. CONCLUSIONS: GB88 significantly reduced osteoclastic alveolar bone loss in mice infected with P. gingivalis, seemingly by preventing the upregulation of several inflammatory cytokines. PAR2 antagonism may be an effective treatment strategy for periodontal disease.

Protease-associated import systems are widespread in Gram-negative bacteria.

Bacteria have evolved sophisticated uptake machineries in order to obtain the nutrients required for growth. Gram-negative plant pathogens of the genus Pectobacterium obtain iron from the protein ferredoxin, which is produced by their plant hosts. This iron-piracy is mediated by the ferredoxin uptake system (Fus), a gene cluster encoding proteins that transport ferredoxin into the bacterial cell and process it proteolytically. In this work we show that gene clusters related to the Fus are widespread in bacterial species. Through structural and biochemical characterisation of the distantly related Fus homologues YddB and PqqL from Escherichia coli, we show that these proteins are analogous to components of the Fus from Pectobacterium. The membrane protein YddB shares common structural features with the outer membrane ferredoxin transporter FusA, including a large extracellular substrate binding site. PqqL is an active protease with an analogous periplasmic localisation and iron-dependent expression to the ferredoxin processing protease FusC. Structural analysis demonstrates that PqqL and FusC share specific features that distinguish them from other members of the M16 protease family. Taken together, these data provide evidence that protease associated import systems analogous to the Fus are widespread in Gram-negative bacteria.

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.

PROSPERous: high-throughput prediction of substrate cleavage sites for 90 proteases with improved accuracy.

Summary: Proteases are enzymes that specifically cleave the peptide backbone of their target proteins. As an important type of irreversible post-translational modification, protein cleavage underlies many key physiological processes. When dysregulated, proteases' actions are associated with numerous diseases. Many proteases are highly specific, cleaving only those target substrates that present certain particular amino acid sequence patterns. Therefore, tools that successfully identify potential target substrates for proteases may also identify previously unknown, physiologically relevant cleavage sites, thus providing insights into biological processes and guiding hypothesis-driven experiments aimed at verifying protease-substrate interaction. In this work, we present PROSPERous, a tool for rapid in silico prediction of protease-specific cleavage sites in substrate sequences. Our tool is based on logistic regression models and uses different scoring functions and their pairwise combinations to subsequently predict potential cleavage sites. PROSPERous represents a state-of-the-art tool that enables fast, accurate and high-throughput prediction of substrate cleavage sites for 90 proteases. Availability and implementation: http://prosperous.erc.monash.edu/. Contact: jiangning.song@monash.edu or geoff.webb@monash.edu or r.pike@latrobe.edu.au. Supplementary information: Supplementary data are available at Bioinformatics online.

Keratinocyte-specific ablation of protease-activated receptor 2 prevents gingival inflammation and bone loss in a mouse model of periodontal disease.

Chronic periodontitis is characterised by gingival inflammation and alveolar bone loss. A major aetiological agent is Porphyromonas gingivalis, which secretes proteases that activate protease-activated receptor 2 (PAR2 ). PAR2 expressed on oral keratinocytes is activated by proteases released by P. gingivalis, inducing secretion of interleukin 6 (IL-6), and global knockout of PAR2 prevents bone loss and inflammation in a periodontal disease model in mice. To test the hypothesis that PAR2 expressed on gingival keratinocytes is required for periodontal disease pathology, keratinocyte-specific PAR2 -null mice were generated using K14-Cre targeted deletion of the PAR2 gene (F2rl1). These mice were subjected to a model of periodontitis involving placement of a ligature around a tooth, combined with P. gingivalis infection ("Lig + Inf"). The intervention caused a significant 44% decrease in alveolar bone volume (assessed by microcomputed tomography) in wildtype (K14-Cre:F2rl1wt/wt ), but not littermate keratinocyte-specific PAR2 -null (K14-Cre:F2rl1fl/fl ) mice. Keratinocyte-specific ablation of PAR2 prevented the significant Lig + Inf-induced increase (2.8-fold) in the number of osteoclasts in alveolar bone and the significant up-regulation (2.4-4-fold) of the inflammatory markers IL-6, IL-1ß, interferon-γ, myeloperoxidase, and CD11b in gingival tissue. These data suggest that PAR2 expressed on oral epithelial cells is a critical regulator of periodontitis-induced bone loss and will help in designing novel therapies with which to treat the disease.

The Structural Basis for Complement Inhibition by Gigastasin, a Protease Inhibitor from the Giant Amazon Leech.

Complement is crucial to the immune response, but dysregulation of the system causes inflammatory disease. Complement is activated by three pathways: classical, lectin, and alternative. The classical and lectin pathways are initiated by the C1r/C1s (classical) and MASP-1/MASP-2 (lectin) proteases. Given the role of complement in disease, there is a requirement for inhibitors to control the initiating proteases. In this article, we show that a novel inhibitor, gigastasin, from the giant Amazon leech, potently inhibits C1s and MASP-2, whereas it is also a good inhibitor of MASP-1. Gigastasin is a poor inhibitor of C1r. The inhibitor blocks the active sites of C1s and MASP-2, as well as the anion-binding exosites of the enzymes via sulfotyrosine residues. Complement deposition assays revealed that gigastasin is an effective inhibitor of complement activation in vivo, especially for activation via the lectin pathway. These data suggest that the cumulative effects of inhibiting both MASP-2 and MASP-1 have a greater effect on the lectin pathway than the more potent inhibition of only C1s of the classical pathway.

Polyphosphate is a novel cofactor for regulation of complement by a serpin, C1 inhibitor.

The complement system plays a key role in innate immunity, inflammation, and coagulation. The system is delicately balanced by negative regulatory mechanisms that modulate the host response to pathogen invasion and injury. The serpin, C1-esterase inhibitor (C1-INH), is the only known plasma inhibitor of C1s, the initiating serine protease of the classical pathway of complement. Like other serpin-protease partners, C1-INH interaction with C1s is accelerated by polyanions such as heparin. Polyphosphate (polyP) is a naturally occurring polyanion with effects on coagulation and complement. We recently found that polyP binds to C1-INH, prompting us to consider whether polyP acts as a cofactor for C1-INH interactions with its target proteases. We show that polyP dampens C1s-mediated activation of the classical pathway in a polymer length- and concentration-dependent manner by accelerating C1-INH neutralization of C1s cleavage of C4 and C2. PolyP significantly increases the rate of interaction between C1s and C1-INH, to an extent comparable to heparin, with an exosite on the serine protease domain of the enzyme playing a major role in this interaction. In a serum-based cell culture system, polyP significantly suppressed C4d deposition on endothelial cells, generated via the classical and lectin pathways. Moreover, polyP and C1-INH colocalize in activated platelets, suggesting that their interactions are physiologically relevant. In summary, like heparin, polyP is a naturally occurring cofactor for the C1s:C1-INH interaction and thus an important regulator of complement activation. The findings may provide novel insights into mechanisms underlying inflammatory diseases and the development of new therapies.

Assembly of the type II secretion system such as found in Vibrio cholerae depends on the novel Pilotin AspS.

The Type II Secretion System (T2SS) is a molecular machine that drives the secretion of fully-folded protein substrates across the bacterial outer membrane. A key element in the machinery is the secretin: an integral, multimeric outer membrane protein that forms the secretion pore. We show that three distinct forms of T2SSs can be distinguished based on the sequence characteristics of their secretin pores. Detailed comparative analysis of two of these, the Klebsiella-type and Vibrio-type, showed them to be further distinguished by the pilotin that mediates their transport and assembly into the outer membrane. We have determined the crystal structure of the novel pilotin AspS from Vibrio cholerae, demonstrating convergent evolution wherein AspS is functionally equivalent and yet structurally unrelated to the pilotins found in Klebsiella and other bacteria. AspS binds to a specific targeting sequence in the Vibrio-type secretins, enhances the kinetics of secretin assembly, and homologs of AspS are found in all species of Vibrio as well those few strains of Escherichia and Shigella that have acquired a Vibrio-type T2SS.

The protease cathepsin L regulates Th17 cell differentiation.

Previously we reported that IL-17(+) T cells, primarily IL-17(+) γδ cells, are increased in mice lacking the protease inhibitor serpinB1 (serpinb1(-/-) mice). Here we show that serpinB1-deficient CD4 cells exhibit a cell-autonomous and selective deficiency in suppressing T helper 17 (Th17) cell differentiation. This suggested an opposing role for one or more protease in promoting Th17 differentiation. We found that several SerpinB1-inhibitable cysteine cathepsins are induced in Th17 cells, most prominently cathepsin L (catL); this was verified by peptidase assays, active site labeling and Western blots. Moreover, Th17 differentiation was suppressed by both broad cathepsin inhibitors and catL selective inhibitors. CatL is present in Th17 cells as single chain (SC)- and two-chain (TC)-forms. Inhibiting asparagine endopeptidase (AEP) blocked conversion of SC-catL to TC-catL and increased generation of serpinb1(-/-) Th17 cells, but not wild-type Th17 cells. These findings suggest that SC-catL is biologically active in promoting Th17 generation and is counter-regulated by serpinB1 and secondarily by AEP. Thus, in addition to regulation by cytokines and transcription factors, differentiation of CD4 cells to Th17 cells is actively regulated by a catL-serpinB1-AEP module. Targeting this protease regulatory module could be an approach to treating Th17 cell-driven autoimmune disorders.

The x-ray crystal structure of mannose-binding lectin-associated serine proteinase-3 reveals the structural basis for enzyme inactivity associated with the Carnevale, Mingarelli, Malpuech, and Michels (3MC) syndrome.

The mannose-binding lectin associated-protease-3 (MASP-3) is a member of the lectin pathway of the complement system, a key component of human innate and active immunity. Mutations in MASP-3 have recently been found to be associated with Carnevale, Mingarelli, Malpuech, and Michels (3MC) syndrome, a severe developmental disorder manifested by cleft palate, intellectual disability, and skeletal abnormalities. However, the molecular basis for MASP-3 function remains to be understood. Here we characterize the substrate specificity of MASP-3 by screening against a combinatorial peptide substrate library. Through this approach, we successfully identified a peptide substrate that was 20-fold more efficiently cleaved than any other identified to date. Furthermore, we demonstrated that mutant forms of the enzyme associated with 3MC syndrome were completely inactive against this substrate. To address the structural basis for this defect, we determined the 2.6-Å structure of the zymogen form of the G666E mutant of MASP-3. These data reveal that the mutation disrupts the active site and perturbs the position of the catalytic serine residue. Together, these insights into the function of MASP-3 reveal how a mutation in this enzyme causes it to be inactive and thus contribute to the 3MC syndrome.

Molecular determinants of the substrate specificity of the complement-initiating protease, C1r.

The serine protease, C1r, initiates activation of the classical pathway of complement, which is a crucial innate defense mechanism against pathogens and altered-self cells. C1r both autoactivates and subsequently cleaves and activates C1s. Because complement is implicated in many inflammatory diseases, an understanding of the interaction between C1r and its target substrates is required for the design of effective inhibitors of complement activation. Examination of the active site specificity of C1r using phage library technology revealed clear specificity for Gln at P2 and Ile at P1', which are found in these positions in physiological substrates of C1r. Removal of one or both of the Gln at P2 and Ile at P1' in the C1s substrate reduced the rate of C1r activation. Substituting a Gln residue into the P2 of the activation site of MASP-3, a protein with similar domain structure to C1s that is not normally cleaved by C1r, enabled efficient activation of this enzyme. Molecular dynamics simulations and structural modeling of the interaction of the C1s activation peptide with the active site of C1r revealed the molecular mechanisms that particularly underpin the specificity of the enzyme for the P2 Gln residue. The complement control protein domains of C1r also made important contributions to efficient activation of C1s by this enzyme, indicating that exosite interactions were also important. These data show that C1r specificity is well suited to its cleavage targets and that efficient cleavage of C1s is achieved through both active site and exosite contributions.

A molecular switch governs the interaction between the human complement protease C1s and its substrate, complement C4.

The complement system is an ancient innate immune defense pathway that plays a front line role in eliminating microbial pathogens. Recognition of foreign targets by antibodies drives sequential activation of two serine proteases, C1r and C1s, which reside within the complement Component 1 (C1) complex. Active C1s propagates the immune response through its ability to bind and cleave the effector molecule complement Component 4 (C4). Currently, the precise structural and biochemical basis for the control of the interaction between C1s and C4 is unclear. Here, using surface plasmon resonance, we show that the transition of the C1s zymogen to the active form is essential for C1s binding to C4. To understand this, we determined the crystal structure of a zymogen C1s construct (comprising two complement control protein (CCP) domains and the serine protease (SP) domain). These data reveal that two loops (492-499 and 573-580) in the zymogen serine protease domain adopt a conformation that would be predicted to sterically abrogate C4 binding. The transition from zymogen to active C1s repositions both loops such that they would be able to interact with sulfotyrosine residues on C4. The structure also shows the junction of the CCP1 and CCP2 domains of C1s for the first time, yielding valuable information about the exosite for C4 binding located at this position. Together, these data provide a structural explanation for the control of the interaction with C1s and C4 and, furthermore, point to alternative strategies for developing therapeutic approaches for controlling activation of the complement cascade.

Identification of a catalytic exosite for complement component C4 on the serine protease domain of C1s.

The classical pathway of complement is crucial to the immune system, but it also contributes to inflammatory diseases when dysregulated. Binding of the C1 complex to ligands activates the pathway by inducing autoactivation of associated C1r, after which C1r activates C1s. C1s cleaves complement component C4 and then C2 to cause full activation of the system. The interaction between C1s and C4 involves active site and exosite-mediated events, but the molecular details are unknown. In this study, we identified four positively charged amino acids on the serine protease domain that appear to form a catalytic exosite that is required for efficient cleavage of C4. These residues are coincidentally involved in coordinating a sulfate ion in the crystal structure of the protease. Together with other evidence, this pointed to the involvement of sulfate ions in the interaction with the C4 substrate, and we showed that the protease interacts with a peptide from C4 containing three sulfotyrosine residues. We present a molecular model for the interaction between C1s and C4 that provides support for the above data and poses questions for future research into this aspect of complement activation.

Total synthesis of homogeneous variants of hirudin P6: a post-translationally modified anti-thrombotic leech-derived protein.

Hirudin P6 is a leech-derived anti-thrombotic protein which possesses two post-translational modifications, O-glycosylation and tyrosine sulfation. In this study we report the ligation-based synthesis of a library of hirudin P6 proteins possessing homogeneous glycosylation and sulfation modifications. The nature of the modifications incorporated was shown to have a drastic effect on inhibition against both the fibrinogenolytic and amidolytic activities of thrombin and thus highlights a potential means for attenuating the biological activity of the protein.

Mannose-binding lectin serine proteases and associated proteins of the lectin pathway of complement: two genes, five proteins and many functions?

The lectin pathway of the complement system is activated following the binding of carbohydrate-based ligands by recognition molecules such as mannose-binding lectin (MBL) or ficolins. Engagement of the recognition molecules causes activation of associated MBL-associated serine proteases or MASPs, which in turn activate downstream complement molecules to activate the system. Two MASP genes are alternatively spliced during expression to yield 5 proteins, including three proteases (MASP-1, -2 and -3) and two truncated proteins, MAp19 and MAp44. Here we discuss what is currently known about these proteins in terms of their structure and function. MASP-2 is autoactivated following the initial binding events of the pathway and is able to subsequently activate the C4 and C2 substrates required to activate the rest of the pathway. MASP-1 is able to augment MASP-2 activation, but also appears to play other roles, although the physiological significance of these is not yet clear. The roles of the truncated Map19 and Map44 proteins and the MASP-3 protease are currently unknown. The proteases form an interesting sub-family of proteins that clearly should be the focus of future research in order to establish their biological roles. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

S1 pocket of a bacterially derived subtilisin-like protease underpins effective tissue destruction.

The ovine footrot pathogen, Dichelobacter nodosus, secretes three subtilisin-like proteases that play an important role in the pathogenesis of footrot through their ability to mediate tissue destruction. Virulent and benign strains of D. nodosus secrete the basic proteases BprV and BprB, respectively, with the catalytic domain of these enzymes having 96% sequence identity. At present, it is not known how sequence variation between these two putative virulence factors influences their respective biological activity. We have determined the high resolution crystal structures of BprV and BprB. These data reveal that that the S1 pocket of BprV is more hydrophobic but smaller than that of BprB. We show that BprV is more effective than BprB in degrading extracellular matrix components of the host tissue. Mutation of two residues around the S1 pocket of BprB to the equivalent residues in BprV dramatically enhanced its proteolytic activity against elastin substrates. Application of a novel approach for profiling substrate specificity, the Rapid Endopeptidase Profiling Library (REPLi) method, revealed that both enzymes prefer cleaving after hydrophobic residues (and in particular P1 leucine) but that BprV has more restricted primary substrate specificity than BprB. Furthermore, for P1 Leu-containing substrates we found that BprV is a significantly more efficient enzyme than BprB. Collectively, these data illuminate how subtle changes in D. nodosus proteases may significantly influence tissue destruction as part of the ovine footrot pathogenesis process.

The subtilisin-like protease AprV2 is required for virulence and uses a novel disulphide-tethered exosite to bind substrates.

Many bacterial pathogens produce extracellular proteases that degrade the extracellular matrix of the host and therefore are involved in disease pathogenesis. Dichelobacter nodosus is the causative agent of ovine footrot, a highly contagious disease that is characterized by the separation of the hoof from the underlying tissue. D. nodosus secretes three subtilisin-like proteases whose analysis forms the basis of diagnostic tests that differentiate between virulent and benign strains and have been postulated to play a role in virulence. We have constructed protease mutants of D. nodosus; their analysis in a sheep virulence model revealed that one of these enzymes, AprV2, was required for virulence. These studies challenge the previous hypothesis that the elastase activity of AprV2 is important for disease progression, since aprV2 mutants were virulent when complemented with aprB2, which encodes a variant that has impaired elastase activity. We have determined the crystal structures of both AprV2 and AprB2 and characterized the biological activity of these enzymes. These data reveal that an unusual extended disulphide-tethered loop functions as an exosite, mediating effective enzyme-substrate interactions. The disulphide bond and Tyr92, which was located at the exposed end of the loop, were functionally important. Bioinformatic analyses suggested that other pathogenic bacteria may have proteases that utilize a similar mechanism. In conclusion, we have used an integrated multidisciplinary combination of bacterial genetics, whole animal virulence trials in the original host, biochemical studies, and comprehensive analysis of crystal structures to provide the first definitive evidence that the extracellular secreted proteases produced by D. nodosus are required for virulence and to elucidate the molecular mechanism by which these proteases bind to their natural substrates. We postulate that this exosite mechanism may be used by proteases produced by other bacterial pathogens of both humans and animals.

GABA production by glutamic acid decarboxylase is regulated by a dynamic catalytic loop.

Gamma-aminobutyric acid (GABA) is synthesized by two isoforms of the pyridoxal 5'-phosphate-dependent enzyme glutamic acid decarboxylase (GAD65 and GAD67). GAD67 is constitutively active and is responsible for basal GABA production. In contrast, GAD65, an autoantigen in type I diabetes, is transiently activated in response to the demand for extra GABA in neurotransmission, and cycles between an active holo form and an inactive apo form. We have determined the crystal structures of N-terminal truncations of both GAD isoforms. The structure of GAD67 shows a tethered loop covering the active site, providing a catalytic environment that sustains GABA production. In contrast, the same catalytic loop is inherently mobile in GAD65. Kinetic studies suggest that mobility in the catalytic loop promotes a side reaction that results in cofactor release and GAD65 autoinactivation. These data reveal the molecular basis for regulation of GABA homeostasis.

Discovery of amino acid motifs for thrombin cleavage and validation using a model substrate.

Understanding the active site preferences of an enzyme is critical to the design of effective inhibitors and to gaining insights into its mechanisms of action on substrates. While the subsite specificity of thrombin is understood, it is not clear whether the enzyme prefers individual amino acids at each subsite in isolation or prefers to cleave combinations of amino acids as a motif. To investigate whether preferred peptide motifs for cleavage could be identified for thrombin, we exposed a phage-displayed peptide library to thrombin. The resulting preferentially cleaved substrates were analyzed using the technique of association rule discovery. The results revealed that thrombin selected for amino acid motifs in cleavage sites. The contribution of these hypothetical motifs to substrate cleavage efficiency was further investigated using the B1 IgG-binding domain of streptococcal protein G as a model substrate. Introduction of a P(2)-P(1)' LRS thrombin cleavage sequence within a major loop of the protein led to cleavage of the protein by thrombin, with the cleavage efficiency increasing with the length of the loop. Introduction of further P(3)-P(1) and P(1)-P(1)'-P(3)' amino acid motifs into the loop region yielded greater cleavage efficiencies, suggesting that the susceptibility of a protein substrate to cleavage by thrombin is influenced by these motifs, perhaps because of cooperative effects between subsites closest to the scissile peptide bond.