Rational design of serine protease inhibitors

Proteases regulate a spectrum of diverse physiological processes, and dysregulation of proteolytic activity drives a plethora of pathological conditions. Understanding protease function is essential to appreciating many aspects of normal physiology and progression of disease. Consequently, development of potent and specific inhibitors of proteolytic enzymes is vital to provide tools for the dissection of protease function in biological systems and for the treatment of diseases linked to aberrant proteolytic activity. The studies in this thesis describe the rational design of potent inhibitors of three proteases that are implicated in disease development. Additionally, key features of the interaction of proteases and their cognate inhibitors or substrates are analysed and a series of rational inhibitor design principles are expounded and tested. Rational design of protease inhibitors relies on a comprehensive understanding of protease structure and biochemistry. Analysis of known protease cleavage sites in proteins and peptides is a commonly used source of such information. However, model peptide substrate and protein sequences have widely differing levels of backbone constraint and hence can adopt highly divergent structures when binding to a protease’s active site. This may result in identical sequences in peptides and proteins having different conformations and diverse spatial distribution of amino acid functionalities. Regardless of this, protein and peptide cleavage sites are often regarded as being equivalent. One of the key findings in the following studies is a definitive demonstration of the lack of equivalence between these two classes of substrate and invalidation of the common practice of using the sequences of model peptide substrates to predict cleavage of proteins in vivo. Another important feature for protease substrate recognition is subsite cooperativity. This type of cooperativity is commonly referred to as protease or substrate binding subsite cooperativity and is distinct from allosteric cooperativity, where binding of a molecule distant from the protease active site affects the binding affinity of a substrate. Subsite cooperativity may be intramolecular where neighbouring residues in substrates are interacting, affecting the scissile bond’s susceptibility to protease cleavage. Subsite cooperativity can also be intermolecular where a particular residue’s contribution to binding affinity changes depending on the identity of neighbouring amino acids. Although numerous studies have identified subsite cooperativity effects, these findings are frequently ignored in investigations probing subsite selectivity by screening against diverse combinatorial libraries of peptides (positional scanning synthetic combinatorial library; PS-SCL). This strategy for determining cleavage specificity relies on the averaged rates of hydrolysis for an uncharacterised ensemble of peptide sequences, as opposed to the defined rate of hydrolysis of a known specific substrate. Further, since PS-SCL screens probe the preference of the various protease subsites independently, this method is inherently unable to detect subsite cooperativity. However, mean hydrolysis rates from PS-SCL screens are often interpreted as being comparable to those produced by single peptide cleavages. Before this study no large systematic evaluation had been made to determine the level of correlation between protease selectivity as predicted by screening against a library of combinatorial peptides and cleavage of individual peptides. This subject is specifically explored in the studies described here. In order to establish whether PS-SCL screens could accurately determine the substrate preferences of proteases, a systematic comparison of data from PS-SCLs with libraries containing individually synthesised peptides (sparse matrix library; SML) was carried out. These SML libraries were designed to include all possible sequence combinations of the residues that were suggested to be preferred by a protease using the PS-SCL method. SML screening against the three serine proteases kallikrein 4 (KLK4), kallikrein 14 (KLK14) and plasmin revealed highly preferred peptide substrates that could not have been deduced by PS-SCL screening alone. Comparing protease subsite preference profiles from screens of the two types of peptide libraries showed that the most preferred substrates were not detected by PS SCL screening as a consequence of intermolecular cooperativity being negated by the very nature of PS SCL screening. Sequences that are highly favoured as result of intermolecular cooperativity achieve optimal protease subsite occupancy, and thereby interact with very specific determinants of the protease. Identifying these substrate sequences is important since they may be used to produce potent and selective inhibitors of protolytic enzymes. This study found that highly favoured substrate sequences that relied on intermolecular cooperativity allowed for the production of potent inhibitors of KLK4, KLK14 and plasmin. Peptide aldehydes based on preferred plasmin sequences produced high affinity transition state analogue inhibitors for this protease. The most potent of these maintained specificity over plasma kallikrein (known to have a very similar substrate preference to plasmin). Furthermore, the efficiency of this inhibitor in blocking fibrinolysis in vitro was comparable to aprotinin, which previously saw clinical use to reduce perioperative bleeding. One substrate sequence particularly favoured by KLK4 was substituted into the 14 amino acid, circular sunflower trypsin inhibitor (SFTI). This resulted in a highly potent and selective inhibitor (SFTI-FCQR) which attenuated protease activated receptor signalling by KLK4 in vitro. Moreover, SFTI-FCQR and paclitaxel synergistically reduced growth of ovarian cancer cells in vitro, making this inhibitor a lead compound for further therapeutic development. Similar incorporation of a preferred KLK14 amino acid sequence into the SFTI scaffold produced a potent inhibitor for this protease. However, the conformationally constrained SFTI backbone enforced a different intramolecular cooperativity, which masked a KLK14 specific determinant. As a consequence, the level of selectivity achievable was lower than that found for the KLK4 inhibitor. Standard mechanism inhibitors such as SFTI rely on a stable acyl-enzyme intermediate for high affinity binding. This is achieved by a conformationally constrained canonical binding loop that allows for reformation of the scissile peptide bond after cleavage. Amino acid substitutions within the inhibitor to target a particular protease may compromise structural determinants that support the rigidity of the binding loop and thereby prevent the engineered inhibitor reaching its full potential. An in silico analysis was carried out to examine the potential for further improvements to the potency and selectivity of the SFTI-based KLK4 and KLK14 inhibitors. Molecular dynamics simulations suggested that the substitutions within SFTI required to target KLK4 and KLK14 had compromised the intramolecular hydrogen bond network of the inhibitor and caused a concomitant loss of binding loop stability. Furthermore in silico amino acid substitution revealed a consistent correlation between a higher frequency of formation and the number of internal hydrogen bonds of SFTI-variants and lower inhibition constants. These predictions allowed for the production of second generation inhibitors with enhanced binding affinity toward both targets and highlight the importance of considering intramolecular cooperativity effects when engineering proteins or circular peptides to target proteases. The findings from this study show that although PS-SCLs are a useful tool for high throughput screening of approximate protease preference, later refinement by SML screening is needed to reveal optimal subsite occupancy due to cooperativity in substrate recognition. This investigation has also demonstrated the importance of maintaining structural determinants of backbone constraint and conformation when engineering standard mechanism inhibitors for new targets. Combined these results show that backbone conformation and amino acid cooperativity have more prominent roles than previously appreciated in determining substrate/inhibitor specificity and binding affinity. The three key inhibitors designed during this investigation are now being developed as lead compounds for cancer chemotherapy, control of fibrinolysis and cosmeceutical applications. These compounds form the basis of a portfolio of intellectual property which will be further developed in the coming years.

Mastering the canonical loop of serine protease inhibitors : enhancing potency by optimising the internal hydrogen bond network

Background Canonical serine protease inhibitors commonly bind to their targets through a rigid loop stabilised by an internal hydrogen bond network and disulfide bond(s). The smallest of these is sunflower trypsin inhibitor (SFTI-1), a potent and broad-range protease inhibitor. Recently, we re-engineered the contact β-sheet of SFTI-1 to produce a selective inhibitor of kallikrein-related peptidase 4 (KLK4), a protease associated with prostate cancer progression. However, modifications in the binding loop to achieve specificity may compromise structural rigidity and prevent re-engineered inhibitors from reaching optimal binding affinity. Methodology/Principal Findings In this study, the effect of amino acid substitutions on the internal hydrogen bonding network of SFTI were investigated using an in silico screen of inhibitor variants in complex with KLK4 or trypsin. Substitutions favouring internal hydrogen bond formation directly correlated with increased potency of inhibition in vitro. This produced a second generation inhibitor (SFTI-FCQR Asn14) which displayed both a 125-fold increased capacity to inhibit KLK4 (Ki = 0.0386±0.0060 nM) and enhanced selectivity over off-target serine proteases. Further, SFTI-FCQR Asn14 was stable in cell culture and bioavailable in mice when administered by intraperitoneal perfusion. Conclusion/Significance These findings highlight the importance of conserving structural rigidity of the binding loop in addition to optimising protease/inhibitor contacts when re-engineering canonical serine protease inhibitors.

Serine protease inhibition and mitochondrial dysfunction associated with cisplatin resistance in human tumor cell lines : targets for therapy

Indicators of mitochondrial function were studied in two different cell culture models of cis-diamminedichloroplatinum-II (CDDP) resistance: the intrinsically resistant human ovarian cancer cell line CI-80-13S, and resistant clones (HeLa-S1a and HeLa-S1b) generated by stable expression of the serine protease inhibitor—plasminogen activator inhibitor type-2 (PAI-2), in the human cervical cancer cell line HeLa. In both models, CDDP resistance was associated with sensitivity to killing by adriamycin, etoposide, auranofin, bis[1,2-bis(diphenylphosphino)ethane]gold(I) chloride {[Au(DPPE)2]Cl}, CdCl2 and the mitochondrial inhibitors rhodamine-123 (Rhl23), dequalinium chloride (DeCH), tetraphenylphosphonium (TPP), and ethidium bromide (EtBr) and with lower constitutive levels of ATP. Unlike the HeLa clones, CI-80-13S cells were additionally sensitive to chloramphenicol, 1-methyl-4-phenylpyridinium ion (MPP+), rotenone, thenoyltrifluoroacetone (TTFA), and antimycin A, and showed poor reduction of 1-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (MTT), suggesting a deficiency in NADH dehydrogenase and/or succinate dehydrogenase activities. Total platinum uptake and DNA-bound platinum were slightly lower in CI-80-13S than in sensitive cells. The HeLa-S1a and HeLa-S1b clones, on the other hand, showed poor reduction of triphenyltetrazolium chloride (TTC), indicative of low cytochrome c oxidase activity. Total platinum uptake by HeLa-S1a was similar to HeLa, but DNA-bound platinum was much lower than for the parent cell line. The mitochondria of CI-80-13S and HeLa-S1a showed altered morphology and were fewer in number than those of JAM and HeLa. In both models, CDDP resistance was associated with less platinum accumulation and with mitochondrial and membrane defects, brought about one case with expression of a protease inhibitor which is implicated in tumor progression. Such markers may identify tumors suitable for treatment with gold phosphine complexes or other mitochondrial inhibitors.

Serine protease inhibition and mitochondrial dysfunction associated with cisplatin resistance in human tumor cell lines: targets for therapy

Indicators of mitochondrial function were studied in two different cell culture models of cis-diamminedichloroplatinum-II (CDDP) resistance: the intrinsically resistant human ovarian cancer cell line CI-80-13S, and resistant clones (HeLa-S1a and HeLa-S1b) generated by stable expression of the serine protease inhibitor—plasminogen activator inhibitor type-2 (PAI-2), in the human cervical cancer cell line HeLa. In both models, CDDP resistance was associated with sensitivity to killing by adriamycin, etoposide, auranofin, bis[1,2-bis(diphenylphosphino)ethane]gold(I) chloride {[Au(DPPE)2]Cl}, CdCl2 and the mitochondrial inhibitors rhodamine-123 (Rhl23), dequalinium chloride (DeCH), tetraphenylphosphonium (TPP), and ethidium bromide (EtBr) and with lower constitutive levels of ATP. Unlike the HeLa clones, CI-80-13S cells were additionally sensitive to chloramphenicol, 1-methyl-4-phenylpyridinium ion (MPP+), rotenone, thenoyltrifluoroacetone (TTFA), and antimycin A, and showed poor reduction of 1-[4,5-dimethylthiazol-2-yl]-2,5-diphenyltetrazolium bromide (MTT), suggesting a deficiency in NADH dehydrogenase and/or succinate dehydrogenase activities. Total platinum uptake and DNA-bound platinum were slightly lower in CI-80-13S than in sensitive cells. The HeLa-S1a and HeLa-S1b clones, on the other hand, showed poor reduction of triphenyltetrazolium chloride (TTC), indicative of low cytochrome c oxidase activity. Total platinum uptake by HeLa-S1a was similar to HeLa, but DNA-bound platinum was much lower than for the parent cell line. The mitochondria of CI-80-13S and HeLa-S1a showed altered morphology and were fewer in number than those of JAM and HeLa. In both models, CDDP resistance was associated with less platinum accumulation and with mitochondrial and membrane defects, brought about one case with expression of a protease inhibitor which is implicated in tumor progression. Such markers may identify tumors suitable for treatment with gold phosphine complexes or other mitochondrial inhibitors.

Identification of a serine protease inhibitor which causes inclusion vacuole reduction and is lethal to Chlamydia trachomatis

The mechanistic details of the pathogenesis of Chlamydia, an obligate intracellular pathogen of global importance, have eluded scientists due to the scarcity of traditional molecular genetic tools to investigate this organism. Here we report a chemical biology strategy that has uncovered the first essential protease for this organism. Identification and application of a unique CtHtrA inhibitor (JO146) to cultures of Chlamydia resulted in a complete loss of viable elementary body formation. JO146 treatment during the replicative phase of development resulted in a loss of Chlamydia cell morphology, diminishing inclusion size, and ultimate loss of inclusions from the host cells. This completely prevented the formation of viable Chlamydia elementary bodies. In addition to its effect on the human C. trachomatis strain, JO146 inhibited the viability of the mouse strain, Chlamydia muridarum, both in vitro and in vivo. Thus, we report a chemical biology approach to establish an essential role for Chlamydia CtHtrA. The function of CtHtrA for Chlamydia appears to be essential for maintenance of cell morphology during replicative the phase and these findings provide proof of concept that proteases can be targetted for anti-microbial therapy for intracellular pathogens.

Murine, but not human, ephrin-B2 can be efficiently cleaved by the serine protease kallikrein-4: Implications for xenograft models of human prostate cancer

Background Ephrin-B2 is the sole physiologically-relevant ligand of the receptor tyrosine kinase EphB4, which is over-expressed in many epithelial cancers, including 66% of prostate cancers, and contributes to cancer cell survival, invasion and migration. Crucially, however, the cancer-promoting EphB4 signalling pathways are independent of interaction with its ligand ephrin-B2, as activation of ligand-dependent signalling causes tumour suppression. Ephrin-B2, however, is often found on the surface of endothelial cells of the tumour vasculature, where it can regulate angiogenesis to support tumour growth. Proteolytic cleavage of endothelial cell ephrin-B2 has previously been suggested as one mechanism whereby the interaction between tumour cell-expressed EphB4 and endothelial cell ephrin-B2 is regulated to support both cancer promotion and angiogenesis. Methods An in silico approach was used to search accessible surfaces of 3D protein models for cleavage sites for the key prostate cancer serine protease, KLK4, and this identified murine ephrin-B2 as a potential KLK4 substrate. Mouse ephrin-B2 was then confirmed as a KLK4 substrate by in vitro incubation of recombinant mouse ephrin-B2 with active recombinant human KLK4. Cleavage products were visualised by SDS-PAGE, silver staining and Western blot and confirmed by N-terminal sequencing. Results At low molar ratios, KLK4 cleaved murine ephrin-B2 but other prostate-specific KLK family members (KLK2 and KLK3/PSA) were less efficient, suggesting cleavage was KLK4-selective. The primary KLK4 cleavage site in murine ephrin-B2 was verified and shown to correspond to one of the in silico predicted sites between extracellular domain residues arginine 178 and asparagine 179. Surprisingly, the highly homologous human ephrin-B2 was poorly cleaved by KLK4 at these low molar ratios, likely due to the 3 amino acid differences at this primary cleavage site. Conclusion These data suggest that in in vivo mouse xenograft models, endogenous mouse ephrin-B2, but not human tumour ephrin-B2, may be a downstream target of cancer cell secreted human KLK4. This is a critical consideration when interpreting data from murine explants of human EphB4+/KLK4+ cancer cells, such as prostate cancer cells, where differential effects may be seen in mouse models as opposed to human clinical situations.

Molecular and functional characterizations of a Kunitz-type serine protease inhibitor FcKuSPI of the shrimp Fenneropenaeus chinensis

Serine proteinase inhibitors play important and diverse roles in biological processes such as coagulation, defense mechanisms, and immune responses. Here, we identified and characterized a Kunitz-type proteinase inhibitor, designated FcKuSPI, of the BPTI/Kunitz family of serine proteinase inhibitors from the hemocyte cDNA library of the shrimp Fenneropenaeus chinensis. The deduced amino acid sequence of FcKuSPI comprises 80 residues with a putative signal peptide of 15 amino acids. The predicted molecular weight of the mature peptide is 7.66 kDa and its predicted isoelectric point is 8.84. FcKuSPI includes a Kunitz domain containing six conserved cysteine residues that are predicted to form three disulfide bonds. FcKuSPI shares 44e53% homology with BPTI/Kunitz family members from other species. FcKuSPI mRNAwas expressed highly in the hemocytes and moderately in muscle in healthy shrimp. Recombinant FcKuSPI protein demonstrated anti-protease activity against trypsin and anticoagulant activity against citrated human plasma in a dose-dependent manner in in vitro assays.

Differential recognition by tick-resistant cattle of the recombinantly expressed Rhipicephalus microplus serine protease inhibitor-3 (RMS-3)

Rhipicephalus micro plus is an important bovine ectoparasite, widely distributed in tropical and subtropical regions of the world causing large economic losses to the cattle industry. Its success as an ectoparasite is associated with its capacity to disarm the antihemostatic and anti-inflammatory reactions of the host. Serpins are protease inhibitors with an important role in the modulation of host-parasite interactions. The cDNA that encodes for a R. microplus serpin was isolated by RACE and subsequently cloned into the pPICZ alpha A vector. Sequence analysis of the cDNA and predicted amino acid showed that this cDNA has a conserved serpin domain. B- and T-cell epitopes were predicted using bioinformatics tools. The recombinant R. microplus serpin (rRMS-3) was secreted into the culture media of Pichia pastoris after methanol induction at 0.2 mg l(-1) qRT-PCR expression analysis of tissues and life cycle stages demonstrated that RMS-3 was mainly expressed in the salivary glands of female adult ticks. Immunological recognition of the rRMS-3 and predicted B-cell epitopes was tested using tick-resistant and susceptible cattle sera. Only sera from tick-resistant bovines recognized the B-cell epitope AHYNPPPPIEFT (Seq7). The recombinant RMS-3 was expressed in P. pastoris, and ELISA screening also showed higher recognition by tick-resistant bovine sera. The results obtained suggest that RMS-3 is highly and specifically secreted into the bite site of R. microplus feeding on tick-resistant bovines. Capillary feeding of semi-engorged ticks with anti-AHYNPPPPIEFT sheep sera led to an 81.16% reduction in the reproduction capacity of R. microplus. Therefore, it is possible to conclude that R. microplus serpin (RMS-3) has an important role in the host-parasite interaction to overcome the immune responses in resistant cattle. (C) 2012 Elsevier GmbH. All rights reserved.

Rhipicephalus microplus serine protease inhibitor family: annotation, expression and functional characterisation assessment

Background: Rhipicephalus (Boophilus) microplus evades the host's haemostatic system through a complex protein array secreted into tick saliva. Serine protease inhibitors (serpins) conform an important component of saliva which are represented by a large protease inhibitor family in Ixodidae. These secreted and non-secreted inhibitors modulate diverse and essential proteases involved in different physiological processes. Methods: The identification of R. microplus serpin sequences was performed through a web-based bioinformatics environment called Yabi. The database search was conducted on BmiGi V1, BmiGi V2.1, five SSH libraries, Australian tick transcriptome libraries and RmiTR V1 using bioinformatics methods. Semi quantitative PCR was carried out using different adult tissues and tick development stages. The cDNA of four identified R. microplus serpins were cloned and expressed in Pichia pastoris in order to determine biological targets of these serpins utilising protease inhibition assays. Results: A total of four out of twenty-two serpins identified in our analysis are new R. microplus serpins which were named as RmS-19 to RmS-22. The analyses of DNA and predicted amino acid sequences showed high conservation of the R. microplus serpin sequences. The expression data suggested ubiquitous expression of RmS except for RmS-6 and RmS-14 that were expressed only in nymphs and adult female ovaries, respectively. RmS-19, and -20 were expressed in all tissues samples analysed showing their important role in both parasitic and non-parasitic stages of R. microplus development. RmS-21 was not detected in ovaries and RmS-22 was not identified in ovary and nymph samples but were expressed in the rest of the samples analysed. A total of four expressed recombinant serpins showed protease specific inhibition for Chymotrypsin (RmS-1 and RmS-6), Chymotrypsin / Elastase (RmS-3) and Thrombin (RmS-15). Conclusion: This study constitutes an important contribution and improvement to the knowledge about the physiologic role of R. microplus serpins during the host-tick interaction.

Source and target enzyme signature in serine protease inhibitor active site sequences

Amino acid sequences of proteinaceous proteinase inhibitors have been extensively analysed for deriving information regarding the molecular evolution and functional relationship of these proteins. These sequences have been grouped into several well defined families. It was found that the phylogeny constructed with the sequences corresponding to the exposed loop responsible for inhibition has several branches that resemble those obtained from comparisons using the entire sequence. The major branches of the unrooted tree corresponded to the families to which the inhibitors belonged. Further branching is related to the enzyme specificity of the inhibitor. Examination of the active site loop sequences of trypsin inhibitors revealed that there are strong preferences for specific amino acids at different positions of the loop. These preferences are inhibitor class specific. Inhibitors active against more than one enzyme occur within a class and confirm to class specific sequence in their loops. Hence, only a few positions in the loop seem to determine the specificity. The ability to inhibit the same enzyme by inhibitors that belong to different classes appears to be a result of convergent evolution

A Kunitz-type serine protease inhibitor from Butea monosperma seed and its influence on developmental physiology of Helicoverpa armigera

A protease inhibitor from the seeds of Butea monosperma (BmPI) was purified, characterized and studied for its influence on developmental physiology of Helicover-pa armigera. BmPI on two-dimensional separations indicated the presence of a 14 kDa protein with an isoelectric point in the acidic region (pl 5.6). Multiple Sequence Analysis data suggested that the BmPI contains a sequence motif which is conserved in various trypsin and chymotrypsin inhibitors of Kunitz-type. The inhibitor exhibited trypsin inhibitory activity in a broad range of pH (4-10) and temperature (10-80 degrees C). The enzyme kinetic studies revealed BmPI as a competitive inhibitor with a K-i value of 1.2 x 10(-9) M. In vitro studies with BmPI indicated measurable inhibitory activity on total gut proteolytic enzymes of H. armigera (IC(50)2.0 mu g/ml) and bovine trypsin. BmPI supplemented artificial diet caused dose dependent mortality and reduction in growth and weight. The fertility and fecundity of H. armigera, declined whereas the larval-pupal duration of the insect life cycle extended. These detrimental effects on H. armigera suggest the usefulness of BmPl in insect pest management of food crops. (C) 2014 Elsevier Ltd. All rights reserved.

A Computational Module Assembled from Different Protease Family Motifs Identifies PI PLC from Bacillus cereus as a Putative Prolyl Peptidase with a Serine Protease Scaffold

Proteolytic enzymes have evolved several mechanisms to cleave peptide bonds. These distinct types have been systematically categorized in the MEROPS database. While a BLAST search on these proteases identifies homologous proteins, sequence alignment methods often fail to identify relationships arising from convergent evolution, exon shuffling, and modular reuse of catalytic units. We have previously established a computational method to detect functions in proteins based on the spatial and electrostatic properties of the catalytic residues (CLASP). CLASP identified a promiscuous serine protease scaffold in alkaline phosphatases (AP) and a scaffold recognizing a beta-lactam (imipenem) in a cold-active Vibrio AP. Subsequently, we defined a methodology to quantify promiscuous activities in a wide range of proteins. Here, we assemble a module which encapsulates the multifarious motifs used by protease families listed in the MEROPS database. Since APs and proteases are an integral component of outer membrane vesicles (OMV), we sought to query other OMV proteins, like phospholipase C (PLC), using this search module. Our analysis indicated that phosphoinositide-specific PLC from Bacillus cereus is a serine protease. This was validated by protease assays, mass spectrometry and by inhibition of the native phospholipase activity of PI-PLC by the well-known serine protease inhibitor AEBSF (IC50 = 0.018 mM). Edman degradation analysis linked the specificity of the protease activity to a proline in the amino terminal, suggesting that the PI-PLC is a prolyl peptidase. Thus, we propose a computational method of extending protein families based on the spatial and electrostatic congruence of active site residues.