The crystal structure of anionic salmon trypsin in complex with bovine pancreatic trypsin inhibitor

Fish trypsins: potential applications in biomedicine and prospects for production

In fishes, trypsins are adapted to different environmental conditions, and the biochemical and kinetic properties of a broad variety of native isoforms have been studied. Proteolytic enzymes remain in high demand in the detergent, food, and feed industries; however, our analysis of the literature showed that, in the last decade, some fish trypsins have been studied for the synthesis of industrial peptides and for specific biomedical uses as antipathogenic agents against viruses and bacteria, which have been recently patented. In addition, innovative strategies of trypsin administration have been studied to ensure that trypsins retain their properties until they exert their action. Biomedical uses require the production of high-quality enzymes. In this context, the production of recombinant trypsins is an alternative. For this purpose, E. coli-based systems have been tested for the production of fish trypsins; however, P. pastoris-based systems also seem to show great potential in the production of fish trypsins with higher production quality. On the other hand, there is a lack of information regarding the specific structures, biochemical and kinetic properties, and characteristics of trypsins produced using heterologous systems. This review describes the potential uses of fish trypsins in biomedicine and the enzymatic and structural properties of native and recombinant fish trypsins obtained to date, outlining some prospects for their study.

Enzyme surface rigidity tunes the temperature dependence of catalytic rates

The structural origin of enzyme adaptation to low temperature, allowing efficient catalysis of chemical reactions even near the freezing point of water, remains a fundamental puzzle in biocatalysis. A remarkable universal fingerprint shared by all cold-active enzymes is a reduction of the activation enthalpy accompanied by a more negative entropy, which alleviates the exponential decrease in chemical reaction rates caused by lowering of the temperature. Herein, we explore the role of protein surface mobility in determining this enthalpy-entropy balance. The effects of modifying surface rigidity in cold- and warm-active trypsins are demonstrated here by calculation of high-precision Arrhenius plots and thermodynamic activation parameters for the peptide hydrolysis reaction, using extensive computer simulations. The protein surface flexibility is systematically varied by applying positional restraints, causing the remarkable effect of turning the cold-active trypsin into a variant with mesophilic characteristics without changing the amino acid sequence. Furthermore, we show that just restraining a key surface loop causes the same effect as a point mutation in that loop between the cold- and warm-active trypsin. Importantly, changes in the activation enthalpy-entropy balance of up to 10 kcal/mol are almost perfectly balanced at room temperature, whereas they yield significantly higher rates at low temperatures for the cold-adapted enzyme.

Qgui: A high-throughput interface for automated setup and analysis of free energy calculations and empirical valence bond simulations in biological systems

Structural information and activity data has increased rapidly for many protein targets during the last decades. In this paper, we present a high-throughput interface (Qgui) for automated free energy and empirical valence bond (EVB) calculations that use molecular dynamics (MD) simulations for conformational sampling. Applications to ligand binding using both the linear interaction energy (LIE) method and the free energy perturbation (FEP) technique are given using the estrogen receptor (ERα) as a model system. Examples of free energy profiles obtained using the EVB method for the rate-limiting step of the enzymatic reaction catalyzed by trypsin are also shown. In addition, we present calculation of high-precision Arrhenius plots to obtain the thermodynamic activation enthalpy and entropy with Qgui from running a large number of EVB simulations.

Three-dimensional Structure of a Kunitz-type Inhibitor in Complex with an Elastase-like Enzyme

Elastase-like enzymes are involved in important diseases such as acute pancreatitis, chronic inflammatory lung diseases, and cancer. Structural insights into their interaction with specific inhibitors will contribute to the development of novel anti-elastase compounds that resist rapid oxidation and proteolysis. Proteinaceous Kunitz-type inhibitors homologous to the bovine pancreatic trypsin inhibitor (BPTI) provide a suitable scaffold, but the structural aspects of their interaction with elastase-like enzymes have not been elucidated. Here, we increased the selectivity of ShPI-1, a versatile serine protease inhibitor from the sea anemone Stichodactyla helianthus with high biomedical and biotechnological potential, toward elastase-like enzymes by substitution of the P1 residue (Lys(13)) with leucine. The variant (rShPI-1/K13L) exhibits a novel anti-porcine pancreatic elastase (PPE) activity together with a significantly improved inhibition of human neuthrophil elastase and chymotrypsin. The crystal structure of the PPE·rShPI-1/K13L complex determined at 2.0 Å resolution provided the first details of the canonical interaction between a BPTI-Kunitz-type domain and elastase-like enzymes. In addition to the essential impact of the variant P1 residue for complex stability, the interface is improved by increased contributions of the primary and secondary binding loop as compared with similar trypsin and chymotrypsin complexes. A comparison of the interaction network with elastase complexes of canonical inhibitors from the chelonian in family supports a key role of the P3 site in ShPI-1 in directing its selectivity against pancreatic and neutrophil elastases. Our results provide the structural basis for site-specific mutagenesis to further improve the binding affinity and/or direct the selectivity of BPTI-Kunitz-type inhibitors toward elastase-like enzymes.

Protein surface softness is the origin of enzyme cold-adaptation of trypsin

Life has effectively colonized most of our planet and extremophilic organisms require specialized enzymes to survive under harsh conditions. Cold-loving organisms (psychrophiles) express heat-labile enzymes that possess a high specific activity and catalytic efficiency at low temperatures. A remarkable universal characteristic of cold-active enzymes is that they show a reduction both in activation enthalpy and entropy, compared to mesophilic orthologs, which makes their reaction rates less sensitive to falling temperature. Despite significant efforts since the early 1970s, the important question of the origin of this effect still largely remains unanswered. Here we use cold- and warm-active trypsins as model systems to investigate the temperature dependence of the reaction rates with extensive molecular dynamics free energy simulations. The calculations quantitatively reproduce the catalytic rates of the two enzymes and further yield high-precision Arrhenius plots, which show the characteristic trends in activation enthalpy and entropy. Detailed structural analysis indicates that the relationship between these parameters and the 3D structure is reflected by significantly different internal protein energy changes during the reaction. The origin of this effect is not localized to the active site, but is found in the outer regions of the protein, where the cold-active enzyme has a higher degree of softness. Several structural mechanisms for softening the protein surface are identified, together with key mutations responsible for this effect. Our simulations further show that single point-mutations can significantly affect the thermodynamic activation parameters, indicating how these can be optimized by evolution.

Elapid snake venom analyses show the specificity of the peptide composition at the level of genera Naja and Notechis

Elapid snake venom is a highly valuable, but till now mainly unexplored, source of pharmacologically important peptides. We analyzed the peptide fractions with molecular masses up to 10 kDa of two elapid snake venoms—that of the African cobra, N. m. mossambica (genus Naja), and the Peninsula tiger snake, N. scutatus, from Kangaroo Island (genus Notechis). A combination of chromatographic methods was used to isolate the peptides, which were characterized by combining complimentary mass spectrometric techniques. Comparative analysis of the peptide compositions of two venoms showed specificity at the genus level. Three-finger (3-F) cytotoxins, bradykinin-potentiating peptides (BPPs) and a bradykinin inhibitor were isolated from the Naja venom. 3-F neurotoxins, Kunitz/basic pancreatic trypsin inhibitor (BPTI)-type inhibitors and a natriuretic peptide were identified in the N. venom. The inhibiting activity of the peptides was confirmed in vitro with a selected array of proteases. Cytotoxin 1 (P01467) from the Naja venom might be involved in the disturbance of cellular processes by inhibiting the cell 20S-proteasome. A high degree of similarity between BPPs from elapid and viperid snake venoms was observed, suggesting that these molecules play a key role in snake venoms and also indicating that these peptides were recruited into the snake venom prior to the evolutionary divergence of the snakes.

Psychrophily and catalysis

Polar and other low temperature environments are characterized by a low content in energy and this factor has a strong incidence on living organisms which populate these rather common habitats. Indeed, low temperatures have a negative effect on ectothermic populations since they can affect their growth, reaction rates of biochemical reactions, membrane permeability, diffusion rates, action potentials, protein folding, nucleic acids dynamics and other temperature-dependent biochemical processes. Since the discovery that these ecosystems, contrary to what was initially expected, sustain a rather high density and broad diversity of living organisms, increasing efforts have been dedicated to the understanding of the molecular mechanisms involved in their successful adaptation to apparently unfavorable physical conditions. The first question that comes to mind is: How do these organisms compensate for the exponential decrease of reaction rate when temperature is lowered? As most of the chemical reactions that occur in living organisms are catalyzed by enzymes, the kinetic and thermodynamic properties of cold-adapted enzymes have been investigated. Presently, many crystallographic structures of these enzymes have been elucidated and allowed for a rather clear view of their adaptation to cold. They are characterized by a high specific activity at low and moderate temperatures and a rather low thermal stability, which induces a high flexibility that prevents the freezing effect of low temperatures on structure dynamics. These enzymes also display a low activation enthalpy that renders them less dependent on temperature fluctuations. This is accompanied by a larger negative value of the activation entropy, thus giving evidence of a more disordered ground state. Appropriate folding kinetics is apparently secured through a large expression of trigger factors and peptidyl–prolyl cis/trans-isomerases.

Structural insights into serine protease inhibition by a marine invertebrate BPTI Kunitz-type inhibitor

Proteins isolated from marine invertebrates are frequently characterized by exceptional structural and functional properties. ShPI-1, a BPTI Kunitz-type inhibitor from the Caribbean Sea anemone Stichodactyla helianthus, displays activity not only against serine-, but also against cysteine-, and aspartate proteases. As an initial step to evaluate the molecular basis of its activities, we describe the crystallographic structure of ShPI-1 in complex with the serine protease bovine pancreatic trypsin at 1.7Å resolution. The overall structure and the important enzyme-inhibitor interactions of this first invertebrate BPTI-like Kunitz-type inhibitor:trypsin complex remained largely conserved compared to mammalian BPTI-Kunitz inhibitor complexes. However, a prominent stabilizing role within the interface was attributed to arginine at position P3. Binding free-energy calculations indicated a 10-fold decrease for the inhibitor affinity against trypsin, if the P3 residue of ShPI-1 is mutated to alanine. Together with the increased role of Arg(11) at P3 position, slightly reduced interactions at the prime side (Pn') of the primary binding loop and at the secondary binding loop of ShPI-1 were detected. In addition, the structure provides important information for site directed mutagenesis to further optimize the activity of rShPI-1A for biotechnological applications.

Isolation and characterization of novel variants of BBI coding genes from the legume Lathyrus sativus

A pool of twelve cDNA sequences coding for Bowman-Birk inhibitors (BBIs) was identified in the legume grass pea (Lathyrus sativus L.). The corresponding amino acid sequences showed a canonical first anti-trypsin domain, predicted according to the identity of the determinant residue P(1). A more variable second binding loop was observed allowing to identify three groups based on the identity of residue P(1): two groups (Ls_BBI_1 and Ls_BBI_2) carried a second reactive site specific for chymotrypsin, while a third group (Ls_BBI_3) was predicted to inhibit elastase. A fourth variant carrying an Asp in the P(1) position of the second reactive site was identified only from genomic DNA. A phylogenetic tree constructed using grass pea BBIs with their homologs from other legume species revealed grouping based on taxonomy and on specificity of the reactive sites. Five BBI sequences, representing five different second reactive sites, were heterologously expressed in the yeast Pichia pastoris. The recombinant proteins demonstrated to be active against trypsin, while three of them were also active against chymotrypsin, and one against human leukocyte elastase. Comparative modeling and protein docking were used to further investigate interactions between two grass pea BBI isoforms and their target proteases. Thus two reliable 3D models have been proposed, representing two potential ternary complexes, each constituted of an inhibitor and its target enzymes.

A family of diverse Kunitz inhibitors from Echinococcus granulosus potentially involved in host-parasite cross-talk

The cestode Echinococcus granulosus, the agent of hydatidosis/echinococcosis, is remarkably well adapted to its definitive host. However, the molecular mechanisms underlying the successful establishment of larval worms (protoscoleces) in the dog duodenum are unknown. With the aim of identifying molecules participating in the E. granulosus-dog cross-talk, we surveyed the transcriptomes of protoscoleces and protoscoleces treated with pepsin at pH 2. This analysis identified a multigene family of secreted monodomain Kunitz proteins associated mostly with pepsin/H(+)-treated worms, suggesting that they play a role at the onset of infection. We present the relevant molecular features of eight members of the E. granulosus Kunitz family (EgKU-1 - EgKU-8). Although diverse, the family includes three pairs of close paralogs (EgKU-1/EgKU-4; EgKU-3/EgKU-8; EgKU-6/EgKU-7), which would be the products of recent gene duplications. In addition, we describe the purification of EgKU-1 and EgKU-8 from larval worms, and provide data indicating that some members of the family (notably, EgKU-3 and EgKU-8) are secreted by protoscoleces. Detailed kinetic studies with native EgKU-1 and EgKU-8 highlighted their functional diversity. Like most monodomain Kunitz proteins, EgKU-8 behaved as a slow, tight-binding inhibitor of serine proteases, with global inhibition constants (K(I) (*)) versus trypsins in the picomolar range. In sharp contrast, EgKU-1 did not inhibit any of the assayed peptidases. Interestingly, molecular modeling revealed structural elements associated with activity in Kunitz cation-channel blockers. We propose that this family of inhibitors has the potential to act at the E. granulosus-dog interface and interfere with host physiological processes at the initial stages of infection.

Stepwise disassembly and apparent nonstepwise reassembly for the oligomeric RbsD protein

Many cellular proteins exist as homo-oligomers. The mechanism of the assembly process of such proteins is still poorly understood. We have previously observed that Hsp16.3, a protein exhibiting chaperone-like activity, undergoes stepwise disassembly and nonstepwise reassembly. Here, the disassembly and reassembly of a nonchaperone protein RbsD, from Escherichia coli, was studied in vitro. The protein was found to mainly exist as decamers with a small portion of apparently larger oligomeric forms, both of which are able to refold/reassemble effectively in a spontaneous way after being completely unfolded. Disassembly RbsD intermediates including pentamers, tetramers, trimers, dimers, and monomers were detected by using urea-containing pore gradient polyacrylamide gel electrophoresis, while only pentamers were detected for its reassembly. The observation of stepwise disassembly and apparent nonstepwise reassembly for both a chaperone protein (Hsp16.3) and a nonchaperone protein (RbsD) strongly suggests that such a feature is most likely general for homo-oligomeric proteins.

Enzyme:substrate hydrogen bond shortening during the acylation phase of serine protease catalysis

Atomic resolution (<or=1.2 A) serine protease intermediate structures revealed that the strength of the hydrogen bonds between the enzyme and the substrate changed during catalysis. The well-conserved hydrogen bonds of antiparallel beta-sheet between the enzyme and the substrate become significantly shorter in the transition from a Michaelis complex analogue (Pontastacus leptodactylus (narrow-fingered crayfish) trypsin (CFT) in complex with Schistocerca gregaria (desert locust) trypsin inhibitor (SGTI) at 1.2 A resolution) to an acyl-enzyme intermediate (N-acetyl-Asn-Pro-Ile acyl-enzyme intermediate of porcine pancreatic elastase at 0.95 A resolution) presumably synchronously with the nucleophilic attack on the carbonyl carbon atom of the scissile peptide bond. This is interpreted as an active mechanism that utilizes the energy released from the stronger hydrogen bonds to overcome the energetic barrier of the nucleophilic attack by the hydroxyl group of the catalytic serine. In the CFT:SGTI complex this hydrogen bond shortening may be hindered by the 27I-32I disulfide bridge and Asn-15I of SGTI. The position of the catalytic histidine changes slightly as it adapts to the different nucleophilic attacker during the transition from the Michaelis complex to the acyl-enzyme state, and simultaneously its interaction with Asp-102 and Ser-214 becomes stronger. The oxyanion hole hydrogen bonds provide additional stabilization for acyl-ester bond in the acyl-enzyme than for scissile peptide bond of the Michaelis complex. Significant deviation from planarity is not observed in the reactive bonds of either the Michaelis complex or the acyl-enzyme. In the Michaelis complex the electron distribution of the carbonyl bond is distorted toward the oxygen atom compared to other peptide bonds in the structure, which indicates the polarization effect of the oxyanion hole.

Conformation and stability of elastase from Atlantic cod, Gadus morhua

Metal binding and conformational stability characteristics of psychrophilic elastase (ACE) from Atlantic cod (Gadus morhua) has been investigated. Chelation to Ca(2+) was found to be important for maintaining the biologically active conformation and for the thermal stability of the enzyme. However, presence of metal ions such as Zn(2+), Fe(3+) and Cu(2+) was found to inhibit its hydrolytic activity and so did the chelating agent EDTA. Both pH and guanidinium chloride induced denaturation of the enzyme was followed by monitoring the changes in the tryptophan fluorescence. ACE exhibited a simple two-state unfolding pattern in both acidic and basic conditions with the midpoint of transition at pH values 4.08 and 10.29, respectively. Guanidinium chloride and heat induced denaturation of the enzyme was investigated at two pH values, 5.50 and 8.00, wherein the enzyme possesses similar tertiary structure but differ in its hydrolytic activity. Guanidinium chloride induced denaturation indicated that the enzyme unfolds with a C(m) of 1.53 M at pH 8.0 and a DeltaG(H2O) of 6.91 kJ mol(-1) (28.65 J mol(-1) residue(-1)) which is the lowest reported for psychrophilic enzymes investigated till-date. However, at pH 5.50, DeltaG(H2O) value is slightly lowered by 0.65 kJ mol(-1) consistent with the observed increase in the apparent quenching constant obtained with acrylamide. On the other hand, increase in T(m) by 38.45 degrees C was observed for the enzyme at acid pH (5.50) in comparison to the heat induced unfolding at pH 8.0. The increase in the apparent T(m) has been attributed to the possible weak intermolecular association of the enzyme molecules at moderately high temperatures that is favoured by the increase in the accessible surface area / dynamics under acidic conditions. The stability characteristics of ACE have been compared with the available data for mesophilic porcine pancreatic elastase and possible mechanism for the low temperature adaptation of ACE has been proposed.

Trypsin specificity as elucidated by LIE calculations, X-ray structures, and association constant measurements

The variation in inhibitor specificity for five different amine inhibitors bound to CST, BT, and the cold-adapted AST has been studied by use of association constant measurements, structural analysis of high-resolution crystal structures, and the LIE method. Experimental data show that AST binds the 1BZA and 2BEA inhibitors 0.8 and 0.5 kcal/mole more strongly than BT. However, structural interactions and orientations of the inhibitors within the S1 site have been found to be virtually identical in the three enzymes studied. For example, the four water molecules in the inhibitor-free structures of AST and BT are channeled into similar positions in the S1 site, and the nitrogen atom(s) of the inhibitors are found in two cationic binding sites denoted Position1 and Position2. The hydrophobic binding contributions for all five inhibitors, estimated by the LIE calculations, are also in the same order (-2.1 +/- 0.2 kcal/mole) for all three enzymes. Our hypothesis is therefore that the observed variation in inhibitor binding arises from different electrostatic interactions originating from residues outside the S1 site. This is well illustrated by AST, in which Asp 150 and Glu 221B, despite some distance from the S1 binding site, lower the electrostatic potential of the S1 site and thus enhance substrate binding. Because the trends in the experimentally determined binding energies were reproduced by the LIE calculations after adding the contribution from long-range interactions, we find this method very suitable for rational studies of protein-substrate interactions.

Atypical Kunitz-type serine proteinase inhibitors produced by the ruminant placenta

Recently, an unusual family of genes was identified with expression confined to the trophoblast of ruminant ungulate species. The members of this family (the trophoblast Kunitz domain proteins, or TKDPs) are characterized by the presence of one or more similar, approximately 80-residue repeat sequences placed ahead of a Kunitz serine proteinase-inhibitor domain. To examine the specificity of the Kunitz moiety, the Kunitz domains of selected TKDPs and a control Kunitz protein, bovine pancreatic trypsin inhibitor (BPTI), were produced as glutathione S-transferase fusions, and their abilities to inhibit six serine proteinases were examined. Circular dichroism spectroscopy confirmed that the Kunitz fold was intact. Three of the TKDPs had unusual residues at their P1 "warhead" (ovine TKDP-1, Asn; bovine TKDP-3, Thr; and bovine TKDP-5, Ile) and exhibited no measurable inhibitory activity toward any of the proteinases. Three (ovine TKDP-3, bovine TKDP-3, and bovine TKDP-4) lacked the conserved cysteines at residues 14 and 38 that form one of the highly conserved disulfide bonds that are structurally important in all known mammalian Kunitz proteins. Ovine TKDP-3 and bovine TKDP-4 had P1 lysines and inhibited trypsin and plasmin with K(i) values only approximately 10-fold higher than that of BPTI. Bovine TKDP-2 had a P1 lysine and the three conserved disulfides, but it possessed an unusual residue (Asp) at P2. It exhibited no inhibitory activity. These data suggest that the function of the TKDP, like certain Kunitz proteins found in snake venoms, may not be in proteinase inhibition.

Identification of substrate binding sites in enzymes by computational solvent mapping

Enzyme structures determined in organic solvents show that most organic molecules cluster in the active site, delineating the binding pocket. We have developed algorithms to perform solvent mapping computationally, rather than experimentally, by placing molecular probes (small molecules or functional groups) on a protein surface, and finding the regions with the most favorable binding free energy. The method then finds the consensus site that binds the highest number of different probes. The probe-protein interactions at this site are compared to the intermolecular interactions seen in the known complexes of the enzyme with various ligands (substrate analogs, products, and inhibitors). We have mapped thermolysin, for which experimental mapping results are also available, and six further enzymes that have no experimental mapping data, but whose binding sites are well characterized. With the exception of haloalkane dehalogenase, which binds very small substrates in a narrow channel, the consensus site found by the mapping is always a major subsite of the substrate-binding site. Furthermore, the probes at this location form hydrogen bonds and non-bonded interactions with the same residues that interact with the specific ligands of the enzyme. Thus, once the structure of an enzyme is known, computational solvent mapping can provide detailed and reliable information on its substrate-binding site. Calculations on ligand-bound and apo structures of enzymes show that the mapping results are not very sensitive to moderate variations in the protein coordinates.

Electrostatic effects play a central role in cold adaptation of trypsin

Organisms that live in constantly cold environments have to adapt their metabolism to low temperatures, but mechanisms of enzymatic adaptation to cold environments are not fully understood. Cold active trypsin catalyses reactions more efficiently and binds ligands more strongly in comparison to warm active trypsin. We have addressed this issue by means of comparative free energy calculations studying the binding of positively charged ligands to two trypsin homologues. Stronger inhibition of the cold active trypsin by benzamidine and positively charged P1-variants of BPTI is caused by rather subtle electrostatic effects. The different affinity of benzamidine originates solely from long range interactions, while the increased binding of P1-Lys and -Arg variants of BPTI is attributed to both long and short range effects that are enhanced in the cold active trypsin compared to the warm active counterpart. Electrostatic interactions thus provide an efficient strategy for cold adaptation of trypsin.

Structural comparison of psychrophilic and mesophilic trypsins. Elucidating the molecular basis of cold-adaptation

Structural rationalizations for differences in catalytic efficiency and stability between mesophilic and cold-adapted trypsins have been suggested from a detailed comparison of eight trypsin structures. Two trypsins, from Antarctic fish and Atlantic cod, have been constructed by homology modeling techniques and compared with six existing X-ray structures of both cold-adapted and mesophilic trypsins. The structural analysis focuses on the cold trypsin residue determinants found in a more extensive comparison of 27 trypsin sequences, and reveals a number of structural features unique to the cold-adapted trypsins. The increased substrate affinity of the psychrophilic trypsins is probably achieved by a lower electrostatic potential of the S1 binding pocket particularly arising from Glu221B, and from the lack of five hydrogen bonds adjacent to the catalytic triad. The reduced stability of the cold trypsins is expected to arise from reduced packing in two distinct core regions, fewer interdomain hydrogen bonds and from a destabilized C-terminal alpha-helix. The helices of the cold trypsins lack four hydrogen bonds and two salt-bridges, and they have poorer van der Waals packing interactions to the body of the molecule, compared to the mesophilic counterparts.

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.

The crystal structures of the complexes between bovine beta-trypsin and ten P1 variants of BPTI

The high-resolution X-ray structures have been determined for ten complexes formed between bovine beta-trypsin and P1 variants (Gly, Asp, Glu, Gln, Thr, Met, Lys, His, Phe, Trp) of bovine pancreatic trypsin inhibitor (BPTI). All the complexes were crystallised from the same conditions. The structures of the P1 variants Asp, Glu, Gln and Thr, are reported here for the first time in complex with any serine proteinase. The resolution of the structures ranged from 1.75 to 2.05 A and the R-factors were about 19-20 %. The association constants of the mutants ranged from 1.5x10(4) to 1.7x10(13) M-1. All the structures could be fitted into well-defined electron density, and all had very similar global conformations. All the P1 mutant side-chains could be accomodated at the primary binding site, but relative to the P1 Lys, there were small local changes within the P1-S1 interaction site. These comprised: (1) changes in the number and dynamics of water molecules inside the pocket; (2) multiple conformations and non-optimal dihedral angles for some of the P1 side-chains, Ser190 and Gln192; and (3) changes in temperature factors of the pocket walls as well as the introduced P1 side-chain. Binding of the cognate P1 Lys is characterised by almost optimal dihedral angles, hydrogen bonding distances and angles, in addition to considerably lower temperature factors. Thus, the trypsin S1 pocket seems to be designed particularly for lysine binding.