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.     

Pathway of oxidative folding of alpha-lactalbumin: a model for illustrating the diversity of disulfide folding pathways

The pathway of oxidative folding of alpha-lactalbumin (alpha LA) (four disulfide bonds) has been characterized by structural and kinetic analysis of the acid-trapped folding intermediates. In the absence of calcium, oxidative folding of alpha LA proceeds through highly heterogeneous species of one-, two-, three-, and four-disulfide (scrambled) intermediates to reach the native structure. In the presence of calcium, the folding intermediates of alpha LA comprise two predominant isomers (alpha LA-IIA and alpha LA-IIIA) adopting exclusively native disulfide bonds, including the two disulfide bonds (Cys(61)-Cys(77) and Cys(73)-Cys(91)) located within the beta-sheet calcium binding domain. alpha LA-IIA is a two-disulfide species consisting of Cys(61)-Cys(77) and Cys(73)-Cys(91) disulfide bonds. alpha LA-IIIA contains Cys(61)-Cys(77), Cys(73)-Cys(91), and Cys(28)-Cys(111) disulfide bonds. The underlying mechanism of the contrasting folding pathways of calcium-bound and calcium-depleted alpha LA is congruent with the cause of diversity of disulfide folding pathways observed among many well-characterized three-disulfide proteins, including bovine pancreatic trypsin inhibitor and hirudin. Our study also reveals novel aspects of the folding mechanism of alpha LA that have not been described previously.     

Balancing oxidative protein folding: The influences of reducing pathways on disulfide bond formation

Oxidative protein folding is confined to few compartments, including the endoplasmic reticulum, the mitochondrial intermembrane space and the bacterial periplasm. Conversely, in compartments in which proteins are translated such as the cytosol, the mitochondrial matrix and the chloroplast stroma proteins are kept reduced by the thioredoxin and glutaredoxin systems that functionally overlap. The highly reducing NADPH pool thereby serves as electron donor that enables glutathione reductase and thioredoxin reductase to keep glutathione pools and thioredoxins in their reduced redox state, respectively. Notably, also compartments containing oxidizing machineries are linked to these reducing pathways. Reducing pathways aid in proofreading of disulfide bond formation by isomerization or they provide reducing equivalents for the reduction of disulfides prior to degradation. In addition, they contribute to the thiol-dependent regulation of protein activities, and they help to counteract oxidative stress. The existence of oxidizing and reducing pathways in the same compartment poses a potential problem as the cell has to avoid futile cycles of oxidation and subsequent reduction reactions. Thus, compartments that contain oxidizing machineries have developed sophisticated ways to spatiotemporally balance and regulate oxidation and reduction. In this review, we discuss oxidizing and reducing pathways in the endoplasmic reticulum, the periplasm and the mitochondrial intermembrane space and highlight the role of glutathione especially in the endoplasmic reticulum and the intermembrane space. This article is part of a Special Issue entitled: Thiol-Based Redox Processes.     

Graph-representation of oxidative folding pathways

Background  

The process of oxidative folding combines the formation of native disulfide bond with conformational folding resulting in the native three-dimensional fold. Oxidative folding pathways can be described in terms of disulfide intermediate species (DIS) which can also be isolated and characterized. Each DIS corresponds to a family of folding states (conformations) that the given DIS can adopt in three dimensions.     

The oxidative folding of proteins by disulfide plus thiol does not correlate with redox potential

The rate and yield of oxidative renaturation of reduced RNase A has been studied as a function of [-S-S-]/[-SH]. The principal conclusion of these studies is that rates and yields of oxidative renaturation are strongly dependent on the low mol. wt disulfide/thiol ratio. The relationships are complex and do not parallel the redox potential of the system. The present results are consistent with earlier findings on other proteins, and lead us to believe that the above conclusion is general. Kinetic studies of oxidative renaturation should recognize and account for the dependence of reaction rate and extent on the disulfide/thiol ratio. This ratio can change substantially over the course of a reaction, either due to stoichiometric transfer of disulfide to protein, and/or adventitious air oxidation of thiols. Failure to account for changes in the disulfide/thiol ratio may compromise the interpretation of such experiments.     

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.     

Diversity in the disulfide folding pathways of cystine knot peptides

Summary The plant cyclotides are a fascinating family of circular proteins that contain a cyclic cystine knot motif (CCK). This unique family was discovered only recently but contains over 50 known sequences to date. Various biological activities are associated with these peptides including antimicrobial and insecticidal activity. The knotted topology and cyclic nature of the cyclotides poses interesting questions about the folding mechanisms and how the knotted arrangement of disulfide bonds is formed. Some studies have been performed on related inhibitor cystine knot (ICK) containing peptides, but little is known about the folding mechanisms of CCK molecules. We have examined the oxidative refolding and reductive unfolding of the prototypic member of the cyclotide family, kalata B1. Analysis of the rates of formation of the intermediates along the reductive unfolding pathway highlights the stability conferred by the cystine knot motif. Significant differences are observed between the folding of kalata B1 and an acyclic cystine knot protein, EETI-II, suggesting that the circular backbone has a significant influence in directing the folding pathway.     

Evidence for the underlying cause of diversity of the disulfide folding pathway

The pathways of oxidative folding of disulfide proteins exhibit a high degree of diversity, which is illustrated by the varied extent of (a) the heterogeneity of folding intermediates, (b) the predominance of intermediates containing native disulfide bonds, and (c) the level of accumulation of fully oxidized scrambled isomers as intermediates. BPTI and hirudin exemplify two extreme cases of such divergent folding pathways. We previously proposed that the underlying cause of this diversity is associated with the degree of stability of protein subdomains. Here we present compelling evidence that substantiates this hypothesis by studying the folding pathway of alphaLA-IIA. alphaLA-IIA is a partially folded intermediate of alpha-lactalbumin (alphaLA). It comprises a structured beta-sheet (calcium-binding) domain linked by two native disulfide bonds (Cys(61)-Cys(77) and Cys(73)-Cys(91)) and a disordered alpha-helical domain with four free cysteines (Cys(6), Cys(28), Cys(111), and Cys(120)). Purified alphaLA-IIA was allowed to refold without and with stabilization of its structured beta-sheet domain by calcium. In the absence of calcium, the folding pathway of alphaLA-IIA resembles that of hirudin, displaying a highly heterogeneous population of folding intermediates, including fully oxidized scrambled species. Upon stabilization of its beta-sheet domain by bound calcium, oxidative folding of alphaLA-IIA undergoes a pathway conspicuously similar to that of BPTI, exhibiting limited species of folding intermediates containing mostly native disulfide bonds.     

Coupling of conformational folding and disulfide-bond reactions in oxidative folding of proteins

The oxidative folding of proteins consists of conformational folding and disulfide-bond reactions. These two processes are coupled significantly in folding-coupled regeneration steps, in which a single chemical reaction (the "forward" reaction) converts a conformationally unstable precursor species into a conformationally stable, disulfide-protected successor species. Two limiting-case mechanisms for folding-coupled regeneration steps are described. In the folded-precursor mechanism, the precursor species is preferentially folded at the moment of the forward reaction. The (transient) native structure increases the effective concentrations of the reactive thiol and disulfide groups, thus favoring the forward reaction. By contrast, in the quasi-stochastic mechanism, the forward reaction occurs quasi-stochastically in an unfolded precursor; i.e., reactive groups encounter each other with a probability determined primarily by loop entropy, albeit modified by conformational biases in the unfolded state. The resulting successor species is initially unfolded, and its folding competes with backward chemical reactions to the unfolded precursors. The folded-precursor and quasi-stochastic mechanisms may be distinguished experimentally by the dependence of their kinetics on factors affecting the rates of thiol--disulfide exchange and conformational (un)folding. Experimental data and structural and biochemical arguments suggest that the quasi-stochastic mechanism is more plausible than the folded-precursor mechanism for most proteins.     

Heat capacity of folding of proteins corrected for disulfide cross-links

The heat capacities (DeltaC(p,f)) for the temperature-induced folding of proteins: barnase, lysozyme T4, papain, trypsin, ribonuclease T1, chymotrypsin, lysozyme and ribonuclease A have been calculated from the change in solvent accessible surface area between the native state and extended polypeptide chain. To visualize the effect of disulfide cross-links on molar heat capacity, loops of varying number of alanine residues and extended alanine chains with terminal cystein are modeled. The difference in DeltaC(p) values between the extended state and the loop conformation of proteins is linearly related to the number of residues in the loop. Corrections to the heat capacity of folding (DeltaC(p,f)) are applied for proteins with cross-links based on this observation. There is good correlation between corrected values of DeltaC(p,f) and experimental values.     

Contribution of individual disulfide bonds to the oxidative folding of ribonuclease A

The eight cysteine residues of ribonuclease A form four disulfide bonds in the native protein. We have analyzed the folding of three double RNase A mutants (C65A/C72A, C58A/C110A, and C26A/C84A, lacking the C65-C72, C58-C110, and C26-C84 disulfide bonds, respectively) and two single mutants (C110A and C26A), in which a single cysteine is replaced with an alanine and the paired cysteine is present in the reduced form. The folding of these mutants was carried out in the presence of oxidized and reduced glutathione, which constitute the main redox agents present within the ER. The use of mass spectrometry in the analysis of the folding processes allowed us (i) to follow the formation of intermediates and thus the pathway of folding of the RNase A mutants, (ii) to quantitate the intermediates that formed, and (iii) to compare the rates of formation of intermediates. By comparison of the folding kinetics of the mutants with that of wild-type RNase A, the contribution of each disulfide bond to the folding process has been evaluated. In particular, we have found that the folding of the C65A/C72A mutant occurs on the same time scale as that of the wild-type protein, thus suggesting that the removal of the C65-C72 disulfide bond has no effect on the kinetics of RNase A folding. Conversely, the C58A/C110A and C26A/C84A mutants fold much more slowly than the wild-type protein. The removal of the C58-C110 and C26-C84 disulfide bonds has a dramatic effect on the kinetics of RNase A folding. Results described in this paper provide specific information about conformational folding events in the regions involving the mutated cysteine residues, thus contributing to a better understanding of the complex mechanism of oxidative folding.     

Local and non-local native topologies reveal the underlying folding landscape of proteins

Due to Plaxco, Simons, Baker and others, it is now well known that the two-state single domain protein folding rate is fairly well predicted from knowledge of the topology of the native structure. Plaxco et al found that the folding rates of two-state proteins correlate with the average degree to which native contacts are 'local' within the chain sequence: fast-folders usually have mostly local structures. Here, we dissected the native topology further by focusing on non-local and local contacts using lower and upper bounds of allowable sequence separation in computing the average contact order. We analyzed non-local and local contacts of 82 two-state proteins whose experimental folding rates span over six orders of magnitude. We observed that both the number of non-local contacts and the average sequence separation of non-local contacts (non-local CO) are both negatively correlated with the folding rate, showing that the non-local contacts dominate the barrier-crossing process. Surprisingly, the local contact orders of the proteins also correlate with the folding rates. However, this correlation shows a strong positive trend indicating the role of a diffusive search in the denatured basin.     

Diselenides as universal oxidative folding catalysts of diverse proteins

Small-molecule diselenides show considerable potential as catalysts of oxidative protein folding. To explore their scope, diselenide-containing redox buffers were used to promote the folding of proteins that varied in properties such as size, overall tertiary structure, number of disulfide bonds, pI value, and difficulty of in vitro folding. Diselenides are able to catalyze the oxidative folding of all proteins tested, providing significant increases in both rate and yield relative to analogous disulfides. Compared to the disulfide-linked dimer of glutathione (the most commonly used oxidant for in vitro protein folding), selenoglutathione provided markedly improved efficiencies in the folding of biotechnologically important proteins such as hirudin, lysozyme, human epidermal growth factor and interferon α-2a. Selenoglutathione also enhances the renaturation of more challenging targets such as bovine serum albumin, whose native state contains 17 disulfide bonds, and the Fab fragment of an antibody. In the latter case, micromolar amounts of selenoglutathione are able to match the modest yield provided by a previously optimized redox buffer, which contains millimolar levels of glutathione. Taken together, the folding reactions of these diverse proteins exemplify the advantages and limitations of diselenide catalysts.     

Efficient traversal of beta-sheet protein folding pathways using ensemble models

Molecular dynamics (MD) simulations can now predict ms-timescale folding processes of small proteins; however, this presently requires hundreds of thousands of CPU hours and is primarily applicable to short peptides with few long-range interactions. Larger and slower-folding proteins, such as many with extended β-sheet structure, would require orders of magnitude more time and computing resources. Furthermore, when the objective is to determine only which folding events are necessary and limiting, atomistic detail MD simulations can prove unnecessary. Here, we introduce the program tFolder as an efficient method for modelling the folding process of large β-sheet proteins using sequence data alone. To do so, we extend existing ensemble β-sheet prediction techniques, which permitted only a fixed anti-parallel β-barrel shape, with a method that predicts arbitrary β-strand/β-strand orientations and strand-order permutations. By accounting for all partial and final structural states, we can then model the transition from random coil to native state as a Markov process, using a master equation to simulate population dynamics of folding over time. Thus, all putative folding pathways can be energetically scored, including which transitions present the greatest barriers. Since correct folding pathway prediction is likely determined by the accuracy of contact prediction, we demonstrate the accuracy of tFolder to be comparable with state-of-the-art methods designed specifically for the contact prediction problem alone. We validate our method for dynamics prediction by applying it to the folding pathway of the well-studied Protein G. With relatively very little computation time, tFolder is able to reveal critical features of the folding pathways which were only previously observed through time-consuming MD simulations and experimental studies. Such a result greatly expands the number of proteins whose folding pathways can be studied, while the algorithmic integration of ensemble prediction with Markovian dynamics can be applied to many other problems.     

Low energy pathways and non-native interactions: the influence of artificial disulfide bridges on the mechanism of folding

Four versions of a beta-sheet protein (CD2.d1) have been made, each with a single artificial disulfide bond inserted into hairpin structures. Folding kinetics of reduced and oxidized forms shows bridge position strongly influences its effect on the folding reaction. Bridging residues 58 and 62 does not affect the rapidly formed intermediate (I) or rate-limiting transition (t) state, whereas bridging 33 and 38, or 31 and 41, lowers the t-state energy, with the latter having the stronger influence. Bridging residues 79 and 90 stabilizes both I- and t-states. To assess additivity in the energetic effects of these bridges, four double-bridge variants have also been made. All show precise additivity of overall stability, with two showing additivity when ground states and the rate-limiting t-state are assessed, i.e. no measurable change in the folding mechanism occurs. However, combining 31-41 and 79-90 bridges produces a molecule that folds through a different pathway, with a much more stable intermediate than expected and a much higher t-state barrier. This is explained by the artificial introduction of stabilizing, non-native contacts in the I-state. More surprisingly, for another double-bridge version (58-62 and 79-90) both I- and t-states are less stable than expected, showing that conformational constraints introduced by the two bridges prevent formation of non-native contacts that would otherwise stabilize the I- and t-states, thereby lowering the energy of the folding landscape in the wild-type (unbridged) molecule. We conclude that the lowest energy path for folding has I- and t-state structures that are stabilized by non-native interactions.     

Catalysis of the oxidative folding of bovine pancreatic ribonuclease A by protein disulfide isomerase

The major oxidative folding pathways of bovine pancreatic ribonuclease A at pH 8.0 and 25 degrees C involve a pre-equilibrium steady state among ensembles of intermediates with zero, one, two, three and four disulfide bonds. The rate-determining steps are the reshuffling of the unstructured three-disulfide ensemble to two native-like three-disulfide species, des-[65-72] and des-[40-95], that convert to the native structure during oxidative formation of the fourth disulfide bond. Under the same regeneration conditions, with oxidized and reduced DTT, used previously for kinetic oxidative-folding studies of this protein, the addition of 4 microM protein disulfide isomerase (PDI) was found to lead to catalysis of each disulfide-formation step, including the rate-limiting rearrangement steps in which the native-like intermediates des-[65-72] and des-[40-95] are formed. The changes in the distribution of intermediates were also determined in the presence and absence of PDI at three different temperatures (with the DTT redox system) as well as at 25 degrees C (with the glutathione redox system). The results indicate that the acceleration of the formation of native protein by PDI, which we observed earlier, is due to PDI catalysis of each of the intermediate steps without changing the overall pathways or folding mechanism.     

Oxidative folding of Amaranthus alpha-amylase inhibitor: disulfide bond formation and conformational folding

Oxidative folding is the fusion of native disulfide bond formation with conformational folding. This complex process is guided by two types of interactions: first, covalent interactions between cysteine residues, which transform into native disulfide bridges, and second, non-covalent interactions giving rise to secondary and tertiary protein structure. The aim of this work is to understand both types of interactions in the oxidative folding of Amaranthus alpha-amylase inhibitor (AAI) by providing information both at the level of individual disulfide species and at the level of amino acid residue conformation. The cystine-knot disulfides of AAI protein are stabilized in an interdependent manner, and the oxidative folding is characterized by a high heterogeneity of one-, two-, and three-disulfide intermediates. The formation of the most abundant species, the main folding intermediate, is favored over other species even in the absence of non-covalent sequential preferences. Time-resolved NMR and photochemically induced dynamic nuclear polarization spectroscopies were used to follow the oxidative folding at the level of amino acid residue conformation. Because this is the first time that a complete oxidative folding process has been monitored with these two techniques, their results were compared with those obtained at the level of an individual disulfide species. The techniques proved to be valuable for the study of conformational developments and aromatic accessibility changes along oxidative folding pathways. A detailed picture of the oxidative folding of AAI provides a model study that combines different biochemical and biophysical techniques for a fuller understanding of a complex process.     

Unfolding and refolding pathways of a major kinetic trap in the oxidative folding of alpha-lactalbumin

Alpha-lactalbumin (alphaLA)-IIIA is a major kinetic intermediate present along the pathways of reductive unfolding and oxidative folding of bovine alpha-lactalbumin (alphaLA). It is a three-disulfide variant of native alphaLA lacking Cys(6)-Cys(120) at the alpha-helical domain. Stability and the unfolding/refolding mechanism of carboxymethylated alphaLA-IIIA have been investigated previously by stop-flow circular dichroism (CD) and fluorescence spectroscopy. A stable intermediate compatible with molten globule was shown to exist along the pathways of unfolding-refolding of alphaLA-IIIA [Ikeguchi et al. (1992) Biochemistry 31, 16695-12700; Horng et al. (2003) Proteins 52, 193-202]. We investigate here the unfolding-refolding pathways and conformational stability of alphaLA-IIIA using the method of disulfide scrambling with the following specific aims: (a) to isolate and characterize the observed stable molten globule, (b) to analyze the heterogeneity of folding-unfolding intermediates, (c) to elucidate the disulfide structure of extensively unfolded isomer of alphaLA-IIIA, and (d) to clarify the relative conformational stability between alphaLA-IIIA and alphaLA. Two scrambled isomers, designated as X-alphaLA-IIIA-c and X-alphaLA-IIIA-a (X stands for scrambled), were isolated under mild and strong denaturing conditions. Their disulfide structures, CD spectra, and manners of refolding to form the native alphaLA-IIIA were analyzed in this report. The results are consistent with the notion that X-alphaLA-IIIA-c and X-alphaLA-IIIA-a represent a partially unfolded and an extensively unfolded isomers of native alphaLA-IIIA, respectively. The unfolding-refolding pathways of alphaLA-IIIA are elaborated and compared with that of intact alphaLA. These results display new insight into one of the most extensively studied molecules in the field of protein folding and unfolding.     

Cadmium causes the oxidative modification of proteins in pea plants

In pea (Pisum sativum L.) leaves from plants grown in the presence of 50 µm CdCl₂ the oxidative production of carbonyl groups in proteins, the rate of protein degradation and the proteolytic activity were investigated. In leaf extracts the content of carbonyl groups measured by derivatization with 2,4‐dinitrophenylhydrazine (DNPH), was two‐fold higher in plants treated with Cd than in control plants. The identification of oxidized proteins was carried out by sodium dodecyl sulphate‐polyacrylamide gel electrophoresis of proteins derivatized with DNPH and immunochemical detection with an antibody against DNPH. The intensity of the reactive bands was higher in plants exposed to Cd than in controls. By using different antibodies some of the oxidized proteins were identified as Rubisco, glutathione reductase, manganese superoxide dismutase, and catalase. The incubation of leaf crude extracts with increasing H₂O₂ concentrations showed a progressive enhancement in carbonyl content and the pattern of oxidized proteins was similar to that found in Cd‐treated plants. Oxidized proteins were more efficiently degraded, and the proteolytic activity increased 20% due to the metal treatment. In peroxisomes purified from pea leaves a rise in the carbonyl content similar to that obtained in crude extracts from Cd‐treated plants was observed, but the functionality of the peroxisomal membrane was not apparently affected by Cd. Results obtained demonstrate the participation of both oxidative stress, probably mediated by H₂O₂, and proteolytic degradation in the mechanism of Cd toxicity in leaves of pea plants, and they appear to be involved in the Cd‐induced senescence previously reported in these plants.