Perturbations of the denatured state ensemble: modeling their effects on protein stability and folding kinetics.

By considering the denatured state of a protein as an ensemble of conformations with varying numbers of sequence-specific interactions, the effects on stability, folding kinetics, and aggregation of perturbing these interactions can be predicted from changes in the molecular partition function. From general considerations, the following conclusions are drawn: (1) A perturbation that enhances a native interaction in denatured state conformations always increases the stability of the native state. (2) A perturbation that promotes a non-native interaction in the denatured state always decreases the stability of the native state. (3) A change in the denatured state ensemble can alter the kinetics of aggregation and folding. (4) The loss (or increase) in stability accompanying two mutations, each of which lowers (or raises) the free energy of the denatured state, will be less than the sum of the effects of the single mutations, except in cases where both mutations affect the same set of partially folded conformations. By modeling the denatured state as the ensemble of all non-native conformations of hydrophobic-polar (HP) chains configured on a square lattice, it can be shown that the stabilization obtained from enhancement of native interactions derives in large measure from the avoidance of non-native interactions in the D state. In addition, the kinetic effects of fixing single native contacts in the denatured state or imposing linear gradients in the HH contact probabilities are found, for some sequences, to significantly enhance the efficiency of folding by a simple hydrophobic zippering algorithm. Again, the dominant mechanism appears to be avoidance of non-native interactions. These results suggest stabilization of native interactions and imposition of gradients in the stability of local structure are two plausible mechanisms involving the denatured state that could play a role in the evolution of protein folding and stability.     

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.     

Molecular basis of cooperativity in protein folding. V. Thermodynamic and structural conditions for the stabilization of compact denatured states

The heat-denatured state of proteins has been usually assumed to be a fully hydrated random coil. It is now evident that under certain solvent conditions or after chemical or genetic modifications, the protein molecule may exhibit a hydrophobic core and residual secondary structure after thermal denaturation. This state of the protein has been called the “compact denatured” or “molten globule” state. Recently is has been shown that α-lactalbumin at pH < 5 denatures into a molten globule state upon increasing the temperature (Griko, Y., Freire, E., Privalov, P. L. Biochemistry 33:1889–1899, 1994). This state has a lower heat capacity and a higher enthalpy at low temperatures than the unfolded state. At those temperatures the stabilization of the molten globule state is of an entropic origin since the enthalpy contributes unfavorably to the Gibbs free energy. Since the molten globule is more structured than the unfolded state and, therefore, is expected to have a lower configurational entropy, the net entropic gain must originate primarily from solvent related entropy arising from the hydrophobic effect, and to a lesser extent from protonation or electrostatic effects. In this work, we have examined a large ensemble of partly folded states derived from the native structure of α-lactalbumin in order to identify those states that satisfy the energetic criteria of the molten globule. It was found that only few states satisfied the experimental constraints and that, furthermore, those states were part of the same structural family. In particular, the regions corresponding to the A, B, and C helices were found to be folded, while the β sheet and the D helix were found to be unfolded. At temperatures below 45°C the states exhibiting those structural characteristics are enthalpically higher than the unfolded state in agreement with the experimental data. Interestingly, those states have a heat capacity close to that observed for the acid pH compact denatured state of α-lactalbumin [980 cal (mol.K)^?l]. In addition, the folded regions of these states include those residues found to be highly protected by NMR hydrogen exchange experiments. This work represents an initial attempt to model the structural origin of the thermodynamic properties of partly folded states. The results suggest a number of structural features that are consistent with experimental data. © 1994 Wiley-Liss, Inc.     

Charge-charge interactions in the denatured state influence the folding kinetics of ribonuclease Sa

Gaining a better understanding of the denatured state ensemble of proteins is important for understanding protein stability and the mechanism of protein folding. We studied the folding kinetics of ribonuclease Sa (RNase Sa) and a charge-reversal variant (D17R). The refolding kinetics are similar, but the unfolding rate constant is 10-fold greater for the variant. This suggests that charge-charge interactions in the denatured state and the transition state ensembles are more favorable in the variant than in RNase Sa, and shows that charge-charge interactions can influence the kinetics and mechanism of protein folding.     

Denatured collapsed states in protein folding: example of apomyoglobin.

Experimental approaches, including circular dichroism, small angle X-ray scattering, steady-state fluorescence, and fluorescence energy transfer, were applied to study the 3D-structure of apomyolgobin in different conformational states. These included the native and molten globules, along with either less ordered conformations induced by the addition of anions or completely unfolded states. The results show that the partially folded forms of apomyoglobin stabilized by KCl and/or Na(2)SO(4) under unfolding conditions (pH 2) exhibit a significant amount of secondary structure (circular dichroism), low packing density of protein molecules (SAXS), and native-like dimensions of the AGH core (fluorescence energy transfer). This finding indicates that a native-like tertiary fold of the polypeptide chain, i.e., the spatial organization of secondary structure elements, most likely emerges prior to the formation of the molten globule state.     

Thermodynamics and kinetics of non-native interactions in protein folding: a single point mutant significantly stabilizes the N-terminal domain of L9 by modulating non-native interactions in the denatured state

Comparatively little is known about the role of non-native interactions in protein folding and their role in both folding and stability is controversial. We demonstrate that non-native electrostatic interactions involving specific residues in the denatured state can have a significant effect upon protein stability and can persist in the transition state for folding. Mutation of a single surface exposed residue, Lys12 to Met, in the N-terminal domain of the ribosomal protein L9 (NTL9), significantly increased the stability of the protein and led to faster folding. Structural and energetic studies of the wild-type and K12M mutant show that the 1.9 kcal mol(-1) increase in stability is not due to native state effects, but rather is caused by modulation of specific non-native electrostatic interactions in the denatured state. pH dependent stability measurements confirm that the increased stability of the K12M is due to the elimination of favorable non-native interactions in the denatured state. Kinetic studies show that the non-native electrostatic interactions involving K12 persist in the transition state. The analysis demonstrates that canonical Phi-values can arise from the disruption of non-native interactions as well as from the development of native interactions.     

Effect of denaturant and protein concentrations upon protein refolding and aggregation: a simple lattice model.

We present a study of the competition between protein refolding and aggregation for simple lattice model proteins. The effect of solvent conditions (i.e., the denaturant concentration and the protein concentration) on the folding and aggregation behavior of a system of simple, two-dimensional lattice protein molecules has been investigated via (dynamic Monte Carlo simulations. The population profiles and aggregation propensities of the nine most populated intermediate configurations exhibit a complex dependence on the solution conditions that can be understood by considering the competition between intra- and interchain interactions. Some of these configurations are not even seen in isolated chain simulations; they are observed to be highly aggregation prone and are stabilized primarily by the aggregation reaction in multiple-chain systems. Aggregation arises from the association of partially folded intermediates rather than from the association of denatured random-coil states. The aggregation reaction dominates over the folding reaction at high protein concentration and low denaturant concentration, resulting in low refolding yields at those conditions. However, optimum folding conditions exist at which the refolding yield is a maximum, in agreement with some experimental observations.     

The protein folding network

The conformation space of a 20 residue antiparallel beta-sheet peptide, sampled by molecular dynamics simulations, is mapped to a network. Snapshots saved along the trajectory are grouped according to secondary structure into nodes of the network and the transitions between them are links. The conformation space network describes the significant free energy minima and their dynamic connectivity without requiring arbitrarily chosen reaction coordinates. As previously found for the Internet and the World-Wide Web as well as for social and biological networks, the conformation space network is scale-free and contains highly connected hubs like the native state which is the most populated free energy basin. Furthermore, the native basin exhibits a hierarchical organization, which is not found for a random heteropolymer lacking a predominant free-energy minimum. The network topology is used to identify conformations in the folding transition state (TS) ensemble, and provides a basis for understanding the heterogeneity of the TS and denatured state ensemble as well as the existence of multiple pathways.     

Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding

Do stable intermediates form very early in the protein folding process? New results and a quantity of literature that bear on this issue are examined here. Results available provide little support for early intermediate accumulation before an initial search-dependent nucleation barrier.     

Hydrophobic clustering in nonnative states of a protein: interpretation of chemical shifts in NMR spectra of denatured states of lysozyme.

Chemical shifts of resonances of specific protons in the 1H NMR spectrum of thermally denatured hen lysozyme have been determined by exchange correlation with assigned native state resonances in 2D NOESY spectra obtained under conditions where the two states are interconverting. There are subtle but widespread deviations of the measured shifts from the values which would be anticipated for a random coil; in the case of side chain protons these are virtually all net upfield shifts and it is shown that this may be the averaged effect of interactions with aromatic rings in a partially collapsed denatured state. In a very few cases, notably that of two sequential tryptophan residues, it is possible to interpret these effects in terms of specific, local interresidue interactions. Generally, however, there is no correlation with either native state shift perturbations or with sequence proximity to aromatic groups. Diminution of most of the residual shift perturbations on reduction of the disulfide cross-links confirms that they are not simply effects of residues adjacent in the sequence. Similar effects of chemical denaturants, with the disulfides intact, demonstrate that the shift perturbations reflect an enhanced tendency to side chain clustering in the thermally denatured state. The temperature dependences of the shift perturbations suggest that this clustering is noncooperative and is driven by small, favorable enthalpy changes. While the extent of conformational averaging is clearly much greater than that observed for a homologous protein, alpha-lactalbumin, in its partially folded "molten globule" state, the results clearly show that thermally denatured lysozyme differs substantially from a random coil, principally in that it is partially hydrophobically collapsed.     

High concentrations of viscogens decrease the protein folding rate constant by prematurely collapsing the coil

In several studies, viscogenic osmolytes have been suggested to decrease the folding rate constant of polypeptides by slowing their motion through the solvent. Here, we show that osmolytes may slow protein folding by prematurely collapsing the coil. At low or moderate concentrations of osmolytes (<30%), folding of the two-state protein CI2 becomes faster with increasing osmolyte concentrations, suggesting that the kinetics are governed by protein stability. However, at higher concentrations of osmolyte, the coil collapses in the dead-time of the refolding experiment, causing a dramatic drop in the folding rate. The collapsed state is non-native and appears to be different for different osmolytes.     

Hydrogen exchange in native and denatured states of hen egg-white lysozyme.

The hydrogen exchange kinetics of 68 individual amide protons in the native state of hen lysozyme have been measured at pH 7.5 and 30 degrees C by 2D NMR methods. These constitute the most protected subset of amides, with exchange half lives some 10(5)-10(7) times longer than anticipated from studies of small model peptides. The observed distribution of rates under these conditions can be rationalized to a large extent in terms of the hydrogen bonding of individual amides and their burial from bulk solvent. Exchange rates have also been measured in a reversibly denatured state of lysozyme; this was made possible under very mild conditions, pH 2.0 35 degrees C, by lowering the stability of the native state through selective cleavage of the Cys-6-Cys-127 disulfide cross-link (CM6-127 lysozyme). In this state the exchange rates for the majority of amides approach, within a factor of 5, the values anticipated from small model peptides. For a few amides, however, there is evidence for significant retardation (up to nearly 20-fold) relative to the predicted rates. The pattern of protection observed under these conditions does not reflect the behavior of the protein under strongly native conditions, suggesting that regions of native-like structure do not persist significantly in the denatured state of CM6-127 lysozyme. The pattern of exchange rates from the native protein at high temperature, pH 3.8 69 degrees C, resembles that of the acid-denatured state, suggesting that under these conditions the exchange kinetics are dominated by transient global unfolding. The rates of folding and unfolding under these conditions were determined independently by magnetization transfer NMR methods, enabling the intrinsic exchange rates from the denatured state to be deduced on the basis of this model, under conditions where the predominant equilibrium species is the native state. Again, in the case of most amides these rates showed only limited deviation from those predicted by a simple random coil model. This reinforces the view that these denatured states of lysozyme have little persistent residual order and contrasts with the behavior found for compact partially folded states of proteins, including an intermediate detected transiently during the refolding of hen lysozyme.     

Molecular basis of cooperativity in protein folding IV. Core: A general cooperative folding model

The cooperative nature of the protein folding process is independent of the characteristic fold and the specific secondary structure attributes of a globular protein. A general folding/unfolding model should, therefore, be based upon structural features that transcend the peculiarities of α-helices, β-sheets, and other structural motifs found in proteins. The studies presented in this paper suggest that a single structural characteristic common to all globular proteins is essential for cooperative folding. The formation of a partly folded state from the native state results in the exposure to solvent of two distinct regions: (1) the portions of the protein that are unfolded; and (2) the “complementary surfaces,” located in the regions of the protein that remain folded. The cooperative character of the folding/unfolding transition is determined largely by the energetics of exposing complementary surface regions to the solvent. By definition, complementary regions are present only in partly folded states; they are absent from the native and unfolded states. An unfavorable free energy lowers the probability of partly folded states and increases the cooperativity of the transition. In this paper we present a mathematical formulation of this behavior and develop a general cooperative folding/unfolding model, termed the “complementary region” (CORE) model. This model successfully reproduces the main properties of folding/unfolding transitions without limiting the number of partly folded states accessible to the protein, thereby permitting a systematic examination of the structural and solvent conditions under which intermediates become populated. It is shown that the CORE model predicts two-state folding/unfolding behavior, even though the two-state character is not assumed in the model. © 1993 Wiley-Liss, Inc.     

Compact dimension of denatured states of staphylococcal nuclease

Fluorescence and circular dichroism stopped-flow have been widely used to determine the kinetics of protein folding including folding rates and possible folding pathways. Yet, these measurements are not able to provide spatial information of protein folding/unfolding. Especially, conformations of denatured states cannot be elaborated in detail. In this study, we apply the method of fluorescence energy transfer with a stopped-flow technique to study global structural changes of the staphylococcal nuclease (SNase) mutant K45C, where lysine 45 is replaced by cysteine, during folding and unfolding. By labeling the thiol group of cysteine with TNB (5,5'-dithiobis-2-nitrobenzoic acid) as an energy acceptor and the tryptophan at position 140 as a donor, distance changes between the acceptor and the donor during folding and unfolding are measured from the efficiency of energy transfer. Results indicate that the denatured states of SNase are highly compact regardless of how the denatured states (pH-induced or GdmCl-induced) are induced. The range of distance changes between two probes is between 25.6 and 25.4 A while it is 20.4 A for the native state. Furthermore, the folding process consists of three kinetic phases while the unfolding process is a single phase. These observations agree with our previous sequential model: N(0) left arrow over right arrow D(1) left arrow over right arrow D(2) left arrow over right arrow D(3) (Chen et al., J Mol Biol 1991;220:771-778). The efficiency of protein folding may be attributed to initiating the folding process from these compact denatured structures.     

Understanding the folding of GFP using biophysical techniques

Green fluorescent protein (GFP) and its many variants are probably the most widely used proteins in medical and biological research, having been extensively engineered to act as markers of gene expression and protein localization, indicators of protein–protein interactions and biosensors. GFP first folds, before it can undergo an autocatalytic cyclization and oxidation reaction to form the chromophore, and in many applications the folding efficiency of GFP is known to limit its use. Here, we review the recent literature on protein engineering studies that have improved the folding properties of GFP. In addition, we discuss in detail the biophysical work on the folding of GFP that is beginning to reveal how this large and complex structure forms.     

Conformation of P22 tailspike folding and aggregation intermediates probed by monoclonal antibodies.

The partitioning of partially folded polypeptide chains between correctly folded native states and off-pathway inclusion bodies is a critical reaction in biotechnology. Multimeric partially folded intermediates, representing early stages of the aggregation pathway for the P22 tailspike protein, have been trapped in the cold and isolated by nondenaturing polyacrylamide gel electrophoresis (PAGE) (speed MA, Wang DIC, King J. 1995. Protein Sci 4:900-908). Monoclonal antibodies against tailspike chains discriminate between folding intermediates and native states (Friguet B, Djavadi-Ohaniance L, King J, Goldberg ME. 1994. J Biol Chem 269:15945-15949). Here we describe a nondenaturing Western blot procedure to probe the conformation of productive folding intermediates and off-pathway aggregation intermediates. The aggregation intermediates displayed epitopes in common with productive folding intermediates but were not recognized by antibodies against native epitopes. The nonnative epitope on the folding and aggregation intermediates was located on the partially folded N-terminus, indicating that the N-terminus remained accessible and nonnative in the aggregated state. Antibodies against native epitopes blocked folding, but the monoclonal directed against the N-terminal epitope did not, indicating that the conformation of the N-terminus is not a key determinant of the productive folding and chain association pathway.     

Manipulation of unfolding-induced protein aggregation by peptides selected for aggregate-binding ability through phage display library screening

A phage-displayed library of peptides (12-mer) was screened for the ability to bind to thermally aggregated bovine carbonic anhydrase (BCA), with a view toward examining whether peptides possessing this ability might bind to partially structured intermediates on the protein's unfolding pathway and, therefore, constitute useful tools for manipulation of the kinetic partitioning of molecules between the unfolded and aggregated states. Two peptides [N-HPSTMGLRTMHP-C and N-TPSAWKTALVKA-C] were identified and tested. While neither showed thermal aggregation autonomously, both peptides individually elicited remarkable increases in the levels of thermal aggregation of BCA. A possible explanation is that both peptides bind to surfaces on molten BCA that are not directly involved in aggregation. Such binding could slow down interconversions between folded and unfolded states and stabilize aggregation-prone intermediate(s) to make them more prone to aggregation, while failing to achieve any steric prevention of aggregation. The approach has the potential of yielding useful aggregation-aiding/inhibiting agents, and may provide clues to whether amorphous aggregates are "immobilized" forms of folding intermediates.     

Unstructured intermediate states in single protein force experiments

Recent single-molecule force measurements on single-domain proteins have highlighted a three-state folding mechanism where a stabilized intermediate state (I) is observed on the folding trajectory between the stretched state and the native state. Here we investigate on-lattice protein-like heteropolymer models that lead to a three-state mechanism and show that force experiments can be useful to determine the structure of I. We have mostly found that I is composed of a core stabilized by a high number of native contacts, plus an unstructured extended chain. The lifetime of I is shown to be sensitive to modifications of the protein that spoil the core. We then propose three types of modifications--point mutations, cuts, and circular permutations--aiming at: (1) confirming the presence of the core and (2) determining its location, within one amino acid accuracy, along the polypeptide chain. We also propose force jump protocols aiming to probe the on/off-pathway nature of I.     

Thermodynamic characterization of an equilibrium folding intermediate of staphylococcal nuclease.

High-sensitivity differential scanning calorimetry and CD spectroscopy have been used to probe the structural stability and measure the folding/unfolding thermodynamics of a Pro117-->Gly variant of staphylococcal nuclease. It is shown that at neutral pH the thermal denaturation of this protein is well accounted for by a 2-state mechanism and that the thermally denatured state is a fully hydrated unfolded polypeptide. At pH 3.5, thermal denaturation results in a compact denatured state in which most, if not all, of the helical structure is missing and the beta subdomain apparently remains largely intact. At pH 3.0, no thermal transition is observed and the molecule exists in the compact denatured state within the 0-100 degrees C temperature interval. At high salt concentration and pH 3.5, the thermal unfolding transition exhibits 2 cooperative peaks in the heat capacity function, the first one corresponding to the transition from the native to the intermediate state and the second one to the transition from the intermediate to the unfolded state. As is the case with other proteins, the enthalpy of the intermediate is higher than that of the unfolded state at low temperatures, indicating that, under those conditions, its stabilization must be of an entropic origin. The folding intermediate has been modeled by structural thermodynamic calculations. Structure-based thermodynamic calculations also predict that the most probable intermediate is one in which the beta subdomain is essentially intact and the rest of the molecule unfolded, in agreement with the experimental data. The structural features of the equilibrium intermediate are similar to those of a kinetic intermediate previously characterized by hydrogen exchange and NMR spectroscopy.