Site-specific conjugation of fibroblast growth factor 2 (FGF2) based on incorporation of alkyne-reactive unnatural amino acid.

Recent advances in site-specific protein modification include the increasingly popular incorporation of unnatural amino acid(s) using amber codon, a method developed by Schultz and coworkers. In this study, we employ this technique to introduce propargyllysine (PrK) in human fibroblast growth factor 2 (FGF2). Owing to an alkyne moiety in its side chain, PrK is compatible with Cu(I)-catalyzed azide-alkyne 1,3-dipolar cycloaddition (CuAAC). We successfully tested CuAAC-mediated conjugation of FGF2 with two compounds - a fluorophore carboxyrhodamine 110 or a cytotoxic drug monomethyl auristatin E (MMAE). In the case of the MMAE conjugate we improved the initial poor conjugation yield to achieve nearly 100% efficiency after extensive optimization. The detergent-based optimization approach may help overcome problems with the CuAAC reaction yield for protein modification with hydrophobic compounds, such as MMAE.

Structural similarities and functional diversity of eukaryotic discoidin-like domains.

The discoidin domain is a approximately 150 amino acid motif common in both eukaryotic and prokaryotic proteins. It is found in a variety of extracellular, intracellular and transmembrane multidomain proteins characterized by a considerable functional diversity, mostly involved in developmental processes. The biological role of the domain depends on its interactions with different molecules, including growth factors, phospholipids and lipids, galactose or its derivatives, and collagen. The conservation of the motif, as well as the serious physiological consequences of discoidin domain disorders underscore the importance of the fold, while the ability to accommodate such an extraordinarily broad range of ligand molecules makes it a fascinating research target. In present review we characterize the distinctive features of discoidin domains and briefly outline the biological role of this module in various eukaryotic proteins.

Structure-function relationships in serine protease-bovine pancreatic trypsin inhibitor interaction.

We report our progress in understanding the structure-function relationships for the interaction between BPTI and serine proteases. We focused on extensive mutagenesis of four crucial positions from the protease binding loop of BPTI. Selected variants were characterized by determination of association constants, stability parameters and structures of protease-inhibitor complexes.

Canonical protein inhibitors of serine proteases.

Serine proteases and their natural protein inhibitors are among the most intensively studied protein complexes. About 20 structurally diverse inhibitor families have been identified, comprising alpha-helical, beta sheet, and alpha/beta proteins, and different folds of small disulfide-rich proteins. Three different types of inhibitors can be distinguished based on their mechanism of action: canonical (standard mechanism) and non-canonical inhibitors, and serpins. The canonical inhibitors bind to the enzyme through an exposed convex binding loop, which is complementary to the active site of the enzyme. The mechanism of inhibition in this group is always very similar and resembles that of an ideal substrate. The non-canonical inhibitors interact through their N-terminal segment. There are also extensive secondary interactions outside the active site, contributing significantly to the strength, speed, and specificity of recognition. Serpins, similarly to the canonical inhibitors, interact with their target proteases in a substrate-like manner; however, cleavage of a single peptide bond in the binding loop leads to dramatic structural changes.

Structure-function relationship of serine protease-protein inhibitor interaction.

We report our progress in understanding the structure-function relationship of the interaction between protein inhibitors and several serine proteases. Recently, we have determined high resolution solution structures of two inhibitors Apis mellifera chymotrypsin inhibitor-1 (AMCI-I) and Linum usitatissimum trypsin inhibitor (LUTI) in the free state and an ultra high resolution X-ray structure of BPTI. All three inhibitors, despite totally different scaffolds, contain a solvent exposed loop of similar conformation which is highly complementary to the enzyme active site. Isothermal calo- rimetry data show that the interaction between wild type BPTI and chymotrypsin is entropy driven and that the enthalpy component opposes complex formation. Our research is focused on extensive mutagenesis of the four positions from the protease binding loop of BPTI: P1, P1', P3, and P4. We mutated these residues to different amino acids and the variants were characterized by determination of the association constants, stability parameters and crystal structures of protease-inhibitor complexes. Accommodation of the P1 residue in the S1 pocket of four proteases: chymotrypsin, trypsin, neutrophil elastase and cathepsin G was probed with 18 P1 variants. High resolution X-ray structures of ten complexes between bovine trypsin and P1 variants of BPTI have been determined and compared with the cognate P1 Lys side chain. Mutations of the wild type Ala16 (P1') to larger side chains always caused a drop of the association constant. According to the crystal structure of the Leu16 BPTI-trypsin complex, introduction of the larger residue at the P1' position leads to steric conflicts in the vicinity of the mutation. Finally, mutations at the P4 site allowed an improvement of the association with several serine proteases involved in blood clotting. Conversely, introduction of Ser, Val, and Phe in place of Gly12 (P4) had invariably a destabilizing effect on the complex with these proteases.

Enhancement of chymotrypsin-inhibitor/substrate interactions by 3 M NaCl.

A series of 16 bovine pancreatic trypsin inhibitor variants mutated at the P(1) position of the binding loop and seven tetrapeptide p-nitroanilide (pNa) substrates of the general formula: suc-Ala-Ala-Pro-Aaa-pNa (where Aaa denotes either: Phe, Arg, Lys, Leu, Met, Nva, Nle) were used to investigate the influence of high salt concentration on the activity of bovine chymotrypsin. The increase of the association constant (K(a)) and the specificity index (k(cat)/K(m)) in the presence of 3 M NaCl highly depends on the chemical nature of the residue at the P(1) position. The highest increase was observed for inhibitors/substrates containing the basic side chains at this site. Surprisingly, for the remaining 13 residues the observed salt effect is not correlated with any side chain properties. In particular, there is a lack of correlation between the accessible non-polar surface area and the magnitude of the salt effect. It suggests that salt-induced increase of the K(a) and k(cat)/K(m) values is not caused by the enhancement of the hydrophobic interactions in chymotrypsin-inhibitor/substrate complex. Moreover, the increase of the K(a) and k(cat)/K(m) values occurs only in the presence of Na(+) ions, while K(+) and Li(+) ions do not change the activity of chymotrypsin. Additionally, the activities of two other proteinases: bovine trypsin and Streptomyces griseus proteinase B were tested in the presence of 3 M NaCl using their specific substrates. The activity of both enzymes was almost not affected by the presence of high NaCl concentration.

Ultrahigh-resolution structure of a BPTI mutant.

The crystal structure of a mutant of bovine pancreatic trypsin inhibitor has been refined to 0.86 A resolution using low-temperature synchrotron data. The variant contains three mutations in the binding loop (Thr11Ala, Pro13Ala, Lys15Arg) and an unrelated Met52Leu substitution. Refinement with anisotropic displacement parameters and with removal of main-chain stereochemical restraints converged with R = 0.1035. The use of full-matrix refinement provided an estimate of the variances in the derived parameters. Some stereochemical parameters, such as the planarity of the peptide group and the value of the N-C(alpha)-C angle, show a wide spread, suggesting that the standard values used as restraints in protein structure refinements may not always be entirely appropriate. Comparison with the recently determined room-temperature structure of the same mutant at 1.42 A resolution confirms the previous observations and provides new details, such as a double conformation of the main chain at Leu29 and at Gly56-Gly57, a high proportion (over 20%) of residues in double conformations, correlation of disorder through lattice contacts and the positions of H atoms, including those in water molecules, and their involvement in C-H...O and N-H...pi hydrogen bonds.

Amino-acid substitutions at the fully exposed P1 site of bovine pancreatic trypsin inhibitor affect its stability.

It is widely accepted that solvent-exposed sites in proteins play only a negligible role in determining protein energetics. In this paper we show that amino acid substitutions at the fully exposed Lys15 in bovine pancreatic trypsin inhibitor (BPTI) influenced the CD- and DSC-monitored stability: The T(den) difference between the least (P1 Trp) and the most stable (P1 His) mutant is 11.2 degrees C at pH 2.0. The DeltaH(den) versus T(den) plot for all the variants at three pH values (2.0, 2.5, 3.0) is linear (DeltaC(p,den) = 0.41 kcal* mole(-1) * K(-1); 1 cal = 4.18 J) leading to a DeltaG(den) difference of 2.1 kcal*mole(-1). Thermal denaturation of the variants monitored by CD signal at pH 2.0 in the presence of 6 M GdmCl again showed differences in their stability, albeit somewhat smaller (DeltaT(den) =7.1 degrees C). Selective reduction of the Cys14-Cys 38 disulfide bond, which is located in the vicinity of the P1 position did not eliminate the stability differences. A correlation analysis of the P1 stability with different properties of amino acids suggests that two mechanisms may be responsible for the observed stability differences: the reverse hydrophobic effect and amino acid propensities to occur in nonoptimal dihedral angles adopted by the P1 position. The former effect operates at the denatured state level and causes a drop in protein stability for hydrophobic side chains, due to their decreased exposure upon denaturation. The latter factor influences the native state energetics and results from intrinsic properties of amino acids in a way similar to those observed for secondary structure propensities. In conclusion, our results suggest that the protein-stability-derived secondary structure propensity scales should be taken with more caution.

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

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Structure of small G proteins and their regulators.

In recent years small G proteins have become an intensively studied group of regulatory GTP hydrolases involved in cell signaling. More than 100 small G proteins have been identified in eucaryotes from protozoan to human. The small G protein superfamily includes Ras, Rho Rab, Rac, Sarl/Arf and Ran homologs, which take part in numerous and diverse cellular processes, such as gene expression, cytoskeleton reorganization, microtubule organization, and vesicular and nuclear transport. These proteins share a common structural core, described as the G domain, and significant sequence similarity. In this paper we review the available data on G domain structure, together with a detailed analysis of the mechanism of action. We also present small G protein regulators: GTPase activating proteins that bind to a catalytic G domain and increase its low intrinsic hydrolase activity, GTPase dissociation inhibitors that stabilize the GDP-bound, inactive state of G proteins, and guanine nucleotide exchange factors that accelerate nucleotide exchange in response to cellular signals. Additionally, in this paper we describe some aspects of small G protein interactions with down-stream effectors.

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.

Amide proton temperature coefficients as hydrogen bond indicators in proteins.

Correlations between amide proton temperature coefficients (deltasigmaHN/deltaT) and hydrogen bonds were investigated for a data set of 793 amides derived from 14 proteins. For amide protons showing temperature gradients more positive than -4.6 ppb/K there is a hydrogen bond predictivity value exceeding 85%. It increases to over 93% for amides within the range between -4 and -1 ppb/K. Detailed analysis shows an inverse proportionality between amide proton temperature coefficients and hydrogen bond lengths. Furthermore, for hydrogen bonds of similar bond lengths, values of temperature gradients in alpha-helices are on average I ppb/K more negative than in beta-sheets. In consequence. a number of amide protons in alpha-helices involved in hydrogen bonds shorter than 2 A show deltasigmaHN/deltaT < -4.6 ppb/K. Due to longer hydrogen bonds, 90% of amides in 3(10) helices and 98% in beta-turns have temperature coefficients more positive than -4.6ppb/K. Ring current effect also significantly influences temperature coefficients of amide protons. In seven out of eight cases non-hydrogen bonded amides strongly deshielded by neighboring aromatic rings show temperature coefficients more positive than -2 ppb/K. In general, amide proton temperature gradients do not change with pH unless they correspond to conformational changes. Three examples of pH dependent equilibrium showing hydrogen bond formation at higher pH were found. In conclusion, amide proton temperature coefficients offer an attractive and simple way to confirm existence of hydrogen bonds in NMR determined structures.

Substitutions at the P(1) position in BPTI strongly affect the association energy with serine proteinases.

The role of the S(1) subsite in trypsin, chymotrypsin and plasmin has been examined by measuring the association with seven different mutants of bovine pancreatic trypsin inhibitor (BPTI); the mutants contain Gly, Ala, Ser, Val, Leu, Arg, and Trp at the P(1) position of the reactive site. The effects of substitutions at the P(1) position on the association constants are very large, comprising seven orders of magnitude for trypsin and plasmin, and over five orders for chymotrypsin. All mutants showed a decrease of the association constant to the three proteinases in the same order: Ala>Gly>Ser>Arg>Val>Leu>Trp. Calorimetric and circular dichroism methods showed that none of the P1 substitutions, except the P1-Val mutant, lead to destabilisation of the binding loop conformation. The X-ray structure of the complex formed between bovine beta-trypsin and P(1)-Leu BPTI showed that the P(1)-Leu sterically conflicts with the side-chain of P(3)-Ile, which thereby is forced to rotate approximately 90 degrees. Ile18 (P(3)) in its new orientation, in turn interacts with the Tyr39 side-chain of trypsin. Introduction of a large side-chain at the P1' position apparently leads to a cascade of small alterations of the trypsin-BPTI interface that seem to destabilise the complex by it adopting a less optimized packing and by tilting the BPTI molecule up to 15 degrees compared to the native trypsin-BPTI complex.

Inhibition of six serine proteinases of the human coagulation system by mutants of bovine pancreatic trypsin inhibitor.

A series of 12 bovine pancreatic trypsin inhibitor variants mutated in the P(4) and P(3) positions of the canonical binding loop containing additional K15R and M52L mutations were used to probe the role of single amino acid substitutions on binding to bovine trypsin and to the following human proteinases involved in blood clotting: plasmin, plasma kallikrein, factors X(a) and XII(a), thrombin, and protein C. The mutants were expressed in Escherichia coli as fusion proteins with the LE1413 hydrophobic polypeptide and purified from inclusion bodies; these steps were followed by CNBr cleavage and oxidative refolding. The mutants inhibited the blood-clotting proteinases with association constants in the range of 10(3)-10(10) m(-)(1). Inhibition of plasma kallikrein, factors X(a) and XII(a), thrombin, and protein C could be improved by up to 2 orders of magnitude by the K15R substitution. The highest increase in the association constant for P(3) mutant was measured for factor XII(a); P13S substitution increased the K(a) value 58-fold. Several other substitutions at P(3) resulted in about 10-fold increase for factor X(a), thrombin, and protein C. The cumulative P(3) and P(1) effects on K(a) values for the strongest mutant compared with the wild type bovine pancreatic trypsin inhibitor were in the range of 2.2- (plasmin) to 4,000-fold (factors XII(a) and X(a)). The substitutions at the P(4) site always caused negative effects (a decrease in the range from over 1,000- to 1.3-fold) on binding to all studied enzymes, including trypsin. Thermal stability studies showed a very large decrease of the denaturation temperature (about 22 degrees C) for all P(4) mutants, suggesting that substitution of the wild type Gly-12 residue leads to a change in the binding loop conformation manifesting itself in non-optimal binding to the proteinase active site.

NMR solution structure of Apis mellifera chymotrypsin/cathepsin G inhibitor-1 (AMCI-1): structural similarity with Ascaris protease inhibitors.

The three-dimensional structure of the 56 residue polypeptide Apis mellifera chymotrypsin/cathepsin G inhibitor 1 (AMCI-1) isolated from honey bee hemolymph was calculated based on 730 experimental NMR restraints. It consists of two approximately perpendicular beta-sheets, several turns, and a long exposed loop that includes the protease binding site. The lack of extensive secondary structure features or hydrophobic core is compensated by the presence of five disulfide bridges that stabilize both the protein scaffold and the binding loop segment. A detailed analysis of the protease binding loop conformation reveals that it is similar to those found in other canonical serine protease inhibitors. The AMCI-1 structure exhibits a common fold with a novel family of inhibitors from the intestinal parasitic worm Ascaris suum. The pH-induced conformational changes in the binding loop region observed in the Ascaris inhibitor ATI are absent in AMCI-1. Similar binding site sequences and structures strongly suggest that the lack of the conformational change can be attributed to a Glu-->Gln substitution at the P1' position in AMCI-1, compared to ATI. Analysis of amide proton temperature coefficients shows very good correlation with the presence of hydrogen bond donors in the calculated AMCI-1 structure.

Inhibition of serine proteinases from human blood clotting system by squash inhibitor mutants.

A series of six CMTI I variants mutated in the P(2)-P(4)' region of the canonical binding loop were used to probe the role of single amino acid substitutions on binding to the following human proteinases involved in blood clotting: plasmin, plasma kallikrein, factors X(a) and XII(a). The mutants were expressed as fusion proteins with the LE1413 hydrophobic polypeptide in Escherichia coli, purified from inclusion bodies, followed by cyanobromide cleavage and refolding. The mutants inhibited the proteinases with the association constants in the range 10(3)-10(9) M(-1). Inhibition of plasma kallikrein and factors X(a) and XII(a) could be improved up to 30-fold by single mutations. In contrast, neither of the introduced mutations increased inhibitory properties of CMTI I against plasmin. Additionally, using two inhibitors of natural origin, CMTI I (P(1) Arg) and CPTI II (P(1) Lys), we determined the effect of Lys-->Arg on binding to four proteinases. With the exception of plasmin (no effect), P(1) Arg resulted in up to 30-fold stronger binding than P(1) Lys.

High-resolution structure of bovine pancreatic trypsin inhibitor with altered binding loop sequence.

A mutant of bovine pancreatic trypsin inhibitor (BPTI) has been constructed and expressed in Escherichia coli in order to probe the kinetic and structural consequences of truncating the binding loop residues to alanine. In addition to two such mutations (Thr11Ala and Pro13Ala), it has a conservative Lys15Arg substitution at position P(1) and an unrelated Met52Leu change. In spite of the binding loop modification, the affinity for trypsin is only 30 times lower than that of the wild-type protein. At pH 7.5 the protein can be crystallized on the time-scale of hours, yielding very stable crystals of a new (tetragonal) form of BPTI. Conventional source X-ray data collected to 1.4 A at room temperature allowed anisotropic structure refinement characterized by R=0.1048. The structure reveals all 58 residues, including the complete C terminus, which is in a salt-bridge contact with the N terminus. The Cys14-Cys38 disulfide bridge is observed in two distinct chiralities. This bridge, together with an internal water molecule, contributes to the stabilization of the binding loop. The Ala mutations have only an insignificant and localized effect on the binding loop, which retains its wild-type conformation (maximum deviation of loop C(alpha) atoms of 0.7 A at Ala13). Four (instead of the typical three) additional water molecules are buried in an internal cleft and connected to the surface via a sulfate anion. Three more SO(4)(2-) anions are seen in the electron density, one of them located on a 2-fold axis. It participates in the formation of a dimeric structure between symmetry-related BPTI molecules, in which electrostatic and hydrogen bonding interactions resulting from the mutated Lys15Arg substitution are of central importance. This dimeric interaction involves direct recognition loop-recognition loop contacts, part of which are hydrophobic interactions of the patches created by the alanine mutations. Another 2-fold symmetric interaction between the BPTI molecules involves the formation of an antiparallel intermolecular beta-sheet that, together with the adjacent intramolecular beta-hairpin loops, creates a four-stranded structure.

Structure-activity relationships for the interaction of bovine pancreatic trypsin inhibitor with an intracellular site on a large conductance Ca(2+)-activated K(+) channel.

Large conductance Ca(2+)-activated K(+) channels (BK(Ca)) contain an intracellular binding site for bovine pancreatic trypsin inhibitor (BPTI), a well-known inhibitor of various serine proteinase (SerP) enzymes. To investigate the structural basis of this interaction, we examined the activity of 11 BPTI mutants using single BK(Ca) channels from rat skeletal muscle incorporated into planar lipid bilayers. All of the mutants induced discrete substate events at the single-channel level. The dwell time of the substate, which is inversely related to the dissociation rate constant of BPTI, exhibited relatively small changes (<9-fold) for the various mutants. However, the apparent association rate constant varied up to 190-fold and exhibited a positive correlation with the net charge of the molecule, suggesting the presence of a negative electrostatic surface potential in the vicinity of the binding site. The substate current level was unaffected by most of the mutations except for substitutions of Lys15. Different residues at this position were found to modulate the apparent conductance of the BPTI-induced substate to 0% (K15G), 10% (K15F), 30% (K15 wild-type), and 55% (K15V) of the open state at +20 mV. Lys15 is located on a loop of BPTI that forms the primary contact region for binding to many SerPs such as trypsin, chymotrypsin, and elastase. The finding that Lys15 is a determinant of the conductance behavior of the BK(Ca) channel when BPTI is bound implies that the same inhibitory loop that contacts SerP's is located close to the protein interface in the BK(Ca) channel complex. This supports the hypothesis that the C-terminal region of the BK(Ca) channel protein contains a domain homologous to SerP's. We propose a domain interaction model for the mechanism of substate production by Kunitz inhibitors based on current ideas for allosteric activation of BK(Ca) channels by voltage and Ca(2+).

Determination of a high precision structure of a novel protein, Linum usitatissimum trypsin inhibitor (LUTI), using computer-aided assignment of NOESY cross-peaks.

The solution structure of a novel 69 residue proteinase inhibitor, Linum usitatissimum trypsin inhibitor (LUTI), was determined using a method based on computer aided assignment of nuclear Overhauser enhancement spectroscopy (NOESY) data. The approach applied uses the program NOAH/DYANA for automatic assignment of NOESY cross-peaks. Calculations were carried out using two unassigned NOESY peak lists and a set of determined dihedral angle restraints. In addition, hydrogen bonds involving amide protons were identified during calculations using geometrical criteria and values of HN temperature coefficients. Stereospecific assignment of beta-methylene protons was carried out using a standard procedure based on nuclear Overhauser enhancement intensities and 3J(alpha)(beta) coupling constants. Further stereospecific assignment of methylene protons and diastereotopic methyl groups were established upon structure-based method available in the program GLOMSA and chemical shift calculations. The applied algorithm allowed us to assign 1968 out of 2164 peaks (91%) derived from NOESY spectra recorded in H2O and 2H2O. The final experimental data input consisted of 1609 interproton distance restraints, 88 restraints for 44 hydrogen bonds, 63 torsion angle restraints and 32 stereospecifically assigned methylene proton pairs and methyl groups. The algorithm allowed the calculation of a high precision protein structure without the laborious manual assignment of NOESY cross-peaks. For the 20 best conformers selected out of 40 refined ones in the program CNS, the calculated average pairwise rmsd values for residues 3 to 69 were 0.38 A (backbone atoms) and 1.02 A (all heavy atoms). The three-dimensional LUTI structure consists of a mixed parallel and antiparallel beta-sheet, a single alpha-helix and shows the fold of the potato 1 family of proteinase inhibitors. Compared to known structures of the family, LUTI contains Arg and Trp residues at positions P6' and P8', respectively, instead of two Arg residues, involved in the proteinase binding loop stabilization. A consequence of the ArgTrp substitution at P8' is a slightly more compact conformation of the loop relative to the protein core.

Interscaffolding additivity: binding of P1 variants of bovine pancreatic trypsin inhibitor to four serine proteases.

Different families of protein inhibitors of serine proteases share similar conformation of the enzyme-binding loop, while their scaffoldings are completely different. In the enzyme-inhibitor complex, the P1position of the loop makes numerous contacts within the S1pocket and significantly influences the energy of the interaction. Here, we determine the association energies (DeltaGavalues) for the interaction of coded P1variants of bovine pancreatic trypsin inhibitor (BPTI) with bovine beta-trypsin (BT), anionic salmon trypsin (AST), bovine alpha-chymotrypsin (BCHYM), and human neutrophil elastase (HNE). The respective DeltaGaranges are 15, 13, 9, and 8 kcal mol-1(1 cal=4.18 J). Next, through interscaffolding additivity cycles, we compare our set of DeltaGavalues determined for BCHYM and HNE with similar data sets available in the literature for three other inhibitor families. The analysis of the cycles shows that 27 to 83 % of cycles fulfil the criteria of additvity. In one particular case (comparison of associations of P1variants of BPTI and OMTKY3 with BCHYM) there is a structural basis for strongly non-additive behaviour. We argue that the interscaffolding additvity depends on sequential and conformational similarities of sites where the mutation(s) are introduced and on the particular substitution. In the second interscaffolding analysis, we compare binding of the same P1mutants to BT and AST. The high correlation coefficient shows that both trypsins recognize with comparable strength the non-cognate side-chains. However, the cognate Arg and Lys side-chains are recognized significantly more strongly by AST.