Relationship between the native-state hydrogen exchange and the folding pathways of barnase

Mechanical Folding and Unfolding of Protein Barnase at the Single-Molecule Level

The unfolding and folding of protein barnase has been extensively investigated in bulk conditions under the effect of denaturant and temperature. These experiments provided information about structural and kinetic features of both the native and the unfolded states of the protein, and debates about the possible existence of an intermediate state in the folding pathway have arisen. Here, we investigate the folding/unfolding reaction of protein barnase under the action of mechanical force at the single-molecule level using optical tweezers. We measure unfolding and folding force-dependent kinetic rates from pulling and passive experiments, respectively, and using Kramers-based theories (e.g., Bell-Evans and Dudko-Hummer-Szabo models), we extract the position of the transition state and the height of the kinetic barrier mediating unfolding and folding transitions, finding good agreement with previous bulk measurements. Measurements of the force-dependent kinetic barrier using the continuous effective barrier analysis show that protein barnase verifies the Leffler-Hammond postulate under applied force and allow us to extract its free energy of folding, ΔG0. The estimated value of ΔG0 is in agreement with our predictions obtained using fluctuation relations and previous bulk studies. To address the possible existence of an intermediate state on the folding pathway, we measure the power spectrum of force fluctuations at high temporal resolution (50 kHz) when the protein is either folded or unfolded and, additionally, we study the folding transition-path time at different forces. The finite bandwidth of our experimental setup sets the lifetime of potential intermediate states upon barnase folding/unfolding in the submillisecond timescale.

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 alpha-helix and two beta-strands (51 folded residues). Using multidimensional NMR, we solved the structure of this intermediate mimic to an atomic resolution (backbone rmsd, 0.51 A). 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.

The folding pathway of T4 lysozyme: an on-pathway hidden folding intermediate

T4 lysozyme has two easily distinguishable but energetically coupled domains: the N and C-terminal domains. In earlier studies, an amide hydrogen/deuterium exchange pulse-labeling experiment detected a stable submillisecond intermediate that accumulates before the rate-limiting transition state. It involves the formation of structures in both the N and C-terminal regions. However, a native-state hydrogen exchange experiment subsequently detected an equilibrium intermediate that only involves the formation of the C-terminal domain. Here, using stopped-flow circular dichroism and fluorescence, amide hydrogen exchange-folding competition, and protein engineering methods, we re-examined the folding pathway of T4-lysozyme. We found no evidence for the existence of a stable folding intermediate before the rate-limiting transition state at neutral pH. In addition, using native-state hydrogen exchange-directed protein engineering, we created a mimic of the equilibrium intermediate. We found that the intermediate mimic folds with the same rate as the wild-type protein, suggesting that the equilibrium intermediate is an on-pathway intermediate that exists after the rate-limiting transition state.

Determination of the folding transition states of barnase by using PhiI-value-restrained simulations validated by double mutant PhiIJ-values

The protein barnase folds from the denatured state into its native conformation via a high-energy intermediate. Using PhiI-values determined experimentally from single-point mutations as restraints in all-atom molecular dynamics simulations, we have determined ensembles of structures corresponding to the transition states for the formation of the folding intermediate and its conversion into the native state. We have also introduced a stringent validation of the approach used to calculate such structures by predicting interaction PhiIJ-values determined experimentally from a series of double-mutant cycles. The ensembles that we have obtained illustrate the heterogeneity in the nucleation-condensation process by which barnase folds. Obligatory steps of this process include the sequential formation of two folding nuclei, which correspond to the two main hydrophobic cores of the protein. Nonobligatory steps include the elongation of the strand beta1 and the formation of the helix alpha2. The results confirm that the use of experimental observables such as PhiI-values as restraints in molecular dynamics simulations is a powerful general strategy to characterize the relatively heterogeneous structural ensembles that populate nonnative regions of the energy landscapes of proteins.

A protein folding pathway with multiple folding intermediates at atomic resolution

Using native-state hydrogen-exchange-directed protein engineering and multidimensional NMR, we determined the high-resolution structure (rms deviation, 1.1 angstroms) for an intermediate of the four-helix bundle protein: Rd-apocytochrome b562. The intermediate has the N-terminal helix and a part of the C-terminal helix unfolded. In earlier studies, we also solved the structures of two other folding intermediates for the same protein: one with the N-terminal helix alone unfolded and the other with a reorganized hydrophobic core. Together, these structures provide a description of a protein folding pathway with multiple intermediates at atomic resolution. The two general features for the intermediates are (i) native-like backbone topology and (ii) nonnative side-chain interactions. These results have implications for important issues in protein folding studies, including large-scale conformation search, -value analysis, and computer simulations.

The Kinetic Pathway of Folding of Barnase

To search for folding intermediates, we have examined the folding and unfolding kinetics of wild-type barnase and four representative mutants under a wide range of conditions that span two-state and multi-state kinetics. The choice of mutants and conditions provided in-built controls for artifacts that might distort the interpretation of kinetics, such as the non-linearity of kinetic and equilibrium data with concentration of denaturant. We measured unfolding rate constants over a complete range of denaturant concentration by using by 1H/2H-exchange kinetics under conditions that favour folding, conventional stopped-flow methods at higher denaturant concentrations and continuous flow. Under conditions that favour multi-state kinetics, plots of the rate constants for unfolding against denaturant concentration fitted quantitatively to the equation for three-state kinetics, with a sigmoid component for a change of rate determining step, as did the refolding kinetics. The position of the transition state on the reaction pathway, as measured by solvent exposure (the Tanford beta value) also moved with denaturant concentration, fitting quantitatively to the same equations with a change of rate determining step. The sigmoid behaviour disappeared under conditions that favoured two-state kinetics. Those data combined with direct structural observations and simulation support a minimal reaction pathway for the folding of barnase that involves two detectable folding intermediates. The first intermediate, I(1), is the denatured state under physiological conditions, D(Phys), which has native-like topology, is lower in energy than the random-flight denatured state U and is suggested by molecular dynamics simulation of unfolding to be on-pathway. The second intermediate, I(2), is high energy, and is proven by the change in rate determining step in the unfolding kinetics to be on-pathway. The change in rate determining step in unfolding with structure or environment reflects the change in partitioning of this intermediate to products or starting materials.

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.

Hidden intermediates and levinthal paradox in the folding of small proteins

It has long been suggested that existence of partially folded intermediates may be essential for proteins to fold in a biologically meaningful time scale. Although partially folded intermediates have been commonly observed in larger proteins, they are generally not detectable in the kinetic folding of smaller proteins (approximately 100 amino acids or less). Recent native-state hydrogen exchange studies suggest that partially folded intermediates may exist behind the rate-limiting transition state in small proteins and evade detection by conventional kinetic methods.

Fast and slow intermediate accumulation and the initial barrier mechanism in protein folding

Do stable intermediates form very early in the protein folding process? New results and a quantity of literature that bear on this issue are examined here. Results available provide little support for early intermediate accumulation before an initial search-dependent nucleation barrier.

Lack of definable nucleation sites in the rate-limiting transition state of barnase under native conditions

It has been shown that the burst-phase (submillisecond) intermediate of barnase, if it exists, can be only marginally more stable than the fully unfolded state at pH 6.3 and 25 degrees C. In the study reported here, no stable burst-phase intermediate could be detected, even in the presence of stabilizing salt (0.4 M Na(2)SO(4)). These results suggest that a burst-phase intermediate with even marginal stability does not exist. The absence of such an intermediate in turn suggests the need for re-examination of the rate-limiting transition state (RLTS) under native conditions, which was previously characterized by using a three-state model with a stable intermediate and protein engineering. Surprisingly, mutations throughout the structure of barnase do not significantly affect the folding rate, suggesting a lack of specific favorable interactions among the side-chains in the RLTS. This RLTS is clearly different from that previously characterized under denaturing conditions, indicating that changes take place in the RLTS under native and denaturing conditions. The occurrence of such changes is further supported by the observation that the unfolding rate constants of barnase and its mutants were divergent or convergent as a function of denaturant concentrations. Consistent with changes in the RLTS, a re-analysis of data from native-state hydrogen exchange studies has shown that the logarithm of the unfolding rate constant inflects down under low concentrations of denaturant. Here, we discuss in detail the question of whether changes in the RLTS involve a kinetically silent intermediate that occurs after the initial RLTS.

Studies of biomolecular conformations and conformational dynamics by mass spectrometry

In the post-genomic era, a wealth of structural information has been amassed for proteins from NMR and crystallography. However, static protein structures alone are not a sufficient description: knowledge of the dynamic nature of proteins is essential to understand their wide range of functions and behavior during the life cycle from synthesis to degradation. Furthermore, few proteins have the ability to act alone in the crowded cellular environment. Assemblies of multiple proteins governed by complex signaling pathways are often required for the tasks of target recognition, binding, transport, and function. Mass spectrometry has emerged over the past several years as a powerful tool to address many of these questions. Recent improvements in "soft" ionization techniques have enabled researchers to study proteins and biomolecular complexes, both directly and indirectly. Likewise, continuous improvements in instrumental design in recent years have resulted in a dramatic expansion of the m/z range and resolution, enabling observation of large multi-protein assemblies whose structures are retained in the gas phase. In this article, we discuss some of the mass spectrometric techniques applied to investigate the nature of the conformations and dynamical properties that govern protein function.

On the nature of conformational openings: native and unfolded-state hydrogen and thiol-disulfide exchange studies of ferric aquomyoglobin

Native-state amide hydrogen exchange (HX) of proteins in the presence of denaturant has provided valuable details on the structures of equilibrium folding intermediates. Here, we extend HX theory to model thiol group exchange (SX) in single cysteine-containing variants of sperm whale ferric aquomyoglobin. SX is complementary to HX in that it monitors conformational opening events that expose side-chains, rather than the main chain, to solvent. A simple two-process model, consisting of EX2-limited local structural fluctuations and EX1-limited global unfolding, adequately accounts for all HX data. SX is described by the same model except at very low denaturant concentrations and when the bulky labeling reagent 5,5'-dithiobis (2-nitrobenzoic acid) (DTNB) is used. Under these conditions SX can occur by a novel denaturant-dependent process. This anomalous behavior is not observed when the smaller labeling reagent methyl methanethiosulfonate is employed, suggesting that it reflects a denaturant-induced increase in the amplitudes of local structural fluctuations. It also is not seen in heme-free apomyoglobin, which may indicate that local openings are sufficiently large in the absence of denaturant to allow DTNB unhindered access. Differences in SX kinetics obtained using the two labeling reagents provide estimates of the sizes of local opening reactions at different sites in the protein. At all sequence positions examined except for position 73, the same opening event appears to facilitate exchange of both backbone amide and side-chain thiol groups. The C73 thiol group is exposed by a low-energy fluctuation that does not expose its amide group to exchange.

A kinetically significant intermediate in the folding of barnase

A series of studies on the small protein barnase in the 1990s established it as a paradigm for protein folding in which there is a kinetically important intermediate. But, a recent study in PNAS claims that there are no stable intermediates on the folding pathway. I summarize the evidence that proves that the folding kinetics of barnase is inconsistent with the absence of a folding intermediate. I reinterpret the major evidence presented against the intermediate (an inflection in the unfolding limb of a chevron plot) and show that the inflection is precisely what is predicted from the energy diagram for a three-state reaction with a kinetically significant on-pathway intermediate. The inflection is indicative of a change of rate determining step from the formation to breakdown of an intermediate on unfolding. Other evidence presented against the intermediate is, in fact, consistent with a kinetically important intermediate. I show how the complexities in the kinetics provide a means for measuring otherwise unobtainable rate constants and provide a strategy for mapping the structure of the early transition state in folding. Rather than refute multistate kinetics, the presence of the inflection in the unfolding plot constitutes a novel type of evidence for on-pathway folding intermediates.

Intermolecular aggregations are responsible for the slow kinetics observed in the folding of cytochrome c at neutral pH

Folding of equine cytochrome c at a low protein concentration (26 microM) eliminated a slow kinetic phase (time constant three seconds) that was observed in the previous hydrogen exchange pulse-labeling experiments at pH 6.2 and 10 degrees C. It was demonstrated that this slow folding phase was caused by intermolecular aggregations. Because heterogeneous kinetics is a very general feature in the folding of proteins characterized by pulsed hydrogen exchange coupled with two-dimensional NMR, our experimental results suggest aggregations might also be responsible for the complex folding kinetics of other proteins. This is possible since these experiments were performed at relatively high protein concentrations.