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

A growth type pathway for improving the folding of BPTI

The in vitro oxidative folding of the protein bovine pancreatic trypsin inhibitor (BPTI) with oxidized dithiothreitol or glutathione has served as a paradigm for protein folding but could take weeks at physiological pH because of the need to escape from kinetic traps via a rearrangement type pathway. The two major kinetic traps are called N' and N* and contain two of the three native disulfide bonds, which occur between residues 5 and 55, 30 and 51, and 14 and 38. N' is missing the disulfide bond between residues 5 and 55 while N* is missing the disulfide bond between residues 30 and 51. By determining rate constant for the reactions of the kinetic traps N* and N' and their mixed disulfides with glutathione and glutathione disulfide, many for the first time, we demonstrate that growth type pathways are feasible and could even be more efficient than rearrangement type pathways. Thus, formally unproductive pathways became productive. Interestingly, under physiological redox conditions both rearrangement and growth type pathways are important highlighting the redundancy of oxidative protein folding. With the new set of rate constants, modeling indicated that in vitro oxidative protein folding of BPTI via a growth type pathway using an oxidation, reduction and oxidation cycle would significantly improve protein folding efficiency, albeit under non-physiological redox conditions. With these changing conditions 91 ± 2% of native BPTI was achieved in 12 h compared to 83% native protein in 24 h using our previous best conditions of 5 mM GSSG and 5 mM GSH. Therefore, changing redox conditions via an oxidation, reduction and oxidation cycle may become an additional methodology for enhancing in vitro protein folding in aqueous solution.

Oxidative folding of conotoxins sharing an identical disulfide bridging framework

Conotoxins are short, disulfide-rich peptide neurotoxins produced in the venom of predatory marine cone snails. It is generally accepted that an estimated 100,000 unique conotoxins fall into only a handful of structural groups, based on their disulfide bridging frameworks. This unique molecular diversity poses a protein folding problem of relationships between hypervariability of amino acid sequences and mechanism(s) of oxidative folding. In this study, we present a comparative analysis of the folding properties of four conotoxins sharing an identical pattern of cysteine residues forming three disulfide bridges, but otherwise differing significantly in their primary amino acid sequence. Oxidative folding properties of M-superfamily conotoxins GIIIA, PIIIA, SmIIIA and RIIIK varied with respect to kinetics and thermodynamics. Based on rates for establishing the steady-state distribution of the folding species, two distinct folding mechanisms could be distinguished: first, rapid-collapse folding characterized by very fast, but low-yield accumulation of the correctly folded form; and second, slow-rearrangement folding resulting in higher accumulation of the properly folded form via the reshuffling of disulfide bonds within folding intermediates. Effects of changing the folding conditions indicated that the rapid-collapse and the slow-rearrangement mechanisms were mainly determined by either repulsive electrostatic or productive noncovalent interactions, respectively. The differences in folding kinetics for these two mechanisms were minimized in the presence of protein disulfide isomerase. Taken together, folding properties of conotoxins from the M-superfamily presented in this work and from the O-superfamily published previously suggest that conotoxin sequence diversity is also reflected in their folding properties, and that sequence information rather than a cysteine pattern determines the in vitro folding mechanisms of conotoxins.

Divergent folding pathways of two homologous proteins, BPTI and tick anticoagulant peptide: Compartmentalization of folding intermediates and identification of kinetic traps

Bovine pancreatic trypsin inhibitor (BPTI) and tick anticoagulant peptide (TAP) are two structurally homologous proteins, which have been shown to exhibit distinct mechanisms of oxidative folding. We demonstrate here differences of their folding properties using the technique of disulfide scrambling. Two extensively unfolded homologous isomers (beads-form) of BPTI (Cys5-Cys14, Cys30-Cys38, Cys51-Cys55) and TAP (Cys5-Cys15, Cys33-Cys39, Cys55-Cys59) were allowed to refold in parallel via disulfide shuffling of 13 potential isomers to form the native structure. Folding intermediates were trapped by acid quenching and analyzed by HPLC. The results reveal the following diversities: (a) there are two predominant folding intermediates of BPTI and seven well-populated folding intermediates of TAP. One of the two predominant BPTI intermediates (Cys5-Cys14, Cys30-Cys51, Cys38-Cys55) contains a native disulfide Cys30-Cys51 and constitutes about 34% of the total BPTI folding intermediates. In contrast, the TAP counterpart (Cys5-Cys15, Cys33-Cys55, Cys39-Cys59) represents only 5% of the total TAP intermediates; (b) three isomers of TAP sharing a stable non-native disulfide bond Cys15-Cys33 are shown to act as kinetic traps of TAP folding. Their counterparts are conspicuously absent in the BPTI folding; and (c) most significantly, folding intermediates of BPTI are found to be energetically compartmentalized, whereas most folding intermediates of TAP are inter-convertible and exist in dynamic equilibrium. Our studies further demonstrate optimal concentrations of denaturant required for destabilization of kinetic traps and acceleration of TAP folding.

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.

Role of kinetic intermediates in the folding of leech carboxypeptidase inhibitor

The oxidative folding and reductive unfolding pathways of leech carboxypeptidase inhibitor (LCI; four disulfides) have been characterized in this work by structural and kinetic analysis of the acid-trapped folding intermediates. The oxidative folding of reduced and denatured LCI proceeds rapidly through a sequential flow of 1-, 2-, 3-, and 4-disulfide (scrambled) species to reach the native form. Folding intermediates of LCI comprise two predominant 3-disulfide species (designated as III-A and III-B) and a heterogeneous population of scrambled isomers that consecutively accumulate along the folding reaction. Our study reveals that forms III-A and III-B exclusively contain native disulfide bonds and correspond to stable and partially structured species that interconvert, reaching an equilibrium prior to the formation of the scrambled isomers. Given that these intermediates act as kinetic traps during the oxidative folding, their accumulation is prevented when they are destabilized, thus leading to a significant acceleration of the folding kinetics. III-A and III-B forms appear to have both native disulfides bonds and free thiols similarly protected from the solvent; major structural rearrangements through the formation of scrambled isomers are required to render native LCI. The reductive unfolding pathway of LCI undergoes an apparent all-or-none mechanism, although low amounts of intermediates III-A and III-B can be detected, suggesting differences in protection against reduction among the disulfide bonds. The characterization of III-A and III-B forms shows that the former intermediate structurally and functionally resembles native LCI, whereas the III-B form bears more resemblance to scrambled isomers.

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.

Mechanism of insulin chain combination. Asymmetric roles of A-chain alpha-helices in disulfide pairing

The A and B chains of insulin combine to form native disulfide bridges without detectable isomers. The fidelity of chain combination thus recapitulates the folding of proinsulin, a precursor protein in which the two chains are tethered by a disordered connecting peptide. We have recently shown that chain combination is blocked by seemingly conservative substitutions in the C-terminal alpha-helix of the A chain. Such analogs, once formed, nevertheless retain high biological activity. By contrast, we demonstrate here that chain combination is robust to non-conservative substitutions in the N-terminal alpha-helix. Introduction of multiple glycine substitutions into the N-terminal segment of the A chain (residues A1-A5) yields analogs that are less stable than native insulin and essentially without biological activity. (1)H NMR studies of a representative analog lacking invariant side chains Ile(A2) and Val(A3) (A chain sequence GGGEQCCTSICSLYQLENYCN; substitutions are italicized and cysteines are underlined) demonstrate local unfolding of the A1-A5 segment in an otherwise native-like structure. That this and related partial folds retain efficient disulfide pairing suggests that the native N-terminal alpha-helix does not participate in the transition state of the reaction. Implications for the hierarchical folding mechanisms of proinsulin and insulin-like growth factors are discussed.

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.

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.

Single amino acid substitutions globally suppress the folding defects of temperature-sensitive folding mutants of phage P22 coat protein

The amino acid sequence of a polypeptide defines both the folding pathway and the final three-dimensional structure of a protein. Eighteen amino acid substitutions have been identified in bacteriophage P22 coat protein that are defective in folding and cause their folding intermediates to be substrates for GroEL and GroES. These temperature-sensitive folding (tsf) substitutions identify amino acids that are critical for directing the folding of coat protein. Additional amino acid residues that are critical to the folding process of P22 coat protein were identified by isolating second site suppressors of the tsf coat proteins. Suppressor substitutions isolated from the phage carrying the tsf coat protein substitutions included global suppressors, which are substitutions capable of alleviating the folding defects of numerous tsf coat protein mutants. In addition, potential global and site-specific suppressors were isolated, as well as a group of same site amino acid substitutions that had a less severe phenotype than the tsf parent. The global suppressors were located at positions 163, 166, and 170 in the coat protein sequence and were 8-190 amino acid residues away from the tsf parent. Although the folding of coat proteins with tsf amino acid substitutions was improved by the global suppressor substitutions, GroEL remained necessary for folding. Therefore, we believe that the global suppressor sites identify a region that is critical to the folding of coat protein.

Using loop length variants to dissect the folding pathway of a four-helix-bundle protein

Rop is a four-helix-bundle protein formed by the association of two helix-loop-helix monomers. The short helix-connecting loop was replaced with a series of polyglycine linkers of increasing length. These mutant proteins all appear to fold via the same general mechanism as that of the wild-type protein, even at the longest loop lengths. Replacement of the wild-type two-residue loop (Asp-Ala) with a (Gly-Gly) linker accelerates both unfolding and refolding rates. These changes in folding and unfolding kinetics likely indicate an alteration in the energy of the transition state. As the length of the glycine linker is further increased, the unfolding rate increases while the refolding rates decrease. The influence of loop length is not limited to these rates, but also impacts upon the stability of the folding intermediate. These dependences underscore the importance of loop closure and help refine the model for Rop's folding, implicating a dimeric intermediate involving hairpin formation. These observations show that loop alteration may be useful as a general technique for dissecting protein folding pathways.

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.

Simulations of the structural and dynamical properties of denatured proteins: the "molten coil" state of bovine pancreatic trypsin inhibitor

The dynamic nature of denatured, unfolded proteins makes it difficult to characterize their structures experimentally. To complement experiment and to obtain more detailed information about the structure and dynamic behavior of the denatured state, we have performed eleven 2.5 ns molecular dynamics simulations of reduced bovine pancreatic trypsin inhibitor (BPTI) at high temperature in water and a control simulation at 298 K, for a total of 30 ns of simulation time. In a neutral pH environment (acidic residues ionized), the unfolded protein structures were compact with an average radius of gyration 9% greater than the native state. The compact conformations resulted from the transient formation of non-native hydrophobic clusters, turns and salt bridges. However, when the acidic residues were protonated, the protein periodically expanded to a radius of gyration of 18 to 20 A. The early steps in unfolding were similar in the different simulations until passing through the major transition state of unfolding. Afterwards, unfolding proceeded through one of two general pathways with respect to secondary structure: loss of the C-terminal helix followed by loss of beta-structure or the opposite. To determine whether the protein preferentially sampled particular conformational substates in the denatured state, pairwise Calpha root-mean-square deviations were measured between all structures, but similar structures were found between only two trajectories. Yet, similar composite properties (secondary structure content, side-chain and water contacts, solvent accessible surface area, etc.) were observed for the structures that unfolded through different pathways. Somewhat surprisingly, the unfolded structures are in agreement with both past experiments suggesting that reduced BPTI is a random coil and more recent experiments providing evidence for non-random structure, demonstrating how ensembles of fluctuating structures can give rise to experimental observables that are seemingly at odds.

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.