Folding Transition State Ensemble Research Articles

Topological Reaction Coordinate Captures the Folding Transition State Ensemble in a Pierced Lasso Protein.

Proteins with a pierced lasso topology (PLT) have a covalent loop created by a disulfide bond, and the backbone circles back to thread the loop. This threaded topology has unique features compared to knotted topologies; notably, the topology is controlled by the chemical environment and the covalent loop remains intact even when denatured. In this work, we use the hormone leptin as our model PLT system and study its folding using molecular dynamics simulations that employ a structure-based (Go̅-like) model. We find that the reduced protein has a two-state folding mechanism with a transition state ensemble (TSE) that can be characterized by the reaction coordinate Q, the fraction of native contacts formed. In contrast, the oxidized protein, which must thread part of the polypeptide chain through a covalent loop, has a folding process that is poorly characterized by Q. Instead, we find that a topological coordinate that monitors the residue crossing the loop can identify the TSE of oxidized leptin. By precisely identifying the predicted TSE, one may now reliably calculate theoretical phi-values for the PLT protein, thereby enabling a comparison with experimental measurements. We find the loop-threading constraint leads to noncanonical phi-values that are uniformly small because this PLT protein has a flat energy landscape through the TSE.

Characterizing the Folding Transition-State Ensembles in the Energy Landscape of an RNA Tetraloop.

Molecular dynamics (MD) simulations have become increasingly powerful and can now describe the folding/unfolding of small biomolecules in atomic detail. However, a major challenge in MD simulations is to represent the complex energy landscape of biomolecules using a small number of reaction coordinates. In this study, we investigate the folding pathways of an RNA tetraloop, gcGCAAgc, using five classical MD simulations with a combined simulation time of approximately 120 μs. Our approach involves analyzing the tetraloop dynamics, including the folding transition state ensembles, using the energy landscape visualization method (ELViM). The ELViM is an approach that uses internal distances to compare any two conformations, allowing for a detailed description of the folding process without requiring root mean square alignment of structures. This method has previously been applied to describe the energy landscape of disordered β-amyloid peptides and other proteins. The ELViM results in a non-linear projection of the multidimensional space, providing a comprehensive representation of the tetraloop's energy landscape. Our results reveal four distinct transition-state regions and establish the paths that lead to the folded tetraloop structure. This detailed analysis of the tetraloop's folding process has important implications for understanding RNA folding, and the ELViM approach can be used to study other biomolecules.

Volumetric Properties of the Transition State Ensemble for Protein Folding.

Understanding how high hydrostatic pressure affects biomacromolecular interaction is important for deciphering the molecular mechanisms by which organisms adapt to live at the bottom of the ocean. The relative effect of hydrostatic pressure on the rates of folding/unfolding reactions is defined by the volumetric properties of the transition state ensemble relative to the folded and unfolded states. All-atom structure-based molecular dynamics simulations combined with quantitative computational protocol to compute volumes from three-dimensional coordinates allow volumetric mapping of protein folding landscape. This, is turn, provides qualitative understanding of the effects of hydrostatic pressure on energy landscape of proteins. The computational results for six different proteins are directly benchmark against experimental data and show an excellent agreement. Both experiments and computation show that the transition-state ensemble volume appears to be in-between the folded and unfolded state volumes, and thus the hydrostatic pressure accelerates protein unfolding.

Mutational Analysis of Protein Folding Transition States: Phi Values.

The analysis of protein folding reactions by monitoring the kinetic effects of specifically designed single-point mutations, the so-termed phi-value analysis, has been a favorite technique to experimentally probe the mechanisms of protein folding. The idea behind phi-value analysis is that the effects that mutations have on the folding and unfolding rate constants report on the energetic/structural features of the folding transition state ensemble (TSE), which is the highest point in the free energy surface connecting the native and unfolded states, and thus the rate limiting step that ultimately defines the folding mechanism. For single-domain, two-state folding proteins, the general procedure to perform the phi-value analysis of protein folding is relatively simple to implement in the lab. Once the mutations have been produced and purified, the researcher needs to follow a few specific guidelines to perform the experiments and to analyze the data so produced. In this chapter, a step-by-step description of how to measure and interpret the effects induced by site-directed mutations on the folding and unfolding rate constants of a protein of interest is provided. Some possible solutions to the most typical problems that arise when performing phi-value analysis in the lab are also provided.

Sampling the Folding Transition State Ensemble in a Tube-like Model of Protein

We used the tube model with Go-like potential for native contacts to study the folding transition of a designed three-helix bundle and a designed protein G-like structure. It is shown that both proteins in this model are two-state folders with a cooperative folding transition coincided with the collapse transition. We defined the transition states as protein conformations in a small region around the saddle point on a free energy surface with the energy and the conformationalroot mean square deviation (rmsd) from the native state as the coordinates. The transition state region on the free energy surface then was sampled by using umbrella sampling technique. We show that the transition state ensemble is broad consisting of different conformations that have different folded and unfolded elements.

Geometrical Frustration in Interleukin-33 Decouples the Dynamics of the Functional Element from the Folding Transition State Ensemble.

Interleukin-33 (IL-33) is currently the focus of multiple investigations into targeting pernicious inflammatory disorders. This mediator of inflammation plays a prevalent role in chronic disorders such as asthma, rheumatoid arthritis, and progressive heart disease. In order to better understand the possible link between the folding free energy landscape and functional regions in IL-33, a combined experimental and theoretical approach was applied. IL-33 is a pseudo- symmetrical protein composed of three distinct structural elements that complicate the folding mechanism due to competition for nucleation on the dominant folding route. Trefoil 1 constitutes the majority of the binding interface with the receptor whereas Trefoils 2 and 3 provide the stable scaffold to anchor Trefoil 1. We identified that IL-33 folds with a three-state mechanism, leading to a rollover in the refolding arm of its chevron plots in strongly native conditions. In addition, there is a second slower refolding phase that exhibits the same rollover suggesting similar limitations in folding along parallel routes. Characterization of the intermediate state and the rate limiting steps required for folding suggests that the rollover is attributable to a moving transition state, shifting from a post- to pre-intermediate transition state as you move from strongly native conditions to the midpoint of the transition. On a structural level, we found that initially, all independent Trefoil units fold equally well until a QCA of 0.35 when Trefoil 1 will backtrack in order to allow Trefoils 2 and 3 to fold in the intermediate state, creating a stable scaffold for Trefoil 1 to fold onto during the final folding transition. The formation of this intermediate state and subsequent moving transition state is a result of balancing the difficulty in folding the functionally important Trefoil 1 onto the remainder of the protein. Taken together our results indicate that the functional element of the protein is geometrically frustrated, requiring the more stable elements to fold first, acting as a scaffold for docking of the functional element to allow productive folding to the native state.

Driving Forces in Pressure-Induced Protein Transitions.

The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when exposed to water in an unfolded conformation. This chapter presents a volumetric analysis of the peptide group and the 20 naturally occurring amino acid side chains in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and in the unfolded conformation modeled by low-molecular analogs of proteins. The transfer of a peptide group from the protein interior to water becomes increasingly favorable as pressure increases. This observation classifies solvation of peptide groups as a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. An inference can be drawn that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures this core, owing to the absence of structural constraints, may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced protein denaturation. Volumetric data presented here have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensembles.

Effects of Energetic Heterogeneity on Protein Folding Dynamics Across Many Non-Homologous Proteins

Native structure based models (SBMs, “Go-models”) are based on energy landscape theory and the principle of minimal frustration [1,2]. Using this framework strongly reduces computational demands and allows us to investigate the influence of the contact-energy weight onto the folding process for many different protein families. A large set (approx. 200) of non-homologous monomeric proteins sized from 50-150 amino acids are simulated on a coarse-grained level, representing each amino acid by a single bead. A fully automatized workflow implemented with the help of eSBMTools [3] allows to quantify typical folding properties like phi-values, folding free energy landscapes and transition state ensembles (TSEs) for the simulated proteins.Conventional SBMs use sequence independent homogeneous energy weights for the native contact energy. We compare these “vanilla” models with sequence dependent heterogeneous “flavored” energy weights as introduced by Miyazawa and Jernigan [4]. Subsequently, the dynamics of forming TSEs are compared by the sequential formation of tertiary interfaces between secondary structure elements in vanilla and flavored simulations. We find that, despite the energetic dissimilarity, the sequential ordering is similar between both types of simulations and appears as a robust property. This suggests that protein folding dynamics are strongly influenced by the native protein topology.1. Onuchic, J.N. and P.G. Wolynes, Theory of protein folding. Curr. Opin. Struct. Bio., 2004. 14(1): p. 70-75.2. Schug, A. and J.N. Onuchic, Current opinion in pharmacology, 2010. 10(6): p. 709-14.3. Lutz, B., C. Sinner, G. Heuermann, A. Verma, and A. Schug, eSBMTools1.0: enhanced native structure-based modeling tools. Bioinformatics, 2013. doi:10.1093/bioinformatics/btt478.4. Miyazawa, S. and R.L. Jernigan, Residue-residue potentials with a favorable contact pair term and an unfavorable high packing density term, for simulation and threading. Journal of Molecular Biology, 1996. 256(3): p. 623-44.

Analyzing Protein Folding by High-throughput Simulations

Molecular Dynamics allows investigating the dynamical properties of biomolecules. Protein folding simulations are computationally challenging when simulating all involved (solvent) atoms considering timescales of ms or slower. Native structure based models (SBM,‘Go-models') reduce computational complexity and have been proven to be a robust and efficient way for exploring the protein folding process (Schug and Onuchic 2010; Thirumalai, O'Brien et al. 2010). They are based on energy landscape theory and the principle of minimal frustration. Using this framework, we simulate protein folding for a large set (∼ 200) of non-homologous monomeric proteins sized from 50-150 amino acids in coarse-grained simulations, representing each amino acid by a single bead. A fully automatized workflow implemented with the help of eSBM Tools (Lutz et al.) guides these simulations. From the simulations, we extract typical folding properties like phi-values, folding free energy landscape and transition state ensembles. We repeat the simulations for a variant SBM with flavored contact strengths pending on amino acids composition. The resulting database estimates the robustness of folding parameters, quantifies the folding behavior, compares the behavior to existing experimental data and can serve as a baseline for comparison to future experiments or simulations of protein folding.Schug, A. and J. N. Onuchic (2010). “From protein folding to protein function and biomolecular binding by energy landscape theory.” Curr Opin Pharmacol 10(6): 709-714.Thirumalai, D., E. P. O'Brien, et al. (2010). “Theoretical perspectives on protein folding.” Annual review of biophysics 39: 159-183.Lutz, B., et al. (2012), in Proceedings of the 11th International Conference on Modeling and Applied Simulation, Wien, 2012, edited by M. Affenzeller, et al., pp. 237.

Equilibrium and kinetic studies of protein cooperativity using urea-induced folding/unfolding of a Ubq–UIM fusion protein

Understanding the origins of cooperativity in proteins remains an important topic in protein folding. This study describes experimental folding/unfolding equilibrium and kinetic studies of the engineered protein Ubq–UIM, consisting of ubiquitin (Ubq) fused to the sequence of the ubiquitin interacting motif (UIM) via a short linker. Urea-induced folding/unfolding profiles of Ubq–UIM were monitored by far-UV circular dichroism and fluorescence spectroscopies and compared to those of the isolated Ubq domain. It was found that the equilibrium data for Ubq–UIM is inconsistent with a two-state model. Analysis of the kinetics of folding shows similarity in the folding transition state ensemble between Ubq and Ubq–UIM, suggesting that formation of Ubq domain is independent of UIM. The major contribution to the stabilization of Ubq–UIM, relative to Ubq, was found to be in the rates of unfolding. Moreover, it was found that the kinetic m-values for Ubq–UIM unfolding, monitored by different probes (far-UV circular dichroism and fluorescence spectroscopies), are different; thereby, further supporting deviations from a two-state behavior. A thermodynamic linkage model that involves four states was found to be applicable to the urea-induced unfolding of Ubq–UIM, which is in agreement with the previous temperature-induced unfolding study. The applicability of the model was further supported by site-directed variants of Ubq–UIM that have altered stabilities of Ubq/UIM interface and/or stabilities of individual Ubq- and UIM-domains. All variants show increased cooperativity and one variant, E43N_Ubq–UIM, appears to behave very close to an equilibrium two-state.

The Effects of pK a Tuning on the Thermodynamics and Kinetics of Folding: Design of a Solvent-Shielded Carboxylate Pair at the a-Position of a Coiled-Coil

The tuning of the pK a of ionizable residues plays a critical role in various protein functions, such as ligand-binding, catalysis, and allostery. Proteins harness the free energy of folding to position ionizable groups in highly specific environments that strongly affect their pK a values. To investigate the interplay among protein folding kinetics, thermodynamics, and pK a modulation, we introduced a pair of Asp residues at neighboring interior positions of a coiled-coil. A single Asp residue was replaced for an Asn side chain at the a-position of the coiled-coil from GCN4, which was also crosslinked at the C-terminus via a flexible disulfide bond. The thermodynamic and kinetic stability of the system was measured by circular dichroism and stopped-flow fluorescence as a function of pH and concentration of guanidine HCl. Both sets of data are consistent with a two-state equilibrium between fully folded and unfolded forms. Distinct pK a values of 6.3 and 5.35 are assigned to the first and second protonation of the Asp pair; together they represent an energetic difference of 5 kcal/mol relative to the protonation of two Asp residues with unperturbed pK a values. Analysis of the rate data as a function of pH and denaturant concentration allowed calculation of the kinetic constants for the conformational transitions of the peptide with the Asp residues in the doubly protonated, singly protonated, and unprotonated forms. The doubly and singly protonated forms fold rapidly, and a ϕ-value analysis shows that their contribution to folding occurs subsequent to the transition state ensemble for folding. By contrast, the doubly charged state shows a reduced rate of folding and a ϕ-value near 0.5 indicative of a repulsive interaction, and possibly also heterogeneity in the transition state ensemble.

Unique Features of the Folding Landscape of a Repeat Protein Revealed by Pressure Perturbation

The volumetric properties of proteins yield information about the changes in packing and hydration between various states along the folding reaction coordinate and are also intimately linked to the energetics and dynamics of these conformations. These volumetric characteristics can be accessed via pressure perturbation methods. In this work, we report high-pressure unfolding studies of the ankyrin domain of the Notch receptor (Nank1–7) using fluorescence, small-angle x-ray scattering, and Fourier transform infrared spectroscopy. Both equilibrium and pressure-jump kinetic fluorescence experiments were consistent with a simple two-state folding/unfolding transition under pressure, with a rather small volume change for unfolding compared to proteins of similar molecular weight. High-pressure fluorescence, Fourier transform infrared spectroscopy, and small-angle x-ray scattering measurements revealed that increasing urea over a very small range leads to a more expanded pressure unfolded state with a significant decrease in helical content. These observations underscore the conformational diversity of the unfolded-state basin. The temperature dependence of pressure-jump fluorescence relaxation measurements demonstrated that at low temperatures, the folding transition state ensemble (TSE) lies close in volume to the folded state, consistent with significant dehydration at the barrier. In contrast, the thermal expansivity of the TSE was found to be equivalent to that of the unfolded state, indicating that the interactions that constrain the folded-state thermal expansivity have not been established at the folding barrier. This behavior reveals a high degree of plasticity of the TSE of Nank1–7.

Origins of Pressure-Induced Protein Transitions

The molecular mechanisms underlying pressure-induced protein denaturation can be analyzed based on the pressure-dependent differences in the apparent volume occupied by amino acids inside the protein and when they are exposed to water in an unfolded conformation. We present here an analysis for the peptide group and the 20 naturally occurring amino acid side chains based on volumetric parameters for the amino acids in the interior of the native state, the micelle-like interior of the pressure-induced denatured state, and the unfolded conformation modeled by N-acetyl amino acid amides. The transfer of peptide groups from the protein interior to water becomes increasingly favorable as pressure increases. Thus, solvation of peptide groups represents a major driving force in pressure-induced protein denaturation. Polar side chains do not appear to exhibit significant pressure-dependent changes in their preference for the protein interior or solvent. The transfer of nonpolar side chains from the protein interior to water becomes more unfavorable as pressure increases. We conclude that a sizeable population of nonpolar side chains remains buried inside a solvent-inaccessible core of the pressure-induced denatured state. At elevated pressures, this core may become packed almost as tightly as the interior of the native state. The presence and partial disappearance of large intraglobular voids is another driving force facilitating pressure-induced denaturation of individual proteins. Our data also have implications for the kinetics of protein folding and shed light on the nature of the folding transition state ensemble.

Desolvation Barrier Effects Are a Likely Contributor to the Remarkable Diversity in the Folding Rates of Small Proteins

The variation in folding rate among single-domain natural proteins is tremendous, but common models with explicit representations of the protein chain are either demonstrably insufficient or unclear as to their capability for rationalizing the experimental diversity in folding rates. In view of the critical role of water exclusion in cooperative folding, we apply native-centric, coarse-grained chain modeling with elementary desolvation barriers to investigate solvation effects on folding rates. For a set of 13 proteins, folding rates simulated with desolvation barriers cover ∼4.6 orders of magnitude, spanning a range essentially identical to that observed experimentally. In contrast, folding rates simulated without desolvation barriers cover only ∼2.2 orders of magnitude. Following a Hammond-like trend, the folding transition-state ensemble (TSE) of a protein model with desolvation barriers generally has a higher average number of native contacts and is structurally more specific, that is, less diffused, than the TSE of the corresponding model without desolvation barriers. Folding is generally significantly slower in models with desolvation barriers because of their higher overall macroscopic folding barriers as well as slower conformational diffusion speeds in the TSE that are ≈1/50 times those in models without desolvation barriers. Nonetheless, the average root-mean-square deviation between the TSE and the native conformation is often similar in the two modeling approaches, a finding suggestive of a more robust structural requirement for the folding rate-limiting step. The increased folding rate diversity in models with desolvation barriers originates from the tendency of these microscopic barriers to cause more heightening of the overall macroscopic folding free-energy barriers for proteins with more nonlocal native contacts than those with fewer such contacts. Thus, the enhancement of folding cooperativity by solvation effects is seen as positively correlated with a protein's native topological complexity.

Identification of the protein folding transition state from molecular dynamics trajectories

The rate of protein folding is governed by the transition state so that a detailed characterization of its structure is essential for understanding the folding process. In vitro experiments have provided a coarse-grained description of the folding transition state ensemble (TSE) of small proteins. Atomistic details could be obtained by molecular dynamics (MD) simulations but it is not straightforward to extract the TSE directly from the MD trajectories, even for small peptides. Here, the structures in the TSE are isolated by the cut-based free-energy profile (cFEP) using the network whose nodes and links are configurations sampled by MD and direct transitions among them, respectively. The cFEP is a barrier-preserving projection that does not require arbitrarily chosen progress variables. First, a simple two-dimensional free-energy surface is used to illustrate the successful determination of the TSE by the cFEP approach and to explain the difficulty in defining boundary conditions of the Markov state model for an entropically stabilized free-energy minimum. The cFEP is then used to extract the TSE of a beta-sheet peptide with a complex free-energy surface containing multiple basins and an entropic region. In contrast, Markov state models with boundary conditions defined by projected variables and conventional histogram-based free-energy profiles are not able to identify the TSE of the beta-sheet peptide.

An Error Analysis for Two-State Protein-Folding Kinetic Parameters and φ-Values: Progress toward Precision by Exploring pH Dependencies on Leffler Plots

The interpretation of φ-values has led to an understanding of the folding transition state ensemble of a variety of proteins. Although the main guidelines and equations for calculating φ are well established, there remains some controversy about the quality of the numerical values obtained. By analyzing a complete set of results from kinetic experiments with the SH3 domain of α-spectrin (Spc-SH3) and applying classical error methods and error-propagation formulas, we evaluated the uncertainties involved in two-state-folding kinetic experimental parameters and the corresponding calculated φ-values. We show that kinetic constants in water and m values can be properly estimated from a judicious weighting of fitting errors and describe some procedures to calculate the errors in Gibbs energies and φ-values from a traditional two-point Leffler analysis. Furthermore, on the basis of general assumptions made with the protein engineering method, we show how to generate multipoint Leffler plots via the analysis of pH dependencies of kinetic parameters. We calculated the definitive φ-values for a collection of single mutations previously designed to characterize the folding transition state of the α-spectrin SH3 domain. The effectiveness of the pH-scanning procedure is also discussed in the context of error analysis. Judging from the magnitudes of the error bars obtained from two-point and multipoint Leffler plots, we conclude that the precision obtained for φ-values should be ∼25%, a reasonable limit that takes into account the propagation of experimental errors.

Vi -Value Analysis: A Pressure-Based Method for Mapping the Folding Transition State Ensemble of Proteins

Here we introduce an approach to mapping the folding transition state ensemble of proteins based on the pressure dependence of protein stability. Previously, we have shown that the activation volume for folding of wild type (WT) SNase is large and positive, and hence that the rate-limiting step in folding involves significant dehydration. In contrast, variants bearing buried ionizable residues at position 66 were shown recently to fold through a highly hydrated transition state ensemble (TSE). We present the effects on the pressure-jump folding kinetics of Lys substitutions in different internal positions throughout the structure. We calculate the Vi value of the variants as the activation volume for folding relative to that of the wild type. We find that the structure of the SNase WT includes part of the β-barrel and part of the first α-helix. The unique advantage of Vi-value analysis is that it conveys direct information about the state of hydration of the TSE, which has been recognized as a key factor in the protein folding transition.

Kinetic Definition of Protein Folding Transition State Ensembles and Reaction Coordinates

Using distributed molecular dynamics simulations we located four distinct folding transitions for a 39-residue ββαβ protein fold. To characterize the nature of each room temperature transition, we calculated the probability of transmission for 500 points along each free energy barrier. We introduced a method for determining transition states by employing the transmission probability, P trans, and determined which conformations were transition state ensemble members ( P trans ≈ 0.5). The transmission probability may be used to characterize the barrier in several ways. For example, we ran simulations at 82°C, determined the change in P trans with temperature for all 2,000 conformations, and quantified Hammond behavior directly using P trans correlation. Additionally, we propose that diffusion along P trans may provide the configurational diffusion rate at the top of the barrier. Specifically, given a transition state conformation x 0 with estimated P trans = 0.5, we selected a large set of subsequent conformations from independent trajectories, each exactly a small time δt after x 0 (250 ps). Calculating P trans for the new trial conformations, we generated the P( P trans| δt = 250 ps) distribution that reflected diffusion. This approach provides a novel perspective on the diffusive nature of a protein folding transition and provides a framework for a quantitative study of activated relaxation kinetics.

φ-Value Analysis of Apo-Azurin Folding: Comparison between Experiment and Theory

Pseudomonas aeruginosa azurin is a 128-residue beta-sandwich metalloprotein; in vitro kinetic experiments have shown that it folds in a two-state reaction. Here, we used a variational free energy functional to calculate the characteristics of the transition state ensemble (TSE) for folding of the apo-form of P. aeruginosa azurin and investigate how it responds to thermal and mutational changes. The variational method directly yields predicted chevron plots for wild-type and mutant apo-forms of azurin. In parallel, we performed in vitro kinetic-folding experiments on the same set of azurin variants using chemical perturbation. Like the wild-type protein, all apo-variants fold in apparent two-state reactions both in calculations and in stopped-flow mixing experiments. Comparisons of phi (phi) values determined from the experimental and theoretical chevron parameters reveal an excellent agreement for most positions, indicating a polarized, highly structured TSE for folding of P. aeruginosa apo-azurin. We also demonstrate that careful analysis of side-chain interactions is necessary for appropriate theoretical description of core mutants.

Hydration of the Folding Transition State Ensemble of a Protein

A complete description of the mechanisms of protein folding requires knowledge of the structural and physical character of the folding transition state ensembles (TSEs). A key question concerning the role of hydration of the hydrophobic core in determining folding mechanisms remains. To address this, we probed the state of hydration of the TSE of staphylococcal nuclease (SNase) by examining the fluorescence-detected pressure-jump relaxation behavior of six SNase variants in which a residue in the hydrophobic core, Val-66, was replaced with polar or ionizable residues (Lys, Arg, His, Asp, Glu, and Asn). Because of a large positive activation volume for folding, the major effect of pressure on the wild-type protein is to decrease the folding rate. By the time wild-type SNase reaches the folding transition state, most water has already been expelled from its hydrophobic core. In contrast, the major effect of pressure on the variant proteins is an increase in the unfolding rate due to a large negative activation volume for unfolding. This results from a significant increase in the level of hydration of the TSE when an internal ionizable group is present. These data confirm that the role of water in the folding reaction can differ from protein to protein and that even a single substitution in a critical position can modulate significantly the properties of the TSE.