Increasing Stability Reduces Conformational Heterogeneity in a Protein Folding Intermediate Ensemble

Differentiating between the sequence of structural events on alternative pathways of folding of a heterodimeric protein

Distinguishing between competing pathways of folding of a protein, on the basis of how they differ in their progress of structure acquisition, remains an important challenge in protein folding studies. A previous study had shown that the heterodimeric protein, double chain monellin (dcMN) switches between alternative folding pathways upon a change in guanidine hydrochloride (GdnHCl) concentration. In the current study, the folding of dcMN has been characterized by the pulsed hydrogen exchange (HX) labeling methodology used in conjunction with mass spectrometry. Quantification of the extent to which folding intermediates accumulate and then disappear with time of folding at both low and high GdnHCl concentrations, where the folding pathways are known to be different, shows that the folding mechanism is describable by a triangular three-state mechanism. Structural characterization of the productive folding intermediates populated on the alternative pathways has enabled the pathways to be differentiated on the basis of the progress of structure acquisition that occurs on them. The intermediates on the two pathways differ in the extent to which the α-helix and the rest of the β-sheet have acquired structure that is protective against HX. The major difference is, however, that β2 has not acquired any protective structure in the intermediate formed on one pathway, but it has acquired significant protective structure in the intermediate formed on the alternative pathway. Hence, the sequence of structural events is different on the two alternative pathways.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Structural Characterization of the Cooperativity of Unfolding of a Heterodimeric Protein using Hydrogen Exchange-Mass Spectrometry

Little is known about how the sequence of structural changes in one chain of a heterodimeric protein is coupled to those in the other chain during protein folding and unfolding reactions, and whether individual secondary structural changes in the two chains occur in one or many coordinated steps. Here, the unfolding mechanism of a small heterodimeric protein, double chain monellin, has been characterized using hydrogen exchange-mass spectrometry. Transient structure opening, which enables HX, was found to be describable by a five state N ↔ I₁ ↔ I₂ ↔ I₃ ↔ U mechanism. Structural changes occur gradually in the first three steps, and cooperatively in the last step. β strands 2, 4 and 5, as well as the α-helix undergo transient unfolding during all three non-cooperative steps, while β1 and the two loops on both sides of the helix undergo transient unfolding during the first two steps. In the absence of GdnHCl, only β3 in chain A of the protein unfolds during the last cooperative step, while in the presence of 1 M GdnHCl, not only β3, but also β2 in chain B unfolds cooperatively. Hence, the extent of cooperative structural change and size of the cooperative unfolding unit increase when the protein is destabilized by denaturant. The naturally evolved two-chain variant of monellin folds and unfolds in a more cooperative manner than does a single chain variant created artificially, suggesting that increasing folding cooperativity, even at the cost of decreasing stability, may be a driving force in the evolution of proteins.

Mapping Distinct Sequences of Structure Formation Differentiating Multiple Folding Pathways of a Small Protein

To determine experimentally how the multiple folding pathways of a protein differ, in the order in which the structural parts are assembled, has been a long-standing challenge. To resolve whether structure formation during folding can progress in multiple ways, the complex folding landscape of monellin has been characterized, structurally and temporally, using the multisite time-resolved FRET methodology. After an initial heterogeneous polypeptide chain collapse, structure formation proceeds on parallel pathways. Kinetic analysis of the population evolution data across various protein segments provides a clear structural distinction between the parallel pathways. The analysis leads to a phenomenological model that describes how and when discrete segments acquire structure independently of each other in different subensembles of protein molecules. When averaged over all molecules, structure formation is seen to progress as α-helix formation, followed by core consolidation, then β-sheet formation, and last end-to-end distance compaction. Parts of the protein that are closer in the primary sequence acquire structure before parts separated by longer sequence.

Snapshotting the transient conformations and tracing the multiple pathways of single peptide folding using a solid-state nanopore

A fundamental question relating to protein folding/unfolding is the time evolution of the folding of a protein into its precisely defined native structure. The proper identification of transition conformations is essential for accurately describing the dynamic protein folding/unfolding pathways. Owing to the rapid transitions and sub-nm conformation differences involved, the acquisition of the transient conformations and dynamics of proteins is difficult due to limited instrumental resolution. Using the electrochemical confinement effect of a solid-state nanopore, we were able to snapshot the transient conformations and trace the multiple transition pathways of a single peptide inside a nanopore. By combining the results with a Markov chain model, this new single-molecule technique is applied to clarify the transition pathways of the β-hairpin peptide, which shows nonequilibrium fluctuations among several blockage current stages. This method enables the high-throughput investigation of transition pathways experimentally to access previously obscure peptide dynamics, which is significant for understanding the folding/unfolding mechanisms and misfolding of peptides or proteins.

A solid-state nanopore based method is described for resolving protein-folding-related problems via snapshotting the folding intermediates and characterizing the kinetics of a single peptide.

Enhanced dynamics of conformationally heterogeneous T7 bacteriophage lysozyme native state attenuates its stability and activity

Proteins are dynamic in nature and exist in a set of equilibrium conformations on various timescale motions. The flexibility of proteins governs various biological functions, and therefore elucidation of such functional dynamics is essential. In this context, we have studied the structure–dynamics–stability–activity relationship of bacteriophage T7 lysozyme/endolysin (T7L) native-state ensemble in the pH range of 6–8. Our studies established that T7L native state is conformationally heterogeneous, as several residues of its C-terminal half are present in two conformations (major and minor) in the slow exchange time scale of nuclear magnetic resonance (NMR). Structural and dynamic studies suggested that the residues belonging to minor conformations do exhibit native-like structural and dynamic features. Furthermore, the NMR relaxation experiments unraveled that the native state is highly dynamic and the dynamic behavior is regulated by the pH, as the pH 6 conformation exhibited enhanced dynamics compared with pH 7 and 8. The stability measurements and cell-based activity studies on T7L indicated that the native protein at pH 6 is ∼2 kcal less stable and is ∼50% less active than those of pH 7 and 8. A comprehensive analysis of the T7L active site, unfolding initiation sites and the residues with altered dynamics outlined that the attenuation of stability and activity is a resultant of its enhanced dynamic properties, which, in turn, can be attributed to the protonation/deprotonation of its partially buried His residues. Our study on T7L structure–dynamics–activity paradigm could assist in engineering novel amidase-based endolysins with enhanced activity and stability over a broad pH range.

Single-molecule studies of flavivirus envelope dynamics: Experiment and computation

Flaviviruses are simple enveloped viruses exhibiting complex structural and functional heterogeneities. Decades of research have provided crucial basic insights, antiviral medication and moderately successful gene therapy trials. The most infectious particle is, however, not always the most abundant one in a population, questioning the utility of classic ensemble-averaging virology approaches. Indeed, viral replication is often not particularly efficient, prone to errors or containing parallel routes. Here, we review different single-molecule sensitive fluorescence methods that are employed to investigate flaviviruses. In particular, we review how (i) time-resolved Förster resonance energy transfer (trFRET) was applied to probe dengue envelope conformations; (ii) FRET-fluorescence correlation spectroscopy to investigate dengue envelope intrinsic dynamics and (iii) single particle tracking to follow the path of dengue viruses in cells. We also discuss how such methods may be supported by molecular dynamics (MD) simulations over a range of spatio-temporal scales, to provide complementary data on the structure and dynamics of flaviviral systems. We describe recent improvements in multiscale MD approaches that allowed the simulation of dengue particle envelopes in near-atomic resolution. We hope this review is an incentive for setting up and applying similar single-molecule studies and combine them with MD simulations to investigate structural dynamics of entire flavivirus particles over the nanosecond-to-millisecond time-scale and follow viruses during infection in cells over milliseconds to minutes.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Microsecond Rearrangements of Hydrophobic Clusters in an Initially Collapsed Globule Prime Structure Formation during the Folding of a Small Protein

Determining how polypeptide chain collapse initiates structure formation during protein folding is a long standing goal. It has been challenging to characterize experimentally the dynamics of the polypeptide chain, which lead to the formation of a compact kinetic molten globule (MG) in about a millisecond. In this study, the sub-millisecond events that occur early during the folding of monellin from the guanidine hydrochloride-unfolded state have been characterized using multiple fluorescence and fluorescence resonance energy transfer probes. The kinetic MG is shown to form in a noncooperative manner from the unfolded (U) state as a result of at least three different processes happening during the first millisecond of folding. Initial chain compaction completes within the first 37μs, and further compaction occurs only after structure formation commences at a few milliseconds of folding. The transient nonnative and native-like hydrophobic clusters with side chains of certain residues buried form during the initial chain collapse and the nonnative clusters quickly disassemble. Subsequently, partial chain desolvation occurs, leading to the formation of a kinetic MG. The initial chain compaction and subsequent chain rearrangement appear to be barrierless processes. The two structural rearrangements within the collapsed globule appear to prime the protein for the actual folding transition.

Resolution of Two Sub-Populations of Conformers and Their Individual Dynamics by Time Resolved Ensemble Level FRET Measurements

Most active biopolymers are dynamic structures; thus, ensembles of such molecules should be characterized by distributions of intra- or intermolecular distances and their fast fluctuations. A method of choice to determine intramolecular distances is based on Förster resonance energy transfer (FRET) measurements. Major advances in such measurements were achieved by single molecule FRET measurements. Here, we show that by global analysis of the decay of the emission of both the donor and the acceptor it is also possible to resolve two sub-populations in a mixture of two ensembles of biopolymers by time resolved FRET (trFRET) measurements at the ensemble level. We show that two individual intramolecular distance distributions can be determined and characterized in terms of their individual means, full width at half maximum (FWHM), and two corresponding diffusion coefficients which reflect the rates of fast ns fluctuations within each sub-population. An important advantage of the ensemble level trFRET measurements is the ability to use low molecular weight small-sized probes and to determine nanosecond fluctuations of the distance between the probes. The limits of the possible resolution were first tested by simulation and then by preparation of mixtures of two model peptides. The first labeled polypeptide was a relatively rigid Pro₇ and the second polypeptide was a flexible molecule consisting of (Gly-Ser)₇ repeats. The end to end distance distributions and the diffusion coefficients of each peptide were determined. Global analysis of trFRET measurements of a series of mixtures of polypeptides recovered two end-to-end distance distributions and associated intramolecular diffusion coefficients, which were very close to those determined from each of the pure samples. This study is a proof of concept study demonstrating the power of ensemble level trFRET based methods in resolution of subpopulations in ensembles of flexible macromolecules.

Unfolding of a small protein proceeds via dry and wet globules and a solvated transition state

Dissecting a protein unfolding process into individual steps can provide valuable information on the forces that maintain the integrity of the folded structure. Solvation of the protein core determines stability, but it is not clear when such solvation occurs during unfolding. In this study, far-UV circular dichroism measurements suggest a simplistic two-state view of the unfolding of barstar, but the use of multiple other probes brings out the complexity of the unfolding reaction. Near-UV circular dichroism measurements show that unfolding commences with the loosening of tertiary interactions in a native-like intermediate, N(∗). Fluorescence resonance energy transfer measurements show that N(∗) then expands rapidly but partially to form an early unfolding intermediate IE. Fluorescence spectral measurements indicate that both N(∗) and IE have retained native-like solvent accessibility of the core, suggesting that they are dry molten globules. Dynamic quenching measurements at the single tryptophan buried in the core suggest that the core becomes solvated only later in a late wet molten globule, IL, which precedes the unfolded form. Fluorescence anisotropy decay measurements show that tight packing around the core tryptophan is lost when IL forms. Of importance, the slowest step is unfolding of the wet molten globule and involves a solvated transition state.

Native state conformational heterogeneity of HP35 revealed by time-resolved FRET

The villin headpiece subdomain (HP35) has become one of the most widely used model systems in protein folding studies, due to its small size and ultrafast folding kinetics. Here, we use HP35 as a test bed to show that the fluorescence decay kinetics of an unnatural amino acid, p-cyanophenylalanine (Phe(CN)), which are modulated by a nearby quencher (e.g., tryptophan or 7-azatryptophan) through the mechanism of fluorescence resonance energy transfer (FRET), can be used to detect protein conformational heterogeneity. This method is based on the notion that protein conformations having different donor-acceptor distances and interconverting slowly compared to the fluorescence lifetime of the donor (Phe(CN)) would exhibit different donor fluorescence lifetimes. Our results provide strong evidence suggesting that the native free energy basin of HP35 is populated with conformations that differ mostly in the position and mean helicity of the C-terminal helix. This finding is consistent with several previous experimental and computational studies. Moreover, this result holds strong implications for computational investigation of the folding mechanism of HP35.

A simple model predicts experimental folding rates and a hub-like topology

A simple model is presented that describes general features of protein folding, in good agreement with experimental results and detailed all-atom simulations. Starting from microscopic physics, and with no free parameters, this model predicts that protein folding occurs remarkably quickly because native-like states are kinetic hubs. A hub-like network arises naturally out of microscopic physical concerns, specifically the kinetic longevity of native contacts during a search of globular conformations. The model predicts folding times scaling as τ(f) ~ e(ξN) in the number of residues, but because the model shows ξ is small, the folding times are much faster than Levinthal's approximation. Importantly, the folding time scale is found to be small due to the topology and structure of the network. We show explicitly how our model agrees with generic experimental features of the folding process, including the scaling of τ(f) with N, two-state thermodynamics, a sharp peak in C(V), and native-state fluctuations.

Probing molecular kinetics with Markov models: metastable states, transition pathways and spectroscopic observables

Markov (state) models (MSMs) have attracted a lot of interest recently as they (1) can probe long-term molecular kinetics based on short-time simulations, (2) offer a way to analyze great amounts of simulation data with relatively little subjectivity of the analyst, (3) provide insight into microscopic quantities such as the ensemble of transition pathways, and (4) allow simulation data to be reconciled with measurement data in a rigorous and explicit way. Here we sketch our current perspective of Markov models and explain in short their theoretical basis and assumptions. We describe transition path theory which allows the entire ensemble of protein folding pathways to be investigated and that combines naturally with Markov models. Experimental observations can be naturally linked to Markov models with the dynamical fingerprint theory, by which experimentally observable timescales can be equipped with an understanding of the structural rearrangement processes that take place at these timescales. The concepts of this paper are illustrated by a simple kinetic model of protein folding.

Dynamics of protein folding: probing the kinetic network of folding-unfolding transitions with experiment and theory

The problem of spontaneous folding of amino acid chains into highly organized, biologically functional three-dimensional protein structures continues to challenge the modern science. Understanding how proteins fold requires characterization of the underlying energy landscapes as well as the dynamics of the polypeptide chains in all stages of the folding process. In recent years, important advances toward these goals have been achieved owing to the rapidly growing interdisciplinary interest and significant progress in both experimental techniques and theoretical methods. Improvements in the experimental time resolution led to determination of the timescales of the important elementary events in folding, such as formation of secondary structure and tertiary contacts. Sensitive single molecule methods made possible probing the distributions of the unfolded and folded states and following the folding reaction of individual protein molecules. Discovery of proteins that fold in microseconds opened the possibility of atomic-level theoretical simulations of folding and their direct comparisons with experimental data, as well as of direct experimental observation of the barrier-less folding transition. The ultra-fast folding also brought new questions, concerning the intrinsic limits of the folding rates and experimental signatures of barrier-less "downhill" folding. These problems will require novel approaches for even more detailed experimental investigations of the folding dynamics as well as for the analysis of the folding kinetic data. For theoretical simulations of folding, a main challenge is how to extract the relevant information from overwhelmingly detailed atomistic trajectories. New theoretical methods have been devised to allow a systematic approach towards a quantitative analysis of the kinetic network of folding-unfolding transitions between various configuration states of a protein, revealing the transition states and the associated folding pathways at multiple levels, from atomistic to coarse-grained representations. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

Direct evidence for a dry molten globule intermediate during the unfolding of a small protein

Little is known about how proteins begin to unfold. In particular, how and when water molecules penetrate into the protein interior during unfolding, thereby enabling the dissolution of specific structure, is poorly understood. The hypothesis that the native state expands initially into a dry molten globule, in which tight packing interactions are broken, but whose hydrophobic core has not expanded sufficiently to be able to absorb water molecules, has very little experimental support. Here, we report our analysis of the earliest observable events during the unfolding of single chain monellin (MNEI), a small plant protein. Far- and near-UV circular dichroism measurements of GdnHCl-induced unfolding indicate that a molten globule intermediate forms initially, before the major slow unfolding reaction commences. Steady-state fluorescence resonance energy transfer measurements show that the C-terminal end of the single helix of MNEI initially moves rapidly away from the single tryptophan residue that is close to the N-terminal end of the helix. The average end-to-end distance of the protein also expands during unfolding to the molten globule intermediate. At this time, water has yet to penetrate the protein core, according to the evidence from intrinsic tryptophan fluorescence and 8-anilino-1-naphthalenesulfonic acid fluorescence-monitored kinetic unfolding measurements. Our results therefore provide direct evidence for a dry molten globule intermediate at the initial stage of unfolding. Our results further suggest that the structural transition between the native and dry molten globule states could be an all-or-none transition, whereas further swelling of the globule appears to occur gradually.

Continuous dissolution of structure during the unfolding of a small protein

The unfolding kinetics of many small proteins appears to be first order, when measured by ensemble-averaging probes such as fluorescence and circular dichroism. For one such protein, monellin, it is shown here that hidden behind this deceptive simplicity is a complexity that becomes evident with the use of experimental probes that are able to discriminate between different conformations in an ensemble of structures. In this study, the unfolding of monellin has been probed by measurement of the changes in the distributions of 4 different intramolecular distances, using a multisite, time-resolved fluorescence resonance energy transfer methodology. During the course of unfolding, the protein molecules are seen to undergo slow and continuous, diffusive swelling. The swelling process can be modeled as the slow diffusive swelling of a Rouse-like chain with some additional noncovalent, intramolecular interactions. Here, we show that specific structure is lost during the swelling process gradually, and not in an all-or-none manner, during unfolding.

A superior drug carrier--aponeocarzinostatin in partially unfolded state fully protects the labile antitumor enediyne

Background

Neocarzinostatin is a potent antitumor drug consisting of an enediyne chromophore and a protein carrier.

Methods

We characterized an intermediate in the equilibrium unfolding pathway of aponeocarzinostatin, using a variety of biophysical techniques including 1-anilino-8-napthalene sulfonate binding studies, size-exclusion fast protein liquid chromatography, intrinsic tryptophan fluorescence, circular dichroism, and ¹H-¹⁵N heteronuclear single quantum coherence spectroscopy.

Results

The partially unfolded protein is in molten globule-like state, in which ~60% and ~20% tertiary and secondary structure is disrupted respectively. Despite lacking a fully coordinated tertiary structure for assembling a functional binding cleft, the protein in molten globule-like state is still able to fully protect the labile chromophore. Titration of chromophore leads the partially denatured apoprotein to fold into its native state.

Conclusion

These findings bring insight into conserving mechanism of neocarzinostatin under harsh environment, where even the partially denatured apoprotein exhibits protective effect, confirming the superiority of the drug carrier.

Multiple routes and structural heterogeneity in protein folding

Experimental studies show that many proteins fold along sequential pathways defined by folding intermediates. An intermediate may not always be a single population of molecules but may consist of subpopulations that differ in their average structure. These subpopulations are likely to fold via independent pathways. Parallel folding and unfolding pathways appear to arise because of structural heterogeneity. For some proteins, the folding pathways can effectively switch either because different subpopulations of an intermediate get populated under different folding conditions, or because intermediates on otherwise hidden pathways get stabilized, leading to their utilization becoming discernible, or because mutations stabilize different substructures. Therefore, the same protein may fold via different pathways in different folding conditions. Multiple folding pathways make folding robust, and evolution is likely to have selected for this robustness to ensure that a protein will fold under the varying conditions prevalent in different cellular contexts.

Barrierless evolution of structure during the submillisecond refolding reaction of a small protein

To determine whether a protein folding reaction can occur in the absence of a dominant barrier is crucial for understanding its complexity. Here direct ultrafast kinetic measurements have been used to study the initial submillisecond (sub-ms) folding reaction of the small protein barstar. The cooperativity of the initial folding reaction has been explored by using two probes: fluorescence resonance energy transfer, through which the contraction of two intramolecular distances is measured, and the binding of 8-anilino-1-naphthalene sulfonic acid, through which the formation of hydrophobic clusters is monitored. A fast chain contraction is shown to precede the formation of hydrophobic clusters, indicating that the sub-ms folding reaction is not cooperative. The observed rate constant of the sub-ms folding reaction monitored by 8-anilino-1-naphthalene sulfonic acid fluorescence has been found to be the same in stabilizing conditions (low urea concentrations), in which specific structure is formed, and in marginally stabilizing conditions (higher urea concentrations), where virtually no structure is formed in the product of the sub-ms folding reaction. The observation that the folding rate is independent of the folding conditions suggests that the initial folding reaction occurs in the absence of a dominant free energy barrier. These results provide kinetic evidence that the formation of specific structure need not be slowed down by any significant free energy barrier during the course of a very fast protein folding reaction.

The protein folding problem

The "protein folding problem" consists of three closely related puzzles: (a) What is the folding code? (b) What is the folding mechanism? (c) Can we predict the native structure of a protein from its amino acid sequence? Once regarded as a grand challenge, protein folding has seen great progress in recent years. Now, foldable proteins and nonbiological polymers are being designed routinely and moving toward successful applications. The structures of small proteins are now often well predicted by computer methods. And, there is now a testable explanation for how a protein can fold so quickly: A protein solves its large global optimization problem as a series of smaller local optimization problems, growing and assembling the native structure from peptide fragments, local structures first.

Characterization of the Folding and Unfolding Reactions of Single-Chain Monellin: Evidence for Multiple Intermediates and Competing Pathways

The mechanisms of folding and unfolding of the small plant protein monellin have been delineated in detail. For this study, a single-chain variant of the natively two-chain monellin, MNEI, was used, in which the C terminus of chain B was connected to the N terminus of chain A by a Gly-Phe linker. Equilibrium guanidine hydrochloride (GdnHCl)-induced unfolding experiments failed to detect any partially folded intermediate that is stable enough to be populated at equilibrium to a significant extent. Kinetic experiments in which the refolding of GdnHCl-unfolded protein was monitored by measurement of the change in the intrinsic tryptophan fluorescence of the protein indicated the accumulation of three transient partially structured folding intermediates. The fluorescence change occurred in three kinetic phases: very fast, fast, and slow. It appears that the fast and slow changes in fluorescence occur on competing folding pathways originating from one unfolded form and that the very fast change in fluorescence occurs on a third parallel pathway originating from a second unfolded form of the protein. Kinetic experiments in which the refolding of alkali-unfolded protein was monitored by the change in the fluorescence of the hydrophobic dye 8-anilino-1-naphthalenesulfonic acid (ANS), consequent to the dye binding to the refolding protein, as well as by the change in intrinsic tryptophan fluorescence, not only confirmed the presence of the three kinetic intermediates but also indicated the accumulation of one or more early intermediates at a few milliseconds of refolding. These experiments also exposed a very slow kinetic phase of refolding, which was silent to any change in the intrinsic tryptophan fluorescence of the protein. Hence, the spectroscopic studies indicated that refolding of single-chain monellin occurs in five distinct kinetic phases. Double-jump, interrupted-folding experiments, in which the accumulation of folding intermediates and native protein during the folding process could be determined quantitatively by an unfolding assay, indicated that the fast phase of fluorescence change corresponds to the accumulation of two intermediates of differing stabilities on competing folding pathways. They also indicated that the very slow kinetic phase of refolding, identified by ANS binding, corresponds to the formation of native protein. Kinetic experiments in which the unfolding of native protein in GdnHCl was monitored by the change in intrinsic tryptophan fluorescence indicated that this change occurs in two kinetic phases. Double-jump, interrupted-unfolding experiments, in which the accumulation of unfolding intermediates and native protein during the unfolding process could be determined quantitatively by a refolding assay, indicated that the fast unfolding phase corresponds to the formation of fully unfolded protein via one unfolding pathway and that the slow unfolding phase corresponds to a separate unfolding pathway populated by partially unfolded intermediates. It is shown that the unfolded form produced by the fast unfolding pathway is the one which gives rise to the very fast folding pathway and that the unfolded form produced by the slower unfolding pathway is the one which gives rise to the slow and fast folding pathways.

Dissecting the Non-specific and Specific Components of the Initial Folding Reaction of Barstar by Multi-site FRET Measurements

Initial polypeptide chain collapse plays a major role in the development of subsequent structure during protein folding, but it has been difficult to elucidate the coupling between its cooperativity and specificity. To better understand this important aspect of protein folding, nine different intramolecular distances in the protein have been measured by fluorescence resonance energy transfer (FRET) in the product(s) of the initial, sub-millisecond collapse reaction during the folding of barstar, under different folding conditions. All nine distances contract in these initial folding products, when the denaturant concentration is reduced. Two of these distances were also measured in peptides corresponding to sequence segments 38-55 and 51-69 of the protein. Surprisingly, both distances do not contract in the peptides which remain fully unfolded when the denaturant concentration is reduced. This suggests that the contraction of at least some segments of the polypeptide chain may be facilitated only by contraction of other segments. In the case of the initial product of folding of the protein, the dependence on denaturant concentration of the relative change in each distance suggests that there are two components to the initial folding reaction. One is a nonspecific component, which appears to be driven by the change in denaturant concentration that is used to initiate refolding. This component corresponds to the collapse of completely unfolded protein (U) to unfolded protein in refolding conditions (U(C)). The extent of nonspecific collapse can be predicted by the response of completely unfolded protein to a change in denaturant concentration. All distances undergo such solvent-induced contraction, but each distance contracts to a different extent. There is also a specific component to initial sub-millisecond folding, in which some distances (but not all) contract more than that predicted by solvent-induced contraction. The observation that only some of the distances undergo contraction over and above solvent-induced contraction, suggest that this specific component is associated with the formation of a specific intermediate (I(E)). FRET efficiency and distance change differently for the different donor-acceptor pairs, with a change in denaturant concentration, indicating that the formation or dissolution of structure in U(C) and I(E) does not happen in a synchronized manner across different regions of the protein molecule. Also, all nine FRET efficiencies and intramolecular distances in the product(s) of sub-ms folding, change continuously with a change in denaturant concentration. Hence, it appears that the transitions from U to U(C) and to I(E) are gradual transformations, and not all-or-none structural transitions. Nevertheless, the product of these gradual transitions, I(E), possesses specific structure.

Comparison of C40/82A and P27A C40/82A barstar mutants using 19F NMR

Barstar, an inhibitor of the enzyme barnase, contains two phenylalanine residues, three tryptophan residues, and two proline residues. After incorporating either 2-19F-Phe, 4-19F-Phe, or 6-19F-Trp, the structural, dynamic, and folding properties of two mutants (C40/82A, a double mutant, and P27A C40/82A, a triple mutant) were studied by 19F NMR. Experiments were performed as a function of temperature and urea with the two mutants. We show that the consequences of the P27A mutation are extensive. The effect of the mutation is transmitted to distant residues (Phe56 and Trp53) as well as to a residue deeply buried in the hydrophobic core (Phe74). By incorporating 2-19F-Phe, it is shown that Phe56 undergoes a slow ring flipping on the NMR time scale in the triple mutant that is not observed in the double mutant. On the other hand, incorporating 4-19F-Phe shows that the P27A mutation has little effect along the Cbeta-Cgamma axis of Phe56. Labeling with 4-19F-Phe shows, from line broadening, that Phe74 experiences more dynamic motion than does Phe56 in both the double and triple mutant. After incorporating 6-19F-Trp, it is found that, in the triple mutant, Trp53 shows conformational heterogeneity at low temperature while Trp44, which is close to the P27A mutation, does not. At 20 degrees C, residual native-like structure was detected around Trp53 at high concentrations of denaturant. Barstar is cold denatured in the presence of urea. For the double mutant at temperatures below 15 degrees C, and in the presence of 2.5-3.5 M urea, the resonance for Phe74 broadens, and two peaks are observed at 5 degrees C indicative of an exchange process. From line-shape analysis, assuming a two-site conformational exchange, the rate constants as a function of temperature can be extracted. An Eyring plot is linear at 0 M urea but deviates from linearity below 20 degrees C in the presence of 2.5 or 3.5 M urea. The data as a function of urea suggest sequential events in the unfolding process.

Diffusional Barrier in the Unfolding of a Small Protein

To determine how the dynamics of the polypeptide chain in a protein molecule are coupled to the bulk solvent viscosity, the unfolding by urea of the small protein barstar was studied in the presence of two viscogens, xylose and glycerol. Thermodynamic studies of unfolding show that both viscogens stabilize barstar by a preferential hydration mechanism, and that viscogen and urea act independently on protein stability. Kinetic studies of unfolding show that while the rate-limiting conformational change during unfolding is dependent on the bulk solvent viscosity, eta, its rate does not show an inverse dependence on eta, as expected by Kramers' theory. Instead, the rate is found to be inversely proportional to an effective viscosity, eta + xi, where xi is an adjustable parameter which needs to be included in the rate equation. xi is found to have a value of -0.7 cP in xylose and -0.5 cP in glycerol, in the case of unfolding, at constant urea concentration as well as under isostability conditions. Hence, the unfolding protein chain does not experience the bulk solvent viscosity, but instead an effective solvent viscosity, which is lower than the bulk solvent viscosity by either 0.7 cP or 0.5 cP. A second important result is the validation of the isostability assumption, commonly used in protein folding studies but hitherto untested, according to which if a certain concentration of urea can nullify the effect of a certain concentration of viscogen on stability, then the same concentrations of urea and viscogen will also not perturb the free energy of activation of the unfolding of the protein.

Criteria for downhill protein folding: Calorimetry, chevron plot, kinetic relaxation, and single-molecule radius of gyration in chain models with subdued degrees of cooperativity

Recent investigations of possible downhill folding of small proteins such as BBL have focused on the thermodynamics of non-two-state, "barrierless" folding/denaturation transitions. Downhill folding is noncooperative and thermodynamically "one-state," a phenomenon underpinned by a unimodal conformational distribution over chain properties such as enthalpy, hydrophobic exposure, and conformational dimension. In contrast, corresponding distributions for cooperative two-state folding are bimodal with well-separated population peaks. Using simplified atomic modeling of a three-helix bundle-in a scheme that accounts for hydrophobic interactions and hydrogen bonding-and coarse-grained C(alpha) models of four real proteins with various degrees of cooperativity, we evaluate the effectiveness of several observables at defining the underlying distribution. Bimodal distributions generally lead to sharper transitions, with a higher heat capacity peak at the transition midpoint, compared with unimodal distributions. However, the observation of a sigmoidal transition is not a reliable criterion for two-state behavior, and the heat capacity baselines, used to determine the van't Hoff and calorimetric enthalpies of the transition, can introduce ambiguity. Interestingly we find that, if the distribution of the single-molecule radius of gyration were available, it would permit discrimination between unimodal and bimodal underlying distributions. We investigate kinetic implications of thermodynamic noncooperativity using Langevin dynamics. Despite substantial chevron rollovers, the relaxation of the models considered is essentially single-exponential over an extended range of native stabilities. Consistent with experiments, significant deviations from single-exponential behavior occur only under strongly folding conditions.

Characterization of Intra-molecular Distances and Site-specific Dynamics in Chemically Unfolded Barstar: Evidence for Denaturant-dependent Non-random Structure

The structure and dynamics of the unfolded form of a protein are expected to play critical roles in determining folding pathways. In this study, the urea and guanidine hydrochloride (GdnHCl)-unfolded forms of the small protein barstar were explored by time-resolved fluorescence techniques. Barstar was labeled specifically with thionitrobenzoate (TNB), by coupling it to the thiol side-chain of a cysteine residue at one of the following positions on the sequence: 14, 25, 40, 42, 62, 82 and 89, in single cysteine-containing mutant proteins. Seven intra-molecular distances (R(DA)) under unfolding conditions were estimated from measurements of time-resolved fluorescence resonance energy transfer between the donor Trp53 and the non-fluorescent acceptor TNB coupled to one of the seven cysteine side-chains. The unfolded protein chain expands with an increase in the concentration of the denaturants. The extent of expansion was found to be non-uniform, with different intra-molecular distances expanding to different extents. In general, shorter distances were found to expand less when compared to longer spans. The extent of expansion was higher in the case of GdnHCl when compared to urea. A comparison of the measured values of R(DA) with those derived from a model based on excluded volume, revealed that while shorter spans showed good agreement, the experimental values of R(DA) of longer spans were smaller when compared to the theoretical values. Sequence-specific flexibility of the polypeptide was determined by time-resolved fluorescence anisotropy decay measurements on acrylodan or 1,5-IAEDANS labeled single cysteine-containing proteins under unfolding conditions. Rotational dynamics derived from these measurements indicated that the level of flexibility increased with increase in the concentration of denaturants and showed a graded increase towards the C-terminal end. Taken together, these results appear to indicate the presence of specific non-random coil structures and show that the deviation from random coil structure is different for the two denaturants.

Characterization of the Formation of Amyloid Protofibrils from Barstar by Mapping Residue-specific Fluorescence Dynamics

The small protein barstar aggregates at low pH to form soluble oligomers, which can be transformed into fibrillar aggregates at an elevated temperature. To characterize structurally, with residue-specific resolution, the process of amyloid formation of barstar, as well as to monitor the increase in size that accompanies the aggregation process, time-resolved fluorescence anisotropy decay measurements have been introduced as a valuable probe. Seven different single-cysteine-containing mutant forms of barstar were made, to each of which a fluorophore was attached at the thiol group. The rotational dynamics of these seven fluorophores, as well as of the sole intrinsic tryptophan residue in the protein, were determined in the amyloid protofibrils formed, as well as in the soluble oligomers from which the protofibrils arise upon heating. Mapping of the fast rotational dynamics onto the sequence of the protein yields dynamic amplitude maps that allowed identification of the segments of the chain that possess local structure in the soluble oligomer and amyloid protofibrils. The patterns of these maps of the soluble oligomer and protofibrils are seen to be similar; and protofibrils display more local structure than do the soluble oligomers, at all residue positions studied. The observation that transformation from soluble oligomers to protofibrils does not perturb local structure significantly at eight different residue positions, suggests that the soluble oligomers transform directly into protofibrils, without undergoing drastic structural rearrangements.

The Study of Protein Folding and Dynamics by Determination of Intramolecular Distance Distributions and Their Fluctuations Using Ensemble and Single-Molecule FRET Measurements

The folding and dynamics of globular proteins is a multidimensional problem. The structures of the heterogeneous population of refolding protein molecules are characterized by multiple distances and time constants. Deciphering the mechanism of folding depends on studies of the processes rather than the folded structures alone. Spectroscopy is indispensable for these sorts of studies. Herein, it is shown that the determination of intramolecular distance distributions by ensemble and single-molecule FRET experiments enable the exploration of partially folded states of refolding protein molecules.

Native and nonnative conformational preferences in the urea-unfolded state of barstar

The refolding of barstar from its urea-unfolded state has been studied extensively using various spectroscopic probes and real-time NMR, which provide global and residue-specific information, respectively, about the folding process. Here, a preliminary structural characterization by NMR of barstar in 8 M urea has been carried out at pH 6.5 and 25 degrees C. Complete backbone resonance assignments of the urea-unfolded protein were obtained using the recently developed three-dimensional NMR techniques of HNN and HN(C)N. The conformational propensities of the polypeptide backbone in the presence of 8 M urea have been estimated by examining deviations of secondary chemical shifts from random coil values. For some residues that belong to helices in native barstar, 13C(alpha) and 13CO secondary shifts show positive deviations in the urea-unfolded state, indicating that these residues have propensities toward helical conformations. These residues are, however, juxtaposed by residues that display negative deviations indicative of propensities toward extended conformations. Thus, segments that are helical in native barstar are unlikely to preferentially populate the helical conformation in the unfolded state. Similarly, residues belonging to beta-strands 1 and 2 of native barstar do not appear to show any conformational preferences in the unfolded state. On the other hand, residues belonging to the beta-strand 3 segment show weak nonnative helical conformational preferences in the unfolded state, indicating that this segment may possess a weak preference for populating a helical conformation in the unfolded state.

Osmolytes Induce Structure in an Early Intermediate on the Folding Pathway of Barstar

Osmolytes stabilize proteins against denaturation, but little is known about how their stabilizing effect might affect a protein folding pathway. Here, we report the effects of the osmolytes, trimethylamine-N-oxide, and sarcosine on the stability of the native state of barstar as well as on the structural heterogeneity of an early intermediate ensemble, IE, on its folding pathway. Both osmolytes increase the stability of the native protein to a similar extent, with stability increasing linearly with osmolyte concentration. Both osmolytes also increase the stability of IE but to different extents. Such stabilization leads to an acceleration in the folding rate. Both osmolytes also alter the structure of IE but do so differentially; the fluorescence and circular dichroism properties of IE differ in the presence of the different osmolytes. Because these properties also differ from those of the unfolded form in refolding conditions, different burst phase changes in the optical signals are seen for folding in the presence of the different osmolytes. An analysis of the urea dependence of the burst phase changes in fluorescence and circular dichroism demonstrates that the formation of IE is itself a multistep process during folding and that the two osmolytes act by stabilizing, differentially, different structural components present in the IE ensemble. Thus, osmolytes can alter the basic nature of a protein folding pathway by discriminating, through differential stabilization, between different members of an early intermediate ensemble, and in doing so, they thereby appear to channel folding along one route when many routes are available.

Effect of Salt on the Urea-Unfolded Form of Barstar Probed bymValue Measurements†

To probe for residual structure present in the urea-unfolded form of the small protein barstar, to determine how salt might modulate such structure, and to determine how such structure might affect the stability of the protein, mutant variants that display m values different from that of the wild-type protein have been studied. The mutant proteins were obtained by site-directed mutagenesis at residue positions located on the surface of the folded protein. The m value, which represents the preferential free energy of interaction of urea with the unfolded form in comparison to that with the folded state, was determined from equilibrium urea-induced unfolding curves. Mutant proteins for which the m values were significantly greater than (m(+) mutant forms), significantly smaller than (m(-) mutant forms), or similar to (m(0) mutant forms) the m value determined for the wild-type protein were studied. The unfolded forms of the m(0), m(+) and m(-) mutant proteins represent different components within the unfolded form ensemble, which differ from each other in their solvent-exposed surface areas. Hence, the m value has been used as a measure of residual structure in the unfolded form. To further understand the nature of structures present in the unfolded form ensemble, the effects of the salt KCl on the stabilities of the wild-type and the mutant proteins, as well as on the structures present in the unfolded form ensemble, were also studied. It was found that the m values of the m(0), m(+) and m(-) mutant proteins all converge to the wild-type m value in the presence of KCl. This result indicates that the salt modulates residual structure in the unfolded form by screening electrostatic interactions that maintain compact and expanded components in the unfolded protein ensemble. The use of free energy cycles has allowed the effect of salt on the structure and free energy of the unfolded protein to be related to the stability of the protein.