Origin of apparent fast and non-exponential kinetics of lysozyme folding measured in pulsed hydrogen exchange experiments

Accurate prediction of protein folding mechanisms by simple structure-based statistical mechanical models

Recent breakthroughs in highly accurate protein structure prediction using deep neural networks have made considerable progress in solving the structure prediction component of the 'protein folding problem'. However, predicting detailed mechanisms of how proteins fold into specific native structures remains challenging, especially for multidomain proteins constituting most of the proteomes. Here, we develop a simple structure-based statistical mechanical model that introduces nonlocal interactions driving the folding of multidomain proteins. Our model successfully predicts protein folding processes consistent with experiments, without the limitations of protein size and shape. Furthermore, slight modifications of the model allow prediction of disulfide-oxidative and disulfide-intact protein folding. These predictions depict details of the folding processes beyond reproducing experimental results and provide a rationale for the folding mechanisms. Thus, our physics-based models enable accurate prediction of protein folding mechanisms with low computational complexity, paving the way for solving the folding process component of the 'protein folding problem'.

Mapping fast protein folding with multiple-site fluorescent probes

Fast protein folding involves complex dynamics in many degrees of freedom, yet microsecond folding experiments provide only low-resolution structural information. We enhance the structural resolution of the five-helix bundle protein λ6-85 by engineering into it three fluorescent tryptophan-tyrosine contact probes. The probes report on distances between three different helix pairs: 1-2, 1-3, and 3-2. Temperature jump relaxation experiments on these three mutants reveal two different kinetic timescales: a slower timescale for 1-3 and a faster one for the two contacts involving helix 2. We hypothesize that these differences arise from a single folding mechanism that forms contacts on different timescales, and not from changes of mechanism due to adding the probes. To test this hypothesis, we analyzed the corresponding three distances in one published single-trajectory all-atom molecular-dynamics simulation of a similar mutant. Autocorrelation analysis of the trajectory reveals the same "slow" and "fast" distance change as does experiment, but on a faster timescale; smoothing the trajectory in time shows that this ordering is robust and persists into the microsecond folding timescale. Structural investigation of the all-atom computational data suggests that helix 2 misfolds to produce a short-lived off-pathway trap, in agreement with the experimental finding that the 1-2 and 3-2 distances involving helix 2 contacts form a kinetic grouping distinct from 1 to 3. Our work demonstrates that comparison between experiment and simulation can be extended to several order parameters, providing a stronger mechanistic test.

Submillisecond protein folding events monitored by rapid mixing and mass spectrometry-based oxidative labeling

Kinetic measurements can provide insights into protein folding mechanisms. However, the initial (submillisecond) stages of folding still represent a formidable analytical challenge. A number of ultrarapid triggering techniques have been available for some time, but coupling of these techniques with detection methods that are capable of providing detailed structural information has proven to be difficult. The current work addresses this issue by combining submillisecond mixing with laser-induced oxidative labeling. Apomyoglobin (aMb) serves as a model system for our measurements. Exposure of the protein to a brief pulse of hydroxyl radical (·OH) at different time points during folding introduces covalent modifications at solvent accessible side chains. The extent of labeling is monitored using mass spectrometry-based peptide mapping, providing spatially resolved measurements of changes in solvent accessibility. The submillisecond mixer used here improves the time resolution by a factor of 50 compared to earlier ·OH labeling experiments from our laboratory. Data obtained in this way indicate that early aMb folding events are driven by both local and sequence-remote docking of hydrophobic side chains. Assembly of a partially formed A(E)G(H) scaffold after 0.2 ms is followed by stepwise consolidation that ultimately yields the native state. Major conformational changes go to completion within 0.1 s. The technique introduced here is capable of providing in-depth structural information on very short time scales that have thus far been dominated by low resolution (global) spectroscopic probes. By employing submillisecond mixing in conjunction with slower mixing techniques, it is possible to observe complete folding pathways, from fractions of a millisecond all the way to minutes.

Pulsed hydrogen/deuterium exchange mass spectrometry for time-resolved membrane protein folding studies

Kinetic folding experiments by pulsed hydrogen/deuterium exchange (HDX) mass spectrometry (MS) are a well-established tool for water-soluble proteins. To the best of our knowledge, the current study is the first that applies this approach to an integral membrane protein. The native state of bacteriorhodopsin (BR) comprises seven transmembrane helices and a covalently bound retinal cofactor. BR exposure to sodium dodecyl sulfate (SDS) induces partial unfolding and retinal loss. We employ a custom-built three-stage mixing device for pulsed-HDX/MS investigations of BR refolding. The reaction is triggered by mixing SDS-denatured protein with bicelles. After a variable folding time (10 ms to 24 h), the protein is exposed to excess D(2) O buffer under rapid exchange conditions. The HDX pulse is terminated by acid quenching after 24 ms. Subsequent off-line analysis is performed by size exclusion chromatography and electrospray MS. These measurements yield the number of protected backbone N-H sites as a function of folding time, reflecting the recovery of secondary structure. Our results indicate that much of the BR secondary structure is formed quite late during the reaction, on a time scale of 10 s and beyond. It is hoped that in the future it will be possible to extend the pulsed-HDX/MS approach employed here to membrane proteins other than BR.

Characterizing short-lived protein folding intermediates by top-down hydrogen exchange mass spectrometry

This work combines pulsed hydrogen/deuterium exchange (HDX) and top-down mass spectrometry for the structural characterization of short-lived protein folding intermediates. A custom-built flow device with three sequential mixing steps is used for (i) triggering protein folding, (ii) pulsed D(2)O labeling, and (iii) acid quenching. The earliest folding time point that can be studied with this system is 10 ms. The mixing device was coupled online to the electrospray source of a Fourier transform mass spectrometer, where intact protein ions are fragmented by electron capture dissociation (ECD). The viability of this experimental strategy is demonstrated by applying it to the refolding of horse apo-myoglobin (aMb), a reaction known to involve a transient intermediate. Cooling of the mixing device to 0 °C reduces the reaction rate such that the folding process occurs within the experimentally accessible time window. Top-down ECD provides an average spatial resolution of ca. 2 residues, surpassing the resolution typically achieved in traditional proteolytic digestion/HDX studies. Amide back exchange is virtually eliminated by the short (∼1 s) duration of the acid quenching step. The aMb folding intermediate exhibits HDX protection in helices G and H, whereas the remainder of the protein is largely unfolded. Marginal protection is seen for helix A. Overall, the top-down ECD approach used here offers insights into the sequence of events leading from the unfolded state to the native conformation, with envisioned future applications in the areas of protein misfolding and aggregation. The time-resolved experiments reported herein represent an extension of our previous work, where HDX/MS with top-down ECD was employed for monitoring "static" protein structures under equilibrium conditions (Pan et al. J. Am. Chem. Soc. 2009, 131, 12801).

Mass spectrometry combined with oxidative labeling for exploring protein structure and folding

This review discusses various mass spectrometry (MS)-based approaches for exploring structural aspects of proteins in solution. Electrospray ionization (ESI)-MS, in particular, has found fascinating applications in this area. For example, when used in conjunction with solution-phase hydrogen/deuterium exchange (HDX), ESI-MS is a highly sensitive tool for probing conformational dynamics. The main focus of this article is a technique that is complementary to HDX, that is, the covalent labeling of proteins by hydroxyl radicals. The reactivity of individual amino acid side chains with *OH is strongly affected by their degree of solvent exposure. Thus, analysis of the oxidative labeling pattern by peptide mapping and tandem mass spectrometry provides detailed structural information. A convenient method for *OH production is the photolysis of H(2)O(2) by a pulsed UV laser, resulting in oxidative labeling on the microsecond time scale. Selected examples demonstrate the use of this technique for structural studies on membrane proteins, and the combination with rapid mixing devices for characterizing the properties of short-lived protein (un)folding intermediates.

Time-dependent changes in side-chain solvent accessibility during cytochrome c folding probed by pulsed oxidative labeling and mass spectrometry

The current work employs a novel approach for characterizing structural changes during the refolding of acid-denatured cytochrome c (cyt c). At various time points (ranging from 10 ms to 5 min) after a pH jump from 2 to 7, the protein is exposed to a microsecond hydroxyl radical (.OH) pulse that induces oxidative labeling of solvent-exposed side chains. Most of the covalent modifications appear as +16-Da adducts that are readily detectable by mass spectrometry. The overall extent of labeling decreases as folding proceeds, reflecting dramatic changes in the accessibility of numerous residues. Peptide mapping and tandem mass spectrometry reveal that the side chains of C14, C17, H33, F46, Y48, W59, M65, Y67, Y74, M80, I81, and Y97 are among the dominant sites of oxidation. Temporal changes in the accessibility of these residues are consistent with docking of the N- and C-terminal helices as early as 10 ms. However, structural reorganization at the helix interface takes place up to at least 1 s. Initial misligation of the heme iron by H33 leads to distal crowding, giving rise to low solvent accessibility of the displaced (native) M80 ligand and the adjacent I81. W59 retains a surprisingly high level of accessibility long into the folding process, indicating the presence of packing defects in the hydrophobically collapsed core. Overall, the results of this work are consistent with previous hydrogen/deuterium exchange studies that proposed a foldon-mediated mechanism. The structural data obtained by .OH labeling monitor the packing and burial of side chains, whereas hydrogen/deuterium exchange primarily monitors the formation of secondary structure elements. Hence, the two approaches yield complementary information. Considering the very short time scale of pulsed oxidative labeling, an extension of the approach used here to sub-millisecond folding studies should be feasible.

Slow formation of aggregation-resistant beta-sheet folding intermediates

Protein folding has been studied extensively for decades, yet our ability to predict how proteins reach their native state from a mechanistic perspective is still rudimentary at best, limiting our understanding of folding-related processes in vivo and our ability to manipulate proteins in vitro. Here, we investigate the in vitro refolding mechanism of a large beta-helix protein, pertactin, which has an extended, elongated shape. At 55 kDa, this single domain, all-beta-sheet protein allows detailed analysis of the formation of beta-sheet structure in larger proteins. Using a combination of fluorescence and far-UV circular dichroism spectroscopy, we show that the pertactin beta-helix refolds remarkably slowly, with multiexponential kinetics. Surprisingly, despite the slow refolding rates, large size, and beta-sheet-rich topology, pertactin refolding is reversible and not complicated by off-pathway aggregation. The slow pertactin refolding rate is not limited by proline isomerization, and 30% of secondary structure formation occurs within the rate-limiting step. Furthermore, site-specific labeling experiments indicate that the beta-helix refolds in a multistep but concerted process involving the entire protein, rather than via initial formation of the stable core substructure observed in equilibrium titrations. Hence pertactin provides a valuable system for studying the refolding properties of larger, beta-sheet-rich proteins, and raises intriguing questions regarding the prevention of aggregation during the prolonged population of partially folded, beta-sheet-rich refolding intermediates. Proteins 2010. (c) 2009 Wiley-Liss, Inc.

Analysis of proton exchange kinetics with time-dependent exchange rate

Mass spectrometry is used to probe the kinetics of hydrogen-deuterium exchange in lysozyme in pH 5, 6 and 7.4. An analysis based on a Verhulst growth model is proposed and effectively applied to the kinetics of the hydrogen exchange. The data are described by a power-like function which is based on a time-dependence of the exchange rate. Experimental data ranging over many time scales is considered and accurate fits of a power-like function are obtained. Results of fittings show correlation between faster hydrogen-deuterium exchange and increase of pH. Furthermore a model is presented that discriminates between easily exchangeable hydrogens (located in close proximity to the protein surface) and those protected from the exchange (located in the protein interior). A possible interpretation of the model and its biological significance are discussed.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Parallel folding pathways in the SH3 domain protein

The transition-state ensemble (TSE) is the set of protein conformations with an equal probability to fold or unfold. Its characterization is crucial for an understanding of the folding process. We determined the TSE of the src-SH3 domain protein by using extensive molecular dynamics simulations of the Go model and computing the folding probability of a generated set of TSE candidate conformations. We found that the TSE possesses a well-defined hydrophobic core with variable enveloping structures resulting from the superposition of three parallel folding pathways. The most preferred pathway agrees with the experimentally determined TSE, while the two least preferred pathways differ significantly. The knowledge of the different pathways allows us to design the interactions between amino acids that guide the protein to fold through the least preferred pathway. This particular design is akin to a circular permutation of the protein. The finding motivates the hypothesis that the different experimentally observed TSEs in homologous proteins and circular permutants may represent potentially available pathways to the wild-type protein.

Folding mechanism of an ankyrin repeat protein: scaffold and active site formation of human CDK inhibitor p19(INK4d)

The p19(INK4d) protein consists of five ankyrin repeats (ANK) and controls the human cell cycle by inhibiting the cyclin D-dependent kinases (CDK) 4 and 6. We investigated the folding of p19(INK4d) by urea-induced unfolding transitions, kinetic analyses of unfolding and refolding, including double-mixing experiments and a special assay for folding intermediates. Folding is a sequential two-step reaction via a hyperfluorescent on-pathway intermediate. This intermediate is present under all conditions, during unfolding, refolding and at equilibrium. The folding mechanism was confirmed by a quantitative global fit of a consistent set of equilibrium and kinetic data revealing the thermodynamics and intrinsic folding rates of the different states. Surprisingly, the N<-->I transition is much faster compared to the I<-->U transition. The urea-dependence of the intrinsic folding rates causes population of the intermediate at equilibrium close to the transition midpoint. NMR detected hydrogen/deuterium exchange and the analysis of truncated variants showed that the C-terminal repeats ANK3-5 are already folded in the on-pathway intermediate, whereas the N-terminal repeats 1 and 2 are not folded. We suggest that during refolding, repeats ANK3-ANK5 first form the scaffold for the subsequent assembly of repeats ANK1 and ANK2. The binding function of p19(INK4d) resides in the latter repeats. We propose that the graded stability and the facile unfolding of repeats 1 and 2 is a prerequisite for the down-regulation of the inhibitory activity of p19(INK4d) during the cell-cycle.

Effects of Protein Concentration on the Extent of γ-Ray-Mediated Oxidative Labeling Studied by Electrospray Mass Spectrometry

Exposure of aqueous protein solutions to gamma-rays results in the formation of *OH radicals that readily react with solvent-exposed amino acid side chains. The incorporation of oxygen leads to peak progressions with a spacing of 16 Da in the mass distribution of the polypeptide chain. Unlike some other oxidative labeling strategies, these radiolysis experiments do not require solution additives that could interfere with the analysis or that might cause secondary oxidation processes. Using myoglobin as a model system, we demonstrate that the level of oxidative labeling depends critically on the protein concentration. If ignored, this effect may lead to ambiguities in the interpretation of experiments employing *OH labeling for monitoring solvent-accessible surface areas, protein folding, and protein-protein interactions. We present a simple analysis that allows oxidation to be interpreted as a process with exponential kinetics, characterized by an apparent rate constant of the form kapp=kRAD/([P]tot+B), where kRAD is the primary rate of hydroxyl radical production, B is a constant, and [P]tot is the total protein concentration. While oxidative labeling may trigger some changes in protein conformation, it is found that the magnitude of this effect is surprisingly small, a result that is consistent with observations previously made by others.

A unified mechanism for protein folding: predetermined pathways with optional errors

There is a fundamental conflict between two different views of how proteins fold. Kinetic experiments and theoretical calculations are often interpreted in terms of different population fractions folding through different intermediates in independent unrelated pathways (IUP model). However, detailed structural information indicates that all of the protein population folds through a sequence of intermediates predetermined by the foldon substructure of the target protein and a sequential stabilization principle. These contrary views can be resolved by a predetermined pathway--optional error (PPOE) hypothesis. The hypothesis is that any pathway intermediate can incorporate a chance misfolding error that blocks folding and must be reversed for productive folding to continue. Different fractions of the protein population will then block at different steps, populate different intermediates, and fold at different rates, giving the appearance of multiple unrelated pathways. A test of the hypothesis matches the two models against extensive kinetic folding results for hen lysozyme which have been widely cited in support of independent parallel pathways. The PPOE model succeeds with fewer fitting constants. The fitted PPOE reaction scheme leads to known folding behavior, whereas the IUP properties are contradicted by experiment. The appearance of a conflict with multipath theoretical models seems to be due to their different focus, namely on multitrack microscopic behavior versus cooperative macroscopic behavior. The integration of three well-documented principles in the PPOE model (cooperative foldons, sequential stabilization, optional errors) provides a unifying explanation for how proteins fold and why they fold in that way.

Folding kinetics of the S100A11 protein dimer studied by time-resolved electrospray mass spectrometry and pulsed hydrogen-deuterium exchange

This study reports the application of electrospray ionization (ESI) mass spectrometry (MS) with on-line rapid mixing for millisecond time-resolved studies of the refolding and assembly of a dimeric protein complex. Acid denaturation of S100A11 disrupts the native homodimeric protein structure. Circular dichroism and HSQC nuclear magnetic resonance measurements reveal that the monomeric subunits unfold to a moderate degree but retain a significant helicity and some tertiary structural elements. Following a rapid change in solution conditions to a slightly basic pH, the native protein reassembles with an effective rate constant of 6 s(-)(1). The ESI charge state distributions measured during the reaction suggest the presence of three kinetic species, namely, a relatively unfolded monomer (M(U)), a more tightly folded monomeric reaction intermediate (M(F)), and dimeric S100A11. These three forms exhibit distinct calcium binding properties, with very low metal loading levels for M(U), up to two calcium ions for M(F), and up to four for the dimer. Surprisingly, on-line pulsed hydrogen-deuterium exchange (HDX) reveals that each of the monomeric forms of the protein comprises two subspecies that can be distinguished on the basis of their isotope exchange levels. As the reaction proceeds, the more extensively labeled species are depleted. The exponential nature of the measured intensity-time profiles implies that the rate-determining step of the overall process is a unimolecular event. The kinetics are consistent with a sequential folding and assembly mechanism involving two increasingly nativelike monomeric intermediates en route to the native S100A11 dimer.

Refolding of a small all beta-sheet protein proceeds with accumulation of kinetic intermediates

The refolding kinetics of Cobrotoxin (CBTX), a small all beta-sheet protein is investigated using a variety of biophysical techniques including quenched-flow hydrogen-deuterium (H/D) exchange in conjunction with two-dimensional NMR spectroscopy. Urea-induced equilibrium unfolding of CBTX follows a two-state mechanism with no distinct intermediates. The protein is observed to fold very rapidly within 250 ms. Both the refolding and the unfolding limbs of the chevron plot of CBTX show a prominent curvature suggesting the accumulation of kinetic intermediates. Quenched-flow H/D exchange data suggest the presence of a broad continuum of kinetic intermediates between the unfolded and native states of the protein. Comparison of the native state hydrogen exchange data and the results of the quenched-flow H/D exchange experiments, reveals that the residues constituting the folding core of CBTX are not a subset of the slow exchange core. To our knowledge, this is the first report wherein the refolding of a small all beta-sheet protein is shown to be a multi-step process involving the accumulation of kinetic intermediates.

Pulsed Hydrogen Exchange and Electrospray Charge-State Distribution as Complementary Probes of Protein Structure in Kinetic Experiments: Implications for Ubiquitin Folding†

The possible involvement of "hidden" kinetic intermediates in the apparent two-state folding of some proteins is currently a matter of debate. This study uses time-resolved electrospray ionization (ESI) mass spectrometry with on-line pulsed hydrogen-deuterium exchange (HDX) for monitoring the refolding of acid/methanol-denatured ubiquitin. It is demonstrated that the ESI charge-state distribution (CSD) and the extent of HDX represent nonredundant probes of the protein structure in solution. When considered in isolation, the data provided by both of these probes are consistent with a two-state behavior, involving only denatured ubiquitin D and refolded protein F. However, a careful comparison of the CSD and HDX kinetics reveals the presence of an additional species, exhibiting a CSD like the folded protein but showing non-native HDX characteristics. This kinetic intermediate, D*, is in rapid equilibrium with D, such that the overall reaction is consistent with the mechanism D <--> D* --> F. The results of this work suggest that the occurrence of transient intermediates may be more widespread than commonly thought, especially in cases where a cursory analysis indicates two-state behavior.

Two-state folding of lysozyme versus multiple-state folding of alpha-lactalbumin illustrated by the technique of disulfide scrambling

The folding of lysozyme and of alpha-lactalbumin exhibits vastly different kinetics and pathways. Existing evidence indicates that folding intermediates of alphaLA form a well-populated equilibrium molten globule state that is absent in the case of hen lysozyme. We demonstrate here such divergent folding mechanisms of lysozyme and alphaLA using the technique of disulfide scrambling. Two extensively unfolded homologous isomers (beads-form) of lysozyme (Cys6-Cys30, Cys64-Cys76, Cys80-Cys94, Cys115-Cys127) and alphaLA (Cys6-Cys28, Cys61-Cys73, Cys77-Cys91, Cys111-Cys120) were allowed to refold in parallel to form the native protein. Folding kinetics was measured by the recovery of the native structure. Folding intermediates, which illustrate the folding pathway, were trapped by quenching disulfide shuffling and were analyzed by reversed-phase high-pressure liquid chromatography. The results revealed that under identical folding conditions, the folding rate of lysozyme is about 30-fold faster than that of alphaLA. Folding intermediates of lysozyme are far less heterogeneous and sparsely populated than those of alphaLA. Numerous predominant on-pathway and off-pathway intermediates observed along the folding pathway of alphaLA are conspicuously absent in the case of lysozyme. The difference is most striking under fast folding conditions performed in the presence of protein disulfide isomerase. Under these conditions, folding of lysozyme undergoes a near two-state mechanism without accumulation of stable folding intermediates.

A flexible approach for understanding protein stability

A distance constraint model (DCM) is presented that identifies flexible regions within protein structure consistent with specified thermodynamic condition. The DCM is based on a rigorous free energy decomposition scheme representing structure as fluctuating constraint topologies. Entropy non-additivity is problematic for naive decompositions, limiting the success of heat capacity predictions. The DCM resolves non-additivity by summing over independent entropic components determined by an efficient network-rigidity algorithm. A minimal 3-parameter DCM is demonstrated to accurately reproduce experimental heat capacity curves. Free energy landscapes and quantitative stability-flexibility relationships are obtained in terms of global flexibility. Several connections to experiment are made.

Lysozyme FoldedIn SilicoAccording to the Limited Conformational Sub-space

The conformational sub-space oriented on early-stage protein folding is applied to lysozyme folding. The part of the Ramachandran map distinguished on the basis of a geometrical model of the polypeptide chain limited to the mutual orientation of the peptide bond planes is shown to deliver the initial structure of the polypeptide for the energy minimization procedure in the ab initio model of protein folding prediction. Two forms of energy minimization and molecular dynamics simulation procedures were applied to the assumed early-stage protein folding of lysozyme. One of them included the disulphide bond system and the other excluded it. The post-energy-minimization and post-dynamics structures were compared using RMS-D and non-bonding contact maps to estimate the degree of approach to the native, target structure of the protein molecule obtained using the limited conformational sub-space for the early stage of folding.

Multiple parallel-pathway folding of proline-free Staphylococcal nuclease

When a protein exhibits complex kinetics of refolding, we often ascribe the complexity to slow isomerization events in the denatured protein, such as cis/trans isomerization of peptidyl prolyl bonds. Does the complex folding kinetics arise only from this well-known reason? Here, we have investigated the refolding of a proline-free variant of staphylococcal nuclease by stopped-flow, double-jump techniques, to examine the folding reactions without the slow prolyl isomerizations. As a result, the protein folds into the native state along at least two accessible parallel pathways, starting from a macroscopically single denatured-state ensemble. The presence of intermediates on the individual folding pathways has revealed the existence of multiple parallel pathways, and is characterized by multi-exponential folding kinetics with a lag phase. Therefore, a "single" amino acid sequence can fold along the multiple parallel pathways. This observation in staphylococcal nuclease suggests that the multiple folding may be more general than we have expected, because the multiple parallel-pathway folding cannot be excluded from proteins that show simpler kinetics.

Conformational dynamics of partially denatured myoglobin studied by time-resolved electrospray mass spectrometry with online hydrogen-deuterium exchange

This study demonstrates the use of electrospray mass spectrometry in conjunction with rapid online mixing ("time-resolved" ESI-MS) for monitoring protein conformational dynamics under equilibrium conditions. The hydrogen/deuterium exchange (HDX) kinetics of mildly denatured myoglobin (Mb) at pD 9.3, in the presence of 27% acetonitrile, were studied with millisecond time resolution. Analytical ultracentrifugation indicates that the average protein compactness under these solvent conditions is similar to that of native holomyoglobin (hMb). The mass spectrum shows protein ions in a wide array of charge and heme binding states, indicating the presence of multiple coexisting conformations. The experimental approach used allows the HDX kinetics of all of these species to be monitored separately. A combination of EX1 and EX2 behavior was observed for hMb ions in charge states 7+ to 9+, which predominantly represent nativelike hMb in solution. The EX1 kinetics are biphasic, indicating the presence of two protein populations that undergo conformational opening events with different rate constants. The EX2 kinetics observed for nativelike hMb are biphasic as well. All other charge and heme binding states represent non-native protein conformations that are involved in rapid interconversion processes, thus leading to monoexponential EX2 kinetics with a common rate constant. Burst phase labeling for these non-native proteins occurs at 125 sites. In contrast, the nativelike protein conformation shows burst phase labeling only for 88 sites. A kinetic model is developed which is based on the assumption of three distinct (un)folding units in Mb. The model implies that the free energy landscape of the protein exhibits a major barrier. The crossing of this barrier is most likely associated with slow, cooperative opening/closing events of the heme binding pocket. Rapid conformational fluctuations on either side of the barrier give rise to the observed EX2 kinetics. Simulated HDX kinetics based on this model are in excellent agreement with the experimental data.

Protein-folding kinetics and mechanisms studied by pulse-labeling and mass spectrometry

The "protein-folding problem" refers to the question of how and why a denatured polypeptide chain can spontaneously fold into a compact and highly ordered conformation. The classical description of this process in terms of reaction pathways has been complemented by models that describe folding as a biased conformational diffusion on a multidimensional energy landscape. The identification and characterization of short-lived intermediates provide important insights into the mechanism of folding. Pulsed hydrogen/deuterium exchange (HDX) methods are among the most powerful tools for studying the properties of kinetic intermediates. Analysis of pulse-labeled proteins by mass spectrometry (MS) provides information that is complementary to that obtained in nuclear magnetic resonance (NMR) studies; NMR data represent an average of entire protein ensembles, whereas MS can detect co-existing protein species. MS-based pulse-labeling experiments can distinguish between folding scenarios that involve parallel pathways, and those where folding is channeled through obligatory intermediates. The proteolytic digestion/MS technique provides spatially resolved information on the HDX pattern of folding intermediates. This method is especially important for proteins that are too large to be studied by NMR. Although traditional pulsed HDX protocols are based on quench-flow techniques, it is also possible to use electrospray (ESI) MS to analyze the reaction mixture on-line and "quasi-instantaneously" after labeling. This approach allows short-lived protein conformations to be studied by their HDX level, their ESI charge-state distribution, and their ligand-binding state. Covalent labeling of free cysteinyl residues provides an alternative approach to pulsed HDX experiments. Another promising development is the use of synchrotron X-rays to induce oxidation at specific sites within a protein for studying their solvent accessibility during folding.

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.

Cascade blue as a donor for resonance energy transfer studies of heme-containing proteins

Cascade Blue acetyl azide is an amine reactive compound with spectral properties ideally suited for fluorescence resonance energy transfer (FRET) studies in which heme prosthetic groups serve as acceptors. To demonstrate utility of the Cascade Blue-heme spectroscopic ruler, cytochrome c was employed as a test case to calibrate distance measurements obtained from FRET analysis. Following modification, stoichiometrically labeled cytochrome c was digested with trypsin and derivatized fragments were analyzed by matrix-assisted laser desorption/ionization and electrospray ionization mass spectrometry to identify Lys25 as the predominant site of covalent modification. FRET analysis on derivatized protein demonstrated nearly complete quenching of Cascade Blue fluorescence, indicating the labeled lysine residue to reside within 30 A of the heme prosthetic group. Sodium dodecyl sulfate (SDS) denaturation resulted in an approximately 28% recovery of fluorescence, demonstrating the utility of this donor-acceptor pair for evaluating distance changes of 30-90 A. Modeling the Cascade Blue donor molecule onto Lys25 of a cytochrome c NMR structure confirmed a distance of < or =30 A from the heme acceptor, as determined by FRET analysis. Further modeling of the SDS-denatured state as an extended chain suggested a maximum separation distance of 45 A, also consistent with results derived from FRET analysis.