Intracellular folding pathway of human chorionic gonadotropin beta subunit.

Few experimental models have been used to investigate how proteins fold inside a cell. Using the formation of disulfide bonds as an index of conformational changes during protein folding, we have developed a unique system to determine the intracellular folding pathway of the beta subunit of human chorionic gonadotropin (hCG). Three folding intermediates of the beta subunit were purified from [35S]cysteine-labeled JAR choriocarcinoma cells by immunoprecipitation and by reverse-phase high performance liquid chromatography (HPLC). To identify unformed disulfide bonds, nonreduced folding intermediates were treated with trypsin to liberate non-disulfide-bound, [35S]cysteine-containing peptides from the disulfide-linked peptides. Released peptides were purified by HPLC and identified by amino acid sequencing. The amount of a peptide that was released indicated the extent of disulfide bond formation involving the cysteine in that peptide. Of the six disulfide bonds in hCG-beta, bonds 34-88 and 38-57 form first. The rate-limiting event of folding involves the formation of the S-S bonds between cysteines 23 and 72 and cysteines 9 and 90. Disulfide bond 93-100, the formation of which appears to be necessary for assembly with the alpha subunit of the hCG heterodimer, forms next. Finally, disulfide bond 26-110 forms after assembly with the alpha subunit, suggesting that completion of folding of the COOH terminus in the beta subunit occurs after assembly with the alpha subunit.     

The NMR structures of the major intermediates of the two-domain tick carboxypeptidase inhibitor reveal symmetry in its folding and unfolding pathways

There is a lack of experimental structural information about folding intermediates of multidomain proteins. Tick carboxypeptidase inhibitor (TCI) is a small, disulfide-rich protein consisting of two domains that fold and unfold autonomously through the formation of two major intermediates, IIIa and IIIb. Each intermediate contains three native disulfide bonds in one domain and six free cysteines in the other domain. Here we have determined the NMR structures of these two intermediates trapped and isolated at acidic pH in which they are stable and compared their structures with that of the native protein analyzed under the same conditions. Both IIIa and IIIb were found to contain a folded region that corresponds to the N- and C-terminal domains of TCI, respectively, with structures very similar to the corresponding regions of the native protein. The remainder of the polypeptide chains of the intermediates was shown to be unfolded in a random coil conformation. Solvent exchange measurements further indicated that the two protein domains are not completely independent, but affect each other in terms of dynamics and stability, in agreement with reported inhibitory activity data. The derived results provide structural evidence for symmetric TCI folding and unfolding mechanisms that converge in IIIa and IIIb and reveal the structural basis that accounts for the strong and simultaneous accumulation of both intermediates. Altogether, this work has important implications for a better understanding of the folding mechanisms of multidomain, disulfide-rich proteins.     

Effect of a single amino acid substitution on the folding of the alpha subunit of tryptophan synthase

The urea-induced unfolding of a missense mutant of the alpha subunit of tryptophan synthase from Escherichia coli involving the replacement of Gly by Glu at position 211 has been monitored by absorbance changes at 286 nm. Like the wild-type protein, the equilibrium unfolding curve demonstrates the presence of one or more stable intermediates. Comparison of these results with those from the wild-type alpha subunit [Matthews, C. R., & Crisanti, M. M. (1981) Biochemistry 20, 784] shows that the transition from the native conformation to the stable intermediates is displaced to higher urea concentration in the mutant alpha subunit; however, the transition from the intermediates to the unfolded form is unaffected. Kinetic studies show that the amino acid replacement slows the rate of unfolding by an order of magnitude. The effect on refolding rates is complex. One phase, previously assigned to proline isomerization [Crisanti, M. M., & Matthews, C. R. (1981) Biochemistry 20, 2700], is unaffected by the substitution. The rate of the second phase, which is urea dependent down to about 1 M urea, is slower than the corresponding phase in the wild-type protein by approximately a factor of 2. Below about 1 M urea, the rate of this phase becomes urea independent and identical with that of the wild-type alpha subunit. This change in urea dependence has been ascribed to a change in the nature of the rate-limiting step for this process from one involving folding to one involving proline isomerization. The results support the folding model for the alpha subunit proposed previously [Matthews, C. R., & Crisanti, M. M. (1981) Biochemistry 20, 784] and clarify the role of proline isomerization in limiting the rate of folding.     

Molecular dynamics simulations of folding processes of a beta-hairpin in an implicit solvent

Computer simulations of beta-hairpin folding are relatively difficult, especially those based on the explicit water model. This greatly limits the complete analysis and understanding of their folding mechanisms. In this paper, we use the generalized Born/solvent accessible implicit solvent model to simulate the folding processes of a nine-residue beta-hairpin. We find that the beta-hairpin can fold into its native structure very easily, even using the traditional molecular dynamics method. This allows us to extract 21 complete folding events and investigate the folding process sufficiently. Our results show that there exist four most stable states on the free energy landscape of the short peptide, one native state and three intermediates. We find that two of the non-native stable states have almost the same potential energy as the native state but with lower entropy. This suggests that the native state can be stabilized entropically. Furthermore, we find that the folding processes of this peptide have common features: to fold into its native state, the peptide undergoes a continuous collapsing-extending-recollapsing process to adjust the positions of the side chains in order to form the native middle inter-strand hydrogen bonds. The formations of these bonds are the key step of the folding process. Once these bonds are formed, the peptide can fold into the native state quickly.     

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.     

Rapid folding and unfolding of Apaf-1 CARD

Caspase recruitment domains (CARDs) are members of the death domain superfamily and contain six antiparallel helices in an alpha-helical Greek key topology. We have examined the equilibrium and kinetic folding of the CARD of Apaf-1 (apoptotic protease activating factor 1), which consists of 97 amino acid residues, at pH 6 and pH 8. The results showed that an apparent two state equilibrium mechanism is not adequate to describe the folding of Apaf-1 CARD at either pH, suggesting the presence of intermediates in equilibrium unfolding. Interestingly, the results showed that the secondary structure is less stable than the tertiary structure, based on the transition mid-points for unfolding. Single mixing and sequential mixing stopped-flow studies showed that Apaf-1 CARD folds and unfolds rapidly and suggest a folding mechanism that contains parallel channels with two unfolded conformations folding to the native conformation. Kinetic simulations show that a slow folding phase is described by a third conformation in the unfolded ensemble that interconverts with one or both unfolded species. Overall, the native ensemble is formed rapidly upon refolding. This is in contrast to other CARDs in which folding appears to be dominated by formation of kinetic traps.     

Effects of the phenylalanine-22----leucine, glutamic acid-49----methionine, glycine-234----aspartic acid, and glycine-234----lysine mutations on the folding and stability of the alpha subunit of tryptophan synthase from Escherichia coli

The effects of four single amino acid replacements on the stability and folding of the alpha subunit of tryptophan synthase from Escherichia coli have been investigated by ultraviolet differences spectroscopy. In previous studies [Miles, E. W., Yutani, K., & Ogasahara, K. (1982) Biochemistry 21, 2586], it had been shown that the urea-induced unfolding at pH 7.8, 25 degrees C, proceeds by the initial unfolding of the less stable carboxyl domain (residues 189-268) followed by the unfolding of the more stable amino domain (residues 1-188). The effects of the Phe-22----Leu, Glu-49----Met, Gly-234----Asp, and Gly-234----Lys mutants on the equilibrium unfolding process can all be understood in terms of the domain unfolding model. With the exception of the Glu-49----Met replacement, the effects on stability are small. In contrast, the effects of three of the four mutations on the kinetics of interconversion of the native form and one of the stable partially folded intermediates are dramatic. The results for the Phe-22----Leu and Gly-234----Asp mutations indicate that these residues play a key role in the rate-limiting step. The Glu-49----Met mutation increases the stability of the native form with respect to that of the intermediate but does not affect the rate-limiting step. The Gly-234----Lys mutation does not affect either the stability or the kinetics of folding for the transition between native and intermediate forms. The changes in stability calculated from the unfolding and refolding rate constants agree quantitatively with those obtained from the equilibrium data. When considered with the results from a previous study on the Gly-211----Glu replacement [Matthews, C. R., Crisanti, M. M., Manz, J. T., & Gepner G. L. (1983) Biochemistry 22, 1445], it can be concluded that the rate-limiting step in the conversion of the intermediate to the native conformation involves either domain association or some other type of molecule-wide phenomenon.     

Fast folding of a ribozyme by stabilizing core interactions: evidence for multiple folding pathways in RNA

Folding of the Tetrahymena ribozyme under physiological conditions in vitro is limited by slow conversion of long-lived intermediates to the active structure. These intermediates arise because the most stable domain of the ribozyme folds 10-50 times more rapidly than the core region containing helix P3. Native gel electrophoresis and time-resolved X-ray-dependent hydroxyl radical cleavage revealed that mutations that weaken peripheral interactions between domains accelerated folding fivefold, while a point mutation that stabilizes P3 enabled 80 % of the mutant RNA to reach the native conformation within 30 seconds at 22 degrees C. The P3 mutation increased the folding rate of the catalytic core as much as 50-fold, so that both domains of the ribozyme were formed at approximately the same rate. The results show that the ribozyme folds rapidly without significantly populating metastable intermediates when native interactions in the ribozyme core are stabilized relative to peripheral structural elements.     

Native hydrogen bonds in a molten globule: the apoflavodoxin thermal intermediate

The structure and energetics of protein-folding intermediates are poorly understood. We have identified, in the thermal unfolding of the apoflavodoxin from Anabaena PCC 7119, an equilibrium intermediate with spectroscopic properties of a molten globule and substantial enthalpy and heat capacity of unfolding. The structure of the intermediate is probed by mutagenesis (and phi analysis) of polar residues involved in surface-exposed hydrogen bonds connecting secondary-structure elements in the native protein. All hydrogen bonds analysed are formed in the molten globule intermediate, either with native strength or debilitated. This suggests the overall intermediate's topology and surface tertiary interactions are close to native, and indicates that hydrogen bonding may contribute significantly to shape the conformation and energetics of folding intermediates.     

Cooperative folding mechanism of a beta-hairpin peptide studied by a multicanonical replica-exchange molecular dynamics simulation

G-peptide is a 16-residue peptide of the C-terminal end of streptococcal protein G B1 domain, which is known to fold into a specific beta-hairpin within 6 micros. Here, we study molecular mechanism on the stability and folding of G-peptide by performing a multicanonical replica-exchange (MUCAREM) molecular dynamics simulation with explicit solvent. Unlike the preceding simulations of the same peptide, the simulation was started from an unfolded conformation without any experimental information on the native conformation. In the 278-ns trajectory, we observed three independent folding events. Thus MUCAREM can be estimated to accelerate the folding reaction more than 60 times than the conventional molecular dynamics simulations. The free-energy landscape of the peptide at room temperature shows that there are three essential subevents in the folding pathway to construct the native-like beta-hairpin conformation: (i) a hydrophobic collapse of the peptide occurs with the side-chain contacts between Tyr45 and Phe52, (ii) then, the native-like turn is formed accompanying with the hydrogen-bonded network around the turn region, and (iii) finally, the rest of the backbone hydrogen bonds are formed. A number of stable native hydrogen bonds are formed cooperatively during the second stage, suggesting the importance of the formation of the specific turn structure. This is also supported by the accumulation of the nonnative conformations only with the hydrophobic cluster around Tyr45 and Phe52. These simulation results are consistent with high phi-values of the turn region observed by experiment.     

NMR Structure of a Monomeric Intermediate on the Evolutionarily Optimized Assembly Pathway of a Small Trimerization Domain

Efficient formation of specific intermolecular interactions is essential for self-assembly of biological structures. The foldon domain is an evolutionarily optimized trimerization module required for assembly of the large, trimeric structural protein fibritin from phage T4. Monomers consisting of the 27 amino acids comprising a single foldon domain subunit spontaneously form a natively folded trimer. During assembly of the foldon domain, a monomeric intermediate is formed on the submillisecond time scale, which provides the basis for two consecutive very fast association reactions. Mutation of an intermolecular salt bridge leads to a monomeric protein that resembles the kinetic intermediate in its spectroscopic properties. NMR spectroscopy revealed essentially native topology of the monomeric intermediate with defined hydrogen bonds and side-chain interactions but largely reduced stability compared to the native trimer. This structural preorganization leads to an asymmetric charge distribution on the surface that can direct rapid subunit recognition. The low stability of the intermediate allows a large free-energy gain upon trimerization, which serves as driving force for rapid assembly. These results indicate different free-energy landscapes for folding of small oligomeric proteins compared to monomeric proteins, which typically avoid the transient population of intermediates.     

Protein disulfide isomerases exploit synergy between catalytic and specific binding domains

Protein disulfide isomerases (PDIs) catalyse the formation of native disulfide bonds in protein folding pathways. The key steps involve disulfide formation and isomerization in compact folding intermediates. The high-resolution structures of the a and b domains of PDI are now known, and the overall domain architecture of PDI and its homologues can be inferred. The isolated a and a′ domains of PDI are good catalysts of simple thiol–disulfide interchange reactions but require additional domains to be effective as catalysts of the rate-limiting disulfide isomerizations in protein folding pathways. The b′ domain of PDI has a specific binding site for peptides and its binding properties differ in specificity between members of the PDI family. A model of PDI function can be deduced in which the domains function synergically: the b′ domain binds unstructured regions of polypeptide, while the a and a′ domains catalyse the chemical isomerization steps.     

Scrambled isomers as key intermediates in the oxidative folding of ligand binding module 5 of the low density lipoprotein receptor

The ligand binding module five (LA5) of the low density lipoprotein receptor is a small, single-domain protein of 40 residues and three disulfide bonds with a calcium binding motif that is essential for its structure and function. Several mutations in LA5 have been reported to cause familial hypercholesterolemia by impairing a proper folding of the module. The current study reports the oxidative folding and reductive unfolding pathways of wild type and mutant LA5 modules through kinetic and structural analysis of the trapped intermediates. Wild type LA5 folding involves an initial phase of nonspecific packing where the sequential oxidation of its cysteines gives rise to complex equilibrated populations of intermediates. In the presence of calcium, the attainment of a coordination-competent conformation becomes the rate-limiting step of folding while binding of the ion funnels both thermodynamically and kinetically the folding reaction toward the native state. In the absence of calcium, a scrambled isomer (termed Xa) constitutes the global free energy minimum of the folding process. Xa and the native form are stable, inter-convertible species whose relative populations at equilibrium appear displaced in disease-linked mutants toward the scrambled form. Because stable scrambled isomers such as Xa avoid the exposition of reactive cysteines in misfolded modules, they might constitute a strategy to prevent wrong interactions with other domains during folding of the receptor. Comparison of the folding pathways of wild type and mutant LA5 provides the molecular basis to understand how LA modules fold into a functional conformation or upon mutation misfold and lead to disease.     

Conformational restrictions on the pathway of folding and unfolding of the pancreatic trypsin inhibitor

The kinetic roles of the partially folded, intermediate protein species with two disulphide bonds in folding and unfolding of the pancreatic trypsin inhibitor have been investigated further. Formation of a second disulphide bond between Cys5 and Cys55 during refolding of the reduced inhibitor, which would yield the species with the 30–51 and 5–55 disulphide bonds and, possibly, the native-like conformation of the protein, is not significant. Instead, three other second disulphide bonds (5–14, 5–38 and 14–38) are formed approximately 10⁵ times more readily, but each of these two-disulphide species then rearranges intramolecularly to the native-like, two-disulphide intermediate. Therefore, the reduced protein does not simply form sequentially the three disulphide bonds of the native state. Unfolding of the native state takes place by the reverse of this process.The kinetic importance for folding and unfolding of this transition between two-disulphide intermediates under the conditions used here was illustrated experimentally by a modified form of the inhibitor in which the thiols of Cys14 and Cys38 were blocked irreversibly. In the folded conformation, this modified protein is more stable to unfolding than normal, but after unfolding cannot readily regain the native-like conformation, because Cys14 or Cys38 are required to be involved in disulphide bonds during the interconversion of the two-disulphide intermediates.Some conception of the conformational transitions that take place at each stage of the folding transition may be inferred from the relative propensities of the six cysteine residues to make or rearrange disulphide bonds. It is concluded that the inhibitor probably does not refold by sequential adoption of the native conformation by the unfolded polypeptide chain. Instead, it appears that essentially all elements of the native conformation are attained simultaneously in the final stage of folding, within an unstable and flexible, yet relatively compact, form of the entire polypeptide chain produced by weak interactions between groups distant in the primary structure.     

Common Structural Transitions in Explicit-Solvent Simulations of Villin Headpiece Folding

Molecular dynamics simulations of protein folding can provide very high-resolution data on the folding process; however, due to computational challenges most studies of protein folding have been limited to small peptides, or made use of approximations such as Gō potentials or implicit solvent models. We have performed a set of molecular dynamics simulations totaling >50 μs on the villin headpiece subdomain, one of the most stable and fastest-folding naturally occurring proteins, in explicit solvent. We find that the wild-type villin headpiece reliably folds to a native conformation on timescales similar to experimentally observed folding, but that a fast folding double-norleucine mutant shows significantly more heterogeneous behavior. Along with other recent simulation studies, we note the occurrence of nonnative structures intermediates, which may yield a nativelike signal in the fluorescence measurements typically used to study villin folding. Based on the wild-type simulations, we propose alternative approaches to measure the formation of the native state.     

Deciphering the Structural Basis That Guides the Oxidative Folding of Leech-derived Tryptase Inhibitor

Protein folding mechanisms have remained elusive mainly because of the transient nature of intermediates. Leech-derived tryptase inhibitor (LDTI) is a Kazal-type serine proteinase inhibitor that is emerging as an attractive model for folding studies. It comprises 46 amino acid residues with three disulfide bonds, with one located inside a small triple-stranded antiparallel β-sheet and with two involved in a cystine-stabilized α-helix, a motif that is widely distributed in bioactive peptides. Here, we analyzed the oxidative folding and reductive unfolding of LDTI by chromatographic and disulfide analyses of acid-trapped intermediates. It folds and unfolds, respectively, via sequential oxidation and reduction of the cysteine residues that give rise to a few 1- and 2-disulfide intermediates. Species containing two native disulfide bonds predominate during LDTI folding (IIa and IIc) and unfolding (IIa and IIb). Stop/go folding experiments demonstrate that only intermediate IIa is productive and oxidizes directly into the native form. The NMR structures of acid-trapped and further isolated IIa, IIb, and IIc reveal global folds similar to that of the native protein, including a native-like canonical inhibitory loop. Enzyme kinetics shows that both IIa and IIc are inhibitory-active, which may substantially reduce proteolysis of LDTI during its folding process. The results reported show that the kinetics of the folding reaction is modulated by the specific structural properties of the intermediates and together provide insights into the interdependence of conformational folding and the assembly of native disulfides during oxidative folding.     

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.     

Molecular dynamics simulations of synthetic peptide folding

Because the time scale of protein folding is much greater than that of the widely used simulations of native structures, a detailed report of molecular dynamics simulations of folding has not been available. In this study, we Included the average solvent effect in the potential functions to simplify the calculation of the solvent effect and carried out long molecular dynamics simulations of the alanine-based synthetic peptides at 274 K. From either an extended or a randomly generated conformation, the simulations approached a helix-coil equilibrium in about 3 ns. The multiple minima problem did not prevent helix folding. The calculated helical ratio of Ac-AAQAAAAQAAAAQAAY-NH₂ was 47%, in good agreement with the circular dichroism measurement (about 50%). A helical segment with frayed ends was the most stable conformation, but the hydrophobic interaction favored the compact, distorted helix-turn-helix conformations. The transition between the two types of conformations occurred in a much larger time scale than helix propagation. The transient hydrogen bonds between the glutamine side chain and the backbone carbonyl group could reduce the free energy barrier of helix folding and unfolding. The substitution of a single alanine residue in the middle of the peptide with valine or glycine decreased the average helical ratio significantly, in agreement with experimental observations. © 1996 Wiley-Liss, Inc.     

Atomistic Picture for the Folding Pathway of a Hybrid-1 Type Human Telomeric DNA G-quadruplex

In this work we studied the folding process of the hybrid-1 type human telomeric DNA G-quadruplex with solvent and ions explicitly modeled. Enabled by the powerful bias-exchange metadynamics and large-scale conventional molecular dynamic simulations, the free energy landscape of this G-DNA was obtained for the first time and four folding intermediates were identified, including a triplex and a basically formed quadruplex. The simulations also provided atomistic pictures for the structures and cation binding patterns of the intermediates. The results showed that the structure formation and cation binding are cooperative and mutually supporting each other. The syn/anti reorientation dynamics of the intermediates was also investigated. It was found that the nucleotides usually take correct syn/anti configurations when they form native and stable hydrogen bonds with the others, while fluctuating between two configurations when they do not. Misfolded intermediates with wrong syn/anti configurations were observed in the early intermediates but not in the later ones. Based on the simulations, we also discussed the roles of the non-native interactions. Besides, the formation process of the parallel conformation in the first two G-repeats and the associated reversal loop were studied. Based on the above results, we proposed a folding pathway for the hybrid-1 type G-quadruplex with atomistic details, which is new and more complete compared with previous ones. The knowledge gained for this type of G-DNA may provide a general insight for the folding of the other G-quadruplexes.     

Kinetic characterization of early intermediates in the folding of E. coli tryptophan-synthase beta 2 subunit

This report describes the use of fluorescence energy transfer between an intrinsic energy donor (tryptophan 177) and two chemically added acceptors to study intermediates in the folding of the beta 2 subunit of E. coli tryptophan-synthase. Two early folding steps are thus identified and characterized. One is very rapid (its rate constant at 12 degrees C is 0.02 sec-1) and corresponds to the folding of the N-terminal domain into a structure whose overall features approximate well those of the native domain. The second step is somewhat slower (its rate constant at 12 degrees C is 0.008 sec-1) and involves a conformational rearrangement of the N-terminal domain brought about by the interactions between the N- and C-terminal domains within a monomeric beta chain. This brings to five the number of intermediates which have been identified and ordered on the folding pathway of the dimeric beta 2 subunit.