Dynamic Monte Carlo simulations of globular protein folding/unfolding pathways. II. Alpha-helical motifs

Dynamic Monte Carlo simulations of the folding pathways of alpha-helical protein motifs have been undertaken in the context of a diamond lattice model of globular proteins. The first question addressed in the nature of the assembly process of an alpha-helical hairpin. While the hairpin could, in principle, be formed via the diffusion-collision-adhesion of isolated performed helices, this is not the dominant mechanism of assembly found in the simulations. Rather, the helices that form native hairpins are constructed on-site, with folding initiating at or near the turn in almost all cases. Next, the folding/unfolding pathways of four-helix bundles having tight bends and one and two long loops in the native state are explored. Once again, an on-site construction mechanism of folding obtains, with a hairpin forming first, followed by the formation of a three-helix bundle, and finally the fourth helix of the native bundle assembles. Unfolding is essentially the reverse of folding. A simplified analytic theory is developed that reproduces the equilibrium folding transitions obtained from the simulations remarkably well and, for the dominant folding pathway, correctly identifies the intermediates seen in the simulations. The analytic theory provides the free energy along the reaction co-ordinate and identifies the transition state for all three motifs as being quite close to the native state, with three of the four helices assembled, and approximately one turn of the fourth helix in place. The transition state is separated from the native conformation by a free-energy barrier of mainly energetic origin and from the denatured state by a barrier of mainly entropic origin. The general features of the folding pathway seen in all variants of the model four-helix bundles are similar to those observed in the folding of beta-barrel, Greek key proteins; this suggests that many of the qualitative aspects of folding are invariant to the particular native state topology and secondary structure.     

Folding of bovine growth hormone is consistent with the molten globule hypothesis

Previous results from equilibrium and kinetic studies of the folding of bovine growth hormone (bGH) have demonstrated that bGH does not follow a simple two-step folding mechanism. These results are summarized and interpreted according to the "molten globule" model. The molten globule state of bGH is characterized as a folding intermediate which is largely alpha-helical, retains a compact hydrodynamic radius, has packing of the aromatic side chains that is similar to the unfolded state, and possesses a solvent-exposed hydrophobic surface along helix 106-127 that readily leads to association.     

Direct kinetic evidence for folding via a highly compact, misfolded state.

The 2 S seed storage protein, sunflower albumin 8 (SFA-8), contains an unusually high proportion of hydrophobic residues including 16 methionines (some of which may form a surface hydrophobic patch) in a disulfide cross-linked, alpha-helical structure. Circular dichroism and fluorescence spectroscopy show that SFA-8 is highly stable to denaturation by heating or chaotropic agents, the latter resulting in a reversible two-state unfolding transition. The small m(U) (-4.7 M(-1) at 10 degrees C) and DeltaC(p) (-0.95 kcal mol(-1) K(-1)) values indicate that relatively little nonpolar surface of the protein is exposed during unfolding. Commensurate with the unusual distribution of hydrophobic residues, stopped-flow fluorescence data show that the folding pathway of SFA-8 is highly atypical, in that the initial product of the rapid collapse phase of folding is a compact nonnative state (or collection of nonnative states) that must unfold before acquiring the native conformation. The inhibited folding reaction of SFA-8, in which the misfolded state (m(M) = -0.95 M(-1) at 10 degrees C) is more compact than the transition state for folding (m(T) = -2.5 M(-1) at 10 degrees C), provides direct kinetic evidence for the transient misfolding of a protein.     

Understanding the mechanism of beta-hairpin folding via phi-value analysis

The folding kinetics of a 16-residue beta-hairpin (trpzip4) and five mutants were studied by a laser-induced temperature-jump infrared method. Our results indicate that mutations which affect the strength of the hydrophobic cluster lead to a decrease in the thermal stability of the beta-hairpin, as a result of increased unfolding rates. For example, the W45Y mutant has a phi-value of approximately zero, implying a folding transition state in which the native contacts involving Trp45 are not yet formed. On the other hand, mutations in the turn or loop region mostly affect the folding rate. In particular, replacing Asp46 with Ala leads to a decrease in the folding rate by roughly 9 times. Accordingly, the phi-value for D46A is determined to be approximately 0.77, suggesting that this residue plays a key role in stabilizing the folding transition state. This is most likely due to the fact that the main chain and side chain of Asp46 form a characteristic hydrogen bond network with other residues in the turn region. Taken together, these results support the folding mechanism we proposed before, which suggests that the turn formation is the rate-limiting step in beta-hairpin folding and, consequently, a stronger turn-promoting sequence increases the stability of a beta-hairpin primarily by increasing its folding rate, whereas a stronger hydrophobic cluster increases the stability of a beta-hairpin primarily by decreasing its unfolding rate. In addition, we have examined the compactness of the thermally denatured and urea-denatured states of another 16-residue beta-hairpin, using the method of fluorescence resonance energy transfer. Our results show that the thermally denatured state of this beta-hairpin is significantly more compact than the urea-denatured state, suggesting that the very first step in beta-hairpin folding, when initiated from an extended conformation, probably corresponds to a process of hydrophobic collapse.     

Dissecting the stability of a beta-hairpin peptide that folds in water: NMR and molecular dynamics analysis of the beta-turn and beta-strand contributions to folding.

NMR studies of the folding and conformational properties of a beta-hairpin peptide, several peptide fragments of the hairpin, and sequence-modified analogues, have enabled the various contributions to beta-hairpin stability in water to be dissected. Temperature and pH-induced unfolding studies indicate that the folding-unfolding equilibrium approximates to a two-state model. The hairpin is highly resistant to denaturation and is still significantly folded in 7 M urea at 298 K. Thermodynamic analysis shows the hairpin to fold in water with a significant change in heat capacity, however, DeltaCp degrees in 7 M urea is reduced. V/Y-->A mutations on one strand of the hairpin reduce folding to <10 %, consistent with a hydrophobic stabilisation model. We show that in a truncated peptide (residues 6-16) lacking the hydrophobic residues on one beta-strand, the type I' Asn-Gly turn in the sequence SINGKK is significantly populated in water in the absence of interstrand hydrophobic contacts. Unrestrained molecular dynamics simulations of unfolding, using an explicit solvation model, show that the conformation of the NG turn persists for longer than the AG analogue, which has a much lower propensity for type I' turn formation from a data base analysis of preferred turns. The origin of the high stability of the Asn-Gly turn is not entirely clear; data base analysis of 66 NG turns, together with molecular dynamics simulations, reveals no participation of the Asn side-chain in turn-stabilising interactions with the peptide backbone. However, hydration analysis of the molecular dynamics simulations reveals a pocket of "high density" water bridging between the Asn side-chain and peptide main-chain that suggests solvent-mediated interactions may play an important role in modulating phi,psi propensities in the NG turn region.     

Non-functional conserved residues in globins and their possible role as a folding nucleus.

Structure-based sequence alignment of 728 sequences of different globin subfamilies shows that in each subfamily there are two clusters of consensually conserved residues. The first is the well-known "functional" cluster which includes six heme-binding conserved residues (Phe CD1, His F8; aliphatic E11, FG5; hydrophobic F4, G5) and seven other conserved residues (Pro C2; aliphatic H19; hydrophobic B10, B13, B14, CD4, E4) that do not bind the heme but belong to its immediate neighborhood. The second cluster revealed here (aliphatic A8, G16, G12; aromatic A12; hydrophobic H8 and possibly H12) is distant from the heme. It is entirely non-polar and includes one turn (i, i+4 positions) from each of helices A, G, and H. It is known that A, G, and H helices formed at the earliest stage of apomyoglobin folding remain relatively stable in the equilibrium molten globule state, and are likely to be tightly packed with each other in this state. We have shown the existence of two similar conserved clusters in c -type cytochromes, heme-binding and distal from the heme. The second cluster in c -cytochromes includes one turn from each of the N and C-terminal alpha-helices. These N and C-terminal helices in cytochrome c are formed at the earliest stage of protein folding, remain relatively stable in the molten globule state, and are tightly packed with each other in this state, similar to the observed behavior of the globins. At least these two large protein families (c -type cytochromes and globins) have a close similarity in the existence and mutual positions of non-functional conserved residues. We assume that non-functional conserved residues are requisite for the fast and correct folding of both of these protein families into their stable 3D structures.     

Design and solution structure of a well-folded stack of two beta-hairpins based on the amino-terminal fragment of human granulin A

Four amino acid substitutions were introduced into a peptide corresponding to the amino-terminal subdomain (30-31 residues) of human granulin A (HGA) in order to assess the contributions of a hydrophobic framework and other interactions to structure stabilization of the stack of two beta-hairpins. The resulting hybrid peptide, HGA 1-31 (D1V, K3H, S9I, Q20P) with four free cysteines, spontaneously formed a uniquely disulfide-bonded isomer with an expected [1-3, 2-4] disulfide pairing pattern. This peptide was characterized in detail by use of NMR and shown to assume a highly stable structure in solution, in contrast to the amino-terminal 1-30 fragment of human granulin A. The prototype peptide, or HGA 1-30 (C17S, C27S), had lower resistance to chemical reduction and proteolysis, broad NH and H(alpha) proton resonances, lower proton resonance dispersion, and no slowly exchanging amide protons. Two other peptides, HGA 1-30 (C17S, Q20P, C27S) and HGA 1-31 (D1V, K3H, S9I, C17S, C27S), with either Pro20 stabilizing a potential reverse turn or with a hydrophobic cluster consisting of Val1, His3, and Ile9, had sharper and slightly better dispersed NH and H(alpha) proton resonances, but still no slowly exchanging amide protons. The solution structure of HGA 1-31 (D1V, K3H, S9I, Q20P) indicates that it adopts a well-folded conformation of a stack of two beta-hairpins, as found for the amino-terminal subdomain of the prototypic carp granulin-1 with representative beta-hairpin stacks. These results highlight the importance of both hydrophobic and turn-stabilizing interactions for the structural integrity of the hairpin stack scaffold. The conformational stability appears to be maintained by a combination of the well-formed second beta-hairpin and two hydrophobic clusters, one located at the interface between the two beta-hairpins and the other on "top" of the first beta-hairpin. The implications of these findings for the design of conformationally stable hairpin stacks are discussed.     

Folding and oxidation of the antibody domain C(H)3

The non-covalent homodimer formed by the C-terminal domains of the IgG1 heavy chains (C(H)3) is the simplest naturally occurring model system for studying immunoglobulin folding and assembly. In the native state, the intrachain disulfide bridge, which connects a three-stranded and a four-stranded beta-sheet is buried in the hydrophobic core of the protein. Here, we show that the disulfide bridge is not required for folding and association, since the reduced C(H)3 domain folds to a dimer with defined secondary and tertiary structure. However, the thermodynamic stability of the reduced C(H)3 dimer is much lower than that of the oxidized state. This allows the formation of disulfide bonds either concomitant with folding (starting from the reduced, denatured state) or after folding (starting from the reduced dimer). The analysis of the two processes revealed that, under all conditions investigated, one of the cysteine residues, Cys 86, reacts preferentially with oxidized glutathione to a mixed disulfide that subsequently interacts with the less-reactive second thiol group of the intra-molecular disulfide bond. For folded C(H)3, the second step in the oxidation process is slow. In contrast, starting from the unfolded and reduced protein, the oxidation reaction is faster. However, the overall folding reaction of C(H)3 during oxidative folding is a slow process. Especially, dimerization is slow, compared to the association starting from the denatured oxidized state. This deceleration may be due to misfolded conformations trapped by the disulfide bridge.     

Predicted Folding of β-Structure in Myelin Basic Protein

Predictions of myelin basic protein secondary structure have not previously considered a major role for beta-structure in the organization of the native molecule because optical rotatory dispersion and circular dichroism studies have provided little, if any, evidence for beta-structure, and because a polycationic protein is generally considered to resist folding into a compact structure. However, the Chou-Fasman, Lim, and Robson algorithms identify a total of five beta-strands in the amino acid sequence. Four of these hydrophobic amino acid sequences (37-45, 87-95, 110-118, and 150-158) could form a hairpin intermediate that initiates folding of a Greek-key-type beta-structure. A second fold on the more hydrophobic side, with the addition of a strand from the N-terminus (residues 13-21), would complete the five-stranded antiparallel beta-sheet. A unique strand alignment can be predicted by phasing the hydrophobic residues. The unusual triproline sequence of myelin basic protein (100-102) is enclosed in the 14-residue hairpin loop. If these prolines are in the trans conformation, models show that a reverse turn could occur at residues 102-105 (Pro-Ser-Gln-Gly). Algorithms do not agree on the prediction of alpha-helices, but each of the two large loops could accommodate an alpha-helix. Myelin basic protein is known to be phosphorylated in vivo on as many as five Ser/Thr residues. Phosphorylation might alter the dynamics of folding if the nascent polypeptide were phosphorylated in the cytoplasm. In particular, phosphorylation of Thr-99 could neutralize cationic residues Lys-106 and Arg-108 within the hairpin loop. In addition, the methylation of Arg-108 might stabilize the hairpin loop structure through hydrophobic interaction with the side chain of Pro-97. The cationic side chains of arginine and lysine residues located on the faces of the beta-sheet (Arg-43, Arg-114, Lys-13, Lys-92, Lys-153, and Lys-156) could provide sites for interaction with phospholipids and other anionic structures on the surface of the myelin lipid bilayer.     

Disulfide bonds and protein folding

The applications of disulfide-bond chemistry to studies of protein folding, structure, and stability are reviewed and illustrated with bovine pancreatic ribonuclease A (RNase A). After surveying the general properties and advantages of disulfide-bond studies, we illustrate the mechanism of reductive unfolding with RNase A, and discuss its application to probing structural fluctuations in folded proteins. The oxidative folding of RNase A is then described, focusing on the role of structure formation in the regeneration of the native disulfide bonds. The development of structure and conformational order in the disulfide intermediates during oxidative folding is characterized. Partially folded disulfide species are not observed, indicating that disulfide-coupled folding is highly cooperative. Contrary to the predictions of "rugged funnel" models of protein folding, misfolded disulfide species are also not observed despite the potentially stabilizing effect of many nonnative disulfide bonds. The mechanism of regenerating the native disulfide bonds suggests an analogous scenario for conformational folding. Finally, engineered covalent cross-links may be used to assay for the association of protein segments in the folding transition state, as illustrated with RNase A.     

Mechanism of protein folding. IV. Forming and breaking of disulfide bonds in bovine pancreatic tripsin inhibitor.

The folding mechanism of bovine pancreatic tripsin inhibitor (BPTI) is explained theoretically on the basis of the island model, where the driving force of folding is hydrophobic interaction. For this purpose, we take a look at the formation and breaking of disulfide bonds during the folding process of BPTI. The intermediate conformations and the native one are successfully obtained, which satisfy the so-called "lampshade" geometrical criterion for the formation of the disulfide bonds. The folding pathway is consistent with the renaturation experiment by Creighton. In addition, an elaborate treatment of side chains of amino acid residues by the software programme CHARMm confirms quantitatively the formation of disulfide bridges.     

Biphasic folding kinetics of RNA pseudoknots and telomerase RNA activity

Using a combined master equation and kinetic cluster approach, we investigate RNA pseudoknot folding and unfolding kinetics. The energetic parameters are computed from a recently developed Vfold model for RNA secondary structure and pseudoknot folding thermodynamics. The folding kinetics theory is based on the complete conformational ensemble, including all the native-like and non-native states. The predicted folding and unfolding pathways, activation barriers, Arrhenius plots, and rate-limiting steps lead to several findings. First, for the PK5 pseudoknot, a misfolded 5' hairpin emerges as a stable kinetic trap in the folding process, and the detrapping from this misfolded state is the rate-limiting step for the overall folding process. The calculated rate constant and activation barrier agree well with the experimental data. Second, as an application of the model, we investigate the kinetic folding pathways for human telomerase RNA (hTR) pseudoknot. The predicted folding and unfolding pathways not only support the proposed role of conformational switch between hairpin and pseudoknot in hTR activity, but also reveal molecular mechanism for the conformational switch. Furthermore, for an experimentally studied hTR mutation, whose hairpin intermediate is destabilized, the model predicts a long-lived transient hairpin structure, and the switch between the transient hairpin intermediate and the native pseudoknot may be responsible for the observed hTR activity. Such finding would help resolve the apparent contradiction between the observed hTR activity and the absence of a stable hairpin.     

The role of the turn in beta-hairpin formation during WW domain folding

The folding of WW domains is rate limited by formation of a beta-hairpin comprising residues from strands 1 and 2. Residues in the turn of this hairpin have reported Phi-values for folding close to 1 and have been proposed to nucleate folding. High Phi-values do not necessarily imply that the energetics of formation are a driving force for initiating folding. We demonstrate by NMR studies and molecular dynamics simulations that the first turn of the hYAP, FBP28, and PIN1 WW domains is structurally dynamic and solvent exposed in the native and folding transition states. It is, therefore, unlikely that the formation of the beta-turn per se provides the energetic driving force for hairpin folding. It is more likely that the turn acts as an easily formed hinge that facilitates the formation of the hairpin; it is a nucleus as defined by the nucleation-condensation mechanism whereby a diffuse nucleus is stabilized by associated interactions.     

Oxidative folding of human lysozyme: effects of the loss of two disulfide bonds and the introduction of a calcium-binding site

Mutant human lysozymes (HLZ) lacking two disulfide bonds were constructed to study the importance of each disulfide bond on oxidative refolding. To avoid destabilization, a calcium-binding site was introduced. Five of the six species of two-disulfide mutants could be obtained with enzymatic activity. Based on the information obtained from refolding and unfolding experiments, the order of importance in oxidative refolding was found to be as follows: SS2(Cys30-Cys116) > SS1(Cys6-Cys128) SS3(Cys65-Cys81) > SS4(Cys77-Cys95). Without SS2, these mutants refolded with low efficiency or did not refold at all. The bond SS2 is located in the interface of B-and D-helices, and a small hydrophobic cluster is formed near SS2. This cluster may play an important role in the folding process and stabilization, and SS2 may act as a stabilizer through its polypeptide linkage. The bond SS2 is the most important disulfide bond for oxidative folding of lysozymes.     

The structural basis of compstatin activity examined by structure-function-based design of peptide analogs and NMR

We have previously identified compstatin, a 13-residue cyclic peptide, that inhibits complement activation by binding to C3 and preventing C3 cleavage to C3a and C3b. The structure of compstatin consists of a disulfide bridge and a type I beta-turn located at opposite sides to each other. The disulfide bridge is part of a hydrophobic cluster, and the beta-turn is part of a polar surface. We present the design of compstatin analogs in which we have introduced a series of perturbations in key structural elements of their parent peptide, compstatin. We have examined the consistency of the structures of the designed analogs compared with compstatin using NMR, and we have used the resulting structural information to make structure-complement inhibitory activity correlations. We propose the following. 1) Even in the absence of the disulfide bridge, a linear analog has a propensity for structure formation consistent with a turn of a 3(10)-helix or a beta-turn. 2) The type I beta-turn is a necessary but not a sufficient condition for activity. 3) Our substitutions outside the type I beta-turn of compstatin have altered the turn population but not the turn structure. 4) Flexibility of the beta-turn is essential for activity. 5) The type I beta-turn introduces reversibility and sufficiently separates the two sides of the peptide, whereas the disulfide bridge prevents the termini from drifting apart, thus aiding in the formation of the hydrophobic cluster. 6) The hydrophobic cluster at the linked termini is involved in binding to C3 and activity but alone is not sufficient for activity. 7) beta-Turn residues Gln(5) (Asn(5))-Asp(6)-Trp(7)(Phe(7))-Gly(8) are specific for the turn formation, but only Gln(5)(Asn(5))-Asp(6)-Trp(7)-Gly(8) residues are specific for activity. 8) Trp(7) is likely to be involved in direct interaction with C3, possibly through the formation of a hydrogen bond. Finally we propose a binding model for the C3-compstatin complex.     

Understanding beta-hairpin formation by molecular dynamics simulations of unfolding

We have studied the mechanism of formation of a 16-residue beta-hairpin from the protein GB1 using molecular dynamics simulations in an aqueous environment. The analysis of unfolding trajectories at high temperatures suggests a refolding pathway consisting of several transient intermediates. The changes in the interaction energies of residues are related with the structural changes during the unfolding of the hairpin. The electrostatic energies of the residues in the turn region are found to be responsible for the transition between the folded state and the hydrophobic core state. The van der Waals interaction energies of the residues in the hydrophobic core reflect the behavior of the radius of gyration of the core region. We have examined the opposing influences of the protein-protein (PP) energy, which favors the native state, and the protein-solvent (PS) energy, which favors unfolding, in the formation of the beta-hairpin structure. It is found that the behavior of the electrostatic components of PP and PS energies reflects the structural changes associated with the loss of backbone hydrogen bonding. Relative changes in the PP and PS van der Waals interactions are related with the disruption of the hydrophobic core of a protein. The results of the simulations support the hydrophobic collapse mechanism of beta-hairpin folding.     

Engineering stabilising beta-sheet interactions into a conformationally flexible region of the folding transition state of ubiquitin

Protein engineering studies suggest that the transition state for the folding of ubiquitin is highly polarised towards the N-terminal part of the sequence and involves a nucleus of residues within the beta-hairpin (residues 1-17) and main alpha-helix (residues 23-34). In contrast, the observation of small phi-values for residues in the C-terminal portion of the sequence (residues 35-76), coupled with a folding topology that results in a much higher contact order, suggests that fast folding of ubiquitin is dependent upon configurational flexibility in the C-terminal part of the polypeptide chain to ensure passage down a relatively smooth folding funnel to the native state. We show that the introduction of a small mini-hairpin motif as an extension of the native 43-50 hairpin stabilises local interactions in the C-terminal part of the sequence, resulting largely in a deceleration of the unfolding kinetics without perturbing the apparent two-state folding mechanism. However, a single-point Leu-->Phe substitution within the engineered hairpin sequence leads to the premature collapse of the denatured ensemble through the stabilisation of non-native interactions and the population of a compact intermediate. Non-linear effects in the kinetic data at low concentrations of denaturant suggest that the collapsed state, which is further stabilised in the presence of cosmotropic salts, may subsequently fold directly to the native state through a "triangular" reaction scheme involving internal rearrangement rather than unfolding and refolding.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Structure and dynamics of an acid-denatured protein G mutant

NMR studies of protein denatured states provide insights into potential initiation sites for folding that may be too transient to be observed kinetically. We have characterized the structure and dynamics of the acid-denatured state of protein G by using a F30H mutant of G(B1) which is on the margin of stability. At 5 degrees C, F30H-G(B1) is greater than 95% folded at pH 7.0 and is greater than 95% unfolded at pH 4.0. This range of stability is useful because the denatured state can be examined under relatively mild conditions which are optimal for folding G(B1). We have assigned almost all backbone (15)N, H(N), and H(alpha) resonances in the acid-denatured state. Chemical shift, coupling constant, and NOE data indicate that the denatured state has considerably more residual structure when studied under these mild conditions than in the presence of chemical denaturants. The acid-denatured state populates nativelike conformations with both alpha-helical and beta-hairpin characteristics. To our knowledge, this is the first example of a denatured state with NOE and coupling constant evidence for beta-hairpin character. A number of non-native turn structures are also detected, particularly in the region corresponding to the beta1-beta2 hairpin of the folded state. Steady-state ?(1)H-(15)N? NOE results demonstrate restricted backbone flexibility in more structured regions of the denatured protein. Overall, our studies suggest that regions of the helix, the beta3-beta4 hairpin, and the beta1-beta2 turn may serve as potential initiation sites for folding of G(B). Furthermore, residual structure in acid-denatured F30H-G(B1) is more extensive than in peptide fragments corresponding to the beta1-beta2, alpha-helix, and beta3-beta4 regions, suggesting additional medium-to-long-range interactions in the full-length polypeptide chain.     

Pathway of oxidative folding of bovine alpha-interferon: predominance of native disulfide-bonded folding intermediates

Bovine alpha-interferon (BoINF-alpha) is a single polypeptide protein containing 166 amino acids, two disulfide bonds (Cys1-Cys99 and Cys29-Cys138), and five stretches of alpha-helical structure. The pathway of oxidative folding of BoINF-alpha has been investigated here. Of the eight possible one- and two-disulfide isomers, only two nativelike one-disulfide isomers, BoINF-alpha (Cys1-Cys99) and BoINF-alpha (Cys29-Cys138), predominate as intermediates along the folding pathway. More strikingly, alpha-helical structures formed almost quantitatively before any detectable formation of a disulfide bond. This is demonstrated by the observation that fully reduced BoINF-alpha (starting material of oxidative folding) and reduced carboxymethylated BoINF-alpha both exhibit alpha-helical structure content indistinguishable form that of native BoINF-alpha. The folding mechanism of BoINF-alpha appears to be compatible with the framework model, in which secondary structures fold first, followed by docking (compaction) of preformed secondary structural elements yielding the native structure.