Diverse pathways of oxidative folding of disulfide proteins: underlying causes and folding models

The pathway of oxidative folding of disulfide proteins exhibits a high degree of diversity, which is manifested mainly by distinct structural heterogeneity and diverse rearrangement pathways of folding intermediates. During the past two decades, the scope of this diversity has widened through studies of more than 30 disulfide-rich proteins by various laboratories. A more comprehensive landscape of the mechanism of protein oxidative folding has emerged. This review will cover three themes. (1) Elaboration of the scope of diversity of disulfide folding pathways, including the two opposite extreme models, represented by bovine pancreatic trypsin inhibitor (BPTI) and hirudin. (2) Demonstration of experimental evidence accounting for the underlying mechanism of the folding diversity. (3) Discussion of the convergence between the extreme models of oxidative folding and models of conventional conformational folding (framework model, hydrophobic collapse model).     

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.     

Predicting protein folding pathways

A structured folding pathway, which is a time ordered sequence of folding events, plays an important role in the protein folding process and hence, in the conformational search. Pathway prediction, thus gives more insight into the folding process and is a valuable guiding tool to search the conformation space. In this paper, we propose a novel 'unfolding' approach to predict the folding pathway. We apply graph-based methods on a weighted secondary structure graph of a protein to predict the sequence of unfolding events. When viewed in reverse this yields the folding pathway. We demonstrate the success of our approach on several proteins whose pathway is partially known.     

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.     

An analysis of protein folding pathways

We have developed a model of the protein folding process based on three primary assumptions: that burying of hydrophobic area is the dominant contribution to the relative free energy of a conformation, that a record of the folding process is largely preserved in the final structure, and that the denatured state is a random coil. Detailed folding pathways are identified for 19 protein structures. The picture of the folding process that emerges from this analysis is one of nucleation by regions of 8-16 residues. Nucleation sites then lead to larger structures by two mechanisms: propagation and diffusion/collision. A Monte Carlo simulation is used to follow the folding pathway when propagation is the dominant mechanism. Because detailed pathways are derived for each protein, the models are susceptible to experimental verification.     

Prediction of protein folding pathways.

Recent 1H nuclear magnetic resonance (n.m.r.) hydrogen exchange experiments on five different proteins have delineated the secondary structures formed in trapped, partially folded intermediates. The early forming structural elements are identifiable through a technique described in this work to predict folding pathways. The method assumes that the sequential selection of structural fragments such as alpha-helices and beta-strands involved in the folding process is founded upon the maximal burial of solvent accessible surface from both the formation of internal structure and substructure association. The substructural elements were defined objectively by major changes in main-chain direction. The predicted folding pathways are in complete correspondence with the n.m.r. results in that the formed structural fragments found in the folding intermediates are those predicted earliest in the pathways. The technique was also applied to proteins of known tertiary structure and with fold similar to one of the five proteins examined by 1H n.m.r. The pathways for these structures also showed general consistency with the n.m.r. observations, suggesting conservation of a secondary structural framework or molten globule about which folding nucleates and proceeds.     

Detergent-assisted oxidative folding of delta-conotoxins.

Conotoxins comprise a diverse group of disulfide-rich peptides found in venoms of predatory Conus species. The native conformation of these peptides is marginally stable in comparison with alternative conformations, often resulting in low folding yields. The oxidative folding of hydrophobic delta-conotoxins was found to produce less than 1% of the native peptide [Bulaj, G. et al. (2001) Biochemistry 40, 13201]. In order to identify factors that might improve folding yields, we screened a number of additives including water-soluble polymers, detergents and osmolytes for their ability to increase steady-state accumulation of the native delta-conotoxin PVIA. The presence of a non-ionic detergent Tween and low temperature appeared to be the most effective factors in improving the oxidative folding. The detergent was also effective in promoting folding of other hydrophobic delta-conotoxins. Based on our findings, we discuss a possible mechanism for detergent-assisted folding and the general applicability of this mechanism to facilitating the proper folding of hydrophobic, cysteine-rich peptides.     

Small-molecule catalysts of oxidative protein folding

Oxidative protein folding occurs both in vivo and in vitro and involves the formation and rearrangement of protein disulfide bonds (SS bonds). In vivo these reactions are catalyzed by enzymes, including the eukaryotic enzyme protein disulfide isomerase (PDI). Using the physical properties of PDI as a guide, several small-molecule catalysts of oxidative protein folding have been designed, synthesized, and tested. These small molecules can improve the folding rate of the model substrate ribonuclease A by a factor of over 10 and improve the yield by up to a factor of 3 over traditional conditions. The molecules have also been demonstrated to significantly improve the in vivo folding of proteins as well.     

The alpha- and beta-tubulin folding pathways

The alpha-beta tubulin heterodimer is the subunit from which microtubules are assembled. The pathway leading to correctly folded alpha- and beta-tubulins is unusually complex: it involves cycles of ATP-dependent interaction of newly synthesized tubulin subunits with cytosolic chaperonin, resulting in the production of quasi-native folding intermediates, which must then be acted upon by additional protein cofactors. These cofactors form a supercomplex containing both alpha- and beta-tubulin polypeptides, from which native heterodimer is released in a GTP-dependent reaction. Here, we discuss the current state of our understanding of the function of cytosolic chaperonin and cofactors in tubulin folding.     

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.     

Role of ascorbate in oxidative protein folding

Both in prokaryotic and eukaryotic cells, disulfide bond formation (oxidation and isomerization steps) are catalyzed exclusively in extracytoplasmic compartments. In eukaryotes, protein folding and disulfide bond formation are coupled processes that occur both co- and posttranslationally in the endoplasmic reticulum (ER), which is the main site of the synthesis and posttranslational modification of secretory and membrane proteins. The formation of a disulfide bond from the thiol groups of two cysteine residues requires the removal of two electrons, consequently, these bonds cannot form spontaneously; an oxidant is needed to accept the electrons. In aerobic conditions the ultimate electron acceptor is usually oxygen; however, oxygen itself is not effective in protein thiol oxidation. Therefore, a small molecular weight membrane permeable compound should be supposed for the transfer of electrons from the ER lumen. The aim of the present study was the investigation of the role of ascorbate/dehydroascorbate redox couple in oxidative folding of proteins. We demonstrated that ascorbate addition or its in situ synthesis from gulonolactone results in protein thiol (and/or glutathione; GSH) oxidation in rat liver microsomes. Since microsomal membrane is hardly permeable to ascorbate, the existence of a transport metabolon was hypothesized. Three components of the system have been described and partially characterized: (i) A microsomal metalloenzyme is responsible for ascorbate oxidation on the outer surface of the ER. Ascorbate oxidation results in ascorbate free radical and dehydroascorbate production. (ii) Facilitated diffusion of dehydroascorbate is present in microsomal vesicles. The transport is presumably mediated by a GLUT-type transporter. On the contrary, the previously hypothesized glutathione disulfide (GSSG) transport is practically absent, while GSH is transported with a moderate velocity. (iii) Protein disulfide isomerase catalyzes the reduction of dehydroascorbate in the ER lumen. Both GSH and protein thiols can be electron donors in the process. Intraluminal dehydroascorbate reduction and the consequent ascorbate accumulation strictly correlate with protein disulfide isomerase activity and protein thiol concentration. The concerted action of the three components of the system results in the intraluminal accumulation of ascorbate, protein disulfide and GSSG. In fact, intraluminal ascorbate and GSSG accumulation could be observed upon dehydroascorbate and GSH uptake. In conclusion, ascorbate is able to promote protein disulfide formation in an in vitro system. Further work is needed to justify its role in intact cellular and in vivo systems, as well as to explore the participation of other antioxidants (e.g. tocopherol, ubiquinone, and vitamin K) in the electron transfer chain responsible for oxidative protein folding in the ER.     

Computing folding pathways between RNA secondary structures

Given an RNA sequence and two designated secondary structures A, B, we describe a new algorithm that computes a nearly optimal folding pathway from A to B. The algorithm, RNAtabupath, employs a tabu semi-greedy heuristic, known to be an effective search strategy in combinatorial optimization. Folding pathways, sometimes called routes or trajectories, are computed by RNAtabupath in a fraction of the time required by the barriers program of Vienna RNA Package. We benchmark RNAtabupath with other algorithms to compute low energy folding pathways between experimentally known structures of several conformational switches. The RNApathfinder web server, source code for algorithms to compute and analyze pathways and supplementary data are available at http://bioinformatics.bc.edu/clotelab/RNApathfinder.     

Cellular pathways controlling integron cassette site folding

By mobilizing small DNA units, integrons have a major function in the dissemination of antibiotic resistance among bacteria. The acquisition of gene cassettes occurs by recombination between the attI and attC sites catalysed by the IntI1 integron integrase. These recombination reactions use an unconventional mechanism involving a folded single‐stranded attC site. We show that cellular bacterial processes delivering ssDNA, such as conjugation and replication, favour proper folding of the attC site. By developing a very sensitive in vivo assay, we also provide evidence that attC sites can recombine as cruciform structures by extrusion from double‐stranded DNA. Moreover, we show an influence of DNA superhelicity on attC site extrusion in vitro and in vivo. We show that the proper folding of the attC site depends on both the propensity to form non‐recombinogenic structures and the length of their variable terminal structures. These results draw the network of cell processes that regulate integron recombination.     

Minimal folding pathways for coarse-grained biopolymer fragments

The minimal folding pathway or trajectory for a biopolymer can be defined as the transformation that minimizes the total distance traveled between a folded and an unfolded structure. This involves generalizing the usual Euclidean distance from points to one-dimensional objects such as a polymer. We apply this distance here to find minimal folding pathways for several candidate protein fragments, including the helix, the β-hairpin, and a nonplanar structure where chain noncrossing is important. Comparing the distances traveled with root mean-squared distance and mean root-squared distance, we show that chain noncrossing can have large effects on the kinetic proximity of apparently similar conformations. Structures that are aligned to the β-hairpin by minimizing mean root-squared distance, a quantity that closely approximates the true distance for long chains, show globally different orientation than structures aligned by minimizing root mean-squared distance.     

Optimum folding pathways for growing protein chains

The folding of a protein is studied as it grows residue by residue from the N-terminus and enters an environment that stabilizes the folded state. This mode of folding of a growing chain is different from refolding where the full chain folds from a disordered initial configuration to the native state. We propose a sequential dynamic optimization method that computes the evolution of optimum folding pathways as amino acid residues are added to the peptide chain one by one. The dynamic optimization formulation is deterministic and uses Newton's equations of motion and a Go-type potential that establishes the native contacts and excluded volume effects. The method predicts the optimal energy-minimizing path among all the alternative feasible pathways. As two examples, the folding of the chicken villin headpiece, a 36-residue protein, and chymotrypsin inhibitor 2 (CI2), a 64-residue protein, are studied. Results on the villin headpiece show significant differences from the refolding of the same chain studied previously. Results on CI2 mostly agree with the results of refolding experiments and computational work.