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.