Search

A congeneric series of benzamidine-type ligands with a central proline moiety and a terminal cycloalkyl group--linked by a secondary amine, ether, or methylene bridge--was synthesized as trypsin inhibitors. This series of inhibitors was investigated by isothermal titration calorimetry, crystal structure analysis in two crystal forms, and molecular dynamics simulations. Even though all of these congeneric ligands exhibited essentially the same affinity for trypsin, their binding profiles at the structural, dynamic, and thermodynamic levels are very distinct. The ligands display a pronounced enthalpy/entropy compensation that results in a nearly unchanged free energy of binding, even though individual enthalpy and entropy terms change significantly across the series. Crystal structures revealed that the secondary amine-linked analogs scatter over two distinct conformational families of binding modes that occupy either the inside or of the outside the protein's S3/S4 specificity pocket. In contrast, the ether-linked and methylene-linked ligands preferentially occupy the hydrophobic specificity pocket. This also explains why the latter ligands could only be crystallized in the conformationally restricting closed crystal form whereas the derivative with the highest residual mobility in the series escaped our attempts to crystallize it in the closed form; instead, a well-resolved structure could only be achieved in the open form with the ligand in disordered orientation. These distinct binding modes are supported by molecular dynamics simulations and correlate with the shifting enthalpic/entropic signatures of ligand binding. The examples demonstrate that, at the molecular level, binding modes and thermodynamic binding signatures can be very different even for closely related ligands. However, deviating binding profiles provide the opportunity to optimally address a given target., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

The large variety of serine protease inhibitors, available from various sources such as tissues, microorganisms, plants, etc., play an important role in regulating the proteolytic enzymes. The analysis of protease-inhibitor complexes helps in understanding the mechanism of action, as well as in designing inhibitors. Vasopressin, an anti-diuretic nonapeptide hormone, is found to be an effective inhibitor of trypsin, with a K(i) value of 5 nM. The crystal structure of the trypsin-vasopressin complex revealed that vasopressin fulfils all the important interactions for an inhibitor, without any break in the scissile peptide bond. The cyclic nature due to a disulfide bridge between Cys1 and Cys6 of vasopressin provides structural rigidity to the peptide hormone. The trypsin-binding site is located at the C terminus, while the neurophysin-binding site is at the N terminus of vasopressin. This study will assist in designing new peptide inhibitors. This study suggests that vasopressin inhibition of trypsin may have unexplored biological implications.

An extensive structural manifold of short hydrogen bond-mediated, active site-directed, serine protease inhibition motifs is revealed in a set of over 300 crystal structures involving a large suite of small molecule inhibitors (2-(2-phenol)-indoles and 2-(2-phenol)-benzimidazoles) determined over a wide range of pH (3.5-11.4). The active site hydrogen-bonding mode was found to vary markedly with pH, with the steric and electronic properties of the inhibitor, and with the type of protease (trypsin, thrombin or urokinase type plasminogen activator (uPA)). The pH dependence of the active site hydrogen-bonding motif is often intricate, constituting a distinct fingerprint of each complex. Isosteric replacements or minor substitutions within the inhibitor that modulate the pK(a) of the phenol hydroxyl involved in short hydrogen bonding, or that affect steric interactions distal to the active site, can significantly shift the pH-dependent structural profile characteristic of the parent scaffold, or produce active site-binding motifs unique to the bound analog. Ionization equilibria at the active site associated with inhibitor binding are probed in a series of the protease-inhibitor complexes through analysis of the pH dependence of the structure and environment of the active site-binding groups involved in short hydrogen bond arrays. Structures determined at high pH (>11), suggest that the pK(a) of His57 is dramatically elevated, to a value as high as approximately 11 in certain complexes. K(i) values involving uPA and trypsin determined as a function of pH for a set of inhibitors show pronounced parabolic pH dependence, the pH for optimal inhibition governed by the pK(a) of the inhibitor phenol involved in short hydrogen bonds. Comparison of structures of trypsin, thrombin and uPA, each bound by the same inhibitor, highlights important structural variations in the S1 and active sites accessible for engineering notable selectivity into remarkably small molecules with low nanomolar K(i) values.

Severe neurodegradative brain diseases, like Alzheimer, are tightly linked with proteolytic activity in the human brain. Proteinases expressed in the brain, such as human trypsin IV, are likely to be involved in the pathomechanism of these diseases. The observation of amyloid formed in the brain of transgenic mice expressing human trypsin IV supports this hypothesis. Human trypsin IV is also resistant towards all studied naturally occurring polypeptide inhibitors. It has been postulated that the substitution of Gly193 to arginine is responsible for this inhibitor resistance. Here we report the X-ray structure of human trypsin IV in complex with the inhibitor benzamidine at 1.7 A resolution. The overall fold of human trypsin IV is similar to human trypsin I, with a root-mean square deviation of only 0.5 A for all C(alpha) positions. The crystal structure reveals the orientation of the side-chain of Arg193, which occupies an extended conformation and fills the S2' subsite. An analysis of surface electrostatic potentials shows an unusually strong clustering of positive charges around the primary specificity pocket, to which the side-chain of Arg193 also contributes. These unique features of the crystal structure provide a structural basis for the enhanced inhibitor resistance, and enhanced substrate restriction, of human trypsin IV., (Copyright 2002 Elsevier Science Ltd.)

The native conformation of proteins in the serpin superfamily is metastable. In order to understand why serpins attain the native state instead of more stable conformations we have begun investigations into the equilibrium-unfolding of alpha(1)-antitrypsin. alpha(1)-Antitrypsin contains two tryptophan residues, Trp194 and Trp238, situated on the A and B beta-sheets, respectively. Site-directed mutagenesis was used to construct two single-tryptophan variants. Both variants were fully active and had similar secondary structure and stabilities to alpha(1)-antitrypsin. The denaturation of alpha(1)-antitrypsin and its variants was extremely similar when followed by far-UV CD, indicating the presence of a single intermediate. Fluorescence analysis of the unfolding behavior of each single tryptophan variant indicated that the sole tryptophan residue reported the structural changes within its immediate environment. These data suggest that the A beta-sheet is expanded in the intermediate state whilst no structural change around the B beta-sheet has occurred. In the urea-induced unfolded state, Trp238 does not become fully solvated, suggesting the persistence of structure around this residue. The implications of these data on the folding, misfolding and function of the serpin superfamily are discussed., (Copyright 2001 Academic Press.)

The binding of a series of low molecular weight ligands towards trypsin and thrombin has been studied by isothermal titration calorimetry and protein crystallography. In a series of congeneric ligands, surprising changes of protonation states occur and are overlaid on the binding process. They result from induced pK(a) shifts depending on the local environment experienced by the ligand and protein functional groups in the complex (induced dielectric fit). They involve additional heat effects that must be corrected before any conclusion on the binding enthalpy (DeltaH) and entropy (DeltaS) can be drawn. After correction, trends in both contributions can be interpreted in structural terms with respect to the hydrogen bond inventory or residual ligand motions. For all inhibitors studied, a strong negative heat capacity change (DeltaC(p)) is detected, thus binding becomes more exothermic and entropically less favourable with increasing temperature. Due to a mutual compensation, Gibbs free energy remains virtually unchanged. The strong negative DeltaC(p) value cannot solely be explained by the removal of hydrophobic surface portions of the protein or ligand from water exposure. Additional contributions must be considered, presumably arising from modulations of the local water structure, changes in vibrational modes or other ordering parameters. For thrombin, smaller negative DeltaC(p) values are observed for ligand binding in the presence of sodium ions compared to the other alkali ions, probably due to stabilising effects on the protein or changes in the bound water structure., (Copyright 2001 Academic Press.)

SFTI-1 is a recently discovered cyclic peptide trypsin inhibitor from sunflower seeds comprising 14 amino acid residues. It is the most potent known Bowman-Birk inhibitor and the only naturally occurring cyclic one. The solution structure of SFTI-1 has been determined by 1H-NMR spectroscopy and compared with a synthetic acyclic permutant. The solution structures of both are remarkably similar. The lowest energy structures from each family of 20 structures of cyclic and acyclic SFTI-1 have an rmsd over the backbone and heavy atoms of 0.29 A and 0.66 A, respectively. The structures consist of two short antiparallel beta-strands joined by an extended loop containing the active site at one end. Cyclic SFTI-1 also has a hairpin turn completing the cycle. Both molecules contain particularly stable arrangements of cross-linking hydrogen bonds between the beta-strands and a single disulfide bridge, making them rigid and well defined in solution. These stable arrangements allow both the cyclic and acyclic variants of SFTI-1 to inhibit trypsin with very high potencies (0.5 nM and 12.1 nM, respectively). The cyclic nature of SFTI-1 appears to have evolved to provide higher trypsin inhibition as well as higher stability. The solution structures are similar to the crystal structure of the cyclic inhibitor in complex with trypsin. The lack of a major conformational change upon binding suggests that the structure of SFTI-1 is rigid and already pre-organized for maximal binding due to minimization of entropic losses compared to a more flexible ligand. These properties make SFTI-1 an ideal platform for the design of small peptidic pharmaceuticals or pesticides., (Copyright 2001 Academic Press.)

We have determined the NMR structure in aqueous solution of a disulphide-cyclised 11-residue peptide that forms a stable beta-hairpin, incorporating a type VIb beta-turn. The structure is found to be extremely well ordered for a short peptide, with the 30 lowest energy simulated annealing structures having an average pairwise r.m.s. deviation of only 0.36 A over the backbone. All but three side-chains adopt distinct conformations, allowing a detailed analysis of their involvement in cross-strand interactions. The peptide sequence analysed originates from a previously reported study, which identified potent inhibitors of human leukocyte elastase from screening a combinatorial peptide library based on the short protein beta-sheet segment that forms the reactive site loop of Bowman-Birk inhibitors. A detailed comparison of the peptide's solution structure with the corresponding region in the whole protein structure reveals a very good correspondence not only for the backbone (r.m.s. deviation approximately 0.7 A) but also for the side-chains. This isolated beta-hairpin retains the biologically active "canonical conformation" typical of small serine proteinase inhibitor proteins, which explains why it retains inhibitory activity. Since the structural integrity is sequence-inherent and does not depend upon the presence of the remaining protein, this beta-hairpin represents an independent structural motif and so provides a useful model of this type of protein architecture and its relation to biological function. The relationship between the conformation of this beta-hairpin and its biological activity is discussed., (Copyright 2001 Academic Press.)

To establish the structural basis underlying the activity of a novel series of metal-chelate trypsin inhibitors, the structures of p-amidinosalicylidene-l-alaninato(aqua)copper(II) (1a), m-amidinosalicylidene-l-alaninato(aqua)copper(II) (1b), bis(p-amidinosalicylidene-l-alaninato)iron(III) (2a), and bis(m-amidinosalicylidene-l-alaninato)iron(III) (2b) bound to bovine beta-trypsin were studied by X-ray crystallography. The amidinium group of the inhibitor donates hydrogen bonds to Asp189, Gly219 and Ser190, as seen before in trypsin-benzamidine complexes. The copper(II) ion of 1a is situated away from trypsin's catalytic triad residues, and is octahedrally coordinated by a Schiff base and three water molecules. In contrast, the copper(II) ion of 1b is situated close to the catalytic triad and adopts a square pyramidal coordination geometry. The iron(III) ion of 2a is octahedrally coordinated by two Schiff base ligands and, like the copper(II) ion of 1a, is situated away from the catalytic triad. The p-amidinophenyl ring of a second Schiff base ligand of 2a is directed toward a hydrophobic groove formed by Trp215 and Leu99. Finally, the iron(III) ion of 2b appears to be replaced by magnesium(II), which is octahedrally coordinated by a Schiff base, Gln192 and two water molecules. One of the Schiff base ligands seen in the trypsin-2a complex or in the unbound form of 2b is replaced by water molecules and Gln192. His57 and Ser195 form water-mediated interactions with the magnesium(II) ion of 2b, and Ser195 also forms a hydrogen bond with the phenolic oxygen atom of the Schiff base ligand. These structures reveal a novel mode of interaction between metal-chelate inhibitors and serine proteases, thus providing a structural basis for the development of more potent inhibitors against a variety of trypsin-like enzymes., (Copyright 2001 Academic Press.)

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.

Ecotin is a dimeric serine protease inhibitor from Escherichia coli which binds proteases to form a hetero-tetramer with three distinct interfaces: an ecotin-ecotin dimer interface, a larger primary ecotin-protease interface, and a smaller secondary ecotin-protease interface. The contributions of these interfaces to binding and inhibition are unequal. To investigate the contribution and adaptability of each interface, we have solved the structure of two mutant ecotin-trypsin complexes and compared them to the structure of the previously determined wild-type ecotin-trypsin complex. Wild-type ecotin has an affinity of 1 nM for trypsin, while the optimized mutant, ecotin Y69F, D70P, which was found using phage display technologies, inhibits rat trypsin with a K(i) value of 0.08 nM. Ecotin 67-70A, M84R which has four alanine substitutions in the ecotin-trypsin secondary binding site, along with the M84R mutation at the primary site, has a K(i) value against rat trypsin of 0.2 nM. The structure of the ecotin Y69F, D70P-trypsin complex shows minor structural changes in the ecotin-trypsin tetramer. The structure of the ecotin 67-70A, M84R mutant bound to trypsin shows large deviations in the tertiary and quaternary structure of the complex. The trypsin structure shows no significant changes, but the conformation of several loop regions of ecotin are altered, resulting in the secondary site releasing its hold on trypsin. The structure of several regions previously considered to be rigid is also significantly modified. The inherent flexibility of ecotin allows it to accommodate these mutations and still maintain tight binding through the compromises of the protein-protein interfaces in the ecotin-trypsin tetramer. A comparison with two recently described ecotin-like genes from other bacteria suggests that these structural and functional features are conserved in otherwise distant bacterial lineages., (Copyright 2000 Academic Press.)

Understanding molecular recognition on a structural basis is an objective with broad academic and applied significance. In the complexes of serine proteases and their proteinaceous inhibitors, recognition is governed mainly by residue P1 in accord with primary serine protease specificity. The bifunctional soybean Bowman-Birk inhibitor (sBBI) should, therefore, interact at LysI16 (subdomain 1) with trypsin and at LeuI43 (subdomain 2) with chymotrypsin. In contrast with this prediction, a 2:1 assembly with trypsin was observed in solution and in the crystal structure of sBBI in complex with trypsin, determined at 2.3 A resolution by molecular replacement. Strikingly, P1LeuI43 of sBBI was fully embedded into the S(1) pocket of trypsin in contrast to primary specificity. The triple-stranded beta-hairpin unique to the BBI-family and the surface loops surrounding the active site of the enzyme formed a protein-protein-interface far extended beyond the primary contact region. Polar residues, hydrophilic bridges and weak hydrophobic contacts were predominant in subdomain 1, interacting specifically with trypsin. However, close hydrophobic contacts across the interface were characteristic of subdomain 2 reacting with both trypsin and chymotrypsin. A Met27Ile replacement shifted the ratio with trypsin to the predicted 1:1 ratio. Thus, the buried salt-bridge responsible for trypsin specificity was stabilised in a polar, and destabilized in a hydrophobic, environment. This may be used for adjusting the specificity of protease inhibitors for applications such as insecticides and cancer chemopreventive agents., (Copyright 2000 Academic Press.)

A previous comprehensive analysis of the pH dependence of native-state amide hydrogen (NH) exchange in turkey ovomucoid third domain (OMTKY3) yielded apparent opening and closing rate constants (k(op) and k(cl)) at 14 NH groups involved in global conformational changes. This analysis has been extended to 18 additional slowly exchanging NH groups. Quench-flow experiments were performed to monitor NH exchange in native OMTKY3 from neutral to very alkaline pH ( approximately 12) conditions. Above pH 10 the mechanism of exchange switched from one governed by a rapid equilibrium preceding the chemistry of exchange (i.e. EX2 exchange), to one where exchange was limited by the rate of opening (i.e. EX1 exchange). Kinetics of solvent exposure are now known for nearly all backbone NH groups in native OMTKY3, yielding rate constants that span five orders of magnitude, 0.004 to 200 s(-1)., (Copyright 2000 Academic Press.)

Proteinaceous serine proteinase inhibitors are widespread throughout the plant kingdom where they play an important role in protection against pests and pathogens. Here, we describe the isolation and characterisation of a novel 14 amino acid residue cyclic peptide from sunflower seeds, which is a potent inhibitor of trypsin (Ki=100 pM). The crystal structure of this peptide in complex with bovine beta-trypsin shows both sequence and conformational similarity with the trypsin-reactive loop of the Bowman-Birk family of serine proteinase inhibitors. This inhibitor, however, is unique in being monofunctional, cyclic and far shorter (14 amino acid residues) than inhibitors belonging to this family (typically 60-70 amino acid residues). The high potency of this peptide is likely to arise from the considerable structural rigidity achieved through its cyclic nature which is further stabilised by a single internal disulphide bond. This study helps delineate the minimal unit required for effective peptide inhibitors of serine proteinases, and will assist in the further design of inhibitors to this widespread class of enzymes., (Copyright 1999 Academic Press.)

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., (Copyright 1999 Academic Press.)

In order to study the effect of internal dynamics on the accuracy of NMR structures in detail, we generated NOE distance data from a long molecular dynamics trajectory of BPTI. Cross-relaxation rates were calculated from the trajectory by analysis of the appropriate proton-proton vector autocorrelation functions. A criterion for the convergence of correlation functions was developed, and the analysis was restricted to those correlation functions that had converged within the simulation time. Effective distances were determined from the calculated cross-relaxation rates. Internal dynamics affected the derived distances in a realistic way, since they were subject both to radial averaging (which increases the cross-relaxation rate) and angular averaging (which decreases the cross-relaxation rate). The comparison of the effective distances with average distance between the protons during the trajectory showed that for most the effects of angular and distance averaging essentially cancel out. For these distances, the effective distance derived from an NOE is therefore a very good estimate of the average distance, or the distance in the average structure. However, for about 10% of the distances, the effective distance was more than 10% larger than the average distance, while for about 5%, it was more than 10% smaller, in some cases by more than 2 A. Little correlation is observed between the effects on cross-relaxation rates to different protons of the same residue. The results of this analysis have implications for the way structures are calculated from NOE distance data. For many distances, the assumption of a rigid structure is valid, and large error bounds would result in the loss of too much information content. On the other hand, the error bounds very often employed are not wide enough for some of the effects seen in our study., (Copyright 1999 Academic Press.)

Ecotin, a dimeric serine protease inhibitor from Escherichia coli, is a novel platform for inhibitor design. An approach using the three-dimensional structure of the ecotin-trypsin complex to guide combinatiorial design efforts was taken to create potent bidentate ecotin inhibitors for trypsin and human urokinase-type plasminogen activator (uPA). The ecotin surface loop that was redesigned is composed of residues 67 to 70 (60 s loop), and binds to the target protease at a region 25 A from the enzyme active site. Two ecotin phage display libraries were constructed to exploit the binding interactions at the 60 s loop. The ecotin 60X4 library, in which residues 67 to 70 of ecotin were randomized, was panned against rat and bovine trypsin in parallel for four rounds. Panning against bovine trypsin resulted in enrichment of ecotin phage but did not yield a consensus sequence. Panning against rat trypsin resulted in enrichment as well as the ecotin consensus sequence WGFP at positions 67 to 70. The variant ecotin encoded by this sequence inhibited rat trypsin at 80 pM, a 12-fold improvement over ecotin wild-type (WT). A second generation library, ecotin M84R+60X4 including an additional methionine to arginine substitution at position 84 in the primary binding site of ecotin, was generated for panning against uPA and rat trypsin. Panning against rat trypsin resulted in enrichment but no consensus sequence. Panning against uPA resulted in enrichment as well as the different ecotin consensus sequence WGYR at positions 67 to 70. Ecotin M84R+D70R bound to uPA at 50 pM, a 56,000-fold increase in binding compared to ecotin WT. Furthermore, ecotin M84R+D70R achieved a 13,680-fold preference of specificity towards uPA versus rat trypsin. The fact that the 60 s loop of ecotin plays different roles in binding to trypsin and uPA suggests this site can be used to introduce specificity and potency for other members of the serine proteases with a chymotrypsin fold., (Copyright 1998 Academic Press Limited.)

Bikunin is a serine protease inhibitor found in the blood serum and urine of humans and other animals. Its sequence shows internal repetition, suggesting that it contains two domains that resemble bovine pancreatic trypsin inhibitor (BPTI). A fragment of bikunin has been crystallised, its structure solved and subsequently refined against 2.5 A data. The two BPTI-like domains pack closely together and are related by an approximate 60 degrees rotation combined with a translation. These domains are very similar to each other and other proteins with this fold. The largest variations occur in the loops responsible for protease recognition. The loops of the first domain are unobstructed by the remaining protein. However, the loops of the second domain are close to the first domain and it is possible that protease binding may be affected or, in some cases, abolished by the presence of the first domain. Thus, cleavage of the two domains could alter the substrate specificity of domain II. Bikunin has a hydrophobic patch close to the N terminus of domain I, which is the most likely site for cell-surface receptor binding. In addition, there is a basic patch at one end of domain II that may be responsible for the inhibition of calcium oxalate crystallization in urine.

Native bovine pancreatic trypsin inhibitor (BPTI) contains three disulfide bonds: Cys5-Cys55, Cys14-Cys38 and Cys30-Cys51. Correct cysteine pairing, native structure, and full anti-proteinase activity can be restored in the process of oxidative refolding of reduced BPTI. Oxidative refolding starts with the formation of single disulfide intermediates. All 15 single-disulfide variants of BPTI (three native and 12 non-native combinations) have been expressed in Escherichia coli. In each variant the remaining four cysteine residues were replaced by alanine. Four of these variants are shown here to inhibit bovine beta-trypsin: three of them contain native and one non-native (Cys5-Cys51) disulfide. All but one (Cys5-Cys55) variant are slowly digested by the enzyme, therefore measurements were performed at pH 4.0, at which trypsin activity is low. Binding constants of these four single disulfide variants were at least two orders of magnitude lower than for native BPTI. Remarkably, in some of the variants the binding constants were found to be higher for the reduced rather than for the oxidized form of the variant. Also for the fully reduced native BPTI, determined here, the binding constant is of considerable value. Two sets of control experiments demonstrated that the binding of reduced native BPTI to trypsin is specific. In the first, mutation of Lys15 (P1 position) in the binding loop abolished binding of the reduced forms to trypsin. In the second, the binding of reduced native BPTI to anhydrotrypsin yielded the expected UV difference spectra. In general, the results obtained indicate that the inhibitor activity can be induced even in the reduced protein. This activity is not a local effect, such as the nature of residues surrounding the binding loop, but rather is induced by residual structure in the unfolded protein. This structure has been shown to consist of a set of hydrophobic residues and the data presented here indicate that reduced cysteine residues provide further stabilization of such a hydrophobic cluster. On the other hand, improper pairing of the cysteine residues in non-native single disulfide variants destabilizes the enzyme-inhibitor complex by inducing deformations of the binding loop region.

Winged bean albumin-1 (WBA) is the main seed albumin of Psophocarpus tetragonolobus, a legume that has excellent potential as a protein-rich food source for humid tropical climates. WBA crystallises in a tetragonal space group and the structure was solved by X-ray crystallography with a combination of multiple isomorphous replacement using four heavy atom derivatives and molecular replacement with a model based on the structure of Erythrina caffra trypsin inhibitor (ETI). Refinement of the structure proceeded to 1.8 A. WBA has a beta-trefoil fold, similar to that found in the STI-Kunitz type trypsin inhibitors. The final structure has an overall R-factor of 19% for 15 to 1.8 A resolution data, all residues in the allowed regions of the Ramachandran plot, and good agreement with ideal geometry. WBA has sequence similarity with the STI-Kunitz trypsin inhibitors, including the apparent conservation of the functional reactive site residue, lysine 64, at the position of the scissile bond (position P1) in the STI-Kunitz type trypsin inhibitors, however, WBA does not inhibit trypsin. The reason for the lack of inhibitory activity against trypsin is clearly evident from the structure. The loop corresponding to the inhibitory loop in the STI-Kunitz trypsin inhibitors does not conform to the canonical conformation of the inhibitory loops of the "small inhibitors". The lysine residue assigned to the P1 position from sequence alignments is instead part of a four amino acid insertion between residues structurally equivalent to residues P1 and P2 of the inhibitors.

The disulphide folding pathway of bovine pancreatic trypsin inhibitor (BPTI) revealed that the native conformation is still stable in each intermediate state with two native disulphide linkages, in the absence of each of the corresponding third disulphide bonds. This is thought to be a consequence of the extreme stability of the native BPTI conformation. The current study addresses the question of whether the native-like conformation would be populated significantly at the two-disulphide stage in disulphide refolding if the final structure is less stable than in the case of BPTI. Dendrotoxin K from black mamba venom provides a good model to test this, since it contains the BPTI fold and was shown to fold predominantly via the same pathway, but its native conformation is stable than that of BPTI. The conformation of a chemically trapped two-disulphide intermediate in the disulphide refolding of dendrotoxin K, with blocking groups on Cys5 and Cys55 and disulphide bonds between Cys30 and Cys51, and Cys14 and Cys38, respectively, has been determined by 1H NMR spectroscopy and compared to those of the native protein and of the corresponding intermediate in BPTI. The analysis reveals that the dendrotoxin K intermediate adopts a partly-folded conformation, in contrast to the quasi-native conformation of the corresponding BPTI intermediate. It is similar to the partly-folded conformation of the BPTI intermediate with just the Cys30-Cys-51 disulphide bond, but with a more fixed conformation in the region of the Cys14-Cys38 disulphide bond. The destabilisation of the fully native conformation of the dendrotoxin K intermediate, relative to BPTI, appears to reduce the cooperativeity of the folding process.

We present a low resolution lattice model for which we can exhaustively generate all possible compact backbone conformations for small proteins. Using simple structural and energetic criteria, for a variety of proteins, we can select for lattice structures that have significant similarities with their known native structures. Our energetic parameters are based on pairwise amino acid contact frequencies in a database of experimentally determined structures. A key step in our method involves the threading of a sequence onto every lattice model, such that a locally optimal pattern of tertiary interactions is formed. We evaluate our results against statistics collected for structures covering all of conformational space, and against statistics collected for permuted sequences. Despite the low resolution of the model, our low energy structures contain many native features. These results indicate that the overall pattern of hydrophobicity of a sequence significantly constrains the range of folds that sequence is likely to adopt.

Erythrina trypsin inhibitor (ETI) from the seeds of Erythrina caffra is a high-affinity inhibitor of trypsin, chymotrypsin and tissue plasminogen activator. Its 172 amino acid polypeptide chain is stabilized in its compact, native state by two disulfide bonds. In spite of their conservation in all trypsin inhibitors of the soybean trypsin inhibitor (STI-Kunitz) family, their state of oxidation is essential only for protein stability but not for inhibitory function. Reduction/reoxidation of ETI in the presence of glutathione reshuffling buffer (GSH/GSSG; pH 8.3) not only allows the inhibitor to be restored in its native structure, but also does not interfere with its binding affinity; carboxymethylation or carboxamidomethylation of the free thiol groups does not affect K1 significantly (for trypsin (KI)ETIox = 2.3 nM, (KI)ETICM = 1.9 nM; for chymotrypsin (KI)ETIox = 30 microM, (KI)ETICM = 25 microM). The two cystine cross-bridges in the native ETI lead to enhanced stability toward pH and chaotropic agents. As taken from intrinsic protein fluorescence at acid pH and varying ionic strength (pH < 4, I = 0.01 to 0.15 M), the oxidized inhibitor retains its spectral properties, whereas reduced and carboxymethylated or carboxamidomethylated ETI undergo at least partial denaturation. At alkaline pH, the oxidized protein is stable up to pH 9.5, whereas the reduced protein undergoes structural alterations at pH > 7, reaching a final plateau at pH 10.0 to 10.5. In the case of urea (U) or guanidinium chloride (GdmCl) denaturation at pH 7.0, structural transitions of the oxidized inhibitor show "hysteresis" with half-concentrations (cU)1/2 approximately 10 M and (cGdmCl)1/2 approximately 4.5 M for denaturation, and (cU)1/2 = 4.7 M and (cGdmCl)1/2 = 1.5 M for renaturation. In contrast, the reduced (and chemically modified) inhibitors exhibit true equilibrium transitions at (cU)1/2 = 0.9 M and (cGdmCl)1/2 = 0.5 M, respectively. Reduction/reoxidation in the absence and in the presence of denaturants (GdmCl) can also be applied to ETI covalently attached to a solid matrix.

We describe what may be the most accurate approach currently available for the calculation of the pKas of ionizable groups in proteins. The accuracy is assessed by comparison of computed pKas with 60 measured pKas in a total of seven proteins. The overall root-mean-square error is 0.89 pKa units. Linear regression analysis of computed versus measured pKas yields a slope of 0.95, y-intercept of -0.02 and a correlation coefficient of 0.96. The proposed approach also picks out many of the shifted pKas of groups in enzyme active sites and special salt bridges. However, it does yield several over-shifted pKas and tends to underestimate pKa shifts which result from desolvation effects. We examine the ability of the new approach to reproduce the dependence of protein stability upon pH, using the ionization polynomial formalism. Overall features of the stability curves are reproduced, but the quantitative agreement is not particularly good. The reasons for the disagreement may have to do both with insufficient accuracy in the theory and with uncertainty in the nature of the unfolded state of proteins. The methodology described here is based upon finite difference solutions of the Poisson-Boltzmann equation. Its success depends upon the use of the rather high protein dielectric constant of 20. However, theoretical considerations and the fact that pKa shifts which result from desolvation are underestimated here imply that the dielectric constant of the protein interior actually is lower than 20. We suggest that the high protein dielectric constant improves the overall agreement with experiment because it accounts approximately for phenomena which tend to mitigate pKa shifts and which are not specifically included in the model. These include conformational relaxation and specific ion-binding. Future models based upon a low protein dielectric constant and treating such phenomena explicitly might yield improved agreement with experiment.

The crystal structure of Bowman-Birk type protease inhibitor A-II from peanut was refined at 2.3 A resolution using a restrained least-squares method. The crystallographic R-factor is 0.196 for 7697 reflections with F > 3 sigma (F) in the range from 6.0 to 2.3 A resolution. Two molecules in an asymmetric unit are independently refined and, their structures are compared with each other. The inhibitor molecule has an elongated shape with two reactive sites, one at both ends of the longest dimension. As a secondary structure, a 4-stranded beta-sheet-like structure is found, in which two water molecules bind two 2-stranded beta-sheets together with six hydrogen bonds. The molecule is constructed by two homologous domains which are related by an intramolecular pseudo 2-fold axis. The structure and atomic B-factors of peptide loops containing a reactive site were compared with that of adzuki bean Bowman-Birk type inhibitor in the complex with bovine beta-trypsin. This comparison shows that no significant structural change occurs in the reactive site of inhibitor at the formation of the inhibitor-protease complex, but structural rigidity around the reactive site seems to increase.

The natural gene for bovine pancreatic trypsin inhibitor (BPTI) was expressed by in vitro transcription/translation systems as the 100-residue pre-proBPTI, with a signal peptide for translocation into the endoplasmic reticulum. Expression in the presence of microsomes defined the site of co-translational cleavage of the signal peptide. The resulting proBPTI in the microsomes consists of the 58 residues of mature BPTI, plus an additional 13 residues at the N terminus, including a cysteine residue at position -10, and seven residues at the C terminus. ProBPTI remained in the unfolded, reduced form within microsomes when synthesized under reducing conditions, but folded and formed disulphide bonds rapidly when the disulphide form of glutathione was added. Complete folding could occur within about one minute, even when residue Cys10 was replaced by Ser. The structure of proBPTI was determined by circular dichroism and two-dimensional NMR and found to be that of mature BPTI with flexible extensions on both termini. Its inhibition of the activity of alpha-chymotrypsin was indistinguishable from that of the mature protein. The extensions of the precursor appeared to play only very minor roles in refolding in vitro under conditions where folding and disulphide bond formation are coupled. Under pH and redox conditions thought to reflect those in vivo, complete folding and disulphide bond formation required several hours. Addition of protein disulphide isomerase to in vitro folding experiments caused substantial and similar increases in the rate of formation of the fully folded state for both mature BPTI and proBPTI; the half time for folding to the native state was reduced to approximately two minutes, which is comparable to that occurring in microsomes. The absence of substantial effects of the N and C-terminal extensions on the protein structure, inhibitor activity and refolding leaves their functional roles to be discovered.

A synthetic peptide corresponding to the 15 N-terminal residues of bovine pancreatic trypsin inhibitor, with serine replacing the two cysteine residues, has been characterized by 1H-nuclear magnetic resonance spectroscopy. This peptide has a very disordered conformation that is essentially the same when it is part of the analogue of the (30-51) one-disulphide intermediate in folding. This confirms the conclusion of a previous paper, that the (30-51) intermediate is partially folded, with the N-terminal segment disordered. Local elements of non-random conformation were observed in the peptide. Especially prominent was an apparently electrostatic interaction between the face of the aromatic ring of Tyr10 and the amide group of Gly12, which caused the latter to have a very anomalous chemical shift. A similar interaction was observed in shorter peptides, especially in tetrapeptides with the sequences Tyr/Phe-X-Gly-Y. The local nature of this interaction indicates that it should be a general feature in peptides and in unfolded proteins with such a sequence.

The crystal structure of the complex formed by porcine beta-trypsin with the MCTI-A inhibitor (Momordica charantia, Linn. Cucurbitaceae) has been determined at 1.6 A resolution using the molecular replacement method. The sequence of MCTI-A was determined by recognizing the electron density, and shows that MCTI-A is a member of the squash family of trypsin inhibitors. We report the first high-resolution structure of porcine beta-trypsin. Detailed comparisons have been made on the overall structure, solvent structure and active-site geometries between this complex and bovine beta-trypsin and its complexes. On the basis of our results, we discuss the interaction patterns between inhibitor and trypsin. Unlike other complex structures formed by bovine trypsin with inhibitors, no out-of-plane distortion around the inhibitor's scissible peptide was observed. The role of the trypsin catalytic triad is also discussed on the basis of this structure.

The solution structures of two proteins (CMTI-I, a trypsin inhibitor from Cucurbita maxima, and hisactophilin, an actin binding protein of 118 amino acids) have been determined based on the NOE data derived solely from the homonuclear 3D NOE-NOE magnetic resonance spectroscopy. Two different approaches for extraction of the structural information from the 3D NOE-NOE experiment were tested. One approach was based on the transformation of the 3D intensities into distance constraints. In the second, and more robust approach, the 3D NOE intensities were used directly in structure calculations, without the need to transform them into distance constraints. A new 2D potential function representing the 3D NOE-NOE intensity was developed and used in the simulated annealing protocol. For CMTI-I, a comparison between structures determined with the 3D NOE-NOE method and various 2D NOE approaches was carried out. The 3D data set allowed better definition of the structures than was previously possible with the 2D NOE procedures that used the isolated two-spin approximation to derive distance information.

An analogue of the bovine pancreatic trypsin inhibitor (BPTI) folding intermediate that contains only the disulphide bond between Cys5 and Cys55 has been prepared in Escherichia coli by protein engineering methods, with the other four Cys residues replaced by Ser. Two-dimensional 1H nuclear magnetic resonance studies of the analogue have resulted in essentially complete resonance assignments of the folded form of the protein. The folded protein has a compact conformation that is structurally very similar to that of native BPTI, although there are subtle differences and the folded conformation is not very stable. Approximately half of the protein molecules are unfolded at 3 degrees C, and this proportion increases at higher temperatures. The folded and unfolded conformations are in slow exchange. The conformational properties of the analogue can explain many aspects of the kinetic role that the normal (5-55) intermediate plays in the folding of BPTI.

An analogue of the BPTI folding intermediate that contains only the disulphide bonds between Cys14 and Cys38 and between Cys30 and Cys51 has been prepared in Escherichia coli by protein engineering methods. The other two Cys residues of native BPTI (at positions 5 and 55) have been replaced by Ser. Essentially complete proton resonance assignments of the analogue were obtained by employing two-dimensional 1H nuclear magnetic resonance techniques. The intermediate has a more extended conformation in the N-terminal (residues 1 to 7) region and there are other differences in the C-terminal (residues 55 to 58) region. The remainder of the protein is substantially identical to native BPTI. The conformational properties of the analogue can explain several aspects of the kinetic role that the normal (14-38, 30-51) intermediate plays in the folding of BPTI.

A bifunctional alpha-amylase/trypsin inhibitor that has two binding sites has been purified from ragi. The inhibitor has been crystallized from its ammonium sulphate solution by the vapour diffusion method. The crystals belong to the orthogonal space group P2(1)2(1)2(1) with unit cell dimensions a = 30.49 A, b = 56.30 A, c = 73.65 A and Z = 4.

The structure of the small squash trypsin inhibitor CMTI-I is refined by directly minimizing the difference between the observed two-dimensional nuclear Overhauser enhancement (NOE) intensities and those calculated by the full relaxation matrix approach. To achieve this, a term proportional to this difference was added to the potential energy function of the molecular dynamics program X-PLOR. Derivatives with respect to atomic co-ordinates are calculated analytically. Spin diffusion effects are thus accounted for fully during the refinement. Initial structures for the refinement were those determined recently by solution nuclear magnetic resonance using the isolated two-spin approximation to derive distance range estimates. The fits to the nuclear magnetic resonance data improve significantly with only small shifts in the refined structures during a few cycles of conjugate gradient minimization. However, larger changes (approximately 1 A) in the conformation occur during simulated annealing, which is accompanied by a further reduction of the difference between experimental and calculated two-dimensional NOE intensities. The refined structures are closer to the X-ray structure of the inhibitor complexed with trypsin than the initial structures. The root-mean-square difference for backbone atoms between the initial structures and the X-ray structure is 0.96 A, and that between the refined structures and the X-ray structure 0.61 A.