Amino Acid Insertion Reveals a Necessary Three-Helical Intermediate in the Folding Pathway of the Colicin E7 Immunity Protein Im7

The small (87-residue) α-helical protein Im7 (an inhibitor protein for colicin E7 that provides immunity to cells producing colicin E7) folds via a three-state mechanism involving an on-pathway intermediate. This kinetic intermediate contains three of four native helices that are oriented in a non-native manner so as to minimise exposed hydrophobic surface area at this point in folding. The short (6-residue) helix III has been shown to be unstructured in the intermediate ensemble and does not dock onto the developing hydrophobic core until after the rate-limiting transition state has been traversed. After helix III has docked, it adopts an α-helical secondary structure, and the side chains of residues within this region provide contacts that are crucial to native-state stability. In order to probe further the role of helix III in the folding mechanism of Im7, we created a variant that contains an eight-amino-acid polyalanine-like helix stabilised by a Glu-Arg salt bridge and an Asn-Pro-Gly capping motif, juxtaposed C-terminal to the natural 6-residue helix III. The effect of this insertion on the structure of the native protein and its folding mechanism were studied using NMR and ?-value analysis, respectively. The results reveal a robust native structure that is not perturbed by the presence of the extended helix III. Mutational analysis performed to probe the folding mechanism of the redesigned protein revealed a conserved mechanism involving the canonical three-helical intermediate. The results suggest that folding via a three-helical species stabilised by both native and non-native interactions is an essential feature of Im7 folding, independent of the helical propensity of helix III.     

Electrostatic interactions in the denatured state and in the transition state for protein folding: effects of denatured state interactions on the analysis of transition state structure

The development of electrostatic interactions during the folding of the N-terminal domain of the ribosomal protein L9 (NTL9) is investigated by pH-dependent rate equilibrium free energy relationships. We show that Asp8, among six acidic residues, is involved in non-native, electrostatic interactions with K12 in the transition state for folding as well as in the denatured state. The perturbed native state pK(a) of D8 (pK(a) = 3.0) appears to be maintained through non-native interactions in both the transition state and the denatured state. Mutational effects on the stability of the transition state for protein (un)folding are often analyzed in respect to change in ground states. Thus, the interpretation of transition state analysis critically depends on an understanding of mutational effects on both the native and denatured state. Increasing evidence for structurally biased denatured states under physiological conditions raises concerns about possible denatured state effects on folding studies. We show that the structural interpretation of transition state analysis can be altered dramatically by denatured state effects.     

Desolvation and Development of Specific Hydrophobic Core Packing during Im7 Folding

Development of a tightly packed hydrophobic core drives the folding of water-soluble globular proteins and is a key determinant of protein stability. Despite this, there remains much to be learnt about how and when the hydrophobic core becomes desolvated and tightly packed during protein folding. We have used the bacterial immunity protein Im7 to examine the specificity of hydrophobic core packing during folding. This small, four-helix protein has previously been shown to fold via a compact three-helical intermediate state. Here, overpacking substitutions, in which residue side-chain size is increased, were used to examine the specificity and malleability of core packing in the folding intermediate and rate-limiting transition state. In parallel, polar groups were introduced into the Im7 hydrophobic core via Val→Thr or Phe→Tyr substitutions and used to determine the solvation status of core residues at different stages of folding. Over 30 Im7 variants were created allowing both series of substitutions to cover all regions of the protein structure. Φ-value analysis demonstrated that the major changes in Im7 core solvation occur prior to the population of the folding intermediate, with key regions involved in docking of the short helix III remaining solvent-exposed until after the rate-limiting transition state has been traversed. In contrast, overpacking core residues revealed that some regions of the native Im7 core are remarkably malleable to increases in side-chain volume. Overpacking residues in other regions of the Im7 core result in substantial (> 2.5 kJ mol^− 1) destabilisation of the native structure or even prevents efficient folding to the native state. This study provides new insights into Im7 folding; demonstrating that whilst desolvation occurs early during folding, adoption of a specifically packed core is achieved only at the very last step in the folding mechanism.     

Sulfate-induced effects in the on-pathway intermediate of the bacterial immunity protein Im7

Intermediates have now been identified in the folding of a number of small, single-domain proteins. Here we describe experiments to determine the effect of Na(2)SO(4) on the properties of the on-pathway intermediate formed early during the folding of the four-helical protein, Im7. This intermediate, studied previously in 0.4 M Na(2)SO(4), contains three of the four native helices and is fascinating in that several residues in helices I, II, and IV make non-native interactions that stabilize this state. Whether these contacts form as a consequence of the presence of Na(2)SO(4), however, remained unresolved. Using kinetic analysis of the effect of Na(2)SO(4) on the unfolding and refolding kinetics of Im7*, combined with detailed analysis of the resulting chevron plots, we show that decreasing the concentration of Na(2)SO(4) from 0.4 to 0 M destabilizes the intermediate and rate-limiting transition (TS2) states by 7 and 10 kJ mol(-)(1), respectively, and has little effect on the relative compactness of these states compared with that of the unfolded ensemble (beta(I) approximately 0.8, beta(TS2) approximately 0.9 in 0 to 0.4 M Na(2)SO(4)). Analysis of 10 variants of the protein in 0.2 M Na(2)SO(4) using Phi-values showed that the structural properties of the intermediate and TS2 are not altered significantly by the concentration of the kosmotrope. The data demonstrate that the rapid formation of a compact intermediate stabilized by non-native interactions during Im7* folding is not induced by high concentrations of the stabilizing salt, but is a generic feature of the folding of this protein.     

The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H

Elucidation of the high-resolution structures of folding intermediates is a necessary but difficult step toward the ultimate understanding of the mechanism of protein folding. Here, using hydrogen-exchange-directed protein engineering, we populated the folding intermediate of the Thermus thermophilus ribonuclease H, which forms before the rate-limiting transition state, by removing the unfolded regions of the intermediate, including an α-helix and two β-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 Å). It has a native-like backbone topology and shows some local deviations from the native structure, revealing that the structure of the folded region of an early folding intermediate can be as well defined as the native structure. The topological parameters calculated from the structures of the intermediate mimic and the native state predict that the intermediate should fold on a millisecond time scale or less and form much faster than the native state. Other factors that may lead to the slow folding of the native state and the accumulation of the intermediate before the rate-limiting transition state are also discussed.     

Nonnative Electrostatic Interactions Can Modulate Protein Folding: Molecular Dynamics with a Grain of Salt

In recent years, a growing number of protein folding studies have focused on the unfolded state, which is now recognized as playing a major role in the folding process. Some of these studies show that interactions occurring in the unfolded state can significantly affect the stability and kinetics of the protein folding reaction. In this study, we modeled the effect of electrostatic interactions, both native and nonnative, on the folding of three protein systems that underwent selective charge neutralization or reversal or complete charge suppression. In the case of the N-terminal L9 protein domain, our results directly attribute the increase in thermodynamic stability to destabilization of the unfolded ensemble, reaffirming the experimental observations. These results provide a deeper structural insight into the ensemble of the unfolded state and predict a new mutation site for increased protein stability. In the second case, charge reversal mutations of RNase Sa affected protein stability, with the destabilizing mutations being less destabilizing at higher salt concentrations, indicating the formation of charge-charge interactions in the unfolded state. In the N-terminal L9 and RNase Sa systems, changes in electrostatic interactions in the unfolded state that cause an increase in free energy had an overall compaction effect that suggests a decrease in entropy. In the third case, in which we compared the β-lactalbumin and hen egg-white lysozyme protein homologues, we successfully eliminated differences between the folding kinetics of the two systems by suppressing electrostatic interactions, supporting previously reported findings. Our coarse-grained molecular dynamics study not only reproduces experimentally reported findings but also provides a detailed molecular understanding of the elusive unfolded-state ensemble and how charge-charge interactions can modulate the biophysical characteristics of folding.     

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.     

Mutational analysis demonstrates that specific electrostatic interactions can play a key role in the denatured state ensemble of proteins

The nature of the denatured state ensemble has been controversial for decades owing, in large part, to the difficulty in characterizing the structure and energetics of denatured state interactions. There is increasing evidence for relatively non-specific hydrophobic clustering in the denatured states of some proteins but other types of interactions are much less well characterized. Here, we report the characterization of highly specific electrostatic interactions in the denatured state of a small alpha-beta protein, the N-terminal domain of the ribosomal protein L9 (NTL9). Mutation of Lys12 to Met has been shown to increase the stability of NTL9 significantly through the disruption of denatured state interactions. Here, we describe the analysis of the pH-dependent stability of 13 mutants designed to probe the nature of the Lys12 denatured state interaction. Lys12 is located in a lysine-rich region of the protein but analysis of a set of Lys to Met mutants shows that it plays a unique role in the denatured state. Analysis of mutants of all of the acidic residues in NTL9 shows that Lys12 forms a specific non-native electrostatic interaction with Asp8 in the denatured state ensemble. Thus the distribution of charge-charge interactions in the denatured state ensemble of NTL9 appears to be biased by few key interactions and is very different from that expected in a random coil. We propose that these interactions are not encoded by local sequence effects but rather reflect interactions among residues more distant in sequence. These results demonstrate that electrostatic as well as hydrophobic interactions can play an important role in the denatured state ensemble.     

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.     

Structural comparison of the two alternative transition states for folding of TI I27

TI I27, a beta-sandwich domain from the human muscle protein titin, has been shown to fold via two alternative pathways, which correspond to a change in the folding mechanism. Under physiological conditions, TI I27 folds by a classical nucleation-condensation mechanism (diffuse transition state), whereas at extreme conditions of temperature and denaturant it switches to having a polarized transition state. We have used experimental Phi-values as restraints in ensemble-averaged molecular dynamics simulations to determine the ensembles of structures representing the two transition states. The comparison of these ensembles indicates that when native interactions are substantially weakened, a protein may still be able to fold if it can access an alternative transition state characterized by a much larger entropic contribution. Analysis of the probability distribution of Phi-values derived from ensemble averaged simulations, enables us to identify residues that form contacts in some members of the ensemble but not in others illustrating that many interactions present in transition states are not strictly required for the successful completion of the folding process.     

The paradoxical behavior of a highly structured misfolded intermediate in RNA folding

Like many structured RNAs, the Tetrahymena group I ribozyme is prone to misfolding. Here we probe a long-lived misfolded species, referred to as M, and uncover paradoxical aspects of its structure and folding. Previous work indicated that a non-native local secondary structure, termed alt P3, led to formation of M during folding in vitro. Surprisingly, hydroxyl radical footprinting, fluorescence measurements with site-specifically incorporated 2-aminopurine, and functional assays indicate that the native P3, not alt P3, is present in the M state. The paradoxical behavior of alt P3 presumably arises because alt P3 biases folding toward M, but, after commitment to this folding pathway and before formation of M, alt P3 is replaced by P3. Further, structural and functional probes demonstrate that the misfolded ribozyme contains extensive native structure, with only local differences between the two states, and the misfolded structure even possesses partial catalytic activity. Despite the similarity of these structures, re-folding of M to the native state is very slow and is strongly accelerated by urea, Na+, and increased temperature and strongly impeded by Mg2+ and the presence of native peripheral contacts. The paradoxical observations of extensive native structure within the misfolded species but slow conversion of this species to the native state are readily reconciled by a model in which the misfolded state is a topological isomer of the native state, and computational results support the feasibility of this model. We speculate that the complex topology of RNA secondary structures and the inherent rigidity of RNA helices render kinetic traps due to topological isomers considerably more common for RNA than for proteins.     

Mesophile versus thermophile: insights into the structural mechanisms of kinetic stability

Obtaining detailed knowledge of folding intermediate and transition state (TS) structures is critical for understanding protein folding mechanisms. Comparisons between proteins adapted to survive extreme temperatures with their mesophilic homologs are likely to provide valuable information on the interactions relevant to the unfolding transition. For kinetically stable proteins such as alpha-lytic protease (alphaLP) and its family members, their large free energy barrier to unfolding is central to their biological function. To gain new insights into the mechanisms that underlie kinetic stability, we have determined the structure and high temperature unfolding kinetics of a thermophilic homolog, Thermobifida fusca protease A (TFPA). These studies led to the identification of a specific structural element bridging the N and C-terminal domains of the protease (the "domain bridge") proposed to be associated with the enhanced high temperature kinetic stability in TFPA. Mutagenesis experiments exchanging the TFPA domain bridge into alphaLP validate this hypothesis and illustrate key structural details that contribute to TFPA's increased kinetic thermostability. These results lead to an updated model for the unfolding transition state structure for this important class of proteases in which domain bridge undocking and unfolding occurs at or before the TS. The domain bridge appears to be a structural element that can modulate the degree of kinetic stability of the different members of this class of proteases.     

Small proteins fold through transition states with native-like topologies

The folding pathway of common-type acyl phosphatase (ctAcP) is characterized using psi-analysis, which identifies specific chain-chain contacts using bi-histidine (biHis) metal-ion binding sites. In the transition state ensemble (TSE), the majority of the protein is structured with a near-native topology, only lacking one beta-strand and an alpha-helix. psi-Values are zero or unity for all sites except one at the amino terminus of helix H2. This fractional psi-value remains unchanged when three metal ions of differing coordination geometries are used, indicating this end of the helix experiences microscopic heterogeneity through fraying in the TSE. Ubiquitin, the other globular protein characterized using psi-analysis, also exhibits a single consensus TSE structure. Hence, the TSE of both proteins have converged to a single configuration, albeit one that contains some fraying at the periphery. Models of the TSE of both proteins are created using all-atom Langevin dynamics simulations using distance constraints derived from the experimental psi-values. For both proteins, the relative contact order of the TS models is approximately 80% of the native value. This shared value viewed in the context of the known correlation between contact order and folding rates, suggests that other proteins will have a similarly high fraction of the native contact order. This constraint greatly limits the range of possible configurations at the rate-limiting step.     

Specificity of the initial collapse in the folding of the cold shock protein

The two-state folding reaction of the cold shock protein from Bacillus caldolyticus (Bc-Csp) is preceded by a rapid chain collapse. A fast shortening of intra-protein distances was revealed by F?rster resonance energy transfer (FRET) measurements with protein variants that carried individual pairs of donor and acceptor chromophores at various positions along the polypeptide chain. Here we investigated the specificity of this rapid compaction. Energy transfer experiments that probed the stretching of strand beta2 and the close approach between the strands beta1 and beta2 revealed that the beta1-beta2 hairpin is barely formed in the collapsed form, although it is native-like in the folding transition state of Bc-Csp. The time course of the collapse could not be resolved by pressure or temperature jump experiments, indicating that the collapsed and extended forms are not separated by an energy barrier. The co-solute (NH4)2SO4 stabilizes both native Bc-Csp and the collapsed form, which suggests that the large hydrated SO4(2-) ions are excluded from the surface of the collapsed form in a similar fashion as they are excluded from folded Bc-Csp. Ethylene glycol increases the stability of proteins because it is excluded preferentially from the backbone, which is accessible in the unfolded state. The collapsed form of Bc-Csp resembles the unfolded form in its interaction with ethylene glycol, suggesting that in the collapsed form the backbone is still accessible to water and small molecules. Our results thus rule out that the collapsed form is a folding intermediate with native-like chain topology. It is better described as a mixture of compact conformations that belong to the unfolded state ensemble. However, some of its structural elements are reminiscent of the native protein.     

Electrostatic interactions in the denatured state ensemble: their effect upon protein folding and protein stability

It is now recognized that the denatured state ensemble (DSE) of proteins can contain significant amounts of structure, particularly under native conditions. Well-studied examples include small units of hydrogen bonded secondary structure, particularly helices or turns as well as hydrophobic clusters. Other types of interactions are less well characterized and it has often been assumed that electrostatic interactions play at most a minor role in the DSE. However, recent studies have shown that both favorable and unfavorable electrostatic interactions can be formed in the DSE. These can include surprisingly specific non-native interactions that can even persist in the transition state for protein folding. DSE electrostatic interactions can be energetically significant and their modulation either by mutation or by varying solution conditions can have a major impact upon protein stability. pH dependent stability studies have shown that electrostatic interactions can contribute up to 4 kcal mol^-1 to the stability of the DSE.     

Probing the folding intermediate of Rd-apocyt b562 by protein engineering and infrared T-jump

Small proteins often fold in an apparent two-state manner with the absence of detectable early-folding intermediates. Recently, using native-state hydrogen exchange, intermediates that exist after the rate-limiting transition state have been identified for several proteins. However, little is known about the folding kinetics from these post-transition intermediates to their corresponding native states. Herein, we have used protein engineering and a laser-induced temperature-jump (T-jump) technique to investigate this issue and have applied it to Rd-apocyt b(562) , a four-helix bundle protein. Previously, it has been shown that Rd-apocyt b(562) folds via an on-pathway hidden intermediate, which has only the N-terminal helix unfolded. In the present study, a double mutation (V16G/I17A) in the N-terminal helix of Rd-apocyt b(562) was made to further increase the relative population of this intermediate state at high temperature by selectively destabilizing the native state. In the circular dichroism thermal melting experiment, this mutant showed apparent two-state folding behavior. However, in the T-jump experiment, two kinetic phases were observed. Therefore, these results are in agreement with the idea that a folding intermediate is populated on the folding pathway of Rd-apocyt b(562) . Moreover, it was found that the exponential growth rate of the native state from this intermediate state is roughly (25 microsec)(-1) at 65 degrees C.     

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.     

Im7 folding mechanism: misfolding on a path to the native state

Many proteins populate collapsed intermediate states during folding. In order to elucidate the nature and importance of these species, we have mapped the structure of the on-pathway intermediate of the four-helix protein, Im7, together with the conformational changes it undergoes as it folds to the native state. Kinetic data for 29 Im7 point mutants show that the intermediate contains three of the four helices found in the native structure, packed around a specific hydrophobic core. However, the intermediate contains many non-native interactions; as a result, hydrophobic interactions become disrupted in the rate-limiting transition state before the final helix docks onto the developing structure. The results of this study support a hierarchical mechanism of protein folding and explain why the misfolding of Im7 occurs. The data also demonstrate that non-native interactions can play a significant role in folding, even for small proteins with simple topologies.