Genetic dissection of pancreatic trypsin inhibitor

Genetic selection for enhanced folding in vivo targets the Cys14-Cys38 disulfide bond in bovine pancreatic trypsin inhibitor

The periplasm provides a strongly oxidizing environment; however, periplasmic expression of proteins with disulfide bonds is often inefficient. Here, we used two different tripartite fusion systems to perform in vivo selections for mutants of the model protein bovine pancreatic trypsin inhibitor (BPTI) with the aim of enhancing its expression in Escherichia coli. This trypsin inhibitor contains three disulfides that contribute to its extreme stability and protease resistance. The mutants we isolated for increased expression appear to act by eliminating or destabilizing the Cys14-Cys38 disulfide in BPTI. In doing so, they are expected to reduce or eliminate kinetic traps that exist within the well characterized in vitro folding pathway of BPTI. These results suggest that elimination or destabilization of a disulfide bond whose formation is problematic in vitro can enhance in vivo protein folding. The use of these in vivo selections may prove a valuable way to identify and eliminate disulfides and other rate-limiting steps in the folding of proteins, including those proteins whose in vitro folding pathways are unknown.

Functional and structural roles of the Cys14-Cys38 disulfide of bovine pancreatic trypsin inhibitor

The disulfide bond between Cys14 and Cys38 of bovine pancreatic trypsin inhibitor lies on the surface of the inhibitor and forms part of the protease-binding region. The functional properties of three variants lacking this disulfide, with one or both of the Cys residues replaced with Ser, were examined, and X-ray crystal structures of the complexes with bovine trypsin were determined and refined to the 1.58-A resolution limit. The crystal structure of the complex formed with the mutant with both Cys residues replaced was nearly identical with that of the complex containing the wild-type protein, with the Ser oxygen atoms positioned to replace the disulfide bond with a hydrogen bond. The two structures of the complexes with single replacements displayed small local perturbations with alternate conformations of the Ser side chains. Despite the absence of the disulfide bond, the crystallographic temperature factors show no evidence of increased flexibility in the complexes with the mutant inhibitors. All three of the variants were cleaved by trypsin more rapidly than the wild-type inhibitor, by as much as 10,000-fold, indicating that the covalent constraint normally imposed by the disulfide contributes to the remarkable resistance to hydrolysis displayed by the wild-type protein. The rates of hydrolysis display an unusual dependence on pH over the range of 3.5-8.0, decreasing at the more alkaline values, as compared with the increased hydrolysis rates for normal substrates under these conditions. These observations can be accounted for by a model for inhibition in which an acyl-enzyme intermediate forms at a significant rate but is rapidly converted back to the enzyme-inhibitor complex by nucleophilic attack by the newly created amino group. The model suggests that a lack of flexibility in the acyl-enzyme intermediate, rather than the enzyme-inhibitor complex, may be a key factor in the ability of bovine pancreatic trypsin inhibitor and similar inhibitors to resist hydrolysis.

Rigidification of a flexible protease inhibitor variant upon binding to trypsin

The Tyr35-->Gly replacement in bovine pancreatic trypsin inhibitor (BPTI) has previously been shown to dramatically enhance the flexibility of the trypsin-binding region of the free inhibitor and to destabilize the interaction with the protease by about 3 kcal/mol. The effects of this replacement on the enzyme-inhibitor interaction were further studied here by X-ray crystallography and isothermal titration calorimetry (ITC). The co-crystal structure of Y35G BPTI bound to trypsin was determined using 1.65 A resolution X-ray diffraction data collected from cryopreserved crystals, and a new structure of the complex with wild-type BPTI under the same conditions was determined using 1.62 A data. These structures reveal that, in contrast to the free protein, Y35G BPTI adopts a conformation nearly identical with that of the wild-type protein, with a water-filled cavity in place of the missing Tyr side-chain. The crystallographic temperature factors for the two complexes indicate that the mutant inhibitor is nearly as rigid as the wild-type protein when bound to trypsin. Calorimetric measurements show that the change in enthalpy upon dissociation of the complex is 2.5 kcal/mol less favorable for the complex containing Y35G BPTI than for the complex with the wild-type inhibitor. Thus, the destabilization of the complex resulting from the Y35G replacement is due to a more favorable change in entropy upon dissociation. The heat capacity changes for dissociation of the mutant and wild-type complexes were very similar, suggesting that the entropic effects probably do not arise from solvation effects, but are more likely due to an increase in protein conformational entropy upon dissociation of the mutant inhibitor. These results define the biophysical role of a highly conserved core residue located outside of a protein-binding interface, demonstrating that Tyr35 has little impact on the trypsin-bound BPTI structure and acts primarily to define the structure of the free protein so as to maximize binding affinity.

Identification of a Residue Critical for Maintaining the Functional Conformation of BPTI

The effects of amino acid replacements on the backbone dynamics of bovine pancreatic trypsin inhibitor (BPTI) were examined using 15N NMR relaxation experiments. Previous studies have shown that backbone amide groups within the trypsin-binding region of the wild-type protein undergo conformational exchange processes on the micros time scale, and that replacement of Tyr35 with Gly greatly increases the number of backbone atoms involved in such motions. In order to determine whether these mutational effects are specific to the replacement of this residue with Gly, six additional replacements were examined in the present study. In two of these, Tyr35 was replaced with either Ala or Leu, and the other four were single replacements of Tyr23, Phe33, Asn43 or Asn44, all of which are highly buried in the native structure and conserved in homologous proteins. The Y35A and Y35L mutants displayed dynamic properties very similar to those of the Y35G mutant, with the backbone segments including residues 10-19 and 32-44 undergoing motions revealed by enhanced 15N transverse relaxation rates. On the other hand, the Y23L, N43G and N44A substitutions caused almost no detectable changes in backbone dynamics, on either the ns-ps or ms-micros time scales, even though each of these replacements significantly destabilizes the native conformation. Replacement of Phe33 with Leu caused intermediate effects, with several residues that have previously been implicated in motions in the wild-type protein displaying enhanced transverse relaxation rates. These results demonstrate that destabilizing amino acid replacements can be accommodated in a native protein with dramatically different effects on conformational dynamics and that Tyr35 plays a particularly important role in defining the conformation of the trypsin-binding site of BPTI.

Mutational analysis of hydrogen bonding residues in the BPTI folding pathway

Nine BPTI variants with replacements that remove one or more hydrogen bonds from the native protein were constructed, and the folding pathways of these proteins were determined by isolating and identifying the disulfide-bonded intermediates that accumulated during unfolding and refolding. The forward and reverse rate constants for the individual steps in the folding pathways for each protein were measured, providing a detailed description of the energetic effects of the substitutions. The native forms of eight of the nine variants were measurably destabilized, by 1-7 kcal/mol (1 cal=4.184 J), with an average effect of 1.6 kcal/mol per hydrogen bond removed. The folding pathways for the variants were found to be similar to that previously described for the wild-type protein, with the kinetically preferred mechanism involving intramolecular rearrangements of intermediates with two disulfide bonds. Some of the substitutions, however, significantly destabilized the major intermediates and broadened the distribution of species with one or two disulfide bonds, thus identifying residues that play important roles in stabilizing the normal intermediates and defining specificity in the folding process. The kinetic data also suggest that one residue, Asn43, may play a distinctive role in defining the BPTI folding mechanism. Replacement of this residue with either Gly or Ala appeared to stabilize the major transition states for folding and unfolding. In the native protein, the side-chain of Asn43 participates directly in the hydrogen bonding pattern of the central beta-sheet, and the kinetic behavior of the Asn43 variants suggests that the major energy barriers in folding and unfolding may be due in part to the steric constraints imposed by this structural element, together with those imposed by the chemical transition states for thiol-disulfide exchange.

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+).

Early events in the disulfide-coupled folding of BPTI

Recent studies of the refolding of reduced bovine pancreatic trypsin inhibitor (BPTI) have shown that a previously unidentified intermediate with a single disulfide is formed much more rapidly than any other one-disulfide species. This intermediate contains a disulfide that is present in the native protein (between Cys14 and 38), but it is thermodynamically less stable than the other two intermediates with single native disulfides. To characterize the role of the [14-38] intermediate and the factors that favor its formation, detailed kinetic and mutational analyses of the early disulfide-formation steps were carried out. The results of these studies indicate that the formation of [14-38] from the fully reduced protein is favored by both local electrostatic effects, which enhance the reactivities of the Cys14 and 38 thiols, and conformational tendencies that are diminished by the addition of urea and are enhanced at lower temperatures. At 25 degrees C and pH 7.3, approximately 35% of the reduced molecules were found to initially form the 14-38 disulfide, but the majority of these molecules then undergo intramolecular rearrangements to generate non-native disulfides, and subsequently the more stable intermediates with native disulfides. Amino acid replacements, other than those involving Cys residues, were generally found to have only small effects on either the rate of forming [14-38] or its thermodynamic stability, even though many of the same substitutions greatly destabilized the native protein and other disulfide-bonded intermediates. In addition, those replacements that did decrease the steady-state concentration of [14-38] did not adversely affect further folding and disulfide formation. These results suggest that the weak and transient interactions that are often detected in unfolded proteins and early folding intermediates may, in some cases, not persist or promote subsequent folding steps.

Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain

15Nitrogen relaxation experiments were used to characterize the backbone dynamics of two modified forms of bovine pancreatic trypsin inhibitor (BPTI). In one form, the disulfide between Cys14 and Cys38 in the wild-type protein was selectively reduced and methylated to generate an analog of the final intermediate in the disulfide-coupled folding pathway. The second form was generated by similarly modifying a mutant protein in which Tyr35 was replaced with Gly (Y35G). For both selectively reduced proteins, the overall conformation of native BPTI was retained, and the relaxation data for these proteins were compared with those obtained previously with the native wild-type and Y35G proteins. Removing the disulfide from either protein had only small effects on the observed longitudinal relaxation rates (R1) or heteronuclear cross relaxation rates (nuclear Overhauser effect), suggesting that the 14-38 disulfide has little influence on the fast (ps to ns) backbone dynamics of either protein. In the wild-type protein, the pattern of residues undergoing slower (micros to ms) internal motions, reflected in unusually large transverse relaxation rates (R2), was also largely unaffected by the removal of this disulfide. It thus appears that the large R2 rates previously observed in native wild-type protein are not a direct consequence of isomerization of the 14-38 disulfide. In contrast with the wild-type protein, reducing the disulfide in Y35G BPTI significantly decreased the number of backbone amides displaying large R2 rates. In addition, the frequencies of the backbone motions in the modified protein, estimated from R2 values measured at multiple refocusing delays, appear to span a wider range than those seen in native Y35G BPTI. Together, these observations suggest that the slow internal motions in Y35G BPTI are more independent in the absence of the 14-38 disulfide and that formation of this bond may lead to a substantial loss of conformational entropy. These effects may account for the previous observation that the Y35G substitution greatly destabilizes the disulfide. The results also demonstrate that the disulfide and the buried side-chain influence the dynamics of the folded protein in a highly cooperative fashion, with the effects of removing either being much greater in the absence of the other.

Mutational analysis of the BPTI folding pathway: II. Effects of aromatic-->leucine substitutions on folding kinetics and thermodynamics

The rates of the individual steps in the disulfide-coupled folding and unfolding of eight BPTI variants, each containing a single aromatic to leucine amino acid replacement, were measured. From this analysis, the contributions of the four phenylalanine and four tyrosine residues to the stabilities of the native protein and the disulfide-bonded folding intermediates were determined. While the substitutions were found to destabilize the native protein by 2 to 7 kcal/mol, they had significantly smaller effects on the intermediates that represent the earlier stages of folding, even when the site of the substitution was located within the ordered regions of the intermediates. These results suggest that stabilizing interactions contribute less to conformational stability in the context of a partially folded intermediate than in a fully folded native protein, perhaps because of decreased cooperativity among the individual interactions. The kinetic analysis also provides new information about the transition states associated with the slowest steps in folding and unfolding, supporting previous suggestions that these transition states are extensively unfolded. Although the substitutions caused large changes in the distribution of folding intermediates and in the rates of some steps in the folding pathway, the kinetically-preferred pathway for all of the variants involved intramolecular disulfide rearrangements, as observed previously for the wild-type protein. These results suggest that the predominance of the rearrangement mechanism reflects conformational constraints present relatively early in the folding pathway.

Thermodynamic genetics of the folding of the B1 immunoglobulin-binding domain from streptococcal protein G

A method has been developed to select proteins that are thermodynamically destabilized yet still folded and functional. The DNA encoding the B1 IgG-binding domain from Group G Streptococcus (Strp G) has been fused to gene III of bacteriophage M13. The resulting fusion protein is displayed on the surface of the phage thus enabling the phage to bind to IgG molecules. In addition, these phage exhibit a small plaque phenotype that is reversed by mutations that destabilize the Strp G domain. By selecting phage with large plaque morphology that retain their IgG-binding function, it is possible to identify mutants that are folded but destabilized compared with wild-type Strp G. Such mutants can be divided into three general categories: 1) those that disrupt packing of hydrophobic side chains in the protein interior; 2) those that destabilize secondary structure; and 3) those that alter specific hydrogen bonds involving amino acid side chains. A number of the mutants have been physically characterized by circular dichroism and nuclear magnetic resonance and have been shown to have structures similar to wild-type Strp G but stabilities that were decreased by 2-5 kcal/mol.

Thermodynamics of BPTI folding

A calorimetric study of the basic pancreatic trypsin inhibitor (BPTI) has been performed using the new generation of the adiabatic scanning microcalorimeters, operating in an extended temperature range of 5-130 degrees C. Precise measurements of the heat capacities of the native and unfolded states of BPTI show that the heat capacity change upon unfolding strongly depends on temperature; its value is maximal at about 50 degrees C and diminishes as the temperature is increased. The temperature dependencies of the enthalpy and entropy changes upon BPTI unfolding were found to be similar to those normally observed for other small globular proteins. The stability of BPTI has been correlated with its structure.

Crevice-forming mutants in the rigid core of bovine pancreatic trypsin inhibitor: crystal structures of F22A, Y23A, N43G, and F45A

Crystal structures of four mutants of bovine pancreatic trypsin inhibitor (F22A, Y23A, N43G, and F45A), engineered to alter their stability properties, have been determined. The mutated residues, which are highly conserved among Kunitz-type inhibitors, are located in the rigid core of the molecule. Replacement of the partially buried bulky residues of the wild-type protein with smaller residues resulted in crevices open to the exterior of the molecule. The overall three-dimensional structure of these mutants is very similar to that of the wild-type protein and only small rearrangements are observed among the atoms lining the crevices.

Probing the determinants of disulfide stability in native pancreatic trypsin inhibitor

The effects of single amino acid replacements on the stability of the 14-38 disulfide bond in the native form of bovine pancreatic trypsin inhibitor (BPTI) were measured. A total of 17 mutant proteins, with substitutions at one of 7 residues located 5-15 A from the disulfide in the native wild-type protein, were examined. The replacements were found to decrease the thermodynamic stability of the disulfide, as measured by exchange with thiol-disulfide reagents, by 0.6-5 kcal/mol, corresponding to a range of nearly 100 mV in redox potentials. The effects of the substitutions on disulfide stability were roughly correlated with the changes in side-chain volume, suggesting that optimal packing is a major factor in determining the stability of the disulfide in the wild-type protein. With only one exception, the substitutions also led to increases, as large as 50-fold, in the rates of disulfide reduction by dithiothreitol. The increased rates of reduction suggest that at least a fraction of the mutational destabilization of the disulfide is due to strain in the native protein that is relieved in the transition state for reduction. The stability of the disulfide in a peptide corresponding to the segments that are linked together by the 14-38 disulfide in native BPTI was found to be about 5 kcal/mol less than that of the disulfide in the intact wild-type protein. Together, the results with the mutant proteins and the peptide indicate that the stability of the disulfide in the native protein depends on both the local environment of the disulfide and on the ability of the rest of the protein to favor a conformation that promotes disulfide formation.

Small effects of amino acid replacements on the reduced and unfolded state of pancreatic trypsin inhibitor

The effects of amino acid replacements on the hydrodynamic volume of reduced and unfolded bovine pancreatic trypsin inhibitor (BPTI) have been examined by gel electrophoresis. The electrophoretic mobilities of the reduced forms of 46 BPTI variants were compared at room temperature in the absence of denaturants. The single substitutions examined include many different types of replacements at sites throughout the polypeptide, and, collectively, alter 22 of the 58 residues of the wild-type protein. The only substitutions found to alter the electrophoretic mobility of the reduced protein by more than approximately 3% are those that change the net charge of the protein. For nine mutants, the rates of disulfide formation in the reduced protein were also examined and found to be very similar to that of the wild-type protein. These results suggest that any structure that may be present in the reduced protein is either relatively insensitive to amino acid replacements or does not greatly influence the averaged properties of the polypeptide chain.