The nature of the transition state ensemble and the mechanisms of protein folding: a review

Probing the Transition State in Enzyme Catalysis by High-Pressure NMR Dynamics

Protein conformational changes are frequently essential for enzyme catalysis, and in several cases, shown to be the limiting factor for overall catalytic speed. However, a structural understanding of corresponding transition states, needed to rationalize the kinetics, remains obscure due to their fleeting nature. Here, we determine the transition-state ensemble of the rate-limiting conformational transition in the enzyme adenylate kinase, by a synergistic approach between experimental high-pressure NMR relaxation during catalysis and molecular dynamics simulations. By comparing homologous kinases evolved under ambient or high pressure in the deep-sea, we detail transition state ensembles that differ in solvation as directly measured by the pressure dependence of catalysis. Capturing transition-state ensembles begins to complete the catalytic energy landscape that is generally characterized by structures of all intermediates and frequencies of transitions among them.

A critical comparison of coarse-grained structure-based approaches and atomic models of protein folding

Structure-based coarse-grained Gō-like models have been used extensively in deciphering protein folding mechanisms because of their simplicity and tractability. Meanwhile, explicit-solvent molecular dynamics (MD) simulations with physics-based all-atom force fields have been applied successfully to simulate folding/unfolding transitions for several small, fast-folding proteins. To explore the degree to which coarse-grained Gō-like models and their extensions to incorporate nonnative interactions are capable of producing folding processes similar to those in all-atom MD simulations, here we systematically compare the computed unfolded states, transition states, and transition paths obtained using coarse-grained models and all-atom explicit-solvent MD simulations. The conformations in the unfolded state in common Gō models are more extended, and are thus more in line with experiment, than those from all-atom MD simulations. Nevertheless, the structural features of transition states obtained by the two types of models are largely similar. In contrast, the folding transition paths are significantly more sensitive to modeling details. In particular, when common Gō-like models are augmented with nonnative interactions, the predicted dimensions of the unfolded conformations become similar to those computed using all-atom MD. With this connection, the large deviations of all-atom MD from simple diffusion theory are likely caused in part by the presence of significant nonnative effects in folding processes modelled by current atomic force fields. The ramifications of our findings to the application of coarse-grained modeling to more complex biomolecular systems are discussed.

Monitoring protein folding through high pressure NMR spectroscopy

High-pressure is a well-known perturbation method used to destabilize globular proteins. It is perfectly reversible, which is essential for a proper thermodynamic characterization of a protein equilibrium. In contrast to other perturbation methods such as heat or chemical denaturant that destabilize protein structures uniformly, pressure exerts local effects on regions or domains of a protein containing internal cavities. When combined with NMR spectroscopy, hydrostatic pressure offers the possibility to monitor at a residue level the structural transitions occurring upon unfolding and to determine the kinetic properties of the process. High-pressure NMR experiments can now be routinely performed, owing to the recent development of commercially available high-pressure sample cells. This review summarizes recent advances and some future directions of high-pressure NMR techniques for the characterization at atomic resolution of the energy landscape of protein folding.

High Pressure ZZ-Exchange NMR Reveals Key Features of Protein Folding Transition States

Understanding protein folding mechanisms and their sequence dependence requires the determination of residue-specific apparent kinetic rate constants for the folding and unfolding reactions. Conventional two-dimensional NMR, such as HSQC experiments, can provide residue-specific information for proteins. However, folding is generally too fast for such experiments. ZZ-exchange NMR spectroscopy allows determination of folding and unfolding rates on much faster time scales, yet even this regime is not fast enough for many protein folding reactions. The application of high hydrostatic pressure slows folding by orders of magnitude due to positive activation volumes for the folding reaction. We combined high pressure perturbation with ZZ-exchange spectroscopy on two autonomously folding protein domains derived from the ribosomal protein, L9. We obtained residue-specific apparent rates at 2500 bar for the N-terminal domain of L9 (NTL9), and rates at atmospheric pressure for a mutant of the C-terminal domain (CTL9) from pressure dependent ZZ-exchange measurements. Our results revealed that NTL9 folding is almost perfectly two-state, while small deviations from two-state behavior were observed for CTL9. Both domains exhibited large positive activation volumes for folding. The volumetric properties of these domains reveal that their transition states contain most of the internal solvent excluded voids that are found in the hydrophobic cores of the respective native states. These results demonstrate that by coupling it with high pressure, ZZ-exchange can be extended to investigate a large number of protein conformational transitions.

Influences of heterogeneous native contact energy and many-body interactions on the prediction of protein folding mechanisms

Combining heterogenous native contact energies and many-body interactions could improve the prediction of Brønsted plots using a structure-based model.

Fast-folding proteins under stress

Proteins are subject to a variety of stresses in biological organisms, including pressure and temperature, which are the easiest stresses to simulate by molecular dynamics. We discuss the effect of pressure and thermal stress on very-fast-folding model proteins, whose in vitro folding can be fully simulated on computers and compared with experiments. We then discuss experiments that can be used to subject proteins to low- and high-temperature unfolding, as well as low- and high-pressure unfolding. Pressure and temperature are prototypical perturbations that illustrate how close many proteins are to instability, a property that cells can exploit to control protein function. We conclude by reviewing some recent in-cell experiments, and progress being made in simulating and measuring protein stability and function inside live cells.

Theoretical volume profiles as a tool for probing transition states: folding kinetics

The mechanism by which conformational changes, particularly folding and unfolding, occur in proteins and other biopolymers has been widely discussed in the literature. Molecular dynamics (MD) simulations of protein folding present a formidable challenge since these conformational changes occur on a time scale much longer than what can be afforded at the current level of computational technology. Transition state (TS) theory offers a more economic description of kinetic properties of a reaction system by relating them to the properties of the TS, or for flexible systems, the TS ensemble (TSE). The application of TS theory to protein folding is limited by ambiguity in the definition of the TSE for this process. We propose to identify the TSE for conformational changes in flexible systems by comparison of its experimentally determined volumetric property, known as the volume of activation, to the structure-specific volume profile of the process calculated using MD. We illustrate this approach by its successful application to unfolding of a model chain system.

How quickly can a β-hairpin fold from its transition state?

Understanding the structural nature of the free energy bottleneck(s) encountered in protein folding is essential to elucidating the underlying dynamics and mechanism. For this reason, several techniques, including Φ-value analysis, have previously been developed to infer the structural characteristics of such high free-energy or transition states. Herein we propose that one (or few) appropriately placed backbone and/or side chain cross-linkers, such as disulfides, could be used to populate a thermodynamically accessible conformational state that mimics the folding transition state. Specifically, we test this hypothesis on a model β-hairpin, Trpzip4, as its folding mechanism has been extensively studied and is well understood. Our results show that cross-linking the two β-strands near the turn region increases the folding rate by an order of magnitude, to about (500 ns)⁻¹, whereas cross-linking the termini results in a hyperstable β-hairpin that has essentially the same folding rate as the uncross-linked peptide. Taken together, these findings suggest that cross-linking is not only a useful strategy to manipulate folding free energy barriers, as shown in other studies, but also, in some cases, it can be used to stabilize a folding transition state analogue and allow for direct assessment of the folding process on the downhill side of the free energy barrier. The calculated free energy landscape of the cross-linked Trpzip4 also supports this picture. An empirical analysis further suggests, when folding of β-hairpins does not involve a significant free energy barrier, the folding time (τ) follows a power law dependence on the number of hydrogen bonds to be formed (n_H), namely, τ = τ₀n_H^α, with τ₀ = 20 ns and α = 2.3.

Network approach to identify the folding transition states of peptides and proteins

Folding transition states and their structures are crucial in understanding protein folding pathways and folding dynamics. As they cannot be detected directly by experiments due to their instability, many computational methods have been proposed to solve this problem. However, each of these methods can give only one part of the transition state ensemble for a peptide or protein. Here we present a folding-network approach to identify the transition states of peptides or proteins and test it on the β-hairpin peptide trpzip2, with the result that we identify all the folding transition states of tripzip2, which may only be determined separately by other methods. This suggests that the network approach can provide more complete information about the folding transition states, at least for peptides.

Stable folding core in the folding transition state of an alpha-helical integral membrane protein

Defining the structural features of a transition state is important in understanding a folding reaction. Here, we use Φ-value and double mutant analyses to probe the folding transition state of the membrane protein bacteriorhodopsin. We focus on the final C-terminal helix, helix G, of this seven transmembrane helical protein. Φ-values could be derived for 12 amino acid residues in helix G, most of which have low or intermediate values, suggesting that native structure is disrupted at these amino acid positions in the transition state. Notably, a cluster of residues between E204 and M209 all have Φ-values close to zero. Disruption of helix G is further confirmed by a low Φ-value of 0.2 between residues T170 on helix F and S226 on helix G, suggesting the absence of a native hydrogen bond between helices F and G. Φ-values for paired mutations involved in four interhelical hydrogen bonds revealed that all but one of these bonds is absent in the transition state. The unstructured helix G contrasts with Φ-values along helix B that are generally high, implying native structure in helix B in the transition state. Thus helix B seems to constitute part of a stable folding nucleus while the consolidation of helix G is a relatively late folding event. Polarization of secondary structure correlates with sequence position, with a structured helix B near the N terminus contrasting with an unstructured C-terminal helix G.

Three-state thermal unfolding of onconase

Onconase is a member of the ribonuclease A superfamily currently in phase IIIb clinical trials as a treatment for malign mesothelioma due to its cytotoxic activity selective against tumor-cells. In this work, we have studied the equilibrium thermal unfolding of onconase using a combination of several structural and biophysical techniques. Our results indicate that at least one significantly populated intermediate, which implies the exposure of hydrophobic surface and significant changes in the environment around Trp3, occurs during the equilibrium unfolding process of this protein. The intermediate begins to populate at about 30° below the global unfolding temperature, reaching a maximum population of nearly 60%, 10° below the global unfolding temperature.

Constrained proper sampling of conformations of transition state ensemble of protein folding

Characterizing the conformations of protein in the transition state ensemble (TSE) is important for studying protein folding. A promising approach pioneered by Vendruscolo et al. [Nature (London) 409, 641 (2001)] to study TSE is to generate conformations that satisfy all constraints imposed by the experimentally measured φ values that provide information about the native likeness of the transition states. Faísca et al. [J. Chem. Phys. 129, 095108 (2008)] generated conformations of TSE based on the criterion that, starting from a TS conformation, the probabilities of folding and unfolding are about equal through Markov Chain Monte Carlo (MCMC) simulations. In this study, we use the technique of constrained sequential Monte Carlo method [Lin et al., J. Chem. Phys. 129, 094101 (2008); Zhang et al. Proteins 66, 61 (2007)] to generate TSE conformations of acylphosphatase of 98 residues that satisfy the φ-value constraints, as well as the criterion that each conformation has a folding probability of 0.5 by Monte Carlo simulations. We adopt a two stage process and first generate 5000 contact maps satisfying the φ-value constraints. Each contact map is then used to generate 1000 properly weighted conformations. After clustering similar conformations, we obtain a set of properly weighted samples of 4185 candidate clusters. Representative conformation of each of these cluster is then selected and 50 runs of Markov chain Monte Carlo (MCMC) simulation are carried using a regrowth move set. We then select a subset of 1501 conformations that have equal probabilities to fold and to unfold as the set of TSE. These 1501 samples characterize well the distribution of transition state ensemble conformations of acylphosphatase. Compared with previous studies, our approach can access much wider conformational space and can objectively generate conformations that satisfy the φ-value constraints and the criterion of 0.5 folding probability without bias. In contrast to previous studies, our results show that transition state conformations are very diverse and are far from nativelike when measured in cartesian root-mean-square deviation (cRMSD): the average cRMSD between TSE conformations and the native structure is 9.4 Å for this short protein, instead of 6 Å reported in previous studies. In addition, we found that the average fraction of native contacts in the TSE is 0.37, with enrichment in native-like β-sheets and a shortage of long range contacts, suggesting such contacts form at a later stage of folding. We further calculate the first passage time of folding of TSE conformations through calculation of physical time associated with the regrowth moves in MCMC simulation through mapping such moves to a Markovian state model, whose transition time was obtained by Langevin dynamics simulations. Our results indicate that despite the large structural diversity of the TSE, they are characterized by similar folding time. Our approach is general and can be used to study TSE in other macromolecules.

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.

Reaching the protein folding speed limit with large, sub-microsecond pressure jumps

Biomolecules are highly pressure-sensitive, but their dynamics upon return to ambient pressure are often too fast to observe with existing approaches. We describe a sample-efficient method capable of large and very fast pressure drops (<1 nanomole, >2,500 atmospheres and <0.7 microseconds). We validated the method by fluorescence-detected refolding of a genetically engineered lambda repressor mutant from its pressure-denatured state. We resolved barrierless structure formation upon return to ambient pressure; we observed a 2.1 +/- 0.7 microsecond refolding time, which is very close to the 'speed limit' for proteins and much faster than the corresponding temperature-jump refolding of the same protein. The ability to experimentally perform a large and very fast pressure drop opens up a new region of the biomolecular energy landscape for atomic-level simulation.

What have we learned from the studies of two-state folders, and what are the unanswered questions about two-state protein folding?

Small proteins with globular structures often fold by simple all-or-none mechanisms, both in an equilibrium and a kinetic sense, despite the very large number of partly folded conformations available. This type of 'two-state' folding will be discussed in terms of experimental tests, underlying molecular mechanisms, and limits to two-state behavior. Factors that appear to be important for two-state folding include topology (sequence distance of contacts in the native structure), molecular cooperativity and local energy distribution. Because their local stability distributions and cooperativities can be dissected and analyzed separately from topological features, recent studies of the folding of symmetric proteins will be discussed as a means to better understand the origins of two-state folding.

Towards a quantitative understanding of protein hydration and volumetric properties

Herein, we probe by pressure perturbation calorimetry (PPC) the coefficient of thermal expansion, the volumetric and the hydration properties of variants of a hyperstable variant of staphylococcal nuclease (SNase), Delta+PHS. The temperature-dependent volumetric properties of the folded and unfolded states of the wild-type protein are calculated with previously published data. The present PPC results are used to interpret the volume diagram and expansivity at a molecular level. We conclude that the expansivity of the unfolded state is, to a first approximation, temperature independent, while that of the folded state decreases with increasing temperature. Our data suggest that at low temperature the defining contribution to DeltaV comes mainly from excluded volume differences and DeltaV for unfolding is negative. In contrast, at high temperatures, differential solvation due to the increased exposed surface area of the unfolded state and, in particular, its larger thermal volume linked to the increased conformational dynamics of the unfolded state ensemble takes over and DeltaV for unfolding eventually becomes positive.

Practical approaches to protein folding and assembly: spectroscopic strategies in thermodynamics and kinetics

We describe here the use of several spectroscopies, such as fluorescence emission, circular dichroism, and differential quenching by acrylamide, in examining the equilibrium and kinetic folding of proteins. The first section regarding equilibrium techniques provides practical information for determining the conformational stability of a protein. In addition, several equilibrium-folding models are discussed, from two-state monomer to four-state homodimer, providing a comprehensive protocol for interpretation of folding curves. The second section focuses on the experimental design and interpretation of kinetic data, such as burst-phase analysis and exponential fits, used in elucidating kinetic folding pathways. In addition, simulation programs are used routinely to support folding models generated by kinetic experiments, and the fundamentals of simulations are covered.