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

Scaffold stability and P14' residue steric hindrance in the differential inhibition of FXIIa by Aedes aegypti trypsin inhibitor versus Infestin-4

Kazal-type protease inhibitors strictly regulate Factor XIIa (FXIIa), a blood-clotting serine protease. However, when negatively-charged surface of prosthetic device come into contact with FXII, it undergoes conformational change and auto-activation, leading to thrombus formation. Some research suggests that Kazal-type protease inhibitor specificity against FXIIa is governed solely by the reactive-site loop sequence, as this sequence makes most-if not all-of the direct contacts with FXIIa. Here, we sought to compare the inhibitory properties of two Kazal-type inhibitors, Infestin-4 (Inf4), a potent inhibitor of FXIIa, and Aedes aegypti trypsin inhibitor (AaTI), which does not inhibit FXIIa, to better understand Kazal-type protease specificity and determine the structural components responsible for inhibition. There are only 3 residue differences in the reactive-site loop between AaTI and Inf4. Through site-directed mutagenesis, we show that the reactive-site loop is only partially responsible for the inhibitory specificity of these proteases. The protein scaffold of AaTI is unstable due to an elongated C5C6 region. Through chimeric study, we show that swapping the protease-binding loop and the C5C6 region from Inf4 with that of AaTI can partially enhance the inhibitory activity of the AaTI_Inf4 chimera. Furthermore, the additional substitution of Asn at the P14' position of AaTI with Gly (Gly27 in Inf4) absolves the steric clashing between AaTI and the surface 140-loop of FXIIa, and increases the inhibition of the chimeric AaTI to match that of wild-type Inf4. Our findings suggest that ancillary regions in addition to the reactive-site loop sequence are important factors driving Kazal-type inhibitor specificity.

Crystal structure of Aedes aegypti trypsin inhibitor in complex with μ-plasmin reveals role for scaffold stability in Kazal-type serine protease inhibitor

Kazal-type protease inhibitor specificity is believed to be determined by sequence of the reactive-site loop that make most, if not all, contacts with the serine protease. Here, we determined the complex crystal structure of Aedes aegypti trypsin inhibitor (AaTI) with μ-plasmin, and compared its reactivities with other Kazal-type inhibitors, infestin-1 and infestin-4. We show that the shortened 99-loop of plasmin creates an S2 pocket, which is filled by phenylalanine at the P2 position of the reactive-site loop of infestin-4. In contrast, AaTI and infestin-1 retain a proline at P2, rendering the S2 pocket unfilled, which leads to lower plasmin inhibitions. Furthermore, the protein scaffold of AaTI is unstable, due to an elongated Cys-V to Cys-VI region leading to a less compact hydrophobic core. Chimeric study shows that the stability of the scaffold can be modified by swapping of this Cys-V to Cys-VI region between AaTI and infestin-4. The scaffold instability causes steric clashing of the bulky P2 residue, leading to significantly reduced inhibition of plasmin by AaTI or infestin-4 chimera. Our findings suggest that surface loops of protease and scaffold stability of Kazal-type inhibitor are both necessary for specific protease inhibition, in addition to reactive site loop sequence. PDB ID code: 7E50.

Identification and Characterization of ShSPI, a Kazal-Type Elastase Inhibitor from the Venom of Scolopendra Hainanum

Elastase is a globular glycoprotein and belongs to the chymotrypsin family. It is involved in several inflammatory cascades on the basis of cleaving the important connective tissue protein elastin, and is strictly regulated to a balance by several endogenous inhibitors. When elastase and its inhibitors are out of balance, severe diseases will develop, especially those involved in the cardiopulmonary system. Much attention has been attracted in seeking innovative elastase inhibitors and various advancements have been taken on clinical trials of these inhibitors. Natural functional peptides from venomous animals have been shown to have anti-protease properties. Here, we identified a kazal-type serine protease inhibitor named ShSPI from the cDNA library of the venom glands of Scolopendra hainanum. ShSPI showed significant inhibitory effects on porcine pancreatic elastase and human neutrophils elastase with Ki values of 225.83 ± 20 nM and 12.61 ± 2 nM, respectively. Together, our results suggest that ShSPI may be an excellent candidate to develop a drug for cardiopulmonary diseases.

Directed Evolution of Canonical Loops and Their Swapping between Unrelated Serine Proteinase Inhibitors Disprove the Interscaffolding Additivity Model

Reversible serine proteinase inhibitors comprise 18 unrelated families. Each family has a distinct representative structure but contains a surface loop that adopts the same, canonical conformation in the enzyme-inhibitor complex. The Laskowski mechanism universally applies for the action of all canonical inhibitors independent of their scaffold, but it has two nontrivial extrapolations. Intrascaffolding additivity states that all enzyme-contacting loop residues act independently of each other, while interscaffolding additivity claims that these residues act independently of the scaffold. These theories have great importance for engineering proteinase inhibitors but have not been comprehensively challenged. Therefore, we tested the interscaffolding additivity theory by hard-randomizing all enzyme-contacting canonical loop positions of a Kazal- and a Pacifastin-scaffold inhibitor, displaying the variants on M13 phage, and selecting the libraries on trypsin and chymotrypsin. Directed evolution delivered different patterns on both scaffolds against both enzymes, which contradicts interscaffolding additivity. To quantitatively assess the extent of non-additivity, we measured the affinities of the optimal binding loop variants and their binding loop-swapped versions. While optimal variants have picomolar affinities, swapping the evolved loops results in up to 200,000-fold affinity loss. To decipher the underlying causes, we characterized the stability, overall structure and dynamics of the inhibitors with differential scanning calorimetry, circular dichroism and NMR spectroscopy and molecular dynamic simulations. These studies revealed that the foreign loop destabilizes the lower-stability Pacifastin scaffold, while the higher-stability Kazal scaffold distorts the foreign loop. Our findings disprove interscaffolding additivity and show that loop and scaffold form one integrated unit that needs to be coevolved to provide high-affinity inhibition.

Characterization of Spink6 in mouse skin: the conserved inhibitor of kallikrein-related peptidases is reduced by barrier injury

The proteolytic regulation of the desquamation process by kallikrein-related peptidases (KLKs) is crucial for epidermal barrier function, and elevated KLK levels have been reported in atopic dermatitis. KLKs are controlled by specific inhibitors of the serine protease inhibitor of Kazal-type (Spink) family. Recently, SPINK6 was shown to be present in human stratum corneum. In order to investigate its role in epidermal barrier function, we studied mouse Spink6. Sequence alignment revealed that the Kazal domain of Spink6 is highly conserved in animals. Recombinant Spink6 efficiently inhibited mouse Klk5 and human KLK2, KLK4, KLK5, KLK6, KLK7, KLK12, KLK13, and KLK14, whereas human KLK1 and KLK8 were not inhibited. Spink6 was expressed in mouse epidermis mainly in the stratum granulosum, and the inner root sheath of hair follicles. Stimulation with flagellin, EGF, and IL-1β did not alter Spink6 expression, whereas stimulation with tumor necrosis factor-α (TNFα)/IFNγ and all-trans retinoic acid resulted in a significant downregulation of Spink6 expression in cultured primary mouse keratinocytes. Mechanically and metabolically induced skin barrier dysfunction resulted both in a downregulation of Spink6 expression. Our study indicates that Spink6 is a potent inhibitor of KLKs and involved in skin barrier function.

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.

Cleavage of peptide bonds bearing ionizable amino acids at P(1) by serine proteases with hydrophobic S(1) pocket

Enzymatic hydrolysis of the synthetic substrate succinyl-Ala-Ala-Pro-Xxx-pNA (where Xxx=Leu, Asp or Lys) catalyzed by bovine chymotrypsin (CHYM) or Streptomyces griseus protease B (SGPB) has been studied at different pH values in the pH range 3-11. The pH optima for substrates having Leu, Asp, and Lys have been found to be 7.5-8.0, 5.5-6.0, and ∼10, respectively. At the normally reported pH optimum (pH 7-8) of CHYM and SGPB, the substrate with Leu at the reactive site is more than 25,000-fold more reactive than that with Asp. However, when fully protonated, Asp is nearly as good a substrate as Leu. The pK values of the side chains of Asp and Lys in the hydrophobic S(1) pocket of CHYM and SGPB have been calculated from pH-dependent hydrolysis data and have been found to be about 9 for Asp and 7.4 and 9.7 for Lys for CHYM and SGPB, respectively. The results presented in this communication suggest a possible application of CHYM like enzymes in cleaving peptide bonds contributed by acidic amino acids between pH 5 and 6.

The structural and energetic basis for high selectivity in a high-affinity protein-protein interaction

High-affinity, high-selectivity protein-protein interactions that are critical for cell survival present an evolutionary paradox: How does selectivity evolve when acquired mutations risk a lethal loss of high-affinity binding? A detailed understanding of selectivity in such complexes requires structural information on weak, noncognate complexes which can be difficult to obtain due to their transient and dynamic nature. Using NMR-based docking as a guide, we deployed a disulfide-trapping strategy on a noncognate complex between the colicin E9 endonuclease (E9 DNase) and immunity protein 2 (Im2), which is seven orders of magnitude weaker binding than the cognate femtomolar E9 DNase-Im9 interaction. The 1.77 A crystal structure of the E9 DNase-Im2 complex reveals an entirely noncovalent interface where the intersubunit disulfide merely supports the crystal lattice. In combination with computational alanine scanning of interfacial residues, the structure reveals that the driving force for binding is so strong that a severely unfavorable specificity contact is tolerated at the interface and as a result the complex becomes weakened through "frustration." As well as rationalizing past mutational and thermodynamic data, comparing our noncognate structure with previous cognate complexes highlights the importance of loop regions in developing selectivity and accentuates the multiple roles of buried water molecules that stabilize, ameliorate, or aggravate interfacial contacts. The study provides direct support for dual-recognition in colicin DNase-Im protein complexes and shows that weakened noncognate complexes are primed for high-affinity binding, which can be achieved by economical mutation of a limited number of residues at the interface.

Structure and function of invertebrate Kazal-type serine proteinase inhibitors

Proteinases and proteinase inhibitors are involved in several biological and physiological processes in all multicellular organisms. The proteinase inhibitors function as modulators for controlling the extent of deleterious proteinase activity. The Kazal-type proteinase inhibitors (KPIs) in family I1 are among the well-known families of proteinase inhibitors, widely found in mammals, avian and a variety of invertebrates. Like those classical KPIs, the invertebrate KPIs can be single or multiple domain proteins containing one or more Kazal inhibitory domains linked together by peptide spacers of variable length. All invertebrate Kazal domains of about 40-60 amino acids in length share a common structure which is dictated by six conserved cysteine residues forming three intra-domain disulfide cross-links despite the variability of amino acid sequences between the half-cystines. Invertebrate KPIs are strong inhibitors as shown by their extremely high association constant of 10(7)-10(13)M(-1). The inhibitory specificity of a Kazal domain varies widely with a different reactive P(1) amino acid. Different invertebrate KPI domains may arise from gene duplication but several KPI proteins can also be derived from alternative splicing. The invertebrate KPIs function as anticoagulants in blood-sucking animals such as leech, mosquitoes and ticks. Several KPIs are likely involved in protecting host from microbial proteinases while some from the parasitic protozoa help protecting the parasites from the host digestive proteinase enzymes. Silk moths produce KPIs to protect their cocoon from predators and microbial destruction.

Systematic mutational analysis of peptide inhibition of the p53-MDM2/MDMX interactions

Inhibition of the interaction between the tumor suppressor protein p53 and its negative regulators MDM2 and MDMX is of great interest in cancer biology and drug design. We previously reported a potent duodecimal peptide inhibitor, termed PMI (TSFAEYWNLLSP), of the p53-MDM2 and -MDMX interactions. PMI competes with p53 for MDM2 and MDMX binding at an affinity roughly 2 orders of magnitude higher than that of (17-28)p53 (ETFSDLWKLLPE) of the same length; both peptides adopt nearly identical alpha-helical conformations in the complexes, where the three highlighted hydrophobic residues Phe, Trp, and Leu dominate PMI or (17-28)p53 binding to MDM2 and MDMX. To elucidate the molecular determinants for PMI activity and specificity, we performed a systematic Ala scanning mutational analysis of PMI and (17-28)p53. The binding affinities for MDM2 and MDMX of a total of 35 peptides including 10 truncation analogs were quantified, affording a complete dissection of energetic contributions of individual residues of PMI and (17-28)p53 to MDM2 and MDMX association. Importantly, the N8A mutation turned PMI into the most potent dual-specific antagonist of MDM2 and MDMX reported to date, registering respective K(d) values of 490 pM and 2.4 nM. The co-crystal structure of N8A-PMI-(25-109)MDM2 was determined at 1.95 A, affirming that high-affinity peptide binding to MDM2/MDMX necessitates, in addition to optimized intermolecular interactions, enhanced helix stability or propensity contributed by non-contact residues. The powerful empirical binding data and crystal structures present a unique opportunity for computational studies of peptide inhibition of the p53-MDM2/MDMX interactions.

Experimental and computational analyses of the energetic basis for dual recognition of immunity proteins by colicin endonucleases

Colicin endonucleases (DNases) are bound and inactivated by immunity (Im) proteins. Im proteins are broadly cross-reactive yet specific inhibitors binding cognate and non-cognate DNases with K(d) values that vary between 10(-4) and 10(-14) M, characteristics that are explained by a 'dual-recognition' mechanism. In this work, we addressed for the first time the energetics of Im protein recognition by colicin DNases through a combination of E9 DNase alanine scanning and double-mutant cycles (DMCs) coupled with kinetic and calorimetric analyses of cognate Im9 and non-cognate Im2 binding, as well as computational analysis of alanine scanning and DMC data. We show that differential DeltaDeltaGs observed for four E9 DNase residues cumulatively distinguish cognate Im9 association from non-cognate Im2 association. E9 DNase Phe86 is the primary specificity hotspot residue in the centre of the interface, which is coordinated by conserved and variable hotspot residues of the cognate Im protein. Experimental DMC analysis reveals that only modest coupling energies to Im9 residues are observed, in agreement with calculated DMCs using the program ROSETTA and consistent with the largely hydrophobic nature of E9 DNase-Im9 specificity contacts. Computed values for the 12 E9 DNase alanine mutants showed reasonable agreement with experimental DeltaDeltaG data, particularly for interactions not mediated by interfacial water molecules. DeltaDeltaG predictions for residues that contact buried water molecules calculated using solvated rotamer models met with mixed success; however, we were able to predict with a high degree of accuracy the location and energetic contribution of one such contact. Our study highlights how colicin DNases are able to utilise both conserved and variable amino acids to distinguish cognate from non-cognate Im proteins, with the energetic contributions of the conserved residues modulated by neighbouring specificity sites.

1.2A-resolution crystal structures reveal the second tetrahedral intermediates of streptogrisin B (SGPB)

Streptogrisin B (SGPB) has served as one of the models for studying the catalytic activities of serine peptidases. Here we report its native crystal structures at pH 4.2 at a resolution of 1.18A, and at pH 7.3 at a resolution of 1.23A. Unexpectedly, outstanding electron density peaks occurred in the active site and the substrate-binding region of SGPB in the computed maps at both pHs. The densities at pH 4.2 were assigned as a tetrapeptide, Asp-Ala-Ile-Tyr, whereas those at pH 7.3 were assigned as a tyrosine molecule and a leucine molecule existing at equal occupancies in both of the SGPB molecules in the asymmetric unit. Refinement with relaxed geometric restraints resulted in molecular structures representing mixtures of the second tetrahedral intermediates and the enzyme-product complexes of SGPB existing in a pH-dependent equilibrium. Structural comparisons with the complexes of SGPB with turkey ovomucoid third domain (OMTKY3) and its variants have shown that, upon the formation of the tetrahedral intermediate, residues Glu192A to Gly193 of SGPB move towards the alpha-carboxylate O of residue P1 of the bound species, and adjustments in the side-chain conformational angles of His57 and Ser195 of SGPB favor the progression of the catalytic mechanism of SGPB.

Structural Insights into the Non-additivity Effects in the Sequence-to-Reactivity Algorithm for Serine Peptidases and their Inhibitors

Sequence-to-reactivity algorithms (SRAs) for proteins have the potential of being broadly applied in molecular design. Recently, Laskowski et al. have reported an additivity-based SRA that accurately predicts most of the standard free energy changes of association for variants of turkey ovomucoid third domain (OMTKY3) with six serine peptidases, one of which is streptogrisin B (commonly known as Streptomyces griseus peptidase B, SGPB). Non-additivity effects for residues 18I and 32I, and for residues 20I and 32I of OMTKY3 occurred when the associations with SGPB were predicted using the SRA. To elucidate precisely the mechanics of these non-additivity effects in structural terms, we have determined the crystal structures of the unbound OMTKY3 (with Gly32I as in the wild-type amino acid sequence) at a resolution of 1.16 A, the unbound Ala32I variant of OMTKY3 at a resolution of 1.23 A, and the SGPB:OMTKY3-Ala32I complex (equilibrium association constant K(a)=7.1x10(9) M(-1) at 21(+/-2) C degrees, pH 8.3) at a resolution of 1.70 A. Extensive comparisons with the crystal structure of the unbound OMTKY3 confirm our understanding of some previously addressed non-additivity effects. Unexpectedly, SGPB and OMTKY3-Ala32I form a 1:2 complex in the crystal. Comparison with the SGPB:OMTKY3 complex shows a conformational change in the SGPB:OMTKY3-Ala32I complex, resulting from a hinged rigid-body rotation of the inhibitor caused by the steric hindrance between the methyl group of Ala32IA of the inhibitor and Pro192BE of the peptidase. This perturbs the interactions among residues 18I, 20I, 32I and 36I of the inhibitor, probably resulting in the above non-additivity effects. This conformational change also introduces residue 10I as an additional hyper-variable contact residue to the SRA.

Despite having a common P1 Leu, eglin C inhibits alpha-lytic proteinase a million-fold more strongly than does turkey ovomucoid third domain

Results of the inhibition of alpha-lytic proteinase by two standard mechanism serine proteinase inhibitors, turkey ovomucoid third domain (OMTKY3) and eglin C, and many of their variants are presented. Despite similarities, including an identical P1 residue (Leu) in their primary contact regions, OMTKY3 and eglin C have vastly different association equilibrium constants toward alpha-lytic proteinase, with Ka values of 1.8 x 10(3) and 1.2 x 10(9) M(-1), respectively. Although 12 of the 13 serine proteinases tested in our laboratory for inhibition by OMTKY3 and eglin C are more strongly inhibited by the latter, the million-fold difference observed here with alpha-lytic proteinase is the largest we have seen. The million-fold stronger inhibition by eglin C is retained when the Ka values of the P1 Gly, Ala, Ser, and Ile variants of OMTKY3 and eglin C are compared. Despite the small size of the S1 pocket in alpha-lytic proteinase, interscaffolding additivity for OMTKY3 and eglin C holds well for the four P1 residues tested here. To better understand this difference, we measured Ka values for other OMTKY3 variants, including some that had residues elsewhere in their contact region that corresponded to those of eglin C. Assuming intrascaffolding additivity and using the Ka values obtained for OMTKY3 variants, we designed an OMTKY3-based inhibitor of alpha-lytic proteinase that was predicted to inhibit 10,000-fold more strongly than wild-type OMTKY3. This variant (K13A/P14E/L18A/R21T/N36D OMTKY3) was prepared, and its Ka value was measured against alpha-lytic proteinase. The measured Ka value was in excellent agreement with the predicted one (1.1 x 10(7) and 2.0 x 10(7) M(-1), respectively). Computational protein docking results are consistent with the view that the backbone conformation of eglin C is not significantly altered in the complex with alpha-lytic proteinase. They also show that the strong binding for eglin C correlates well with more favorable atomic contact energy and desolvation energy contributions as compared to OMTKY3.

Analysis of sequence-reactivity space for protein-protein interactions

Sequence-reactivity space is defined by the relationships between amino acid type choices at some residue positions in a protein and the reactivities of the resulting variants. We are studying Kazal superfamily serine proteinase inhibitors, under substitution of any combination of residue types at 10 binding-region positions. Reactivities are defined by the standard free energy of association for an inhibitor against an enzyme, and we are interested in both the strength (the free energy value) and specificity (relative free energy values for one inhibitor against different enzymes). Characterizing the structure of such a space poses several interesting questions: (1) How many sequences achieve particular strength and specificity characteristics? (2) What is the best such sequence? (3) What are some nearly-as-good alternatives? (4) What are their common residue type characteristics (e.g., conservation and correlation)? Although these problems are all highly combinatorial in nature, this article develops an efficient, integrated mechanism to address them under a data-driven model that predicts reactivity for given sequences. We employ sampling and a novel deterministic distribution propagation algorithm, in order to determine both the reactivity distribution and sequence composition statistics; integer programming and a novel branch-and-bound search algorithm, in order to optimize sequences and enumerate near-optimal sequences; and correlation-based sequence decomposition, in order to identify sequence motifs. We demonstrate the value of our mechanism in analyzing the Kazal superfamily sequence-reactivity space, providing insights into the underlying biochemistry and suggesting hypotheses for further experimental consideration. In general, our mechanism offers a valuable tool for investigating the available degrees of freedom in protein design within a combined computational-experimental framework.

Functional evolution within a protein superfamily

The ability to predict and characterize distributions of reactivities over families and even superfamilies of proteins opens the door to an array of analyses regarding functional evolution. In this article, insights into functional evolution in the Kazal inhibitor superfamily are gained by analyzing and comparing predicted association free energy distributions against six serine proteinases, over a number of groups of inhibitors: all possible Kazal inhibitors, natural avian ovomucoid first and third domains, and sets of Kazal inhibitors with statistically weighted combinations of residues. The results indicate that, despite the great hypervariability of residues in the 10 proteinase-binding positions, avian ovomucoid third domains evolved to inhibit enzymes similar to the six enzymes selected, whereas the orthologous first domains are not inhibitors of these enzymes on purpose. Hypervariability arises because of similarity in energetic contribution from multiple residue types; conservation is in terms of functionality, with "good" residues, which make positive or less deleterious contributions to the binding, selected more frequently, and yielding overall the same distributional characteristics. Further analysis of the distributions indicates that while nature did optimize inhibitor strength, the objective may not have been the strongest possible inhibitor against one enzyme but rather an inhibitor that is relatively strong against a number of 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.

The conserved P1' Ser of Bowman-Birk-type proteinase inhibitors is not essential for the integrity of the reactive site loop

The isolated reactive site beta-hairpin loop of Bowman-Birk-type proteinase inhibitors has become a widely studied proteinomimetic because it retains the three-dimensional structure and much of the inhibitory potency of the corresponding region of the complete protein. Here we analyse the role of the P1' Ser residue which is highly conserved and intramolecularly hydrogen bonded in the complete proteins. A combined kinetic and structural analysis of variant proteinomimetic peptides demonstrates that the hydrogen-bond potential of the side-chain oxygen atom of the P1' Ser is not essential for the integrity of the reactive site loop and that it provides only a small contribution to the trypsin affinity and no apparent contribution to the stability against tryptic turnover. We conclude that the potential of the P1' side chain to engineer improved inhibition and selectivity for serine proteinases is best explored further in concert with the side chains of the P2 and P5' residues which may interact or compete for the same space.

Testing of the additivity-based protein sequence to reactivity algorithm

The standard free energies of association (or equilibrium constants) are predicted for 11 multiple variants of the turkey ovomucoid third domain, a member of the Kazal family of protein inhibitors, each interacting with six selected enzymes. The equilibrium constants for 38 of 66 possible interactions are strong enough to measure, and for these, the predicted and measured free energies are compared, thus providing an additional test of the additivity-based sequence to reactivity algorithm. The test appears to be unbiased as the 11 variants were designed a decade ago to study furin inhibition and the specificity of furin differs greatly from the specificities of our six target enzymes. As the contact regions of these inhibitors are highly positive, nonadditivity was expected. Of the 11 variants, one does not satisfy the restriction that either P(2) Thr or P(1)' Glu should be present and all three measurable results on it, as expected, are nonadditive. For the remaining 35 measurements, 22 are additive, 12 are partially additive, and only one is (slightly) nonadditive. These results are comparable to those obtained for a set of 398 equilibrium constants for natural variants of ovomucoid third domains. The expectation that clustering of charges would be nonadditive is modified to the expectation that major nonadditivity will be observed only if the combining sites in both associating proteins involve large charge clusters of the opposite sign. It is also shown here that an analysis of a small variant set can be accomplished with a smaller subset, in this case 13 variants, rather than the whole set of 191 members used for the complete algorithm.