Correlation between disulfide reduction and conformational unfolding in bovine pancreatic trypsin inhibitor

BPTI folding revisited: switching a disulfide into methylene thioacetal reveals a previously hidden path

The folding mechanism of the model protein bovine pancreatic trypsin inhibitor was revisited. By switching the solvent exposed disulfide bond with methylene thioacetal we uncovered a hidden pathway in its folding mechanism. In addition, this moiety enhanced protein stability while fully maintaining the protein structure and biological function.

Bovine pancreatic trypsin inhibitor (BPTI) is a 58-residue protein that is stabilized by three disulfide bonds at positions 5–55, 14–38 and 30–51. Widely studied for about 50 years, BPTI represents a folding model for many disulfide-rich proteins. In the study described below, we replaced the solvent exposed 14–38 disulfide bond with a methylene thioacetal bridge in an attempt to arrest the folding pathway of the protein at its two well-known intermediates, N′ and N*. The modified protein was expected to be unable to undergo the rate-determining step in the widely accepted BPTI folding mechanism: the opening of the 14–38 disulfide bond followed by rearrangements that leads to the native state, N. Surprisingly, instead of halting BPTI folding at N′ and N*, we uncovered a hidden pathway involving a direct reaction between the N* intermediate and the oxidizing reagent glutathione (GSSG) to form the disulfide-mixed intermediate N*–SG, which spontaneously folds into N. On the other hand, N′ was unable to fold into N. In addition, we found that the methylene thioacetal bridge enhances BPTI stability while fully maintaining its structure and biological function. These findings suggest a general strategy for enhancing protein stability without compromising on function or structure, suggesting potential applications for future therapeutic protein production.

Impact of an easily reducible disulfide bond on the oxidative folding rate of multi-disulfide-containing proteins

The burial of native disulfide bonds, formed within stable structure in the regeneration of multi-disulfide-containing proteins from their fully reduced states, is a key step in the folding process, as the burial greatly accelerates the oxidative folding rate of the protein by sequestering the native disulfide bonds from thiol-disulfide exchange reactions. Nevertheless, several proteins retain solvent-exposed disulfide bonds in their native structures. Here, we have examined the impact of an easily reducible native disulfide bond on the oxidative folding rate of a protein. Our studies reveal that the susceptibility of the (40-95) disulfide bond of Y92G bovine pancreatic ribonuclease A (RNase A) to reduction results in a reduced rate of oxidative regeneration, compared with wild-type RNase A. In the native state of RNase A, Tyr 92 lies atop its (40-95) disulfide bond, effectively shielding this bond from the reducing agent, thereby promoting protein oxidative regeneration. Our work sheds light on the unique contribution of a local structural element in promoting the oxidative folding of a multi-disulfide-containing protein.

Reductive unfolding of serum albumins uncovered by Raman spectroscopy

The reductive unfolding of bovine serum albumin (BSA) and human serum albumin (HSA) induced by dithiothreitol (DTT) is investigated using Raman spectroscopy. The resolution of the S-S Raman band into both protein and oxidized DTT contributions provides a reliable basis for directly monitoring the S-S bridge exchange reaction. The related changes in the protein secondary structure are identified by analyzing the protein amide I Raman band. For the reduction of one S-S bridge of BSA, a mean Gibbs free energy of -7 kJ mol(-1) is derived by studying the reaction equilibrium. The corresponding value for the HSA S-S bridge reduction is -2 kJ mol(-1). The reaction kinetics observed via the S-S or amide I Raman bands are identical giving a reaction rate constant of (1.02 +/- 0.11) M(-1) s(-1) for BSA. The contribution of the conformational Gibbs free energy to the overall Gibbs free energy of reaction is further estimated by combining experimental data with ab initio calculations.

Oxidative folding and N-terminal cyclization of onconase

Cyclization of the N-terminal glutamine residue to pyroglutamic acid in onconase, an anti-cancer chemotherapeutic agent, increases the activity and stability of the protein. Here, we examine the correlated effects of the folding/unfolding process and the formation of this N-terminal pyroglutamic acid. The results in this study indicate that cyclization of the N-terminal glutamine has no significant effect on the rate of either reductive unfolding or oxidative folding of the protein. Both the cyclized and uncyclized proteins seem to follow the same oxidative folding pathways; however, cyclization altered the relative flux of the protein in these two pathways by increasing the rate of formation of a kinetically trapped intermediate. Glutaminyl cyclase (QC) catalyzed the cyclization of the unfolded, reduced protein but had no effect on the disulfide-intact, uncyclized, folded protein. The structured intermediates of uncyclized onconase were also resistant to QC catalysis, consistent with their having a native-like fold. These observations suggest that, in vivo, cyclization takes place during the initial stages of oxidative folding, specifically, before the formation of structured intermediates. The competition between oxidative folding and QC-mediated cyclization suggests that QC-catalyzed cyclization of the N-terminal glutamine in onconase occurs in the endoplasmic reticulum, probably co-translationally.

A localized specific interaction alters the unfolding pathways of structural homologues

Reductive unfolding studies of proteins are designed to provide information about intramolecular interactions that govern the formation (and stabilization) of the native state and about folding/unfolding pathways. By mutating Tyr92 to G, A, or L in the model protein, bovine pancreatic ribonuclease A, and through analysis of temperature factors and molecular dynamics simulations of the crystal structures of these mutants, it is demonstrated that the markedly different reductive unfolding rates and pathways of ribonuclease A and its structural homologue onconase can be attributed to a single, localized, ring-stacking interaction between Tyr92 and Pro93 in the bovine variant. The fortuitous location of this specific stabilizing interaction in a disulfide-bond-containing loop region of ribonuclease A results in the localized modulation of protein dynamics that, in turn, enhances the susceptibility of the disulfide bond to reduction leading to an alteration in the reductive unfolding behavior of the homologues. These results have important implications for folding studies involving topological determinants to obtain folding/unfolding rates and pathways, for protein structure-function prediction through fold recognition, and for predicting proteolytic cleavage sites.

Dissimilarity in the Reductive Unfolding Pathways of Two Ribonuclease Homologues

Using DTT(red) as the reducing agent, the kinetics of the reductive unfolding of onconase, a frog ribonuclease, has been examined. An intermediate containing three disulfides, Ir, that is formed rapidly in the reductive pathway, is more resistant to further reduction than the parent molecule, indicating that the remaining disulfides in onconase are less accessible to DTT(red). Disulfide-bond mapping of Ir indicated that it is a single species lacking the (30-75) disulfide bond. The reductive unfolding pattern of onconase is consistent with an analysis of the exposed surface area of the cysteine sulfur atoms in the (30-75) disulfide bond, which reveals that these atoms are about four- and sevenfold, respectively, more exposed than those in the next two maximally exposed disulfides. By contrast, in the reductive unfolding of the homologue, RNase A, there are two intermediates, arising from the reduction of the (40-95) and (65-72) disulfide bonds, which takes place in parallel, and on a much longer time-scale, compared to the initial reduction of onconase; this behavior is consistent with the almost equally exposed surface areas of the cysteine sulfur atoms that form the (40-95) and (65-72) disulfide bonds in RNase A and the fourfold more exposed cysteine sulfur atoms of the (30-75) disulfide bond in onconase. Analysis and in silico mutation of the residues around the (40-95) disulfide bond in RNase A, which is analogous to the (30-75) disulfide bond of onconase, reveal that the side-chain of tyrosine 92 of RNase A, a highly conserved residue among mammalian pancreatic ribonucleases, lies atop the (40-95) disulfide bond, resulting in a shielding of the corresponding sulfur atoms from the solvent; such burial of the (30-75) sulfur atoms is absent from onconase, due to the replacement of Tyr92 by Arg73, which is situated away from the (30-75) disulfide bond and into the solvent, resulting in the large exposed surface-area of the cysteine sulfur atoms forming this bond. Removal of Tyr92 from RNase A resulted in the relatively rapid reduction of the mutant to form a single intermediate (des [40-95] Y92A), i.e. it resulted in an onconase-like reductive unfolding behavior. The reduction of the P93A mutant of RNase A proceeds through a single intermediate, the des [40-95] P93A species, as in onconase. Although mutation of Pro93 to Ala does not increase the exposed surface area of the (40-95) cysteine sulfur atoms, structural analysis of the mutant reveals that there is greater flexibility in the (40-95) disulfide bond compared to the (65-72) disulfide bond that may make the (40-95) disulfide bond much easier to expose, consistent with the reductive unfolding pathway and kinetics of P93A. Mutation of Tyr92 to Phe92 in RNase A has no effect on its reductive unfolding pathway, suggesting that the hydrogen bond between the hydroxyl group of Tyr92 and the carbonyl group of Lys37 has no impact on the local unfolding free energy required to expose the (40-95) disulfide bond. Thus, these data shed light on the differences between the reductive unfolding pathways of the two homologous proteins and provide a structural basis for the origin of this difference.

Simultaneous Characterization of the Reductive Unfolding Pathways of RNase B Isoforms by Top-Down Mass Spectrometry

A novel method for characterization of the simultaneous reductive unfolding pathways of five isoforms of bovine pancreatic ribonuclease B (RNase B) is demonstrated. The results indicate that each isoform unfolds reductively through two three-disulfide-containing structured intermediates before proceeding to the fully reduced form, as in the reductive unfolding pathways of the A variant lacking the carbohydrate chain. The rates of reduction of bovine pancreatic ribonuclease A (RNase A) and RNase B and the formation and consumption of their reductive intermediates are identical, indicating that the unfolding events necessary to expose disulfide bonds for reduction are not affected by the oligosaccharide. The method utilizes top-down mass spectrometry and a naturally occurring tag on the protein, viz. the carbohydrate moiety, to obtain unfolding information of an ensemble of protein isoforms and is a generally applicable methodological advance for conducting folding studies on mixtures of different proteins.

An electrospray ionization mass spectrometry investigation of 1-anilino-8-naphthalene-sulfonate (ANS) binding to proteins

The binding of 1-anilino-8-naphthalene-sulfonic acid (ANS) to various globular proteins at acidic pH has been investigated by electrospray ionization mass spectrometry (ESI-MS). Maximal ANS binding is observed in the pH range 3-5. As many as seven species of dye-bound complexes are detected for myoglobin. Similar studies were carried out with cytochrome c, carbonic anhydrase, triosephosphate isomerase, lysozyme, alpha-lactalbumin, and bovine pancreatic trypsin inhibitor (BPTI). Strong ANS binding was observed wherever molten globule states were postulated in solution. ANS binding is not observed for lysozyme and BPTI, which have tightly folded structures in the native form. Alpha-lactalbumin, which is structurally related to lysozyme but forms a molten globule at acidic pH, exhibited ANS binding. Reduction of disulfide bonds in these proteins leads to the detection of ANS binding even at neutral pH. Binding was suppressed at very low pH (<2.5), presumably due to neutralization of the charge on the sulfonate moiety. The distribution of the relative intensities of the protein bound ANS species varies with the charge state, suggesting heterogeneity of gas phase conformations. The binding strength of these complexes was qualitatively estimated by dissociating them using enhanced nozzle skimmer potentials. The skimmer voltages also affected the lower and higher charge states of these complexes in a different manner.

Separation of cytochrome c-dependent caspase activation from thiol-disulfide redox change in cells lacking mitochondrial DNA

Release of mitochondrial cytochrome c (cyt c) is an early and common event during apoptosis. Previous studies showed that the loss of cyt c triggered superoxide production by mitochondria and contributed to the oxidation of cellular thiol-disulfide redox state. In this study, we tested whether loss of the functional electron transport chain due to depleting mitochondrial DNA (mtDNA) would affect this redox-signaling mechanism during apoptosis. Results showed that cyt c release and caspase activation in response to staurosporine treatment were preserved in cells lacking mitochondrial DNA (rho0 cells). However, unlike the case with rho+ cells, in which a dramatic oxidation of intracellular glutathione (GSH) occurred after mitochondrial cyt c release, the thiol-disulfide redox state in apoptotic rho0 cells remained largely unchanged. Thus, mitochondrial signaling of caspase activation can be separated from the bioenergetic function, and mitochondrial respiratory chain is the principal source of ROS generation in staurosporine-induced apoptosis.

Solubility of recombinant human tissue factor pathway inhibitor

Study of recombinant human tissue factor pathway inhibitor (rhTFPI) solubility shows (1) an inverted bell-shaped pH-solubility profile with a broad solubility minimum between pH 5 and 10 such that the solubility minimum midpoint is 2-3 pH units away from its isoelectric point; (2) a negative temperature-solubility coefficient; (3) a strong dependence of solubility on the valence of electrolytes, with both multivalent cations and anions enhancing this effect; and (4) a significant increase of solubility in the presence of charged polymers. At pH 6-7, rhTFPI solubility-salt profiles display typical salting-in and salting-out biphasic effects. At a slightly lower pH (pH 5), a third phase in addition to the salting-in and salting-out phases was observed at low ionic strength conditions (5 to 50 mM) where rhTFPI solubility increased as salt concentration decreased. The salting-out constant for rhTFPI in NaCl is 1.04 M(-1) and is independent of the pH of the solution. Resolubilization of rhTFPI precipitates revealed that "insolubility precipitates" (seen during buffer exchanges) resulted from protein solute saturation and could be redissolved by "native" solvent conditions. On the other hand, "instability precipitates" (typically seen after exposure to elevated temperatures or extended storage periods) were caused by insoluble protein aggregate formation and required strongly denaturing conditions to redissolve.

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.

Protein folding stability can determine the efficiency of escape from endoplasmic reticulum quality control

A fraction of each secreted protein is retained and degraded by the endoplasmic reticulum (ER) quality control apparatus that restricts export to correctly folded proteins. The intrinsic biophysical attributes that determine efficiency of escape from this proofreading process have been examined by expressing mutants of bovine pancreatic trypsin inhibitor (BPTI) in yeast. Secretion efficiency is strongly correlated with thermodynamic stability for a series of six point mutations of BPTI. No correlation of secretion efficiency with either oxidative folding or refolding rates in vitro is found; both the rapidly folded Y35L BPTI mutant and the slowly unfolded G36D BPTI mutant exhibit low secretion efficiency. Elimination of cysteines 14 and 38 by mutagenesis does not increase secretion efficiency, indicating that intramolecular thiol/disulfide rearrangements are not primarily responsible for retention and degradation of destabilized BPTI variants. Mutant yeast strains with diminished ER-associated degradation do not secrete BPTI more efficiently, indicating that retention and degradation are separable processes. These data support a model for ER quality control, wherein protein folding is functionally reversible and the relative rates of folding, unfolding, vesicular export, and retention determine secretion efficiency.