Low folding cooperativity of HP35 revealed by single-molecule force spectroscopy and molecular dynamics simulation

Some small proteins, such as HP35, fold at submicrosecond timescale with low folding cooperativity. Although these proteins have been extensively investigated, still relatively little is known about their folding mechanism. Here, using single-molecule force spectroscopy and steered molecule dynamics simulation, we study the unfolding of HP35 under external force. Our results show that HP35 unfolds at extremely low forces without a well-defined unfolding transition state. Subsequently, we probe the structure of unfolded HP35 using the persistence length obtained in the force spectroscopy. We found that the persistence length of unfolded HP35 is around 0.72 nm, >40% longer than typical unstructured proteins, suggesting that there are a significant amount of residual secondary structures in the unfolded HP35. Molecular dynamics simulation further confirmed this finding and revealed that many native contacts are preserved in HP35, even its two ends have been extended up to 8 nm. Our results therefore suggest that retaining a significant amount of secondary structures in the unfolded state of HP35 may be an efficient way to reduce the entropic cost for the formation of tertiary structure and increase the folding speed, although the folding cooperativity is compromised. Moreover, we anticipate that the methods we used in this work can be extended to the study of other proteins with complex folding behaviors and even intrinsically disordered ones.     

Folding/Unfolding Kinetics of Apomyoglobin

The equilibrium and kinetic folding/unfolding of apomyoglobin (ApoMb) were studied at pH 6.2, 11 °C by recording tryptophan fluorescence. The equilibrium unfolding of ApoMb in the presence of urea was shown to involve accumulation of an intermediate state, which had a higher fluorescence intensity as compared with the native and unfolded states. The folding proceeded through two kinetic phases, a rapid transition from the unfolded to the intermediate state and a slow transition from the intermediate to the native state. The accumulation of the kinetic intermediate state was observed in a wide range of urea concentrations. The intermediate was detected even in the region corresponding to the unfolding limb of the chevron plot. Urea concentration dependence was obtained for the observed folding/unfolding rate. The shape of the dependence was compared with that of two-state proteins characterized by a direct transition from the unfolded to the native state.     

The mechanism of folding of globular proteins. Equilibria and kinetics of conformational transitions of penicillinase from Staphylococcus aureus involving a state of intermediate conformation.

1. The thermodynamically reversible unfolding and refolding of penicillinase between the native and fully unfolded states were followed by using guanidinium chloride as denaturant. 2. The equilibria, studied by optical rotation, u.v. absorption, viscosity and enzyme activity, show the presence of a state of intermediate conformation, termed state H, which is stable at 20 degrees C in 0.8 M-guanidinium chloride. 3. The physical properties of this state show that it is slightly expanded with an intrinsic viscosity of 8 ml-g-1, that the 13 tyrosine residues, which are distributed through the primary sequence, are maximally exposed to the solvent and that the helix content is the same as that of the native state. 4. The kinetics of the transition between the native state, state H and the fully unfolded state were followed by u.v. absorption and by optical rotation. They are interpreted as showing that state H lies on the folding pathway between the native and fully unfolded states. 5. The transition between the native state and state H exhibits monophasic unfolding kinetics and biphasic refolding kinetics. This indicates that there must be at least two intermediate states in this process, at least one of which lies on the folding pathway which may also involve cul-de-sac paths. 6. The results are discussed in terms of a mechanism involving rapid stabilization of nucleation regions in a moderately compact but internally solvated structure, with 'native format' [Anfinsen (1973) Science 181, 233-230] secondary structure stabilized by tertiary interaction. The final and rate-limiting step in refolding involves shuffling of these structural elements into the native state. 7. This model is discussed in relation to folding in vivo.     

Topology and sequence in the folding of a TIM barrel protein: global analysis highlights partitioning between transient off-pathway and stable on-pathway folding intermediates in the complex folding mechanism of a (betaalpha)8 barrel of unknown function from B. subtilis

The relative contributions of chain topology and amino acid sequence in directing the folding of a (betaalpha)(8) TIM barrel protein of unknown function encoded by the Bacillus subtilis iolI gene (IOLI) were assessed by reversible urea denaturation and a combination of circular dichroism, fluorescence and time-resolved fluorescence anisotropy spectroscopy. The equilibrium reaction for IOLI involves, in addition to the native and unfolded species, a stable intermediate with significant secondary structure and stability and self-associated forms of both the native and intermediate states. Global kinetic analysis revealed that the unfolded state partitions between an off-pathway refolding intermediate and the on-pathway equilibrium intermediate early in folding. Comparisons with the folding mechanisms of two other TIM barrel proteins, indole-3-glycerol phosphate synthase from the thermophile Sulfolobus solfataricus (sIGPS) and the alpha subunit of Escherichia coli tryptophan synthase (alphaTS), reveal striking similarities that argue for a dominant role of the topology in both early and late events in folding. Sequence-specific effects are apparent in the magnitudes of the relaxation times and relative stabilities, in the presence of additional monomeric folding intermediates for alphaTS and sIGPS and in rate-limiting proline isomerization reactions for alphaTS.     

Essential dynamics of reversible peptide folding: memory-free conformational dynamics governed by internal hydrogen bonds

A principal component analysis has been applied on equilibrium simulations of a beta-heptapeptide that shows reversible folding in a methanol solution. The analysis shows that the configurational space contains only three dense sub-states. These states of relatively low free energy correspond to the "native" left-handed helix, a partly helical intermediate, and a hairpin-like structure. The collection of unfolded conformations form a relatively diffuse cloud with little substructure. Internal hydrogen-bonding energies were found to correlate well with the degree of folding. The native helical structure folds from the N terminus; the transition from the major folding intermediate to the native helical structure involves the formation of the two most C-terminal backbone hydrogen bonds. A four-state Markov model was found to describe transition frequencies between the conformational states within error limits, indicating that memory-effects are negligible beyond the nanosecond time-scale. The dominant native state fluctuations were found to be very similar to unfolding motions, suggesting that unfolding pathways can be inferred from fluctuations in the native state. The low-dimensional essential subspace, describing 69% of the collective atomic fluctuations, was found to converge at time-scales of the order of one nanosecond at all temperatures investigated, whereas folding/unfolding takes place at significantly longer time-scales, even above the melting temperature.     

Tryptophan-tryptophan energy migration as a tool to follow apoflavodoxin folding

Submolecular details of Azotobacter vinelandii apoflavodoxin (apoFD) (un)folding are revealed by time-resolved fluorescence anisotropy using wild-type protein and variants lacking one or two of apoFD's three tryptophans. ApoFD equilibrium (un)folding by guanidine hydrochloride follows a three-state model: native ↔ unfolded ↔ intermediate. In native protein, W128 is a sink for Förster resonance energy transfer (FRET). Consequently, unidirectional FRET with a 50-ps transfer correlation time occurs from W167 to W128. FRET from W74 to W167 is much slower (6.9 ns). In the intermediate, W128 and W167 have native-like geometry because the 50-ps transfer time is observed. However, non-native structure exists between W74 and W167 because instead of 6.9 ns the transfer correlation time is 2.0 ns. In unfolded apoFD this 2.0-ns transfer correlation time is also detected. This decrease in transfer correlation time is a result of W74 and W167 becoming solvent accessible and randomly oriented toward one another. Apparently W74 and W167 are near-natively separated in the folding intermediate and in unfolded apoFD. Both tryptophans may actually be slightly closer in space than in the native state, even though apoFD's radius increases substantially upon unfolding. In unfolded apoFD the 50-ps transfer time observed for native and intermediate folding states becomes 200 ps as W128 and W167 are marginally further separated than in the native state. Apparently, apoFD's unfolded state is not a featureless statistical coil but contains well-defined substructures. The approach presented is a powerful tool to study protein folding.     

Analysis on folding of misgurin using two-dimensional HP model

Misgurin is an antimicrobial peptide from the loach, while the hydrophobic-polar (HP) model is a way to study the folding conformations and native states in peptide and protein although several amino acids cannot be classified either hydrophobic or polar. Practically, the HP model requires extremely intensive computations, thus it has yet to be used widely. In this study, we use the two-dimensional HP model to analyze all possible folding conformations and native states of misgurin with conversion of natural amino acids according to the normalized amino acid hydrophobicity index as well as the shortest benchmark HP sequence. The results show that the conversion of misgurin into HP sequence with glycine as hydrophobic amino acid at pH 2 has 1212 folding conformations with the same native state of minimal energy -6; the conversion of glycine as polar amino acid at pH 2 has 13,386 folding conformations with three native states of minimal energy -5; the conversion of glycine as hydrophobic amino acid at pH 7 has 2538 folding conformations with three native states of minimal energy -5; and the conversion of glycine as polar amino acid at pH 7 has 12,852 folding conformations with three native states of minimal energy -4. Those native states can be ranked according to the normalized amino acid hydrophobicity index. The detailed discussions suggest two ways to modify misgurin.     

Cofactor effects on the protein folding reaction: acceleration of alpha-lactalbumin refolding by metal ions

About 30% of proteins require cofactors for their proper folding. The effects of cofactors on the folding reaction have been investigated with alpha-lactalbumin as a model protein and metal ions as cofactors. Metal ions accelerate the refolding of alpha-lactalbumin by lessening the energy barrier between the molten globule state and the transition state, mainly by decreasing the difference of entropy between the two states. These effects are linked to metal ion binding to the protein in the native state. Hence, relationships between the metal affinities for the intermediate states and those for the native state are observed. Some residual specificity for the calcium ion is still observed in the molten globule state, this specificity getting closer in the transition state to that of the native state. The comparison between kinetic and steady-state data in association with the Phi value method indicates the binding of the metal ions on the unfolded state of alpha-lactalbumin. Altogether, these results provide insight into cofactor effects on protein folding. They also suggest new possibilities to investigate the presence of residual native structures in the unfolded state of protein and the effects of such structures on the protein folding reaction and on protein stability.     

Investigation of folding/unfolding kinetics of apomyoglobin

Apomyoglobin kinetic and equilibrium unfolding and folding processes were studied at pH 6.2, 11 degrees C by stopped-flow tryptophan fluorescence. There are two distinct consecutive processes in apomyoglobin folding process, namely, the protein fast transition between the unfolded (U) and an intermediate (I) states (U <----> I) and slow transition between the intermediate and the native (N) states (I <----> N). Accumulation of the intermediate state was observed in the wide range of urea concentrations. The presence of the intermediate state was shown even beyond the middle transition on the unfolding limb. The dependence of observed folding/unfolding rates on urea concentration (chevron plot) was obtained. The shape of this dependence was compared with that of two-state proteins, folding from the U to N state.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Understanding how small helical proteins fold: conformational dynamics of Im proteins relevant to their folding landscapes

Understanding the mechanism of folding of small proteins requires characterization of their starting unfolded states and any partially unfolded states populated during folding. Here, we review what is known from NMR about these states of Im7, a 4-helix bundle protein that folds via an on-pathway intermediate, and show that there is an alignment of non-native structure in urea-unfolded Im7 with the helices of native Im7 that is a consequence of hydrophobic helix-promoting residues also promoting cluster-formation in the unfolded protein. We suggest that this kind of alignment is present in other proteins and is relevant to how native state topology determines folding rates.     

A folding model of α-lactalbumin deduced from the three-state denaturation mechanism

It has been shown that α-lactalbumin undergoes a three-state denaturation, involving a helical intermediate state, on treatment with guanidine hydrochloride. The unfolding of the protein and the characteristics of the intermediate state are examined by means of circular dichroism, difference spectra and pH-jump measurements to investigate the temperature dependence and kinetic properties of the unfolding and refolding, the pH dependence of the transition between the intermediate and the fully unfolded states, and the effect of disulphide bond reduction on the stabilization of the intermediate.The results show that the long-range specific interactions such as specific electrostatic interactions and disulphide linkages are not important for stabilizing the intermediate, and that the transition between the intermediate and the fully unfolded states is extremely rapid (a relaxation time of less than one millisecond) and may correspond to the helix-coil transition of a polypeptide backbone. On the other hand, the activation parameters of the transition between the native and the intermediate states have suggested that the final stabilization by charge-pair interactions is preceded by hydrophobic interactions in the process of going from the intermediate to the native state.The mechanism of folding of the protein is discussed, and the folding process from the fully unfolded to the native state is apparently divided into at least three main steps: (1) the formation of incipient helical structures dictated by local interactions; (2) the packing of the helical segments accompanied with hydrophobic interactions; (3) the final stabilization by the electrostatic interactions. The relevance to the current theoretical results on protein folding is also discussed.     

Test for cooperativity in the early kinetic intermediate in lysozyme folding

During the folding of many proteins, collapsed globular states are formed prior to the native structure. The role of these states for the folding process has been widely discussed. Comparison with properties of synthetic homo and heteropolymers had suggested that the initial collapse represented a shift of the ensemble of unfolded conformations to more compact states without major energy barriers. We investigated the folding/unfolding transition of a collapsed state, which transiently populates early in lysozyme folding. This state forms within the dead-time of stopped-flow mixing and it has been shown to be significantly more compact and globular than the denaturant-induced unfolded state. We used the GdmCl-dependence of the dead-time signal change to characterize the unfolding transition of the burst phase intermediate. Fluorescence and far-UV CD give identical unfolding curves, arguing for a cooperative two-state folding/unfolding transition between unfolded and collapsed lysozyme. These results show that collapse leads to a distinct state in the folding process, which is separated from the ensemble of unfolded molecules by a significant energy barrier. NMR, fluorescence and small angle X-ray scattering data further show that some local interactions in unfolded lysozyme exist at denaturant concentrations above the coil-collapse transition. These interactions might play a crucial role in the kinetic partitioning between fast and slow folding pathways.     

NMR characterization of a peptide model provides evidence for significant structure in the unfolded state of the villin headpiece helical subdomain

The villin headpiece subdomain (HP36) is the smallest naturally occurring protein that folds cooperatively. The protein folds on a microsecond time scale. Its small size and very rapid folding have made it a popular target for biophysical studies of protein folding. Temperature-dependent one-dimensional (1D) NMR studies of the full-length protein together with CD and 1D NMR studies of the 21-residue peptide fragment (HP21) derived from HP36 have shown that there is significant structure in the unfolded state of HP36 and have demonstrated that HP21 is a good model of these interactions. Here, we characterized the model peptide HP21 in detail by two-dimensional NMR. Strongly upfield shifted C(alpha) protons, the magnitude of the 3J(NH,alpha) coupling constants, and the pattern of backbone-backbone and backbone-side chain NOEs indicate that the ensemble of structures populated by HP21 contains alpha-helical structure and native as well as non-native hydrophobic contacts. The hydrogen-bonded secondary structure inferred from the NOEs is, however, not sufficient to confer significant protection against amide H-D exchange. These studies indicate that there is significant secondary structure and hydrophobic clustering in the unfolded state of HP36. The implications for the folding of HP36 are discussed.     

Equilibrium unfolding of mutant apomyoglobins carrying substitutions of conserved nonfunctional residues with alanine

Protein aggregation or misfolding in the cell is connected with many genetic diseases and can result from substitutions in proteins. Substitutions can influence the protein stability and folding rates in both intermediate and native states. The equilibrium urea-induced unfolding was studied for mutant apomyoglobins carrying substitutions of the conserved nonfunctional residues Val10, Trp14, Ile111, Leu115, Met131, and Leu135 with Ala. Conformational transitions were monitored by intrinsic Trp fluorescence and far-UV circular dichroism. Free energy changes upon transition from the native to the intermediate state and from the intermediate to the unfolded state were determined. All substitutions considerably decreased the stability of native apomyoglobin, whereas the effect on the stability of the intermediate state was essentially smaller.     

Mapping the folding free energy surface for metal-free human Cu,Zn superoxide dismutase

Mutations at many different sites in the gene encoding human Cu,Zn superoxide dismutase (SOD) are known to be causative agents in amyotrophic lateral sclerosis (ALS). One explanation for the molecular basis of this pathology is the aggregation of marginally soluble, partially structured states whose populations are enhanced in the protein variants. As a benchmark for testing this hypothesis, the equilibrium and kinetic properties of the reversible folding reaction of a metal-free variant of SOD were investigated. Reversibility was achieved by replacing the two non-essential cysteine residues with non-oxidizable analogs, C6A/C111S, to produce apo-AS-SOD. The metal-free pseudo-wild-type protein is folded and dimeric in the absence of chemical denaturants, and its equilibrium folding behavior is well described by an apparent two-state mechanism involving the unfolded monomer and the native dimer. The apparent free energy of folding in the absence of denaturant and at standard state is -20.37(+/- 1.04) kcal (mol dimer)(-1). A global analysis of circular dichroism kinetic traces for both unfolding and refolding reactions, combined with results from small angle X-ray scattering and time-resolved fluorescence anisotropy measurements, supports a sequential mechanism involving the unfolded monomer, a folded monomeric intermediate, and the native dimer. The rate-limiting monomer folding reaction is followed by a near diffusion-limited self-association reaction to form the native dimer. The relative population of the folded monomeric intermediate is predicted not to exceed 0.5% at micromolar concentrations of protein under equilibrium and both strongly unfolding and refolding conditions for metal-free pseudo-wild-type SOD.     

Native protein sequences are designed to destabilize folding intermediates

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.     

Trehalose favors a cutinase compact intermediate off-folding pathway

The folding of cutinase, an enzyme displaying lipolytic activity, has been studied in the presence of trehalose. Equilibrium unfolding data show that trehalose increases the free energy change between folded and unfolded states. Unfolding kinetics reveal the presence of an intermediate which is ca. 60% folded in terms of solvent exposure. Trehalose stabilizes this intermediate relative to the folded state. In contrast, the intermediate revealed by folding kinetics is more compact than the transition state, as shown by the positive slope observed at low denaturant concentration in the chevron plot, as well as the decrease in the observable rate constant for folding with the increase in trehalose concentration. This intermediate displays more than 50% of area buried from the solvent (relative to the native state) compared to around 40% for the transition state for folding and therefore appears to be off the folding pathway. Trehalose stabilizes and guanidine hydrochloride destabilizes this compact intermediate. Both unfolding and folding kinetics show that compact conformational states are stabilized by trehalose, in agreement with current models on the effect of compatible solutes. This effect occurs even for compact states that decelerate the folding as in the case of the intermediate revealed by folding kinetics.