Kinetic analysis of the folding and unfolding of a mutant form of bovine pancreatic trypsin inhibitor lacking the cysteine-14 and -38 thiols

The kinetics of the disulfide-coupled unfolding-refolding transition of a mutant form of bovine pancreatic trypsin inhibitor (BPTI) lacking Cys-14 and -38 were measured and compared to previous results for the wild-type protein and other modified forms. The altered cysteines, which were changed to serine in the mutant protein, are normally paired in a disulfide in the native protein but from disulfides with Cys-5 in two-disulfide kinetic intermediates during folding. Although the mutant protein could fold efficiently, the kinetics of both folding and unfolding were altered, reflecting the roles of these cysteines in the two-disulfide intermediates with "wrong" disulfides. The intramolecular rate constant for the formation of the second disulfide of the native mutant protein was more than 10(3)-fold lower than that for the formation of a second disulfide during the refolding of the wild-type protein. The observed rate of unfolding of the mutant protein was also lower than that of the wild-type protein, demonstrating that the altered cysteines are involved in the intramolecular rearrangements that are the rate-determining step in the unfolding of the wild-type protein. These results confirm the previous conclusion [Creighton, T.E. (1977) J. Mol. Biol. 113, 275-293] that the energetically preferred pathway for folding and unfolding of BPTI includes intramolecular rearrangements of intermediates in which Cys-14 and -38 are paired in disulfides not present in the native protein. The present results are also consistent with other, less detailed, studies with similar mutants lacking Cys-14 and -38 [Marks, C.B., Naderi, H., Kosen, P.A., Kuntz, I.D., & Anderson, S. (1987) Science (Washington, D.C.) 235, 1370-1371].     

Mutational analysis of the BPTI folding pathway: I. Effects of aromatic-->leucine substitutions on the distribution of folding intermediates.

The roles of aromatic residues in determining the folding pathway of bovine pancreatic trypsin inhibitor (BPTI) were analyzed mutationally by examining the distribution of disulfide-bonded intermediates that accumulated during the refolding of protein variants in which tyrosine or phenylalanine residues were individually replaced with leucine. The eight substitutions examined all caused significant changes in the intermediate distribution. In some cases, the major effect was to decrease the accumulation of intermediates containing two of the three disulfides found in the native protein, without affecting the distribution of earlier intermediates. Other substitutions, however, led to much more random distributions of the intermediates containing only one disulfide. These results indicate that the individual residues making up the hydrophobic core of the native protein make clearly distinguishable contributions to conformation and stability early in folding: The early distribution of intermediates does not appear to be determined by a general hydrophobic collapse. The effects of the substitutions were generally consistent with the structures of the major intermediates determined by NMR studies of analogs, confirming that the distribution of disulfide-bonded species is determined by stabilizing interactions within the ordered regions of the intermediates. The plasticity of the BPTI folding pathway implied by these results can be described using conformational funnels to illustrate the degree to which conformational entropy is lost at different stages in the folding of the wild-type and mutant proteins.     

Effect of phospholipids on conformational structure of bovine pancreatic trypsin inhibitor (BPTI) and its thermolabile mutants

The effects of negatively charged phosphatidylserine-prepared membranes (PS) and neutral phosphatidylcholine-prepared membranes (PC) on the structure of wild-type and mutant bovine pancreatic trypsin inhibitor (BPTI) at neutral pH were investigated. The presence of PC did not have any effect on the protein structure while PS induced a non-native structure in three mutant BPTI proteins. However, the negatively charged membrane did not have any effect on wild-type BPTI. The findings revealed that (i) elimination of some disulphide bonds results in dramatic change in protein structure, and, (ii) that this biochemical interaction is surface-driven and electrostatic interactions may play a very strong role in influencing the fore-stated changes in protein structure. Of further interest were the results obtained from investigating the possible role of PS fluidity and concentration in altering mutant. When the value of Gibbs free-energy change of unfolding (DeltaG(U)) was positive, various non-native structures were formed in a concentration-dependent manner. However, when the value of DeltaG(U) was negative, only two types of non-native structures were formed: one with high beta structure content at low PS fluidity state, and the other with a high alpha-helical content at high PS fluidity state. Our study reveals how particular combinations of phospholipid:protein interactions can induce a protein conformation transition from a native to a non-native one at neutral pH, especially when the native structure is predestabilized by amino acid substitutions. This revelation may open up opportunities to explore alternative ways in which phospholipids may play a role in protein mis-folding and the related pathologies.     

Alteration of the disulfide-coupled folding pathway of BPTI by circular permutation

The kinetics of disulfide-coupled folding and unfolding of four circularly permuted forms of bovine pancreatic trypsin inhibitor (BPTI) were studied and compared with previously published results for both wild-type BPTI and a cyclized form. Each of the permuted proteins was found to be less stable than either the wild-type or circular proteins, by 3-8 kcal/mole. These stability differences were used to estimate effective concentrations of the chain termini in the native proteins, which were 1 mM for the wild-type protein and 2.5 to 4000 M for the permuted forms. The circular permutations increased the rates of unfolding and caused a variety of effects on the kinetics of refolding. For two of the proteins, the rates of a direct disulfide-formation pathway were dramatically increased, making this process as fast or faster than the competing disulfide rearrangement mechanism that predominates in the folding of the wild-type protein. These two permutations break the covalent connectivity among the beta-strands of the native protein, and removal of these constraints appears to facilitate direct formation and reduction of nearby disulfides that are buried in the folded structure. The effects on folding kinetics and mechanism do not appear to be correlated with relative contact order, a measure of overall topological complexity. These observations are consistent with the results of other recent experimental and computational studies suggesting that circular permutation may generally influence folding mechanisms by favoring or disfavoring specific interactions that promote alternative pathways, rather than through effects on the overall topology of the native protein.     

Detection of local structures in reduced unfolded bovine pancreatic trypsin inhibitor.

The structure of BPTI and reduced BPTI in concentrated guanidinium HCl (GUHCl) in the presence of glycerol has been probed by measurements of dynamic nonradiative excitation energy transfer between probes attached to its amino groups. Interprobe distance distributions were obtained from analysis of donor fluorescence decay curves and used to characterize local structures in unordered states of the protein. Site specifically fluorescently labeled BPTI derivatives (1-n)BPTI (n = 15, 20, 41, 46) were used, each carrying a 2-methoxy-naphthyl-1-methylenyl group (MNA) at the N-terminal amino group of arg1 and 7-(dimethylamino)-coumarin-4-yl-acetyl residue (DA-coum) at one of its epsilon-NH2 groups of the lysine side chains. Analysis of donor fluorescence decay kinetics gave the interprobe distance distributions in the native and denatured states. The N-terminal-segment, residues 1-15, is in an extended conformation (with an average interprobe distance of 34 +/- 2 A) in the native state. Upon unfolding by reduction with DTT or beta-mercapto ethanol in 6 M GUHCl/glycerol mixture, the conformation of this segment relaxed to a state characterized by a reduced average interprobe distance and a larger width of the distances distribution. The average distance between residues 1 and 26, i.e., between the N-terminus and the turn of the twisted beta sheet element (residues 18-35), increased upon unfolding. At -30 degrees C in the above solvent, the distribution between these two sites was probably composed of two conformational subpopulations. About 45 +/- 20% of the molecules were characterized by a short interprobe distance (like the native state) representing a compact conformation, and 55 +/- 20% of the molecules showed large interprobe distances representing an expanded (unfolded) conformation. Thus local structures seem to exist in reduced denatured BPTI even under denaturing conditions in 6 M GUHCl/glycerol mixtures. Some of those structures are unstable in guanidinium isothiocyanate (GUSCN). The method introduced here is suitable for probing local structures and very long range interactions in unfolded proteins and for search for folding initiation sites (FISs) and early folding intermediates.     

Secretion incompetence of bovine pancreatic trypsin inhibitor expressed in Escherichia coli.

There is increasing evidence that protein folding and protein export are competing processes in prokaryotic cells. Virtually all secretion studies reported to date, however, have employed proteins that are relatively uncharacterized in terms of their folding behavior and three-dimensional structure. In contrast, the structural and biochemical parameters governing the folding of bovine pancreatic trypsin inhibitor (BPTI) and several of its mutants have been studied intensively. We therefore undertook a study of the secretion behavior in Escherichia coli of recombinant BPTI and its mutants. Wild-type BPTI and two well-characterized folding mutants (C14A, C38A)BPTI and (C30A, C51A)BPTI (missing the 14-38 and 30-51 disulfide bonds, respectively), were investigated by analyzing their expression fused to an E. coli signal sequence or to two synthetic IgG-binding domains of staphylococcal protein A. Both disulfide mutants are destabilized relative to wild-type BPTI and exhibit markedly altered folding kinetics: one (C14A, C38A) folds more slowly than wild-type BPTI and the other (C30A, C51A) unfolds more rapidly. Both mutants were observed to be exported 3-10 times more efficiently than the wild-type molecule. Moreover, the levels of unprocessed preprotein in the cytoplasm were severalfold higher for the wild-type fusion than for the fusion to the two folding mutants. Intracellular degradation of the BPTI moiety was also observed. These results are consistent with traffic of intracellular BPTI preproteins on at least three routes along the secretory pathway: (a) facile secretion of unfolded material, (b) intracellular folding leading to secretion blockage, and (c) degradation followed by export of truncated molecules. A novel feature of these findings is the implication that disulfide bonds can form in the bacterial cytoplasm and lead to secretion incompetence.     

On the role of some conserved and nonconserved amino acid residues in the transitional state and intermediate of apomyoglobin folding

The contributions of some amino acid residues in the A, B, G, and H helices to the formation of the folding nucleus and folding intermediate of apomyoglobin were estimated. The effects of point substitutions of Ala for hydrophobic amino acid residues on the structural stability of the native (N) protein and its folding intermediate (I), as well as on the folding/unfolding rates for four mutant apomyoglobin forms, were studied. The equilibrium and kinetic studies of the folding/unfolding rates of these mutant proteins in a wide range of urea concentrations demonstrated that their native state was considerably destabilized as compared with the wild-type protein, whereas the stability of the intermediate state changed moderately. It was shown that the amino acid residues in the A, G, and H helices contributed insignificantly to the stabilization of the apomyoglobin folding nucleus in the rate-limiting I ? N transition, taking place after the formation of the intermediate, whereas the residue of the B helix was of great importance in the formation of the folding nucleus in this transition.     

Rapid unfolding of a domain populates an aggregation-prone intermediate that can be recognized by GroEL

Some amino acid substitutions in phage P22 coat protein cause a temperature-sensitive folding (tsf) phenotype. In vivo, these tsf amino acid substitutions cause coat protein to aggregate and form intracellular inclusion bodies when folded at high temperatures, but at low temperatures the proteins fold properly. Here the effects of tsf amino acid substitutions on folding and unfolding kinetics and the stability of coat protein in vitro have been investigated to determine how the substitutions change the ability of coat protein to fold properly. The equilibrium unfolding transitions of the tsf variants were best fit to a three-state model, N if I if U, where all species concerned were monomeric, a result confirmed by velocity sedimentation analytical ultracentrifugation. The primary effect of the tsf amino acid substitutions on the equilibrium unfolding pathway was to decrease the stability (DeltaG) and the solvent accessibility (m-value) of the N if I transition. The kinetics of folding and unfolding of the tsf coat proteins were investigated using tryptophan fluorescence and circular dichroism (CD) at 222 nm. The tsf amino acid substitutions increased the rate of unfolding by 8-14-fold, with little effect on the rate of folding, when monitored by tryptophan fluorescence. In contrast, when folding or unfolding reactions were monitored by CD, the reactions were too fast to be observed. The tsf coat proteins are natural substrates for the molecular chaperones, GroEL/S. When native tsf coat protein monomers were incubated with GroEL, they bound efficiently, indicating that a folding intermediate was significantly populated even without denaturant. Thus, the tsf coat proteins aggregate in vivo because of an increased propensity to populate this unfolding intermediate.     

An approach for the assessment of the order of disruption of the elements of protein structure upon protein unfolding: A study of carbonic anhydrase B

An experimental approach named μ-analysis has been developed in order to elucidate the sequence of the loss of ordered structure by elements of a protein during the denaturation of the molecule. This approach is applicable for the analysis of proteins that fold (unfold) in a multistep process that involve the formation (destruction) of a range of intermediate states. The concept of the approach consists in systematic analysis of mutagenized forms of the protein with point substitutions of hydrophobic amino-acid residues and additional cysteine bridges. Importantly, the substitutions of the amino-acid residues must be localized to the same structural elements of the protein. Point substitutions of hydrophobic amino-acid residues mainly provide information on the structural elements of the protein that are disrupted at the final stages of protein denaturation. The addition of cysteine bridges to the surface of the protein molecule allows investigation of structural elements of the protein that are the first to unfold upon protein denaturation. Calorimetric studies of non-equilibrium melting of bovine carbonic anhydrase B yielded information on the rate constants of the unfolding of ten mutant forms of the protein. The analysis of the effects of mutations on the rates of different stages of protein unfolding allowed for elucidation of the order of disruption of structural elements of carbonic anhydrase B upon thermal denaturation.     

Basic residues in azurocidin/HBP contribute to both heparin binding and antimicrobial activity

Azurocidin/CAP37/HBP is an antimicrobial and chemotactic protein that is part of the innate defenses of human neutrophils. In addition, azurocidin is an inactive serine protease homolog with binding sites for diverse ligands including heparin and the bovine pancreatic trypsin inhibitor (BPTI). The structure of the protein reveals a highly cationic domain concentrated on one side of the molecule and responsible for its strong polarity. To investigate the role of this highly basic region, we produced three recombinant azurocidin mutant proteins that were altered in either one or both of two clusters of 4 basic residues located symmetrically on each side of a central cleft in the cationic domain. Two of the mutant proteins (Loop 3: R5Q, K6Q, R8Q, and R10Q; Loop 4: R61Q, R62Q, R63Q, and R65Q) exhibited little or no change in heparin and BPTI binding or in antimicrobial function. In contrast, the Loop 3/Loop 4 mutant (R5Q, K6Q, R8Q, R10Q, R61Q, R62Q, R63Q, and R65Q) in which all 8 basic residues were replaced showed greatly decreased ability to bind heparin and to kill Escherichia coli and Candida albicans. Thus, we report that the 8 basic residues that were altered in the Loop 3/Loop 4 mutant contribute to the ability of the wild-type azurocidin molecule to bind heparin and to kill E. coli and C. albicans. Because BPTI binding was comparable in wild-type and Loop 3/Loop 4 mutant protein, we conclude that the same 8 basic residues are not involved in the binding of BPTI to azurocidin, supporting the notion that the binding site for BPTI is distinct from the site involved in heparin binding and antimicrobial activity. Finally, we show that removal of all 4 positively charged amino acids in the 20-44 azurocidin sequence (DMC1: R23Q,H24S,H32S,R34Q), a region previously thought to contain an antimicrobial domain, does not affect the activity of the protein against E. coli, Streptococcus faecalis, and C. albicans.     

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.     

NMR studies of structure, hydrogen exchange, and main-chain dynamics in a disrupted-core mutant of thioredoxin.

Core-packing mutants of proteins often approach molten globule states, and hence may have attributes of folding intermediates. We have studied a core-packing mutant of thioredoxin, L78K, in which a leucine residue is substituted by lysine, using 15N heteronuclear two- and three-dimensional NMR. Chemical shift differences between the mutant and wild-type main-chain resonances reveal that structural changes caused by the mutation are localized within 12 A of the altered side chain. The majority of resonances are unchanged, as are many 1H-1H NOEs indicative of the main-chain fold, suggesting that the structure of L78K is largely similar to wild type. Hydrogen exchange studies reveal that residues comprising the central beta-sheet of both mutant and wild-type proteins constitute a local unfolding unit, but with the unfolding/folding equilibrium approximately 12 times larger in L78K. The dynamics of main-chain NH bonds in L78K were studied by 15N spin relaxation and compared with a previous study of wild type. Order parameters for angular motion of NH bonds in the mutant are on average lower than in wild type, suggesting greater spatial freedom on a rapid time scale, but may also be related to different rotational correlation times in the two proteins. There is also evidence of greater conformational exchange in the mutant. Differences between mutant and wild type in hydrogen exchange and main-chain dynamics are not confined to the vicinity of the mutation. We infer that mispacking of the protein core in one location affects local dynamics and stability throughout.     

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.     

Two-dimensional 1H nuclear magnetic resonance study of the (5-55) single-disulphide folding intermediate of bovine pancreatic trypsin inhibitor

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.     

Nonlocal interactions stabilize compact folding intermediates in reduced unfolded bovine pancreatic trypsin inhibitor.

To further our understanding of the protein folding process, it is desirable to examine the structural intermediates (equilibrium and kinetic) that are populated between the statistical coil state and the folded molecule. X-ray crystallography and NMR structural studies are unable to determine long-range distances in proteins under denaturing solution conditions. Nonradiative (F?rster) energy transfer, however, has been shown to be a spectroscopic ruler for the measurement of distance distributions and diffusion between selected sites in proteins under a range of different solution conditions. The distributions of distances between a donor probe at the N-terminal residue and an acceptor attached to one of the four lysine residues (15, 26, 41, 46) of reduced and unfolded (in 6 M guanidine hydrochloride and 20 mM dithiothreitol) bovine pancreatic trypsin inhibitor (BPTI) were measured as a function of temperature. Even in strong denaturant and reducing agent, BPTI does not exist as a statistical coil polypeptide. It appears that nonlocal (long-range) interactions are already beginning to "fold" the protein toward a more compact, native conformation. As the temperature is increased under these conditions, hydrophobic interactions lead to an even more compact structure consistent with the predictions of phase diagrams for globular proteins.     

Phi-values for BPTI folding intermediates and implications for transition state analysis

Amino acid replacements were used to probe the roles of 14 sites in two well-characterized intermediates in the folding pathway of bovine pancreatic trypsin inhibitor (BPTI). One of these intermediates contains one of the three disulfides found in the native protein (30--51). NMR studies have shown that approximately two-thirds of this polypeptide has a native-like conformation. The other intermediate contains two native disulfides (30--51 and 5--55) and has a fully folded conformation. The phi-values for a majority of residues were <1, indicating that the native protein was significantly more destabilized than either intermediate even when the altered residue was located in a well-ordered region of the intermediate. These observations suggest that folding intermediates and transition states may generally be more structured than indicated by phi-values alone.     

Structural stability of disulfide mutants of basic pancreatic trypsin inhibitor: A molecular dynamics study

The structure and folding of basic pancreatic trypsin inhibitor (BPTI) has been studied extensively by experimental means. We report a computer simulation study of the structural stability of various disulfide mutants of BPTI, involving eight 250-psec molecular dynamics simulations of the proteins in water, with and without a phosphate counterion. The presence of the latter alters the relative stability of the single disulfide species [5–55] and [30–51]. This conclusion can explain results of mutational studies and the conservation of residues in homologues of BPTI, and suggests a possible role of ions in stabilizing one intermediate over another in unfolding or folding processes. © 1996 Wiley-Liss, Inc.     

The underlying mechanism for the diversity of disulfide folding pathways

The disulfide folding pathway of bovine pancreatic trypsin inhibitor (BPTI) is characterized by the predominance of folding intermediates with native-like structures. Our laboratory has recently analyzed the folding pathway(s) of four 3-disulfide-containing proteins, including hirudin, potato carboxypeptidase inhibitor, epidermal growth factor, and tick anticoagulant peptide. Their folding mechanism(s) differ from that of BPTI by 1) a higher degree of heterogeneity of 1- and 2-disulfide intermediates and 2) the presence of 3-disulfide scrambled isomers as folding intermediates. To search for the underlying causes of these diversities, we conducted kinetic analyses of the reductive unfolding of these five proteins. The experiment of reductive unfolding was designed to evaluate the relative stability and interdependence of disulfide bonds in the native protein. It is demonstrated here that among these five proteins, there exists a striking correlation between the mechanism(s) of reductive unfolding and that of oxidative folding. Those proteins with their native disulfide bonds reduced in a collective and simultaneous manner exhibit both a high degree of heterogeneity of folding intermediates and the accumulation of scrambled isomers along the folding pathway. A sequential reduction of the native disulfide bonds is associated with the presence of predominant intermediates with native- like structures.