Native disulfide bonds greatly accelerate secondary structure formation in the folding of lysozyme.

To assess the respective roles of local and long-range interactions during protein folding, the influence of the native disulfide bonds on the early formation of secondary structure was investigated using continuous-flow circular dichroism. Within the first 4 ms of folding, lysozyme with intact disulfide bonds already had a far-UV CD spectrum reflecting large amounts of secondary structure. Conversely, reduced lysozyme remained essentially unfolded at this early folding time. Thus, native disulfide bonds not only stabilize the cfinal conformation of lysozyme but also provide, in early folding intermediates, the necessary stabilization that favors the formation of secondary structure.     

Cooperative folding of the isolated alpha-helical domain of hen egg-white lysozyme.

Proteins in the alpha-lactalbumin and c-type lysozyme family have been studied extensively as model systems in protein folding. Early formation of the alpha-helical domain is observed in both alpha-lactalbumin and c-type lysozyme; however, the details of the kinetic folding pathways are significantly different. The major folding intermediate of hen egg-white lysozyme has a cooperatively formed tertiary structure, whereas the intermediate of alpha-lactalbumin exhibits the characteristics of a molten globule. In this study, we have designed and constructed an isolated alpha-helical domain of hen egg-white lysozyme, called Lyso-alpha, as a model of the lysozyme folding intermediate that is stable at equilibrium. Disulfide-exchange studies show that under native conditions, the cysteine residues in Lyso-alpha prefer to form the same set of disulfide bonds as in the alpha-helical domain of full-length lysozyme. Under denaturing conditions, formation of the nearest-neighbor disulfide bonds is strongly preferred. In contrast to the isolated alpha-helical domain of alpha-lactalbumin, Lyso-alpha with two native disulfide bonds exhibits a well-defined tertiary structure, as indicated by cooperative thermal unfolding and a well-dispersed NMR spectrum. Thus, the determinants for formation of the cooperative side-chain interactions are located mainly in the alpha-helical domain. Our studies suggest that the difference in kinetic folding pathways between alpha-lactalbumin and lysozyme can be explained by the difference in packing density between secondary structural elements and support the hypothesis that the structured regions in a protein folding intermediate may correspond to regions that can fold independently.     

Evidence for identity between the equilibrium unfolding intermediate and a transient folding intermediate: a comparative study of the folding reactions of alpha-lactalbumin and lysozyme

The refolding kinetics of alpha-lactalbumin at different concentrations of guanidine hydrochloride have been investigated by means of kinetic circular dichroism and stopped-flow absorption measurements. The refolding reaction consists of at least two stages, the instantaneous accumulation of the transient intermediate that has peptide secondary structure and the subsequent slow process associated with formation of tertiary structure. The transient intermediate is compared with the well-characterized equilibrium intermediate observed during the denaturant-induced unfolding. Stabilities of the secondary structures against the denaturant, affinities for Ca2+, and tryptophan absorption properties of the transient and equilibrium intermediates were investigated. In all of these respects, the transient intermediate is identical with the equilibrium one, demonstrating the validity of the use of the equilibrium intermediate as a model of the folding intermediate. Essentially the same transient intermediate was also detected in the folding of lysozyme, the protein known to be homologous to alpha-lactalbumin but whose equilibrium unfolding is represented as a two-state reaction. The stability and cooperativity of the secondary structure of the intermediate of lysozyme are compared with those of alpha-lactalbumin. The results show that the protein folding occurring via the intermediate is not limited to the proteins that show equilibrium intermediates. Although the unfolding equilibria of most proteins are well approximated as a two-state reaction, the two-state hypothesis may not be applicable to the folding reaction under the native condition. Two models of protein folding, intermediate-controlled folding model and multiple-pathway folding model, which are different in view of the role of the intermediate in determining the pathway of folding, are also discussed.     

Identification of the folding inhibitors of hen-egg lysozyme: gathering the right tools

The unfolded state of proteins displays a surprisingly rich amount of local native structure, which appears to be critical for driving the protein to its native state. Peptides with the same sequence of the corresponding structured segments can be used to interfere with the correct folding of the protein. Using model simulations, we investigate the folding of hen-egg lysozyme, identifying its key segments. Activity assays, NMR and circular dichroism experiments are used to screen the peptides which are able to inhibit the folding of lysozyme. Few peptides, corresponding to the segments of the protein which are structured in the unfolded state, are identified to have significant inhibitory effects.     

Native-like tertiary structure formation in the alpha-domain of a hen lysozyme two-disulfide variant.

Structure formation in two species of the two-disulfide variant of hen lysozyme was investigated by means of CD spectroscopy, disulfide exchange measurement, and 1H-NMR spectroscopy. One species, 2SS [6-127, 30-115], which contained the two disulfide bonds found in the alpha-domain of authentic lysozyme, had amounts of secondary and tertiary structures, and bacteriolytic activity comparable to those of authentic lysozyme, and showed a cooperative thermal unfolding. By contrast, the other species, 2SS [64-80, 76-94], which contained the beta-domain disulfide bond as well as the inter-domain one, had a limited amount of secondary structure and little tertiary structure. Disulfide-exchange did not occur for 2SS [6-127, 30-115], whereas it occurred for 2SS [64-80, 76-94], indicating that the protein main-chain fold coupled with the formation of two disulfide bonds is relatively stable for the former variant, while unstable for the latter. 1H-NMR spectra of 2SS [6-127, 30-115] showed that native-like local environment is present within the region that corresponds to the alpha-domain, while it is absent within the region that corresponds to the beta or inter-domain. These results indicate that the alpha-domain of hen lysozyme can be an independent folding domain at equilibrium. Although the bipartite nature in the structure formation of hen lysozyme is similar to that reported for alpha-lactalbumin, differences exist between the disulfide-intermediates of the two proteins in terms of the structural domain that accomplishes tertiary structure.     

Effect of the structure of the denatured state of lysozyme on the aggregation reaction at the early stages of folding from the reduced form

We previously demonstrated that the hydrophobic clusters present in hen lysozyme under denaturing conditions were disrupted by the mutation of Trp62 to Gly (W62G). In order to examine the effects of the structure of the denatured state of W62G lysozyme on folding, we analyzed the early events in the folding of reduced W62G lysozyme in detail. From the exchange measurements of disulfide bonds using the variants containing a pair of cysteine residues (1SS), it was found that the formation of disulfide bond in the W62G1SS lysozyme was not accompanied by a prominent interaction between amino acid residues, indicating that the disruption of the hydrophobic core led to the random folding at the early stages in the process of folding of the reduced lysozyme. On the other hand, analyses of the oxidative-renaturation of reduced W62G lysozymes, as well as measurements of the extent of aggregation of the reduced and carboxy amido methylated W62G lysozyme, indicated that the formation of an aggregate is more prominent in the reduced W62G lysozyme than in the reduced wild-type lysozyme. Moreover, a lag phase was detected in the oxidative-renaturation of reduced W62G lysozyme, as based on observations of the recovery of activity. The simulation of the folding process indicated that intermediates were present at the early stages in the folding of the reduced W62G lysozyme. These results suggest that the presence of the intermediates was derived from the random folding at the early stages in the folding process of reduced W62G lysozyme due to the disruption of the structure of the denatured state. Folding thus appears to have been kinetically delayed by these processes, which then led to the significant aggregation of reduced lysozyme. Moreover, from the analysis of amyloid aggregation of the reduced lysozymes, it was suggested that the disruption of the residual structure in denatured state by W62G mutation deterred the formation of the amyloid fibrils of lysozyme.     

Effect of an alternative disulfide bond on the structure, stability, and folding of human lysozyme

Human lysozyme has four disulfide bonds, one of which, Cys65-Cys81, is included in a long loop of the beta-domain. A cysteine-scanning mutagenesis in which the position of Cys65 was shifted within a continuous segment from positions 61 to 67, with fixed Cys81, has previously shown that only the mutant W64CC65A, which has a nonnative Cys64-Cys81 disulfide, can be correctly folded and secreted by yeast. Here, using the W64CC65A mutant, we investigated the effects of an alternative disulfide bond on the structure, stability, and folding of human lysozyme using circular dichroism (CD) and fluorescence spectroscopy combined with a stopped-flow technique. Although the mutant is expected to have a different main-chain structure from that of the wild-type protein around the loop region, far- and near-UV CD spectra show that the native state of the mutant has tightly packed side chains and secondary structure similar to that of the wild-type. Guanidine hydrochloride-induced equilibrium unfolding transition of the mutant is reversible, showing high stability and cooperativity of folding. In the kinetic folding reaction, both proteins accumulate a similar burst-phase intermediate having pronounced secondary structure within the dead time of the measurement and fold into the native structure by means of a similar folding mechanism. Both the kinetic refolding and unfolding reactions of the mutant protein are faster than those of the wild-type, but the increase in the unfolding rate is larger than that of the refolding rate. The Gibbs' free-energy diagrams obtained from the kinetic analysis suggest that the structure around the loop region in the beta-domain of human lysozyme is formed after the transition state of folding, and thus, the effect of the alternative disulfide bond on the structure, stability, and folding of human lysozyme appears mainly in the native state.     

Multi-State Proteins: Approach Allowing Experimental Determination of the Formation Order of Structure Elements in the Green Fluorescent Protein

The most complex problem in studying multi-state protein folding is the determination of the sequence of formation of protein intermediate states. A far more complex issue is to determine at what stages of protein folding its various parts (secondary structure elements) develop. The structure and properties of different intermediate states depend in particular on these parts. An experimental approach, named μ-analysis, which allows understanding the order of formation of structural elements upon folding of a multi-state protein was used in this study. In this approach the same elements of the protein secondary structure are “tested” by substitutions of single hydrophobic amino acids and by incorporation of cysteine bridges. Single substitutions of hydrophobic amino acids contribute to yielding information on the late stages of protein folding while incorporation of ss-bridges allows obtaining data on the initial stages of folding. As a result of such an μ-analysis, we have determined the order of formation of beta-hairpins upon folding of the green fluorescent protein.     

Stepwise formation of alpha-helices during cytochrome c folding

Two models have been proposed to describe the folding pathways of proteins. The framework model assumes the initial formation of the secondary structures whereas the hydrophobic collapse model supposes their formation after the collapse of backbone structures. To differentiate between these models for real proteins, we have developed a novel CD spectrometer that enables us to observe the submillisecond time frame of protein folding and have characterized the timing of secondary structure formation in the folding process of cytochrome c (cyt c). We found that approximately 20% of the native helical content was organized in the first phase of folding, which is completed within milliseconds. Furthermore, we suggest the presence of a second intermediate, which has alpha-helical content resembling that of the molten globule state. Our results indicate that many of the alpha-helices are organized after collapse in the folding mechanism of cyt c.     

Early events in protein folding

Recent advances have significantly increased the time and spectroscopic resolution of protein folding experiments. We can now study the timescale and nature of polypeptide collapse, and how this correlates with secondary and tertiary structure formation. Studies on ultrafast folding proteins and peptides provide experimental benchmarks on a timescale that overlaps directly with that of molecular dynamics simulations. This makes possible direct tests of both simulations and current models of protein folding.     

Folding path of P5abc RNA involves direct coupling of secondary and tertiary structures

Folding mechanisms in which secondary structures are stabilized through the formation of tertiary interactions are well documented in protein folding but challenge the folding hierarchy normally assumed for RNA. However, it is increasingly clear that RNA could fold by a similar mechanism. P5abc, a small independently folding tertiary domain of the Tetrahymena thermophila group I ribozyme, is known to fold by a secondary structure rearrangement involving helix P5c. However, the extent of this rearrangement and the precise stage of folding that triggers it are unknown. We use experiments and simulations to show that the P5c helix switches to the native secondary structure late in the folding pathway and is directly coupled to the formation of tertiary interactions in the A-rich bulge. P5c mutations show that the switch in P5c is not rate-determining and suggest that non-native interactions in P5c aid folding rather than impede it. Our study illustrates that despite significant differences in the building blocks of proteins and RNA, there may be common ways in which they self-assemble.     

Mimicking the action of folding chaperones in molecular dynamics simulations: Application to the refinement of homology-based protein structures

A novel method for the refinement of misfolded protein structures is proposed in which the properties of the solvent environment are oscillated in order to mimic some aspects of the role of molecular chaperones play in protein folding in vivo. Specifically, the hydrophobicity of the solvent is cycled by repetitively altering the partial charges on solvent molecules (water) during a molecular dynamics simulation. During periods when the hydrophobicity of the solvent is increased, intramolecular hydrogen bonding and secondary structure formation are promoted. During periods of increased solvent polarity, poorly packed regions of secondary structures are destabilized, promoting structural rearrangement. By cycling between these two extremes, the aim is to minimize the formation of long-lived intermediates. The approach has been applied to the refinement of structural models of three proteins generated by using the ROSETTA procedure for ab initio structure prediction. A significant improvement in the deviation of the model structures from the corresponding experimental structures was observed. Although preliminary, the results indicate computationally mimicking some functions of molecular chaperones in molecular dynamics simulations can promote the correct formation of secondary structure and thus be of general use in protein folding simulations and in the refinement of structural models of small- to medium-size proteins.     

The unfolded state of the villin headpiece helical subdomain: computational studies of the role of locally stabilized structure

The 36 residue villin headpiece helical subdomain (HP36) is one of the fastest cooperatively folding proteins, folding on the microsecond timescale. HP36's simple three helix topology, fast folding and small size have made it an attractive model system for computational and experimental studies of protein folding. Recent experimental studies have explored the denatured state of HP36 using fragment analysis coupled with relatively low-resolution spectroscopic techniques. These studies have shown that there is apparently only a small tendency to form locally stabilized secondary structure. Here, we complement the experimental studies by using replica exchange molecular dynamics with explicit solvent to investigate the structural features of these peptide models of unfolded HP36. To ensure convergence, two sets of simulations for each fragment were performed with different initial structures, and simulations were continued until these generated very similar final ensembles. These simulations reveal low populations of native-like structure and early folding events that cannot be resolved by experiment. For each fragment, calculated J-coupling constants and helical propensities are in good agreement with experimental trends. HP-1, corresponding to residues 41 to 53 and including the first alpha-helix, contains the highest helical population. HP-3, corresponding to residues 62 through 75 and including the third alpha-helix, contains a small population of helical turn residing at the N terminus while HP-2, corresponding to residues 52 through 61 and including the second alpha-helix, formed little to no structure in isolation. Overall, HP-1 was the only fragment to adopt a native-like conformation, but the low population suggests that formation of significant structure only occurs after formation of specific tertiary interactions.     

Analysis of the early stage of the folding process of reduced lysozyme using all lysozyme variants containing a pair of cysteines

Twenty-eight hen lysozyme variants that contained a pair of cysteines were constructed to examine the formation of the individual native and nonnative disulfide bonds. We analyzed the extent of the formation of a disulfide bond in each lysozyme variant using a redox buffer (pH 8) containing 1.0 mM reduced and 0.1 mM oxidized glutathione in the absence or presence of 6 M guanidine hydrochloride. In the presence of 6 M guanidine hydrochloride, the extent of the formation of the disulfide bond in each lysozyme variant was proportional to the distance between cysteine residues, indicating that reduced hen lysozyme under a highly denaturing condition adopted a randomly coiled structure. In aqueous solution, the formations of all disulfide bonds occurred much more easily than under a denatured condition. This finding indicated that reduced lysozyme had a somewhat compact structure. Moreover, the scattering data for the extents of the formation of the disulfide bonds among all lysozyme variants were observed. These results suggested that the nonrandom folding occurred in the early stage of the folding of reduced lysozyme, which should provide new insight into the early-stage events in the folding process of reduced lysozyme.     

Gatekeepers in the ribosomal protein s6: thermodynamics, kinetics, and folding pathways revealed by a minimalist protein model

We investigate the effect of structural gatekeepers on the folding of the ribosomal protein S6. Folding thermodynamics and early refolding kinetics are studied for this system utilizing computer simulations of a minimalist protein model. When gatekeepers are eliminated, the thermodynamic signature of a folding intermediate emerges, and a marked decrease in folding efficiency is observed. We explain the prerequisites that determine the "strength" of a given gatekeeper. The investigated gatekeepers are found to have distinct functions, and to guide the folding and time-dependent packing of non-overlapping secondary structure elements in the protein. Gatekeepers avoid kinetic traps during folding by favoring the formation of "productive topologies" on the way to the native state. The trends in folding rates in the presence/absence of gatekeepers observed for our minimalist model of S6 are in very good agreement with experimental data on this protein.     

Protein conformational transitions coupled to binding in molecular recognition of unstructured proteins: deciphering the effect of intermolecular interactions on computational structure prediction of the p27Kip1 protein bound to the cyclin A-cyclin-dependent kinase 2 complex

The relationship between folding mechanism coupled to binding and structure prediction of the tertiary complexes is studied for the p27(Kip) (1) protein which has an intrinsically disordered unbound form and undergoes a functional folding transition during complex formation with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) binary complex. Hierarchy of p27(Kip1) structural loss determined in our earlier studies from temperature-induced Monte Carlo simulations and subsequent characterization of the transition state ensemble (TSE) for the folding reaction have shown that simultaneous ordering of the p27(Kip1) native intermolecular interface for the beta-hairpin and beta-strand secondary structure elements is critical for nucleating a rapid kinetic transition to the native tertiary complex. In the present study, we investigate the effect of forming specific intermolecular interactions on structure prediction of the p27(Kip1) tertiary complex. By constraining different secondary structure elements of p27(Kip1) in their native bound conformations and conducting multiple simulated annealing simulations, we analyze differences in the success rate of predicting the native structure of p27(Kip1) in the tertiary complex. In accordance with the nucleation-condensation mechanism, we have found that further stabilization of the native intermolecular interface for the beta-hairpin and beta-strand elements of p27(Kip1), that become ordered in the TSE, but are hardly populated in the unbound state, results in a consistent acquisition of the native bound structure. Conversely, the excessive stablization of the local secondary structure elements, which are rarely detected in the TSE, has a detrimental effect on convergence to the native bound structure.     

The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.     

Lysozyme folded in silico according to the limited conformational sub-space

The conformational sub-space oriented on early-stage protein folding is applied to lysozyme folding. The part of the Ramachandran map distinguished on the basis of a geometrical model of the polypeptide chain limited to the mutual orientation of the peptide bond planes is shown to deliver the initial structure of the polypeptide for the energy minimization procedure in the ab initio model of protein folding prediction. Two forms of energy minimization and molecular dynamics simulation procedures were applied to the assumed early-stage protein folding of lysozyme. One of them included the disulphide bond system and the other excluded it. The post-energy-minimization and post-dynamics structures were compared using RMS-D and non-bonding contact maps to estimate the degree of approach to the native, target structure of the protein molecule obtained using the limited conformational sub-space for the early stage of folding.     

Independent nucleation and heterogeneous assembly of structure during folding of equine lysozyme.

The refolding of equine lysozyme from guanidinium chloride has been studied using hydrogen exchange pulse labelling in conjunction with NMR spectroscopy and stopped flow optical methods. The stopped flow optical experiments indicate that extensive hydrophobic collapse occurs rapidly after the initiation of refolding. Pulse labelling experiments monitoring nearly 50 sites within the protein have enabled the subsequent formation of native-like structure to be followed in considerable detail. They reveal that an intermediate having persistent structure within three of the four helices of the alpha-domain of the protein is formed for the whole population of molecules within 4 ms. Subsequent to this event, however, the hydrogen exchange protection kinetics are complex and highly heterogeneous. Analysis of the results by fitting to stretched exponential functions shows that a series of other intermediates is formed as a consequence of the stepwise assembly of independently nucleated local regions of structure. In some molecules the next step in folding involves the stabilisation of the remaining helix in the alpha-domain, whilst in others persistent structure begins to form in the beta-domain. The formation of native-like structure throughout the beta-domain is itself heterogeneous, involving at least three kinetically distinguishable steps. Residues in loop regions throughout the protein attain persistent structure more slowly than regions of secondary structure. There is in addition evidence for locally misfolded regions of structure that reorganise on much longer timescales. The results reveal that the native state of the protein is generated by the heterogeneous assembly of a series of locally cooperative regions of structure. This observation has many features in common with the findings of recent theoretical simulations of protein folding.     

Fast in silico protein folding by introduction of alternating hydrogen bond potentials

We accelerate protein folding in all-atom molecular dynamics simulations by introducing alternating hydrogen bond potentials as a supplement to the force field. The alternating hydrogen bond potentials result in accelerated hydrogen bond reordering, which leads to rapid formation of secondary structure elements. The method does not require knowledge of the native state but generates the potentials based on the development of the tertiary structure in the simulation. In protein folding, the formation of secondary structure elements, especially α-helix and β-sheet, is very important, and we show that our method can fold both efficiently and with great speed.