Evidence for Initial Non-specific Polypeptide Chain Collapse During the Refolding of the SH3 Domain of PI3 Kinase

Refolding of the SH3 domain of PI3 kinase from the guanidine hydrochloride (GdnHCl)-unfolded state has been probed with millisecond (stopped flow) and sub-millisecond (continuous flow) measurements of the change in fluorescence, circular dichroism, ANS fluorescence and three-site fluorescence resonance energy transfer (FRET) efficiency. Fluorescence measurements are unable to detect structural changes preceding the rate-limiting step of folding, whereas measurements of changes in ANS fluorescence and FRET efficiency indicate that polypeptide chain collapse precedes the major structural transition. The initial chain collapse reaction is complete within 150 μs. The collapsed form at this time possesses hydrophobic clusters to which ANS binds. Each of the three measured intra-molecular distances has contracted to an extent predicted by the dependence of the FRET signal in completely unfolded protein on denaturant concentration, indicating that contraction is non-specific. The extent of contraction of each intra-molecular distance in the collapsed product of sub-millisecond folding increases continuously with a decrease in [GdnHCl]. The gradual contraction is continuous with the gradual contraction seen in completely unfolded protein, and its dependence on [GdnHCl] is not indicative of an all-or-none collapse reaction. The dependence of the extent of contraction on [GdnHCl] was similar for the three distances, indicating that chain collapse occurs in a synchronous manner across different segments of the polypeptide chain. The sub-millisecond measurements of folding in GdnHCl were unable to determine whether hydrophobic cluster formation, probed by ANS fluorescence measurement, precedes chain contraction probed by FRET. To determine whether hydrogen bonding plays a role in initial chain collapse, folding was initiated by dilution of the urea-unfolded state. The extent of contraction of at least one intra-molecular distance in the collapsed product of sub-millisecond folding in urea is similar to that seen in GdnHCl, and the initial contraction in urea too appears to be gradual.     

Microsecond hydrophobic collapse in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein

Using small-angle X-ray scattering combined with a continuous-flow mixing device, we monitored the microsecond compaction dynamics in the folding of Escherichia coli dihydrofolate reductase, an alpha/beta-type protein. A significant collapse of the radius of gyration from 30 A to 23.2 A occurs within 300 micros after the initiation of refolding by a urea dilution jump. The subsequent folding after the major chain collapse occurs on a considerably longer time-scale. The protein folding trajectories constructed by comparing the development of the compactness and the secondary structure suggest that the specific hydrophobic collapse model rather than the framework model better explains the experimental observations. The folding trajectory of this alpha/beta-type protein is located between those of alpha-helical and beta-sheet proteins, suggesting that native structure determines the folding landscape.     

Pressure-jump small-angle x-ray scattering detected kinetics of staphylococcal nuclease folding

The kinetics of chain disruption and collapse of staphylococcal nuclease after positive or negative pressure jumps was monitored by real-time small-angle x-ray scattering under pressure. We used this method to probe the overall conformation of the protein by measuring its radius of gyration and pair-distance-distribution function p(r) which are sensitive to the spatial extent and shape of the particle. At all pressures and temperatures tested, the relaxation profiles were well described by a single exponential function. No fast collapse was observed, indicating that the rate limiting step for chain collapse is the same as that for secondary and tertiary structure formation. Whereas refolding at low pressures occurred in a few seconds, at high pressures the relaxation was quite slow, approximately 1 h, due to a large positive activation volume for the rate-limiting step for chain collapse. A large increase in the system volume upon folding implies significant dehydration of the transition state and a high degree of similarity in terms of the packing density between the native and transition states in this system. This study of the time-dependence of the tertiary structure in pressure-induced folding/unfolding reactions demonstrates that novel information about the nature of protein folding transitions and transition states can be obtained from a combination of small-angle x-ray scattering using high intensity synchrotron radiation with the high pressure perturbation technique.     

Population of on-pathway intermediates in the folding of ubiquitin

The role that intermediate states play in protein folding is the subject of intense investigation and in the case of ubiquitin has been controversial. We present fluorescence-detected kinetic data derived from single and double mixing stopped-flow experiments to show that the F45W mutant of ubiquitin (WT*), a well-studied single-domain protein and most recently regarded as a simple two-state system, folds via on-pathway intermediates. To account for the discrepancy we observe between equilibrium and kinetic stabilities and m-values, we show that the polypeptide chain undergoes rapid collapse to an intermediate whose presence we infer from a fast lag phase in interrupted refolding experiments. Double-jump kinetic experiments identify two direct folding phases that are not associated with slow isomerisation reactions in the unfolded state. These two phases are explained by kinetic partitioning which allows molecules to reach the native state from the collapsed state via two possible competing routes, which we further examine using two destabilised ubiquitin mutants. Interrupted refolding experiments allow us to observe the formation and decay of an intermediate along one of these pathways. A plausible model for the folding pathway of ubiquitin is presented that demonstrates that obligatory intermediates and/or chain collapse are important events in restricting the conformational search for the native state of ubiquitin.     

Simulation of the thermal denaturation of hen egg white lysozyme: trapping the molten globule state.

In the study of protein folding, much attention has focused on the characterization of folding intermediates. We report here molecular dynamics simulations in which the initial stages of the thermal denaturation of hen egg white lysozyme in aqueous solution are examined in detail. It is found that lysozyme unfolds in a two-stage process with the initial formation a quasi-stable state in which significant rearrangement of the secondary structure takes place. No evidence for distinct folding domains was found. The simulations suggest that the formation of well-defined secondary structure occurs after the initial collapse of the peptide chain and thus tend against the framework model of protein folding.     

Observation of Continuous Contraction and a Metastable Misfolded State during the Collapse and Folding of a Small Protein

To obtain proper insight into how structure develops during a protein folding reaction, it is necessary to understand the nature and mechanism of the polypeptide chain collapse reaction, which marks the initiation of folding. Here, the time-resolved fluorescence resonance energy transfer technique, in which the decay of the fluorescence light intensity with time is used to determine the time evolution of the distribution of intra-molecular distances, has been utilized to study the folding of the small protein, monellin. It is seen that when folding begins, about one-third of the protein molecules collapse into a molten globule state (I_MG), from which they relax by continuous further contraction to transit to the native state. The larger fraction gets trapped into a metastable misfolded state. Exit from this metastable state occurs via collapse to the lower free energy I_MG state. This exit is slow, on a time scale of seconds, because of activation energy barriers. The trapped misfolded molecules as well as the I_MG molecules contract continuously and slowly as structure develops. A phenomenological model of Markovian evolution of the polymer chain undergoing folding, incorporating these features, has been developed, which fits well the experimentally observed time evolution of distance distributions. The observation that the “wrong turn” to the misfolded state has not been eliminated by evolution belies the common belief that the folding of functional protein sequences is very different from that of a random heteropolymer, and the former have been selected by evolution to fold quickly.     

Relationship between chain collapse and secondary structure formation in a partially folded protein

Chain collapse and secondary structure formation are frequently observed during the early stages of protein folding. Is the chain collapse brought about by interactions between secondary structure units or is it due to polymer behavior in a poor solvent (coil‐globule transition)? To answer this question, we measured small‐angle X‐ray scattering for a series of β‐lactoglobulin mutants under conditions in which they assume a partially folded state analogous to the folding intermediates. Mutants that were designed to disrupt the secondary structure units showed the gyration radii similar to that of the wild type protein, indicating that chain collapse is due to coil‐globule transitions. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 651–658, 2014.     

Dissecting the non-specific and specific components of the initial folding reaction of barstar by multi-site FRET measurements

Initial polypeptide chain collapse plays a major role in the development of subsequent structure during protein folding, but it has been difficult to elucidate the coupling between its cooperativity and specificity. To better understand this important aspect of protein folding, nine different intramolecular distances in the protein have been measured by fluorescence resonance energy transfer (FRET) in the product(s) of the initial, sub-millisecond collapse reaction during the folding of barstar, under different folding conditions. All nine distances contract in these initial folding products, when the denaturant concentration is reduced. Two of these distances were also measured in peptides corresponding to sequence segments 38-55 and 51-69 of the protein. Surprisingly, both distances do not contract in the peptides which remain fully unfolded when the denaturant concentration is reduced. This suggests that the contraction of at least some segments of the polypeptide chain may be facilitated only by contraction of other segments. In the case of the initial product of folding of the protein, the dependence on denaturant concentration of the relative change in each distance suggests that there are two components to the initial folding reaction. One is a nonspecific component, which appears to be driven by the change in denaturant concentration that is used to initiate refolding. This component corresponds to the collapse of completely unfolded protein (U) to unfolded protein in refolding conditions (U(C)). The extent of nonspecific collapse can be predicted by the response of completely unfolded protein to a change in denaturant concentration. All distances undergo such solvent-induced contraction, but each distance contracts to a different extent. There is also a specific component to initial sub-millisecond folding, in which some distances (but not all) contract more than that predicted by solvent-induced contraction. The observation that only some of the distances undergo contraction over and above solvent-induced contraction, suggest that this specific component is associated with the formation of a specific intermediate (I(E)). FRET efficiency and distance change differently for the different donor-acceptor pairs, with a change in denaturant concentration, indicating that the formation or dissolution of structure in U(C) and I(E) does not happen in a synchronized manner across different regions of the protein molecule. Also, all nine FRET efficiencies and intramolecular distances in the product(s) of sub-ms folding, change continuously with a change in denaturant concentration. Hence, it appears that the transitions from U to U(C) and to I(E) are gradual transformations, and not all-or-none structural transitions. Nevertheless, the product of these gradual transitions, I(E), possesses specific structure.     

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.     

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.     

Principles of protein folding--a perspective from simple exact models.

General principles of protein structure, stability, and folding kinetics have recently been explored in computer simulations of simple exact lattice models. These models represent protein chains at a rudimentary level, but they involve few parameters, approximations, or implicit biases, and they allow complete explorations of conformational and sequence spaces. Such simulations have resulted in testable predictions that are sometimes unanticipated: The folding code is mainly binary and delocalized throughout the amino acid sequence. The secondary and tertiary structures of a protein are specified mainly by the sequence of polar and nonpolar monomers. More specific interactions may refine the structure, rather than dominate the folding code. Simple exact models can account for the properties that characterize protein folding: two-state cooperativity, secondary and tertiary structures, and multistage folding kinetics--fast hydrophobic collapse followed by slower annealing. These studies suggest the possibility of creating "foldable" chain molecules other than proteins. The encoding of a unique compact chain conformation may not require amino acids; it may require only the ability to synthesize specific monomer sequences in which at least one monomer type is solvent-averse.     

Protein folding and chain collapse

It is currently assumed that the folded structure of a globular protein is controlled in a highly deterministic way by its amino acid sequence. We show here that a very different although not necessarily contrasting viewpoint can be adopted. From statistical treatment of x-ray results, we suggest that the folding pattern essentially follows the collapse to be expected on statistical-mechanical grounds for an ideal chain effectively experiencing self-attraction and comprising identical units, whose conformational properties are obtained as an average over the actual amino acid units. The local details of folding of each protein, obviously dictated byits amino acid sequence, can be regarded as statistical fluctuations. We consider 31 globular fragments belonging to 21 different water-soluble, nonmembrane proteins. By the theory of chain collapse proposed by two of us [G. Allegra and F. Ganazzoli (1985) J. Chem. Phys. 83 , 397], all the average intramolecular distances may be obtained. Accordingly, first an average plot of the mean-square distances between kth neighboring amino acid units is constructed, starting from the observed crystallographic coordinates. Then the plot is basically reproduced with a wormlike chain model undergoing collapse as a result of intramolecular attractive forces. Agreement is especially good for short amino acid sequences (k ? 30), in which case the statistical sampling is more accurate, enabling us to determine the model parameters. The resulting mean-square radius of gyration is also in good agreement with the experimental average, whereas the unperturbed characteristic ratio is roughly consistent with results from conformational calculations by W. L. Mattice [(1977) Macromolecules 10 , 516], based on the rotational isomeric state approach.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

Early Closure of a Long Loop in the Refolding of Adenylate Kinase: A Possible Key Role of Non-Local Interactions in the Initial Folding Steps

Most globular protein chains, when transferred from high to low denaturant concentrations, collapse instantly before they refold to their native state. The initial compaction of the protein molecule is assumed to have a key effect on the folding pathway, but it is not known whether the earliest structures formed during or instantly after collapse are defined by local or by non-local interactions—that is, by secondary structural elements or by loop closure of long segments of the protein chain. Stable closure of one or several long loops can reduce the chain entropy at a very early stage and can prevent the protein from following non-productive pathways whose number grows exponentially with the length of the protein chain. In Escherichia coli adenylate kinase (AK), about seven long loops define the topology of the native structure. We selected four loop-forming sections of the chain and probed the time course of loop formation during refolding of AK. We labeled the termini of the loop segments with tryptophan and cysteine-5-amidosalicylic acid. This donor-acceptor pair of probes used with fluorescence resonance excitation energy transfer spectroscopy (FRET) is suitable for detecting very short distances and thus is able to distinguish between random and specific compactions. Refolding of AK was initiated by stopped-flow mixing, followed simultaneously by donor and acceptor fluorescence, and analyzed in terms of energy transfer efficiency and distance. In the collapsed state of AK, observed after the 5-ms dead time of the instrument, one of the selected segments shows a native-like separation of its termini; it forms a loop already in the collapsed state. A second segment that includes the first but is longer by 15 residues shows an almost native-like separation of its termini. In contrast, a segment that is shorter but part of the second segment shows a distance separation of its termini as high as a segment that spans almost the whole protein chain. We conclude that a specific network of non-local interactions, the closure of one or several loops, can play an important role in determining the protein folding pathway at its early phases.     

Ubiquitin folds via a flip-twist-lock mechanism

To perform specific functional activities, the majority of proteins should fold into their distinct three-dimensional conformations. However, the biologically active conformation of a protein is generally found to be marginally stable than the other conformations that the chain can adopt. How a protein finds its native conformation from its post-synthesis unfolded structure in a complex conformational landscape is the unsolved question that still drives the protein folding community. Here, we report the folding mechanism of a globular protein, ubiquitin, from its chemically denatured state using all-atom molecular dynamics simulations. From the kinetic analysis of the simulated trajectories we show that the folding process can be described by the hydrophobic collapse mechanism, initiated by the “dewetting transition”, and subsequently assisted by the origination of an N-terminal folding nucleus, and finally supported by a native salt-bridge interaction between K11 and E34. We show that ubiquitin folds via an intermediate. Finally, we confirm the presence of “biological water” and explain its role to the folding process.     

An all-atom structure-based potential for proteins: bridging minimal models with all-atom empirical forcefields

Protein dynamics take place on many time and length scales. Coarse-grained structure-based (Go) models utilize the funneled energy landscape theory of protein folding to provide an understanding of both long time and long length scale dynamics. All-atom empirical forcefields with explicit solvent can elucidate our understanding of short time dynamics with high energetic and structural resolution. Thus, structure-based models with atomic details included can be used to bridge our understanding between these two approaches. We report on the robustness of folding mechanisms in one such all-atom model. Results for the B domain of Protein A, the SH3 domain of C-Src Kinase, and Chymotrypsin Inhibitor 2 are reported. The interplay between side chain packing and backbone folding is explored. We also compare this model to a C(alpha) structure-based model and an all-atom empirical forcefield. Key findings include: (1) backbone collapse is accompanied by partial side chain packing in a cooperative transition and residual side chain packing occurs gradually with decreasing temperature, (2) folding mechanisms are robust to variations of the energetic parameters, (3) protein folding free-energy barriers can be manipulated through parametric modifications, (4) the global folding mechanisms in a C(alpha) model and the all-atom model agree, although differences can be attributed to energetic heterogeneity in the all-atom model, and (5) proline residues have significant effects on folding mechanisms, independent of isomerization effects. Because this structure-based model has atomic resolution, this work lays the foundation for future studies to probe the contributions of specific energetic factors on protein folding and function.     

Monte Carlo simulation of protein folding with orientation-dependent monomer–monomer interactions

We present the results of lattice Monte Carlo simulations of protein folding in the framework of a model taking into account (i) the dependence of the energy of interaction of amino-acid residues on their orientation and (ii) the rigidity of the polypeptide chain with respect to the formation of kinks. If the chain is flexible, the final protein structures are predicted to be compact. Increasing the energy cost of creation of kinks is found to favor the formation of flat structures mimicking an ideal antiparallel β sheet. For compact structures, the kinetics of folding exhibit the standard two-phase regime (a rapid collapse to one of the metastable state, followed by slow reconfiguration of the chain to the native structure). For flat structures, the transition to the native state is often gradual. Proteins 29:508–516, 1997. © 1997 Wiley-Liss, Inc.     

Commitment and nucleation in the protein G transition state

An accurate characterization of the transition state ensemble (TSE) is central to furthering our understanding of the protein folding reaction. We have extensively tested a recently reported method for studying a protein's TSE, utilizing phi-value data from protein engineering experiments and computational studies as restraints in all-atom Monte Carlo (MC) simulations. The validity of interpreting experimental phi-values as the fraction of native contacts made by a residue in the TSE was explored, revealing that this definition is unable to uniquely specify a TSE. The identification of protein G's second hairpin, in both pre and post-transition conformations demonstrates that high experimental phi-values do not guarantee a residue's importance in the TSE. An analysis of simulations based on structures restrained by experimental phi-values is necessary to yield this result, which is not obvious from a simplistic interpretation of individual phi-values. The TSE that we obtain corresponds to a single, specific nucleation event, characterized by six residues common to all three observed, convergent folding pathways. The same specific nucleus was independently identified from computational and experimental data, and "Conservation of Conservation" analysis in the protein G fold. When associated strictly with complete nucleus formation and concomitant chain collapse, folding is a well-defined two state event. Once the nucleus has formed, the folding reaction enters a slow relaxation process associated with side-chain packing and small, local backbone rearrangements. A detailed analysis of phi-values and their relationship to the transition state ensemble allows us to construct a unified theoretical model of protein G folding.     

Kinetics of cytochrome C folding: atomically detailed simulations

The vast range of time scales (from nanoseconds to seconds) during protein folding is a challenge for experiments and computations. To make concrete predictions on folding mechanisms, atomically detailed simulations of protein folding, using potentials derived from chemical physics principles, are desired. However, due to their computational complexity, straightforward molecular dynamics simulations of protein folding are impossible today. An alternative algorithm is used that makes it possible to compute approximate atomically detailed long time trajectories (the Stochastic Difference Equation in Length). This algorithm is used to compute 26 atomically detailed folding trajectories of cytochrome c (a millisecond process). The early collapse of the protein chain (with marginal formation of secondary structure), and the earlier formation of the N and C helices (compare to the 60's helix) are consistent with the experiment. The existence of an energy barrier upon entry to the molten globule is examined as well. In addition to (favorable) comparison to experiments, we show that non-native contacts drive the formation of the molten globule. In contrast to popular folding models, the non-native contacts do not form off-pathway kinetic traps in cytochrome c.     

BiP and PDI cooperate in the oxidative folding of antibodies in vitro

Immunoglobulin heavy chain binding protein (BiP), a member of the Hsp70 chaperone family, and the oxidoreductase protein-disulfide isomerase (PDI) play an important role in the folding and oxidation of proteins in the endoplasmic reticulum. However, it was not clear whether both cooperate in this process. We show here that BiP and PDI act synergistically in the in vitro folding of the denatured and reduced Fab fragment. Several ATP-dependent cycles of binding, release, and rebinding of the unfolded antibody chains by BiP are required for efficient reactivation. Our data suggest that in the absence of BiP unfolded antibody chains collapse rapidly upon refolding, rendering cysteine side chains inaccessible for PDI. BiP binds the unfolded polypeptide chains and keeps them in a conformation in which the cysteine residues are accessible for PDI. These findings support the idea of a network of folding helper proteins in the endoplasmic reticulum, which makes this organelle a dedicated protein-processing compartment.