Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

SPINKs in Tumors: Potential Therapeutic Targets

The serine protease inhibitor Kazal type (SPINK) family includes SPINK1-14 and is the largest branch in the serine protease inhibitor family. SPINKs play an important role in pancreatic physiology and disease, sperm maturation and capacitation, Nager syndrome, inflammation and the skin barrier. Evidence shows that the unregulated expression of SPINK1, 2, 4, 5, 6, 7, and 13 is closely related to human tumors. Different SPINKs exhibit various regulatory modes in different tumors and can be used as tumor prognostic markers. This article reviews the role of SPINK1, 2, 4, 5, 6, 7, and 13 in different human cancer processes and helps to identify new cancer treatment targets.

Climbing Up and Down Binding Landscapes through Deep Mutational Scanning of Three Homologous Protein-Protein Complexes

Protein-protein interactions (PPIs) have evolved to display binding affinities that can support their function. As such, cognate and noncognate PPIs could be highly similar structurally but exhibit huge differences in binding affinities. To understand this phenomenon, we study three homologous protease-inhibitor PPIs that span 9 orders of magnitude in binding affinity. Using state-of-the-art methodology that combines protein randomization, affinity sorting, deep sequencing, and data normalization, we report quantitative binding landscapes consisting of ΔΔGbind values for the three PPIs, gleaned from tens of thousands of single and double mutations. We show that binding landscapes of the three complexes are strikingly different and depend on the PPI evolutionary optimality. We observe different patterns of couplings between mutations for the three PPIs with negative and positive epistasis appearing most frequently at hot-spot and cold-spot positions, respectively. The evolutionary trends observed here are likely to be universal to other biological complexes in the cell.

Arginine substitution by alanine at the P1 position increases the selectivity of CmPI-II, a non-classical Kazal inhibitor

CmPI-II is a Kazal-type tight-binding inhibitor isolated from the Caribbean snail Cenchritis muricatus. This inhibitor has an unusual specificity in the Kazal family, as it can inhibit subtilisin A (SUBTA), elastases and trypsin. An alanine in CmPI-II P1 site could avoid trypsin inhibition while improving/maintaining SUBTA and elastases inhibition. Thus, an alanine mutant of this position (rCmPI-II R12A) was obtained by site-directed mutagenesis. The gene cmpiR12A was expressed in P. pastoris KM71H yeast. The recombinant protein (rCmPI-II R12A) was purified by the combination of two ionic exchange chromatography (1:cationic, 2 anionic) followed by and size exclusion chromatography. The N-terminal sequence obtained as well as the experimental molecular weight allowed verifying the identity of the recombinant protein, while the correct folding was confirmed by CD experiments. rCmPI-II R12A shows a slightly increase in potency against SUBTA and elastases. The alanine substitution at P1 site on CmPI-II abolishes the trypsin inhibition, confirming the relevance of an arginine residue at P1 site in CmPI-II for trypsin inhibition and leading to a molecule with more potentialities in biomedicine.

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

Quantitative mapping of binding specificity landscapes for homologous targets by using a high-throughput method

To facilitate investigations of protein-protein interactions (PPIs), we developed a novel platform for quantitative mapping of protein binding specificity landscapes, which combines the multi-target screening of a mutagenesis library into high- and low-affinity populations with sophisticated next-generation sequencing analysis. Importantly, this method generates accurate models to predict affinity and specificity values for any mutation within a protein complex, and requires only a few experimental binding affinity measurements using purified proteins for calibration. We demonstrated the utility of the approach by mapping quantitative landscapes for interactions between the N-terminal domain of the tissue inhibitor of metalloproteinase 2 (N-TIMP2) and three matrix metalloproteinases (MMPs) having homologous structures but different affinities (MMP-1, MMP-3, and MMP-14). The binding landscapes for N-TIMP2/MMP-1 and N-TIMP2/MMP-3 showed the PPIs to be almost fully optimized, with most single mutations giving a loss of affinity. In contrast, the non-optimized PPI for N-TIMP2/MMP-14 was reflected in a wide range of binding affinities, where single mutations exhibited a far more attenuated effect on the PPI. Our new platform reliably and comprehensively identified not only hot- and cold-spot residues, but also specificity-switch mutations that shape target affinity and specificity. Thus, our approach provides a methodology giving an unprecedentedly rich quantitative analysis of the binding specificity landscape, which will broaden the understanding of the mechanisms and evolutionary origins of specific PPIs and facilitate the rational design of specific inhibitors for structurally similar target proteins.

SKEMPI 2.0: an updated benchmark of changes in protein-protein binding energy, kinetics and thermodynamics upon mutation

Motivation

Understanding the relationship between the sequence, structure, binding energy, binding kinetics and binding thermodynamics of protein–protein interactions is crucial to understanding cellular signaling, the assembly and regulation of molecular complexes, the mechanisms through which mutations lead to disease, and protein engineering.

Results

We present SKEMPI 2.0, a major update to our database of binding free energy changes upon mutation for structurally resolved protein–protein interactions. This version now contains manually curated binding data for 7085 mutations, an increase of 133%, including changes in kinetics for 1844 mutations, enthalpy and entropy changes for 443 mutations, and 440 mutations, which abolish detectable binding.

Availability and implementation

The database is available as supplementary data and at https://life.bsc.es/pid/skempi2/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Isolation and Characterization of Poecistasin, an Anti-Thrombotic Antistasin-Type Serine Protease Inhibitor from Leech Poecilobdella manillensis

Antistasin, first identified as a potent inhibitor of the blood coagulation factor Xa, is a novel family of serine protease inhibitors. In this study, we purified a novel antistasin-type inhibitor from leech Poecilobdella manillensis called poecistasin. Amino acid sequencing of this 48-amino-acid protein revealed that poecistasin was an antistasin-type inhibitor known to consist of only one domain. Poecistasin inhibited factor XIIa, kallikrein, trypsin, and elastase, but had no inhibitory effect on factor Xa and thrombin. Poecistasin showed anticoagulant activities. It prolonged the activated partial thromboplastin time and inhibited FeCl₃-induced carotid artery thrombus formation, implying its potent function in helping Poecilobdella manillensis to take a blood meal from the host by inhibiting coagulation. Poecistasin also suppressed ischemic stroke symptoms in transient middle cerebral artery occlusion mice model. Our results suggest that poecistasin from the leech Poecilobdella manillensis plays a crucial role in blood-sucking and may be an excellent candidate for the development of clinical anti-thrombosis and anti-ischemic stroke medicines.

Thrombospondin-1 protects against pathogen-induced lung injury by limiting extracellular matrix proteolysis

Acute lung injury is characterized by excessive extracellular matrix proteolysis and neutrophilic inflammation. A major risk factor for lung injury is bacterial pneumonia. However, host factors that protect against pathogen-induced and host-sustained proteolytic injury following infection are poorly understood. Pseudomonas aeruginosa (PA) is a major cause of nosocomial pneumonia and secretes proteases to amplify tissue injury. We show that thrombospondin-1 (TSP-1), a matricellular glycoprotein released during inflammation, dose-dependently inhibits PA metalloendoprotease LasB, a virulence factor. TSP-1-deficient (Thbs1-/-) mice show reduced survival, impaired host defense, and increased lung permeability with exaggerated neutrophil activation following acute intrapulmonary PA infection. Administration of TSP-1 from platelets corrects the impaired host defense and aberrant injury in Thbs1-/- mice. Although TSP-1 is cleaved into 2 fragments by PA, TSP-1 substantially inhibits Pseudomonas elastolytic activity. Administration of LasB inhibitor, genetic disabling of the PA type II secretion system, or functional deletion of LasB improves host defense and neutrophilic inflammation in mice. Moreover, TSP-1 provides an additional line of defense by directly subduing host-derived proteolysis, with dose-dependent inhibition of neutrophil elastase from airway neutrophils of mechanically ventilated critically ill patients. Thus, a host matricellular protein provides dual levels of protection against pathogen-initiated and host-sustained proteolytic injury following microbial trigger.

Thrombospondin-1 restrains neutrophil granule serine protease function and regulates the innate immune response during Klebsiella pneumoniae infection

Neutrophil elastase (NE) and cathepsin G (CG) contribute to intracellular microbial killing but, if left unchecked and released extracellularly, promote tissue damage. Conversely, mechanisms that constrain neutrophil serine protease activity protect against tissue damage but may have the untoward effect of disabling the microbial killing arsenal. The host elaborates thrombospondin-1 (TSP-1), a matricellular protein released during inflammation, but its role during neutrophil activation following microbial pathogen challenge remains uncertain. Mice deficient in TSP-1 (thbs1(-/-)) showed enhanced lung bacterial clearance, reduced splenic dissemination, and increased survival compared with wild-type (WT) controls during intrapulmonary Klebsiella pneumoniae infection. More effective pathogen containment was associated with reduced burden of inflammation in thbs1(-/-) mouse lungs compared with WT controls. Lung NE activity was increased in thbs1(-/-) mice following K. pneumoniae challenge, and thbs1(-/-) neutrophils showed enhanced intracellular microbial killing that was abrogated with recombinant TSP-1 administration or WT serum. Thbs1(-/-) neutrophils exhibited enhanced NE and CG enzymatic activity, and a peptide corresponding to amino-acid residues 793-801 within the type-III repeat domain of TSP-1 bridled neutrophil proteolytic function and microbial killing in vitro. Thus, TSP-1 restrains proteolytic action during neutrophilic inflammation elicited by K. pneumoniae, providing a mechanism that may regulate the microbial killing arsenal.

Learning sequence determinants of protein:protein interaction specificity with sparse graphical models

In studying the strength and specificity of interaction between members of two protein families, key questions center on which pairs of possible partners actually interact, how well they interact, and why they interact while others do not. The advent of large-scale experimental studies of interactions between members of a target family and a diverse set of possible interaction partners offers the opportunity to address these questions. We develop here a method, DgSpi (data-driven graphical models of specificity in protein:protein interactions), for learning and using graphical models that explicitly represent the amino acid basis for interaction specificity (why) and extend earlier classification-oriented approaches (which) to predict the ΔG of binding (how well). We demonstrate the effectiveness of our approach in analyzing and predicting interactions between a set of 82 PDZ recognition modules against a panel of 217 possible peptide partners, based on data from MacBeath and colleagues. Our predicted ΔG values are highly predictive of the experimentally measured ones, reaching correlation coefficients of 0.69 in 10-fold cross-validation and 0.63 in leave-one-PDZ-out cross-validation. Furthermore, the model serves as a compact representation of amino acid constraints underlying the interactions, enabling protein-level ΔG predictions to be naturally understood in terms of residue-level constraints. Finally, the model DgSpi readily enables the design of new interacting partners, and we demonstrate that designed ligands are novel and diverse.

Plasmodium falciparum infection induces expression of a mosquito salivary protein (Agaphelin) that targets neutrophil function and inhibits thrombosis without impairing hemostasis

Background

Invasion of mosquito salivary glands (SGs) by Plasmodium falciparum sporozoites is an essential step in the malaria life cycle. How infection modulates gene expression, and affects hematophagy remains unclear.

Principal Findings

Using Affimetrix chip microarray, we found that at least 43 genes are differentially expressed in the glands of Plasmodium falciparum-infected Anopheles gambiae mosquitoes. Among the upregulated genes, one codes for Agaphelin, a 58-amino acid protein containing a single Kazal domain with a Leu in the P1 position. Agaphelin displays high homology to orthologs present in Aedes sp and Culex sp salivary glands, indicating an evolutionarily expanded family. Kinetics and surface plasmon resonance experiments determined that chemically synthesized Agaphelin behaves as a slow and tight inhibitor of neutrophil elastase (K_D∼10 nM), but does not affect other enzymes, nor promotes vasodilation, or exhibit antimicrobial activity. TAXIscan chamber assay revealed that Agaphelin inhibits neutrophil chemotaxis toward fMLP, affecting several parameter associated with cell migration. In addition, Agaphelin reduces paw edema formation and accumulation of tissue myeloperoxidase triggered by injection of carrageenan in mice. Agaphelin also blocks elastase/cathepsin-mediated platelet aggregation, abrogates elastase-mediated cleavage of tissue factor pathway inhibitor, and attenuates neutrophil-induced coagulation. Notably, Agaphelin inhibits neutrophil extracellular traps (NETs) formation and prevents FeCl₃-induced arterial thrombosis, without impairing hemostasis.

Conclusions

Blockade of neutrophil elastase emerges as a novel antihemostatic mechanism in hematophagy; it also supports the notion that neutrophils and the innate immune response are targets for antithrombotic therapy. In addition, Agaphelin is the first antihemostatic whose expression is induced by Plasmodium sp infection. These results suggest that an important interplay takes place in parasite-vector-host interactions.

Malaria is transmitted by Plasmodium falciparum-infected Anopheles gambiae mosquitoes. Salivary gland contributes to the development of the parasite by creating a favorable environment for the infection and facilitating blood feeding and reproduction of the vector. However, the molecular mechanism by which the vector salivary gland modulates parasite/host interactions is not understood. We discovered that infection of the mosquito salivary gland upregulates several genes; among them, one codes for a protease inhibitor named Agaphelin. Notably, Agaphelin was found to exhibit multiple antihemostatic functions by targeting elastase. As a result, it inhibits platelet function which is required for blood to clot, and it prevents cleavage of TFPI, an anticoagulant that has recently been found to play a crucial role in thrombus formation in vivo. Agaphelin also attenuates neutrophils chemotaxis and the release of Neutrophil Extracellular Traps. These results provide evidence that neutrophils serve as a link between coagulation and the innate immune response. Agaphelin also exhibits anti-inflammatory and antithrombotic effects in vivo. Furthermore, Agaphelin did not promote bleeding, suggesting that targeting neutrophil exhibits potential therapeutic value. Altogether, these results highlight that the interplay between parasite, vector and host is a dynamic process that contributes and sustains the interface among Plasmodium, Anopheles and humans.

Learning Sequence Determinants of Protein:protein Interaction Specificity with Sparse Graphical Models

In studying the strength and specificity of interaction between members of two protein families, key questions center on which pairs of possible partners actually interact, how well they interact, and why they interact while others do not. The advent of large-scale experimental studies of interactions between members of a target family and a diverse set of possible interaction partners offers the opportunity to address these questions. We develop here a method, DgSpi (Data-driven Graphical models of Specificity in Protein:protein Interactions), for learning and using graphical models that explicitly represent the amino acid basis for interaction specificity (why) and extend earlier classification-oriented approaches (which) to predict the ΔG of binding (how well). We demonstrate the effectiveness of our approach in analyzing and predicting interactions between a set of 82 PDZ recognition modules, against a panel of 217 possible peptide partners, based on data from MacBeath and colleagues. Our predicted ΔG values are highly predictive of the experimentally measured ones, reaching correlation coefficients of 0.69 in 10-fold cross-validation and 0.63 in leave-one-PDZ-out cross-validation. Furthermore, the model serves as a compact representation of amino acid constraints underlying the interactions, enabling protein-level ΔG predictions to be naturally understood in terms of residue-level constraints. Finally, as a generative model, DgSpi readily enables the design of new interacting partners, and we demonstrate that designed ligands are novel and diverse.

Additivity-based design of the strongest possible turkey ovomucoid third domain inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB)

We describe here successful designs of strong inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB). For each enzyme two inhibitor variants were designed. In one, the reactive site residue (position 18) was retained and the best residues were substituted at contact positions 13, 14, and 15. In the other variant the best residues were substituted at all contact positions except the reactive site where a Gly was substituted. The four designed variants were: for PPE, T(13)E(14)Y(15)-OMTKY3 and T(13)E(14)Y(15)G(18)M(21)P(32)V(36)-OMTKY3, and for SGPB, S(13)D(14)Y(15)-OMTKY3 and S(13)D(14)Y(15)G(18)I(19)K(21)-OMTKY3. The free energies of association (ΔG(0)) of expressed variants have been measured with the proteases for which they were designed as well as with five other serine proteases and the results are discussed.

Spink13, an epididymis-specific gene of the Kazal-type serine protease inhibitor (SPINK) family, is essential for the acrosomal integrity and male fertility

Sperm maturation involves numerous surface modifications by a variety of secreted proteins from epididymal epithelia. The sperm surface architecture depends on correct localization of its components and highlights the importance of the sequence of the proteolytic processing of the sperm surface in the epididymal duct. The presence of several protease inhibitors from different families is consistent with the hypothesis that correctly timed epididymal protein processing is essential for proper sperm maturation. Here we show that the rat (Rattus norvegicus) epididymis-specific gene Spink13, an androgen-responsive serine protease inhibitor, could bind to the sperm acrosome region. Furthermore, knockdown of Spink13 in vivo dramatically enhanced the acrosomal exocytosis during the process of capacitation and thus led to a significant reduction in male fertility, indicating that Spink13 was essential for sperm maturation. We conclude that blockade of SPINK13 may provide a new putative target for post-testicular male contraceptives.

A new type of Kazal proteinase inhibitor related to shrimp Penaeus (Litopenaeus) vannamei immunity

A clone encoding a four-Kazal domain-containing protein was isolated from the hemostats of a Penaeus vannamei cDNA library. The full-length cDNA sequence is 975 bp in length and encodes a 24.4 kDa protein (228 residues). Four Kazal domains, each 43-46 residues in length, were detected in the deduced primary structure. The first, third and fourth domains have the CPLREELPVC, CPAVYDPVC and CPLYVDPVC motifs, respectively, suggesting that they are able to inhibit chymotrypsin and elastase. The mRNA levels of the Kazal protein were modified after the injection of Vibrio alginolyticus, indicating the probable role of this protein in the immune response. All these characteristics are similar to previously reported shrimp Kazal, however, based on both domain architecture and expression profile following Vibrio stimulation, this protein represents a new type of Kazal inhibitor associated with shrimp immunity.

A divide-and-conquer approach to determine the Pareto frontier for optimization of protein engineering experiments

In developing improved protein variants by site-directed mutagenesis or recombination, there are often competing objectives that must be considered in designing an experiment (selecting mutations or breakpoints): stability versus novelty, affinity versus specificity, activity versus immunogenicity, and so forth. Pareto optimal experimental designs make the best trade-offs between competing objectives. Such designs are not "dominated"; that is, no other design is better than a Pareto optimal design for one objective without being worse for another objective. Our goal is to produce all the Pareto optimal designs (the Pareto frontier), to characterize the trade-offs and suggest designs most worth considering, but to avoid explicitly considering the large number of dominated designs. To do so, we develop a divide-and-conquer algorithm, Protein Engineering Pareto FRontier (PEPFR), that hierarchically subdivides the objective space, using appropriate dynamic programming or integer programming methods to optimize designs in different regions. This divide-and-conquer approach is efficient in that the number of divisions (and thus calls to the optimizer) is directly proportional to the number of Pareto optimal designs. We demonstrate PEPFR with three protein engineering case studies: site-directed recombination for stability and diversity via dynamic programming, site-directed mutagenesis of interacting proteins for affinity and specificity via integer programming, and site-directed mutagenesis of a therapeutic protein for activity and immunogenicity via integer programming. We show that PEPFR is able to effectively produce all the Pareto optimal designs, discovering many more designs than previous methods. The characterization of the Pareto frontier provides additional insights into the local stability of design choices as well as global trends leading to trade-offs between competing criteria.

Two Kazal-type protease inhibitors from Macrobrachium nipponense and Eriocheir sinensis: comparative analysis of structure and activities

Kazal-type inhibitors (KPIs) play important roles in many biological and physiological processes, such as blood clotting, the immune response and reproduction. In the present study, two male reproductive tract KPIs, termed Man-KPI and Ers-KPI, were identified in Macrobrachium nipponense and Eriocheir sinensis, respectively. The inhibitory activities of recombinant Man-KPI and Ers-KPI against chymotrypsin, elastase, trypsin and thrombin were determined. The results showed that both of them strongly inhibit chymotrypsin and elastase. Kinetic studies were performed to elucidate their inhibition mechanism. Furthermore, individual domains were also expressed to learn further which domain contributes to the inhibitory activities of intact KPIs. Only Man-KPI_domain3 is active in the inhibition of chymotrypsin and elastase. Meanwhile, Ers-KPI_domain2 and 3 are responsible for inhibition of chymotrypsin, and Ers-KPI_domains2, 3 and 4 are responsible for the inhibition of elastase. Meanwhile, the inhibitory activities of these two KPIs toward Macrobrachium rosenbergii, M. nipponense and E. sinensis sperm were compared with that of the Kazal-type peptidase inhibitor (MRPINK) characterized from the M. rosenbergii reproductive tract in a previous study. The results demonstrated that KPIs can completely inhibit the gelatinolytic activities of sperm proteases from their own species, while different levels of cross-inhibition were observed between KPI and proteases from different species. These results may provide new perspective to further clarify the mechanism of KPI-proteases interaction in the male reproductive system.

Determinants of affinity and proteolytic stability in interactions of Kunitz family protease inhibitors with mesotrypsin

An important functional property of protein protease inhibitors is their stability to proteolysis. Mesotrypsin is a human trypsin that has been implicated in the proteolytic inactivation of several protein protease inhibitors. We have found that bovine pancreatic trypsin inhibitor (BPTI), a Kunitz protease inhibitor, inhibits mesotrypsin very weakly and is slowly proteolyzed, whereas, despite close sequence and structural homology, the Kunitz protease inhibitor domain of the amyloid precursor protein (APPI) binds to mesotrypsin 100 times more tightly and is cleaved 300 times more rapidly. To define features responsible for these differences, we have assessed the binding and cleavage by mesotrypsin of APPI and BPTI reciprocally mutated at two nonidentical residues that make direct contact with the enzyme. We find that Arg at P(1) (versus Lys) favors both tighter binding and more rapid cleavage, whereas Met (versus Arg) at P'(2) favors tighter binding but has minimal effect on cleavage. Surprisingly, we find that the APPI scaffold greatly enhances proteolytic cleavage rates, independently of the binding loop. We draw thermodynamic additivity cycles analyzing the interdependence of P(1) and P'(2) substitutions and scaffold differences, finding multiple instances in which the contributions of these features are nonadditive. We also report the crystal structure of the mesotrypsin·APPI complex, in which we find that the binding loop of APPI displays evidence of increased mobility compared with BPTI. Our data suggest that the enhanced vulnerability of APPI to mesotrypsin cleavage may derive from sequence differences in the scaffold that propagate increased flexibility and mobility to the binding loop.

Trp-26 imparts functional versatility to human alpha-defensin HNP1

We performed a comprehensive alanine scan of human alpha-defensin HNP1 and tested the ability of the resulting analogs to kill Staphylococcus aureus, inhibit anthrax lethal factor, and bind human immunodeficiency virus-1 gp120. By far, the most deleterious mutation for all of these functions was W26A. The activities lost by W26A-HNP1 were restored progressively by replacing W26 with non-coded, straight-chain aliphatic amino acids of increasing chain length. The hydrophobicity of residue 26 also correlated with the ability of the analogs to bind immobilized wild type HNP1 and to undergo further self-association. Thus, the hydrophobicity of residue 26 is not only a key determinant of the direct interactions of HNP1 with target molecules, but it also governs the ability of this peptide to form dimers and more complex quaternary structures at micromolar concentrations. Although all defensin peptides are cationic, their amphipathicity is at least as important as their positive charge in enabling them to participate in innate host defense.

Domain inhibitory and bacteriostatic activities of the five-domain Kazal-type serine proteinase inhibitor from black tiger shrimp Penaeus monodon

Serine proteinase inhibitors (SPIs) in multi-cellular organisms are important modulators of proteinase activities in various biological processes. A five-domain Kazal-type SPI SPIPm2 from the black tiger shrimp Penaeus monodon is presumably involved in innate immune response. The SPIPm2 with the domain P1 residues T, A, E, K and E was isolated from the hemocyte cDNA libraries and found to strongly inhibit subtilisin and elastase, and weakly inhibit trypsin. To unravel further the inhibitory activity of each domain, we subcloned, over-expressed and purified each individual SPI domain. Their inhibitory specificities against trypsin, subtilisin and elastase were determined. Domain 1 was found to be inactive. Domains 2, 3 and 5 inhibited subtilisin. Domain 2 inhibited also elastase. Domain 4 weakly inhibited subtilisin and trypsin. The intact SPIPm2 inhibitor was found to possess bacteriostatic activity against the Bacillus subtilis but not the Bacillus megaterium, Staphylococcus aureus, Vibrio harveyi 639 and Escherichia coli JM109. Domains 2, 4 and 5 contributed to this bacteriostatic activity.

Inhibition mechanism and the effects of structure on activity of male reproduction-related peptidase inhibitor Kazal-type (MRPINK) of Macrobrachium rosenbergii

In our previous reports, a Kazal family serine protease inhibitor, male reproduction-related peptidase inhibitor Kazal-type (MRPINK) has been identified from the prawn, Macrobrachium rosenbergii, and discovered having an inhibitory effect on the sperm gelatinolytic activity. MRPINK was predicated to inhibit chymotrypsin since it contains leucine and proline at P(1) positions of the two domains, respectively. In this report, recombinant MRPINK was as expected found to specifically inhibit chymotrypsin, but no inhibition was detected against trypsin or thrombin. By the analysis of kinetic tests, the inhibition mechanism of MRPINK was determined to be typical competitive model with K (i) of 354 nM. To elucidate the effects of structure on activity of MRPINK, the mutants (domain-1 only, domain-2 only, MRPINK(P88I), MRPINK(L37K), MRPINK(L37A), and MRPINK(L37G)) were prepared and their inhibitory activities assayed. The results showed that domain-2 was the key contributor to the inhibition of chymotrypsin (K (i) of 416 nM) and P(1) Pro was crucial for the activity. Nevertheless, whether the P(1) amino acid residue was Leu, or even if it was replaced by Lys, Ala, or Gly, domain-1 was ineffective to the activity.

Enzyme optimization: moving from blind evolution to statistical exploration of sequence-function space

Directed evolution is a powerful tool for the creation of commercially useful enzymes, particularly those approaches that are based on in vitro recombination methods, such as DNA shuffling. Although these types of search algorithms are extraordinarily efficient compared with purely random methods, they do not explicitly represent or interrogate the genotype-phenotype relationship and are essentially blind in nature. Recently, however, researchers have begun to apply multivariate statistical techniques to model protein sequence-function relationships and guide the evolutionary process by rapidly identifying beneficial diversity for recombination. In conjunction with state-of-the-art library generation methods, the statistical approach to sequence optimization is now being used routinely to create enzymes efficiently for industrial applications.

When the Surface Tells What Lies Beneath: Combinatorial Phage-display Mutagenesis Reveals Complex Networks of Surface–Core Interactions in the Pacifastin Protease Inhibitor Family

Pacifastin protease inhibitors are small cysteine-rich motifs of approximately 35 residues that were discovered in arthropods. The family is divided into two related groups on the basis of the composition of their minimalist inner core. In group I, the core is governed by a Lys10-Trp26 interaction, while in group II it is organized around Phe10. Group I inhibitors exhibit intriguing taxon specificity: potent arthropod-trypsin inhibitors from this group are almost inactive against vertebrate enzymes. The group I member SGPI-1 and the group II member SGPI-2 are extensively studied inhibitors. SGPI-1 is taxon-selective, while SGPI-2 is not. Individual mutations failed to explain the causes underlying this difference. We deciphered this phenomenon using comprehensive combinatorial mutagenesis and phage display. We produced a complete chimeric SGPI-1 / SGPI-2 inhibitor-phage library, in which the two sequences were shuffled at the highest possible resolution of individual residues. The library was selected for binding to bovine trypsin and crayfish trypsin. Sequence analysis of the selectants revealed that taxon specificity is due to an intra-molecular functional coupling between a surface loop and the Lys10-Trp26 core. Five SGPI-2 surface residues transplanted into SGPI-1 resulted in a variant that retained the "taxon-specific" core, but potently inhibited both vertebrate and arthropod enzymes. An additional rational point mutation resulted in a picomolar inhibitor of both trypsins. Our results challenge the generally accepted view that surface residues are the exclusive source of selectivity for canonical inhibitors. Moreover, we provide important insights into general principles underlying the structure-function properties of small disulfide-rich polypeptides, molecules that exist at the borderline between peptides and proteins.

Expression, purification and characterization of a three-domain Kazal-type inhibitor from silkworm pupae (Bombyx mori)

Serine protease inhibitors are essential for host physiological and immunological activities in insects. Analyzing the amino-acid sequence of a cDNA coding for a serine protease inhibitor in Bombyx mori (BmSPI), we found that BmSPI contained three homologous domains with a conserved sequence of C-X(3)-C-X(9)-C-X(6)-Y-X(7)-C-X(3)-C-X(11)-C similar to that of Kazal-type serine protease inhibitors, suggesting BmSPI as a new member of the Kazal-type serine protease inhibitor family. To characterize the three-domain Kazal-type inhibitor from silkworm pupae, the recombinant protein was expressed in Escherichia coli BL21 (DE3) Star. After purification with affinity and reversed-phase chromatographies, the recombinant BmSPI with a molecular mass of 33.642 Da was shown to be a specific subtilisin A inhibitor. Further studies indicated that the K(i) value of the recombinant BmSPI was 3.35 nM and the inhibitor seemed to form a 1:1 complex with subtilisin A. This is a first description of the structure and characterization of Kazal-type inhibitor with three domains cloned from silkworm pupae, B. mori.

Ring-Toss: Capping highly exposed tyrosyl or tryptophyl residues in proteins with beta-cyclodextrin

We have used UV difference spectroscopy and fluorescence spectroscopy to study the perturbation by beta-cyclodextrin of tyrosyl or tryptophyl residues located at each of the 10 variable consensus contact positions in the third domain of turkey ovomucoid. The goal was to monitor the accessibility of the side chain rings of these residues when located at these positions. The results indicated that the tyrosyl or tryptophyl rings are most highly exposed when located in the P1 position followed by the P4 position. It was possible to determine the association constants for beta-cyclodextrin binding at these positions. When located at the P2, P5, P6 and P3' positions, the rings of the tyrosyl or tryptophyl residues were exposed but less so than at the P1 or P4 positions. By contrast, when located at the P1', P2', P14' and P18' positions, the tyrosyl or tryptophyl residues were insufficiently exposed to be perturbed by beta-cyclodextrin, although they reacted positively to dimethyl sulfoxide solvent perturbation. These findings indicate that beta-cyclodextrin perturbation provides a convenient way to detect highly exposed tyrosyls or tryptophyls in proteins. Furthermore, we evaluated the ability of beta-cyclodextrin to inhibit the interaction of turkey ovomucoid third domain variants with different P1 residues. The results showed that the presence of beta-cyclodextrin had little effect on the association constant when the P1 residue was a glycyl residue, but greatly decreased the association constant when the P1 residue was a tyrosyl or tryptophyl residue. Thus, beta-cyclodextrin may be used to selectively modulate the interaction between proteinase inhibitors and their cognate enzymes.

A two disulfide bridge Kazal domain from Phytophthora exhibits stable inhibitory activity against serine proteases of the subtilisin family

BACKGROUND

Kazal-like serine protease inhibitors are defined by a conserved sequence motif. A typical Kazal domain contains six cysteine residues leading to three disulfide bonds with a 1-5/2-4/3-6 pattern. Most Kazal domains described so far belong to this class. However, a novel class of Kazal domains with two disulfide bridges resulting from the absence of the third and sixth cysteines have been found in biologically important molecules, such as human LEKTI, a 15-domain inhibitor associated with the severe congenital disease Netherton syndrome. These domains are referred to as atypical Kazal domains. Previously, EPI1, a Kazal-like protease inhibitor from the oomycete plant pathogen Phytophthora infestans, was shown to be a tight-binding inhibitor of subtilisin A. EPI1 also inhibits and interacts with the pathogenesis-related P69B subtilase of the host plant tomato, suggesting a role in virulence. EPI1 is composed of two Kazal domains, the four-cysteine atypical domain EPI1a and the typical domain EPI1b.

RESULTS

In this study, we predicted the inhibition constants of EPI1a and EPI1b to subtilisin A using the additivity-based sequence to reactivity algorithm (Laskowski algorithm). The atypical domain EPI1a, but not the typical domain EPI1b, was predicted to have strong inhibitory activity against subtilisin A. Inhibition assays and coimmunoprecipitation experiments showed that recombinant domain EPI1a exhibited stable inhibitory activity against subilisin A and was solely responsible for inhibition and interaction with tomato P69B subtilase.

CONCLUSION

The finding that the two disulfide bridge atypical Kazal domain EPI1a is a stable inhibitor indicates that the missing two cysteines and their corresponding disulfide bond are not essential for inhibitor reactivity and stability. This report also suggests that the Laskowski algorithm originally developed and validated with typical Kazal domains might operate accurately for atypical Kazal domains.

Expressed sequence tags from the oomycete fish pathogen Saprolegnia parasitica reveal putative virulence factors

Background

The oomycete Saprolegnia parasitica is one of the most economically important fish pathogens. There is a dramatic recrudescence of Saprolegnia infections in aquaculture since the use of the toxic organic dye malachite green was banned in 2002. Little is known about the molecular mechanisms underlying pathogenicity in S. parasitica and other animal pathogenic oomycetes. In this study we used a genomics approach to gain a first insight into the transcriptome of S. parasitica.

Results

We generated 1510 expressed sequence tags (ESTs) from a mycelial cDNA library of S. parasitica. A total of 1279 consensus sequences corresponding to 525944 base pairs were assembled. About half of the unigenes showed similarities to known protein sequences or motifs. The S. parasitica sequences tended to be relatively divergent from Phytophthora sequences. Based on the sequence alignments of 18 conserved proteins, the average amino acid identity between S. parasitica and three Phytophthora species was 77% compared to 93% within Phytophthora. Several S. parasitica cDNAs, such as those with similarity to fungal type I cellulose binding domain proteins, PAN/Apple module proteins, glycosyl hydrolases, proteases, as well as serine and cysteine protease inhibitors, were predicted to encode secreted proteins that could function in virulence. Some of these cDNAs were more similar to fungal proteins than to other eukaryotic proteins confirming that oomycetes and fungi share some virulence components despite their evolutionary distance

Conclusion

We provide a first glimpse into the gene content of S. parasitica, a reemerging oomycete fish pathogen. These resources will greatly accelerate research on this important pathogen. The data is available online through the Oomycete Genomics Database [1].

A Second Kazal-like protease inhibitor from Phytophthora infestans inhibits and interacts with the apoplastic pathogenesis-related protease P69B of tomato

The plant apoplast forms a protease-rich environment in which proteases are integral components of the plant defense response. Plant pathogenic oomycetes, such as the potato (Solanum tuberosum) and tomato (Lycopersicon esculentum) pathogen Phytophthora infestans, secrete a diverse family of serine protease inhibitors of the Kazal family. Among these, the two-domain EPI1 protein was shown to inhibit and interact with the pathogenesis-related protein P69B subtilase of tomato and was implicated in counter-defense. Here, we describe and functionally characterize a second extracellular protease inhibitor, EPI10, from P. infestans. EPI10 contains three Kazal-like domains, one of which was predicted to be an efficient inhibitor of subtilisin A by an additivity-based sequence to reactivity algorithm (Laskowski algorithm). The epi10 gene was up-regulated during infection of tomato, suggesting a potential role during pathogenesis. Recombinant EPI10 specifically inhibited subtilisin A among the major serine proteases, and inhibited and interacted with P69B subtilase of tomato. The finding that P. infestans evolved two distinct and structurally divergent protease inhibitors to target the same plant protease suggests that inhibition of P69B could be an important infection mechanism for this pathogen.

Energetics and Cooperativity of the Hydrogen Bonding and Anchor Interactions that Bind Peptides to MHC Class II Protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).

Directed molecular evolution by machine learning and the influence of nonlinear interactions

Alternative search strategies for the directed evolution of proteins are presented and compared with each other. In particular, two different machine learning strategies based on partial least-squares regression are developed: the first contains only linear terms that represent a given residue's independent contribution to fitness, the second contains additional nonlinear terms to account for potential epistatic coupling between residues. The nonlinear modeling strategy is further divided into two types, one that contains all possible nonlinear terms and another that makes use of a genetic algorithm to select a subset of important interaction terms. The performance of each modeling type as a function of training set size is analysed. Simulated molecular evolution on a synthetic protein landscape shows the use of machine learning techniques to guide library design can be a powerful addition to library generation methods such as DNA shuffling.

Intramolecular cooperativity in a protein binding site assessed by combinatorial shotgun scanning mutagenesis

Combinatorial shotgun alanine-scanning was used to assess intramolecular cooperativity in the high affinity site (site 1) of human growth hormone (hGH) for binding to its receptor. A total of 19 side-chains were analyzed and statistically significant data were obtained for 145 of the 171 side-chain pairs. The analysis revealed that 90% of the side-chain pairs exhibited no statistically significant pair interactions, and the remaining 10% of side-chain pairs exhibited only small interactions corresponding to cooperative interaction energies with magnitudes less than 0.4 kcal/mol. The statistical predictions were tested by measuring affinities for purified mutant proteins and were found to be accurate for five of six side-chain pairs tested. The results reveal that hGH site 1 behaves in a highly additive manner and suggest that shotgun scanning should be useful for assessing cooperative effects in other protein-protein interactions.

The many faces of protease-protein inhibitor interaction

Proteases and their natural protein inhibitors are among the most intensively studied protein-protein complexes. There are about 30 structurally distinct inhibitor families that are able to block serine, cysteine, metallo- and aspartyl proteases. The mechanisms of inhibition can be related to the catalytic mechanism of protease action or include a mechanism-unrelated steric blockage of the active site or its neighborhood. The structural elements that are responsible for the inhibition most often include the N- or the C-terminus or exposed loop(s) either separately or in combination of several such elements. During complex formation, no major conformational changes are usually observed, but sometimes structural transitions of the inhibitor and enzyme occur. In many cases, convergent evolution, with respect to the inhibitors' parts that are responsible for the inhibition, can be inferred from comparisons of their structures or sequences, strongly suggesting that there are only limited ways to inhibit proteases by proteins.

Effects of single amino acid substitution on the collision-induced dissociation of intact protein ions: Turkey ovomucoid third domain

Expanded understanding of the factors that direct polypeptide ion fragmentation can lead to improved specificity in the use of tandem mass spectrometry for the identification and characterization of proteins. Like the fragmentation of peptide cations, the dissociation of whole protein cations shows several preferred cleavages, the likelihood for which is parent ion charge dependent. While such cleavages are often observed, they are far from universally observed, despite the presence of the residues known to promote them. Furthermore, cleavages at residues not noted to be common in a variety of proteins can be dominant for a particular protein or protein ion charge state. Motivated by the ability to study a small protein, turkey ovomucoid third domain, for which a variety of single amino acid variants are available, the effects of changing the identity of one amino acid in the protein sequence on its dissociation behavior were examined. In particular, changes in amino acids associated with C-terminal aspartic acid cleavage and N-terminal proline cleavage were emphasized. Consistent with previous studies, the product ion spectra were found to be dependent upon the parent ion charge state. Furthermore, the fraction of possible C-terminal aspartic acid cleavages observed to occur for this protein was significantly larger than the fraction of possible N-terminal proline cleavages. In fact, very little N-terminal proline cleavage was noted for the wild-type protein despite the presence of three proline residues in the protein. The addition/removal of proline and aspartic acids was studied along with changes in selected residues adjacent to proline residues. Evidence for inhibition of proline cleavage by the presence of nearby basic residues was noted, particularly if the basic residue was likely to be protonated.

Structure and energetics of protein-protein interactions: the role of conformational heterogeneity in OMTKY3 binding to serine proteases

Proteins with flexible binding surfaces can interact with numerous binding partners. However, this promiscuity is more difficult to understand in "rigid-body" proteins, whose binding results in little, or no, change in the position of backbone atoms. The binding of Kazal inhibitors to serine proteases is considered a classic case of rigid-body binding, although they bind to a wide range of proteases. We have studied the thermodynamics of binding of the Kazal serine protease inhibitor, turkey ovomucoid third domain (OMTKY3), to the serine protease subtilisin Carlsberg using isothermal titration calorimetry and have determined the crystal structure of the complex at very high resolution (1.1A). Comparison of the binding energetics and structure to other OMTKY3 interactions demonstrates that small changes in the position of side-chains can make significant contributions to the binding thermodynamics, including the enthalpy of binding. These effects emphasize that small, "rigid-body" proteins are still dynamic structures, and these dynamics make contributions to both the enthalpy and entropy of binding interactions.

Putting engineering back into protein engineering: bioinformatic approaches to catalyst design

Complex multivariate engineering problems are commonplace and not unique to protein engineering. Mathematical and data-mining tools developed in other fields of engineering have now been applied to analyze sequence-activity relationships of peptides and proteins and to assist in the design of proteins and peptides with specified properties. Decreasing costs of DNA sequencing in conjunction with methods to quickly synthesize statistically representative sets of proteins allow modern heuristic statistics to be applied to protein engineering. This provides an alternative approach to expensive assays or unreliable high-throughput surrogate screens.

Additivity-based prediction of equilibrium constants for some protein-protein associations

For many protein families, such as serine proteinases or serine proteinase inhibitors, the family assignment predicts reactivity only in general terms. Both detailed specificity and quantitative reactivity are lacking. We believe that, for many such protein families, algorithms can be devised by defining the subset of n functionally important sequence positions, making the 19n possible single mutants and measuring their reactivity. Given the assumption that the contributions of the n positions are additive, the reactivities of the 20(n) variants can be predicted. This is illustrated by an almost complete algorithm for the Kazal family of protein inhibitors of serine proteinases.

Expression, refolding, and activation of the catalytic domain of human blood coagulation factor XII

Human blood coagulation factor XII (FXII; 80 kDa) contains a C-terminal serine protease zymogen domain, which becomes activated upon contacting a negative surface. Activated FXII (alphaFXIIa) brings about reciprocal activation of FXII and kallikrein that by further hydrolysis produces the free catalytic domain (betaFXIIa; 28 kDa). Increased levels of alphaFXIIa are associated with coronary heart disease, sepsis, and diabetes. Biophysical investigation of the structural basis of activation, substrate specificity, and regulation of FXII requires an efficient bacterial system for producing the wild-type and mutant recombinant proteins. Here, the cDNA of the zymogen domain of FXII (betaFXII) was cloned into the pET-28a(+) vector and the plasmid was transformed into Escherichia coli strain BL21 (DE3) and overexpressed. The multi-disulfide, recombinant protein, His(6)-betaFXII (rbetaFXII), expressed as an inclusion body, was purified by means of a Ni(2+)-charged resin. The matrix-bound rbetaFXII was subjected to refolding with the glutathione redox system and activated by the in vivo activator, kallikrein. The active form, rbetaFXIIa, obtained in milligram quantities, exhibited similar structural and comparable functional properties relative to human betaFXIIa, as indicated by circular dichroism spectroscopy and kinetics of substrate hydrolysis. Thermodynamics of enzyme:inhibitor complex formation, including the expected 1:1 stoichiometry, was determined for rbetaFXIIa by isothermal calorimetric titration with a specific recombinant protein inhibitor, Cucurbita maxima trypsin inhibitor-V (rCMTI-V; 7kDa).

A theoretical analysis of the thermodynamic contributions for the adsorption of individual protein residues on functionalized surfaces

Although the denaturing of adsorbed proteins on biomaterials surfaces is believed to lead to adverse tissue reactions to implanted materials, very little is currently known of the actual mechanisms involved. These mechanisms must be understood if surfaces are to be proactively designed to control protein adsorption behavior. Concepts widely employed in rational drug design and in protein and RNA folding predictions provide a means to approach this problem. Accordingly, a theoretical analysis has been conducted to estimate the thermodynamic contributions (changes in enthalpy, entropy, and Gibbs free energy) for the adsorption of selected individual mid-chain protein residues to functionalized surfaces. Enthalpic contributions from residue-surface interactions were calculated using semi-empirical quantum mechanical-based computational chemistry methods in a simulated aqueous environment (MOPAC/COSMO), and enthalpic and entropic contributions due to water restructuring effects assumed to occur during adsorption were estimated from experimental data for functional group wetting and calculated changes in solvent accessible surface area as each protein residue approached each surface. When combined with intraprotein residue-residue interactions, the understanding of residue-surface adsorption energy relationships provides a means to begin to predict protein adsorption behavior as a function of biomaterials surface chemistry. It is recognized that several assumptions have been made in this approach that could be challenged, and that truncations necessary due to programming limitations have been applied that may neglect potentially important interactions. Therefore, it must be understood that the modeling predictions may not be directly applicable to biomaterials for surface design under actual physiologic conditions at this stage. However, this attempt at modeling fundamental components of protein adsorption is presented as an initial approach to understanding these complex events.

Protein-protein docking with multiple residue conformations and residue substitutions

The protein docking problem has two major aspects: sampling conformations and orientations, and scoring them for fit. To investigate the extent to which the protein docking problem may be attributed to the sampling of ligand side-chain conformations, multiple conformations of multiple residues were calculated for the uncomplexed (unbound) structures of protein ligands. These ligand conformations were docked into both the complexed (bound) and unbound conformations of the cognate receptors, and their energies were evaluated using an atomistic potential function. The following questions were considered: (1) does the ensemble of precalculated ligand conformations contain a structure similar to the bound form of the ligand? (2) Can the large number of conformations that are calculated be efficiently docked into the receptors? (3) Can near-native complexes be distinguished from non-native complexes? Results from seven test systems suggest that the precalculated ensembles do include side-chain conformations similar to those adopted in the experimental complexes. By assuming additivity among the side chains, the ensemble can be docked in less than 12 h on a desktop computer. These multiconformer dockings produce near-native complexes and also non-native complexes. When docked against the bound conformations of the receptors, the near-native complexes of the unbound ligand were always distinguishable from the non-native complexes. When docked against the unbound conformations of the receptors, the near-native dockings could usually, but not always, be distinguished from the non-native complexes. In every case, docking the unbound ligands with flexible side chains led to better energies and a better distinction between near-native and non-native fits. An extension of this algorithm allowed for docking multiple residue substitutions (mutants) in addition to multiple conformations. The rankings of the docked mutant proteins correlated with experimental binding affinities. These results suggest that sampling multiple residue conformations and residue substitutions of the unbound ligand contributes to, but does not fully provide, a solution to the protein docking problem. Conformational sampling allows a classical atomistic scoring function to be used; such a function may contribute to better selectivity between near-native and non-native complexes. Allowing for receptor flexibility may further extend these results.

Analysis of serine proteinase-inhibitor interaction by alanine shaving

We analyzed the energetic importance of residues surrounding the hot spot (the P(1) position) of bovine pancreatic trypsin inhibitor (BPTI) in interaction with two proteinases, trypsin and chymotrypsin, by a procedure called molecular shaving. One to eight residues of the structural epitope, composed of two extended and exposed loops, were mutated to alanine(s). Although truncation of the side chains of residues surrounding the P(1) position to methyl groups caused a decrease in Delta G(den) values up to 6.4 kcal mole(-1), it did not influence the overall conformation of the inhibitor. We found that the replacement of up to six residues with alanines was fully additive at the level of protein stability. To analyze the influence of the structural epitope on the association energy, we determined association constants for BPTI variants and both enzymes and applied the additivity analysis. Shaving of two binding loops led to a progressive drop in the association energy, more pronounced for trypsin (decrease up to 9.6 kcal mole(-1)) than chymotrypsin (decrease up to 3.5 kcal mole(-1)). In the case of extensively mutated variants interacting with chymotrypsin, the association energies agreed very well with the values calculated from single mutational effects. However, when P(1)-neighboring residues were shaved to alanine(s), their contribution to the association energy was not fully removed because of the presence of methyl groups and main chain-main chain intermolecular hydrogen bonds. Moreover, the hot spot had a different contribution to the complex stability in the fully shaved BPTI variant compared with the wild type, which was caused by perturbations of the P(1)-S(1) electrostatic interaction.

Kazal-type chymotrypsin inhibitor from duck pancreas

A chymotrypsin inhibitor of the Kazal-type has been isolated from duck pancreas, by affinity chromatography on immobilized chymotrypsin, gel filtration on Bio-Gel P-10 and reverse phase (RP)-HPLC. It inhibits bovine chymotrypsin Aalpha with an association constant (K(a)) of 2.06x10(7) M(-1). The complete amino acid sequence was determined after digestion of pyridylethylated inhibitor with Staphylococcus aureus V8 protease and chemical cleavage with CNBr. Duck pancreatic chymotrypsin inhibitor (DPCI) was found to be a single polypeptide chain composed of 65 amino acid residues, corresponding to a molecular mass of 7191 Da.

Selection of potent chymotrypsin and elastase inhibitors from M13 phage library of basic pancreatic trypsin inhibitor (BPTI)

The combinatorial approach offered by phage display has proved to be powerful in obtaining novel variants of canonical inhibitors of serine proteinases that show new binding patterns. We applied this strategy to search for variants of basic pancreatic trypsin inhibitor (BPTI) that would be strong inhibitors of two serine proteinases: bovine alpha-chymotrypsin and porcine pancreatic elastase. BPTI only moderately inhibits the first and does not inhibit the second enzyme. A representative library of 3.2 x 10(4) BPTI variants, randomized at P(1), P(1)', P(2)' and P(3)' positions of the proteinase binding loop, was displayed on the surface of phage M13. After four to five rounds of selection on the target proteinase consensus sequences of the inhibitor binding loop were obtained. In both cases, the variants selected differed from BPTI at two to four positions, with a strong preference for selection of hydrophobic residues. Nevertheless, five of these variants expressed in a free form appeared to be correctly folded, stable proteins, and did not aggregate during thermal denaturation. The midpoints of the thermal unfolding curves of these variants were lowered by 5-20 degrees C as compared to BPTI. The expressed variants proved to be new potent inhibitors of the target enzymes with association constants up to 6.9 x 10(9) M(-1) and 3.7 x 10(10) M(-1) for elastase and chymotrypsin, respectively. Thus, the inhibitory properties of BPTI were improved by as much as 7 x 10(6)-fold towards elastase and 420-fold towards chymotrypsin.