Transition State Movement Research Articles

Force-dependent switch in protein unfolding pathways and transition-state movements.

Although it is known that single-domain proteins fold and unfold by parallel pathways, demonstration of this expectation has been difficult to establish in experiments. Unfolding rate, [Formula: see text], as a function of force f, obtained in single-molecule pulling experiments on src SH3 domain, exhibits upward curvature on a [Formula: see text] plot. Similar observations were reported for other proteins for the unfolding rate [Formula: see text]. These findings imply unfolding in these single-domain proteins involves a switch in the pathway as f or [Formula: see text] is increased from a low to a high value. We provide a unified theory demonstrating that if [Formula: see text] as a function of a perturbation (f or [Formula: see text]) exhibits upward curvature then the underlying energy landscape must be strongly multidimensional. Using molecular simulations we provide a structural basis for the switch in the pathways and dramatic shifts in the transition-state ensemble (TSE) in src SH3 domain as f is increased. We show that a single-point mutation shifts the upward curvature in [Formula: see text] to a lower force, thus establishing the malleability of the underlying folding landscape. Our theory, applicable to any perturbation that affects the free energy of the protein linearly, readily explains movement in the TSE in a β-sandwich (I27) protein and single-chain monellin as the denaturant concentration is varied. We predict that in the force range accessible in laser optical tweezer experiments there should be a switch in the unfolding pathways in I27 or its mutants.

Kinetic network models of tryptophan mutations in β-hairpins reveal the importance of non-native interactions.

We present an analysis of the most extensive explicit-solvent simulations of β-hairpins to date (9.4 ms in aggregate), with the aim of probing the effects of tryptophan mutations on folding. From molecular simulations of GB1 hairpin, trpzip4, trpzip5, and trpzip6 performed on Folding@home, Markov State Models (MSMs) were constructed using a unified set of metastable states, enabling objective comparison of folding mechanisms. MSM models display quantitative agreement with experimental structural observables and folding kinetics, and predict multimodal kinetics due to specific non-native kinetic traps, which be identified as on- or off-pathway from the network topology. We quantify kinetic frustration by several means, including the s-ensemble method to evaluate glasslike behavior. Free-energy profiles and transition state movement clearly show stabilization of non-native states as Trp mutations are introduced. Remarkably, we find that "β-capped" sequences (trpzip4 and trpzip5) are able to overcome this frustration and remain cooperative two-state folders with a large time-scale gap. These results suggest that, while β-capping motifs are robust, fold stabilization by tryptophan generally may require overcoming significant non-native kinetic traps, perhaps explaining their under-representation in natural proteins.

The transition state for integral membrane protein folding

Biology relies on the precise self-assembly of its molecular components. Generic principles of protein folding have emerged from extensive studies on small, water-soluble proteins, but it is unclear how these ideas are translated into more complex situations. In particular, the one-third of cellular proteins that reside in biological membranes will not fold like water-soluble proteins because membrane proteins need to expose, not hide, their hydrophobic surfaces. Here, we apply the powerful protein engineering method of Phi-value analysis to investigate the folding transition state of the alpha-helical membrane protein, bacteriorhodopsin, from a partially unfolded state. Our results imply that much of helix B of the seven-transmembrane helical protein is structured in the transition state with single-point alanine mutations in helix B giving Phi values >0.8. However, residues Y43 and T46 give lower Phi values of 0.3 and 0.5, respectively, suggesting a possible reduction in native structure in this region of the helix. Destabilizing mutations also increase the activation energy of folding, which is accompanied by an apparent movement of the transition state toward the partially unfolded state. This apparent transition state movement is most likely due to destabilization of the structured, unfolded state. These results contrast with the Hammond effect seen for several water-soluble proteins in which destabilizing mutations cause the transition state to move toward, and become closer in energy to, the folded state. We thus introduce a classic folding analysis method to membrane proteins, providing critical insight into the folding transition state.

Consequences of Domain Insertion on the Stability and Folding Mechanism of a Protein

SlyD, the sensitive-to-lysis protein from Escherichia coli, consists of two domains. They are not arranged successively along the protein chain, but one domain, the “insert-in-flap” (IF) domain, is inserted internally as a guest into a surface loop of the host domain, which is a prolyl isomerase of the FK506 binding protein (FKBP) type. We used SlyD as a model to elucidate how such a domain insertion affects the stability and folding mechanism of the host and the guest domain. For these studies, the two-domain protein was compared with a single-domain variant SlyDΔIF, SlyD* without the chaperone domain (residues 1–69 and 130–165) in which the IF domain was removed and replaced by a short loop, as present in human FKBP12. Equilibrium unfolding and folding kinetics followed an apparent two-state mechanism in the absence and in the presence of the IF domain. The inserted domain decreased, however, the stability of the host domain in the transition region and decelerated its refolding reaction by about 10-fold. This originates from the interruption of the chain connectivity by the IF domain and its inherent instability. To monitor folding processes in this domain selectively, a Trp residue was introduced as fluorescent probe. Kinetic double-mixing experiments revealed that, in intact SlyD, the IF domain folds and unfolds about 1000-fold more rapidly than the FKBP domain, and that it is strongly stabilized when linked with the folded FKBP domain. The unfolding limbs of the kinetic chevrons of SlyD show a strong downward curvature. This deviation from linearity is not caused by a transition-state movement, as often assumed, but by the accumulation of a silent unfolding intermediate at high denaturant concentrations. In this kinetic intermediate, the FKBP domain is still folded, whereas the IF domain is already unfolded.

Temperature-Dependent Hammond Behavior in a Protein-Folding Reaction: Analysis of Transition-State Movement and Ground-State Effects

Characterization of the transition-state ensemble and the nature of the free-energy barrier for protein folding are areas of intense activity and some controversy. A key issue that has emerged in recent years is the width of the free-energy barrier and the susceptibility of the transition state to movement. Here we report denaturant-induced and temperature-dependent folding studies of a small mixed α–β protein, the N-terminal domain of L9 (NTL9). The folding of NTL9 was determined using fluorescence-detected stopped-flow fluorescence measurements conducted at seven different temperatures between 11 and 40 °C. Plots of the log of the observed first-order rate constant versus denaturant concentration, “chevron plots,” displayed the characteristic V shape expected for two-state folding. There was no hint of deviation from linearity even at the lowest denaturant concentrations. The relative position of the transition state, as judged by the Tanford β parameter, β T, shifts towards the native state as the temperature is increased. Analysis of the temperature dependence of the kinetic and equilibrium m values indicates that the effect is due to significant movement of the transition state and also includes a contribution from temperature-dependent ground-state effects. Analysis of the Leffler plots, plots of Δ G ‡ versus Δ G°, and their cross-interaction parameters confirms the transition-state movement. Since the protein is destabilized at high temperature, the shift represents a temperature-dependent Hammond effect. This provides independent confirmation of a recent theoretical prediction. The magnitude of the temperature-denaturant cross-interaction parameter is larger for NTL9 than has been reported for the few other cases studied. The implications for temperature-dependent studies of protein folding are discussed.

An On-pathway Intermediate in the Folding of a PDZ Domain

The folding pathways of some proteins include the population of partially structured species en route to the native state. Identification and characterization of these folding intermediates are particularly difficult as they are often only transiently populated and play different mechanistic roles, being either on-pathway productive species or off-pathway kinetic traps. To define the role of folding intermediates, a quantitative analysis of the folding and unfolding rate constants over a wide range of denaturant concentration is often required. Such a task is further complicated by the reversible nature of the folding reaction, which implies the observed kinetics to be governed by a complex combination of different microscopic rate constants. Here, we tackled this problem by measuring directly the folding rate constant under highly denaturing conditions, namely by inducing the folding of a PDZ domain through a quasi-irreversible binding reaction with a specific peptide. In analogy with previous works based on hydrogen exchange experiments, we present evidence that the folding pathway of the PDZ domain involves the formation of an obligatory on-pathway intermediate. The results presented exemplify a novel type of kinetic test to detect an on-pathway folding intermediate.

Correlation between Protein Stability Cores and Protein Folding Kinetics: A Case Study on Pseudomonas aeruginosa Apo-Azurin

This paper reports a combined computational and experimental study of the correlation between protein stability cores and folding kinetics. An empirical potential function was developed, and it was used for analyzing interaction energies among secondary structure elements. Studies on a beta sandwich protein, Pseudomonas aeruginosa azurin, showed that the computationally identified substructure with the strongest interactions in the native state is identical to the "interlocked pair" of beta strands, an invariant motif found in most sandwich-like proteins. Moreover, previous and new in vitro folding results revealed that the identified substructure harbors most residues that form native-like interactions in the folding transition state. These observations demonstrate that the potential function is effective in revealing the relative strength of interactions among various protein parts; they also strengthen the suggestion that the most stable regions in native proteins favor stable interactions early during folding.

Forced-Unfolding and Force-Quench Refolding of RNA Hairpins

Nanomanipulation of individual RNA molecules, using laser optical tweezers, has made it possible to infer the major features of their energy landscape. Time-dependent mechanical unfolding trajectories, measured at a constant stretching force (fS) of simple RNA structures (hairpins and three-helix junctions) sandwiched between RNA/DNA hybrid handles show that they unfold in a reversible all-or-none manner. To provide a molecular interpretation of the experiments we use a general coarse-grained off-lattice Gō-like model, in which each nucleotide is represented using three interaction sites. Using the coarse-grained model we have explored forced-unfolding of RNA hairpin as a function of fS and the loading rate (rf). The simulations and theoretical analysis have been done both with and without the handles that are explicitly modeled by semiflexible polymer chains. The mechanisms and timescales for denaturation by temperature jump and mechanical unfolding are vastly different. The directed perturbation of the native state by fS results in a sequential unfolding of the hairpin starting from their ends, whereas thermal denaturation occurs stochastically. From the dependence of the unfolding rates on rf and fS we show that the position of the unfolding transition state is not a constant but moves dramatically as either rf or fS is changed. The transition-state movements are interpreted by adopting the Hammond postulate for forced-unfolding. Forced-unfolding simulations of RNA, with handles attached to the two ends, show that the value of the unfolding force increases (especially at high pulling speeds) as the length of the handles increases. The pathways for refolding of RNA from stretched initial conformation, upon quenching fS to the quench force fQ, are highly heterogeneous. The refolding times, upon force-quench, are at least an order-of-magnitude greater than those obtained by temperature-quench. The long fQ-dependent refolding times starting from fully stretched states are analyzed using a model that accounts for the microscopic steps in the rate-limiting step, which involves the trans to gauche transitions of the dihedral angles in the GAAA tetraloop. The simulations with explicit molecular model for the handles show that the dynamics of force-quench refolding is strongly dependent on the interplay of their contour length and persistence length and the RNA persistence length. Using the generality of our results, we also make a number of precise experimentally testable predictions.

Shape of the Free Energy Barriers for Protein Folding Probed by Multiple Perturbation Analysis

The characterization of the free energy barriers has been a major goal in studies on the mechanism of protein folding. Testing the effect of mutations or denaturants on protein folding reactions revealed that transition state movement is rare, suggesting that folding barriers are robust and narrow maxima on the free energy landscape. Here we demonstrate that the application of multiple perturbations allows the observation of small transition state movements that escape detection in single perturbation experiments. We used tendamistat as a model protein to test the broadness of the free energy barriers. Tendamistat folds over two consecutive transition states and through a high-energy intermediate. Measuring the combined effect of temperature and denaturant on the position of the transition state in the wild-type protein and in several mutants revealed that the early transition state shows significant transition state movement. Its accessible surface area state becomes more native-like with destabilization of the native state by temperature. To the same extent, the entropy of the early transition state becomes more native-like with increasing denaturant concentration, in accordance with Hammond behavior. The position of the late transition state, in contrast, is much less sensitive to the applied perturbations. These results suggest that the barriers in protein folding become increasingly narrow as the folding polypeptide chain approaches the native state.

Correspondence between anomalous m‐ and Δ Cp‐values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.

Scattered Hammond plots reveal second level of site-specific information in protein folding: phi' (beta++).

Site-specific information about structural heterogeneities of the protein-folding transition-state ensemble is commonly derived from the scatter of the Brønsted plot through the individual values of = Delta logk(f)/Delta logK(D-N). Here, we provide a second level of site-specific detail in the transition-state analysis by demonstrating that the scatter of the Hammond plot is related to heterogeneities in the -value growth. That is, the extent of transition-state movement (Delta beta(++)) is proportional to the free-energy gradient of the mutational perturbation across the top of the activation barrier, '(beta(++)) proportional, variant Delta logK(D-N). The analysis is applied to the two-state protein L23 where the site-specific free-energy gradients are used to identify the interactions that show the highest degree of consolidation after crossing the barrier top. These interactions are distributed as a shell around the high- initiation point and denote the side-chain contacts that add criticality to the folding nucleus.

Approaching the speed limit for Greek Key β-barrel formation: transition-state movement tunes folding rate of zinc-substituted azurin

Azurin is a blue-copper protein with a β-barrel structure of Greek Key topology. In vitro, copper can be substituted with zinc without change in protein structure. We here analyze the kinetic folding behavior of zinc-substituted Pseudomonas aeruginosa azurin. Our findings can be summarized in three key conclusions: First, zinc remains strongly bound to the polypeptide upon unfolding, suggesting that the cofactor may bind to the protein before polypeptide folding in vivo. Second, the semi-logarithmic plot of folding and unfolding rates for zinc-substituted azurin as a function of denaturant concentration exhibits curvature due to a changing transition-state structure. Third, the extrapolated folding speed in water for zinc-substituted azurin is similar to that of other proteins with the same topology, implying that there is a speed limit that can be modulated by stability-driven transition-state movement for formation of β-barrel structures with Greek Key topology.

Hammond Behavior versus Ground State Effects in Protein Folding: Evidence for Narrow Free Energy Barriers and Residual Structure in Unfolded States

Apparent transition state movement upon mutation or changes in solvent conditions is frequently observed in protein folding and is often interpreted in terms of Hammond behavior. This led to the conclusion that barrier regions in protein folding are broad maxima on the free energy landscape. Here, we use the concept of self-interaction and cross-interaction parameters to test experimental data of 21 well-characterized proteins for Hammond behavior. This allows us to characterize the origin of transition state movements along different reaction coordinates. Only one of the 21 proteins shows a small but coherent transition state movement in agreement with the Hammond postulate. In most proteins the structure of the transition state is insensitive to changes in protein stability. The apparent change in the position of the transition state upon mutation, which is frequently observed in φ-value analysis, is in most cases due to ground-state effects caused by structural changes in the unfolded state. This argues for significant residual structure in unfolded polypeptide chains of many proteins. Disruption of these residual interactions by mutation often leads to decreased folding rates, which implies that these interactions are still present in the transition state. The failure to detect Hammond behavior shows that the free energy barriers encountered by a folding polypeptide chain are generally rather narrow and robust maxima for all experimentally explorable reaction coordinates.

Weak Cooperativity in the Core Causes a Switch in Folding Mechanism Between Two Proteins of the cks Family

The human protein ckshs1 (cks1) is a 79 residue α/β protein with low thermodynamic and kinetic stability. Its folding mechanism was probed by mutation at sites throughout the structure. Many of the mutations caused changes in the slope of the unfolding arm of the chevron plot. The effects can be rationalised in terms of either transition-state movement or native-state “breathing”, and in either case, the magnitude of the effect enables the sequence of events in the folding reaction to be determined. Those sites that fold early exhibit a small perturbation, whilst those sites that fold late exhibit a large perturbation. The results show that cks1 folds sequential pairs of β-strands first; β1/β2 and β3/β4. Subsequently, these pairs pack against each other and onto the α-helical region to form the core. The folding process of cks1 contrasts with that of the homologue, suc1. The 113 residue suc1 has the same β-sheet core structure but, additionally, two large insertions that confer much greater thermodynamic and kinetic stability. The more extensive network of tertiary interactions in suc1 provides sufficient enthalpic gain to overcome the entropic cost of forming the core and thus tips the balance in favour of non-local interactions: the non-local, central β-strand pair, β2/β4, forms first and the periphery strands pack on later. Moreover, the greater cooperativity of the core of suc1 protects its folding from perturbation and consequently the slope of the unfolding arm of the chevron plot is much less sensitive to mutation.

Non-linear rate-equilibrium free energy relationships and Hammond behavior in protein folding

Non-linear rate-equilibrium relationships upon mutation or changes in solvent conditions are frequently observed in protein folding reactions and are usually interpreted in terms of Hammond behavior. Here we first give a general overview over the concept of transition state movements in chemical reactions and discuss its application to protein folding. We then show examples for genuine Hammond behavior and for apparent transition state movements caused by other effects like changes in the rate-limiting step of the folding reaction or ground state effects, i.e. structural changes in either the native state or the unfolded state. These examples show that apparent transition state movements can easily be mistaken for Hammond behavior. We describe experimental tests using self- and cross-interaction parameters to distinguish between structural changes in a single transition state following Hammond behavior and apparent transition state movements caused by other effects.

Mimic of Photocycle by a Protein Folding Reaction in Photoactive Yellow Protein

The blue light receptor photoactive yellow protein (PYP) displays rhodopsin-like photochemistry based on the trans to cis photoisomerization of its p-coumaric acid chromophore. Here, we report that protein refolding from the acid-denatured state of PYP mimics the last photocycle transition in PYP. This implies a direct link between transient protein unfolding and photosensory signal transduction. We utilize this link to study general issues in protein folding. Chromophore trans to cis photoisomerization in the acid-denatured state strongly decelerates refolding, and converts the pH dependence of the barrier for refolding from linear to nonlinear. We propose transition state movement to explain this phenomenon. The cis chromophore significantly stabilizes the acid-denatured state, but acidification of PYP results in the accumulation of the acid-denatured state containing a trans chromophore. This provides a clear example of kinetic control in a protein unfolding reaction. These results demonstrate the power of PYP as a light-triggered model system to study protein folding.

Structural changes in the transition state of protein folding: alternative interpretations of curved chevron plots.

The interpretation of folding rates is often rationalized within the context of transition state theory. This means that the reaction rate is linked to an activation barrier, the height of which is determined by the free energy difference between a ground state (the starting point) and an apparent transition state. Changes in the folding kinetics are thus caused by effects on either the ground state, the transition state, or both. However, structural changes of the transition state are rarely discussed in connection with experimental data, and kinetic anomalies are commonly ascribed to ground state effects alone, e.g., depletion or accumulation of structural intermediates upon addition of denaturant. In this study, we present kinetic data which are best described by transition state changes. We also show that ground state effects and transition state effects are in general difficult to distinguish kinetically. The analysis is based on the structurally homologous proteins U1A and S6. Both proteins display two-state behavior, but there is a marked difference in their kinetics. S6 exhibits a classical V-shaped chevron plot (log observed rate constant vs denaturant concentration), whereas U1A's chevron plot is symmetrically curved, like an inverted bell curve. However, S6 is readily mutated to display U1A-like kinetics. The seemingly drastic effects of these mutations are readily ascribed to transition state movements where large kinetic differences result from relatively small alterations of a common free energy profile and broad activation barriers.