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

Understanding the heterogeneity intrinsic to protein folding

Relating the native fold of a protein to its amino acid sequence remains a fundamental problem in biology. While computer algorithms have demonstrated recently their prowess in predicting what structure a particular amino acid sequence will fold to, an understanding of how and why a specific protein fold is achieved remains elusive. A major challenge is to define the role of conformational heterogeneity during protein folding. Recent experimental studies, utilizing time-resolved FRET, hydrogen-exchange coupled to mass spectrometry, and single-molecule force spectroscopy, often in conjunction with simulation, have begun to reveal how conformational heterogeneity evolves during folding, and whether an intermediate ensemble of defined free energy consists of different sub-populations of molecules that may differ significantly in conformation, energy and entropy.

Structural insights and aggregation propensity of a super-stable monellin mutant: A new potential building block for protein-based nanostructured materials

Protein fibrillation is commonly associated with pathologic amyloidosis. However, under appropriate conditions several proteins form fibrillar structures in vitro that can be used for biotechnological applications. MNEI and its variants, firstly designed as single chain derivatives of the sweet protein monellin, are also useful models for protein fibrillary aggregation studies. In this work, we have drawn attention to a protein dubbed Mut9, already characterized as a "super stable" MNEI variant. Comparative analysis of the respective X-ray structures revealed how the substitutions present in Mut9 eliminate several unfavorable interactions and stabilize the global structure. Molecular dynamic predictions confirmed the presence of a hydrogen-bonds network in Mut9 which increases its stability, especially at neutral pH. Thioflavin-T (ThT) binding assays and Fourier transform infrared (FTIR) spectroscopy indicated that the aggregation process occurs both at acidic and neutral pH, with and without addition of NaCl, even if with a different kinetics. Accordingly, Transmission Electron Microscopy (TEM) showed a fibrillar organization of the aggregates in all the tested conditions, albeit with some differences in the quantity and in the morphology of the fibrils. Our data underline the great potential of Mut9, which combines great stability in solution with the versatile conversion into nanostructured biomaterials.

Replacement of the native cis prolines by alanine leads to simplification of the complex folding mechanism of a small globular protein

The folding mechanism of MNEI, a single-chain variant of naturally occurring double-chain monellin, is complex, with multiple parallel refolding channels. To determine whether its folding energy landscape could be simplified, the two native cis-prolines, Pro41 and Pro93, were mutated, singly and together, to Ala. The stability of P93A was the same as that of the wild-type protein, pWT; however, P41A and P41AP93A were destabilized by ∼0.9 kcal mol-1. The effects of the mutations on the very fast, fast, slow, and very slow phases of folding were studied. They showed that heterogeneity in the unfolded state arises due to cis to trans isomerization of the Gly92-Pro93 peptide bond. The Pro41 to Ala mutation abolished the very slow phase of folding, whereas surprisingly, the Pro93 to Ala mutation abolished the very fast phase of folding. Double-jump, interrupted folding experiments indicated that two sequential trans to cis proline isomerization steps, of the Gly92-Pro93 peptide bond followed by the Arg40-Pro41 peptide bond, lead to the formation of the native state. They also revealed the accumulation of a late native-like intermediate, N∗, which differs from the native state in the isomeric status of the Arg40-Pro41 bond, as well as in a few tertiary contacts as monitored by near-UV CD measurements. The Pro to Ala mutations not only eliminated the cis to trans Pro isomerization reaction in the unfolded state, but also the two trans to cis Pro isomerization reactions during folding. By doing so, and by differentially affecting the relative stabilities of folding intermediates, the mutations resulted in a simplification of the folding mechanism. The two Pro to Ala mutations together accelerate folding to such an extent that the native state forms more than 1000-fold faster than in the case of pWT.

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.

Understanding the self-assembly pathways of a single chain variant of monellin: A first step towards the design of sweet nanomaterials

Peptides and proteins possess an inherent tendency to self-assemble, prompting the formation of amyloid aggregates from their soluble and functional states. Amyloids are linked to many devastating diseases, but self-assembling proteins can also represent formidable tools to produce new and sustainable biomaterials for biomedical and biotechnological applications. The mechanism of fibrillar aggregation, which influences the morphology and the properties of the protein aggregates, depend on factors such as pH, ionic strength, temperature, agitation, and protein concentration. We have here used intensive mechanical agitation, with or without beads, to prompt the aggregation of the single-chain derivative of the plant protein monellin, named MNEI, which is a well characterized sweet protein. Transmission electron microscopy confirmed the formation of fibrils several micrometers long, morphologically different from the previously characterized fibers of MNEI. Changes in the protein secondary structures during the aggregation process were monitored by Fourier transform infrared spectroscopy, which detected differences in the conformation of the final aggregates obtained under mechanical agitation. Moreover, soluble oligomers could be detected in the early phases of aggregation by polyacrylamide gel-electrophoresis. These findings emphasize the existence of multiple pathways of fibrillar aggregation for MNEI, which could be exploited for the design of innovative protein-based biomaterials.

PFDB: A standardized protein folding database with temperature correction

We constructed a standardized protein folding kinetics database (PFDB) in which the logarithmic rate constants of all listed proteins are calculated at the standard temperature (25 °C). A temperature correction based on the Eyring–Kramers equation was introduced for proteins whose folding kinetics were originally measured at temperatures other than 25 °C. We verified the temperature correction by comparing the logarithmic rate constants predicted and experimentally observed at 25 °C for 14 different proteins, and the results demonstrated improvement of the quality of the database. PFDB consists of 141 (89 two-state and 52 non-two-state) single-domain globular proteins, which has the largest number among the currently available databases of protein folding kinetics. PFDB is thus intended to be used as a standard for developing and testing future predictive and theoretical studies of protein folding. PFDB can be accessed from the following link: http://lee.kias.re.kr/~bala/PFDB.

Functional principles and regulation of molecular chaperones

To be able to perform their biological function, a protein needs to be correctly folded into its three dimensional structure. The protein folding process is spontaneous and does not require the input of energy. However, in the crowded cellular environment where there is high risk of inter-molecular interactions that may lead to protein molecules sticking to each other, hence forming aggregates, protein folding is assisted. Cells have evolved robust machinery called molecular chaperones to deal with the protein folding problem and to maintain proteins in their functional state. Molecular chaperones promote efficient folding of newly synthesized proteins, prevent their aggregation and ensure protein homeostasis in cells. There are different classes of molecular chaperones functioning in a complex interplay. In this review, we discuss the principal characteristics of different classes of molecular chaperones, their structure-function relationships, their mode of regulation and their involvement in human disorders.

Intrinsic Disorder in a Well-Folded Globular Protein

The folded structure of the heterodimeric sweet protein monellin mimics single-chain proteins with topology β1-α1-β2-β3-β4-β5 (chain A: β3-β4-β5; chain B: β1-α1-β2). Furthermore, like naturally occurring single-chain proteins of a similar size, monellin folds cooperatively with no detectable intermediates. However, the two monellin chains, A and B, are marginally structured in isolation and fold only upon binding to each other. Thus, monellin presents a unique opportunity to understand the design of intrinsically disordered proteins that fold upon binding. Here, we study the folding of a single-chain variant of monellin (scMn) using simulations of an all heavy-atom structure-based model. These simulations can explain mechanistic details derived from scMn experiments performed using several different structural probes. scMn folds cooperatively in our structure-based simulations, as is also seen in experiments. We find that structure formation near the transition-state ensemble of scMn is not uniformly distributed but is localized to a hairpin-like structure which contains one strand from each chain (β2, β3). Thus, the sequence and the underlying energetics of heterodimeric monellin promote the early formation of the interchain interface (β2-β3). By studying computational scMn mutants whose "interchain" interactions are deleted, we infer that this energy distribution allows the two protein chains to remain largely disordered when this interface is not folded. From these results, we suggest that cutting the protein backbone of a globular protein between residues which lie within its folding nucleus may be one way to construct two disordered fragments which fold upon binding.

Site-specific time-resolved FRET reveals local variations in the unfolding mechanism in an apparently two-state protein unfolding transition

Using multi-site time-resolved FRET, it is shown that equilibrium unfolding of monellin is not only heterogeneous, but that the degree of non-cooperativity differs between the sole α-helix and different parts of the β-sheet.

Understanding protein domain-swapping using structure-based models of protein folding

In domain-swapping, two or more identical protein monomers exchange structural elements and fold into dimers or multimers whose units are structurally similar to the original monomer. Domain-swapping is of biotechnological interest because inhibiting domain-swapping can reduce disease-causing fibrillar protein aggregation. To achieve such inhibition, it is important to understand both the energetics that stabilize the domain-swapped structure and the protein dynamics that enable the swapping. Structure-based models (SBMs) encode the folded structure of the protein in their potential energy functions. SBMs have been successfully used to understand diverse aspects of monomer folding. Symmetrized SBMs model interactions between two identical protein chains using only intra-monomer interactions. Molecular dynamics simulations of such symmetrized SBMs have been used to correctly predict the domain-swapped structure and to understand the mechanism of domain-swapping. Here, we review such models and illustrate that monomer topology determines key aspects of domain-swapping. However, in some proteins, specifics of local energetic interactions modulate domain-swapping and these need to be added to the symmetrized SBMs. We then summarize some general principles of the mechanism of domain-swapping that emerge from the symmetrized SBM simulations. Finally, using our own results, we explore how symmetrized SBMs could be used to design domain-swapping in proteins.

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.

Protein Domain-Swapping Can Be a Consequence of Functional Residues

Monomer topology has been implicated in domain-swapping, a potential first step on the route to disease-causing protein aggregation. Despite having the same topology (β1-α1-β2-β3-β4-β5), the cysteine protease inhibitor stefin-B domain swaps more readily than a single-chain variant of the heterodimeric sweet protein monellin (scMn). Here, we computationally study the folding of stefin-B and scMn in order to understand the molecular basis for the difference in their domain-swapping propensities. In agreement with experiments, our structure-based simulations show that scMn folds cooperatively without the population of an intermediate while stefin-B populates an equilibrium intermediate state. Since the simulation intermediate has only one domain structured (β3-β4-β5), it can directly lead to domain-swapping. Using computational variants of stefin-B, we show that the population of this intermediate is caused by regions of stefin-B that have been implicated in protease inhibition. We also find that the protease-binding regions are located on two structural elements and localized in space. In contrast, the residues that contribute to the sweetness of monellin are not localized to a few structural elements but are distributed over the protein fold. We conclude that the distributed functional residues of monellin do not induce large local perturbations in the protein structure, eliminating the formation of folding intermediates and in turn domain-swapping. On the other hand, the localized protease-binding regions of stefin-B promote the formation of a folding intermediate which can lead to domain-swapping. Thus, domain-swapping can be a direct consequence of the constraints that function imposes on the protein structure.

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.

Direct observation of parallel folding pathways revealed using a symmetric repeat protein system

Although progress has been made to determine the native fold of a polypeptide from its primary structure, the diversity of pathways that connect the unfolded and folded states has not been adequately explored. Theoretical and computational studies predict that proteins fold through parallel pathways on funneled energy landscapes, although experimental detection of pathway diversity has been challenging. Here, we exploit the high translational symmetry and the direct length variation afforded by linear repeat proteins to directly detect folding through parallel pathways. By comparing folding rates of consensus ankyrin repeat proteins (CARPs), we find a clear increase in folding rates with increasing size and repeat number, although the size of the transition states (estimated from denaturant sensitivity) remains unchanged. The increase in folding rate with chain length, as opposed to a decrease expected from typical models for globular proteins, is a clear demonstration of parallel pathways. This conclusion is not dependent on extensive curve-fitting or structural perturbation of protein structure. By globally fitting a simple parallel-Ising pathway model, we have directly measured nucleation and propagation rates in protein folding, and have quantified the fluxes along each path, providing a detailed energy landscape for folding. This finding of parallel pathways differs from results from kinetic studies of repeat-proteins composed of sequence-variable repeats, where modest repeat-to-repeat energy variation coalesces folding into a single, dominant channel. Thus, for globular proteins, which have much higher variation in local structure and topology, parallel pathways are expected to be the exception rather than the rule.

Equilibrium unfolding studies of monellin: the double-chain variant appears to be more stable than the single-chain variant

To improve our understanding of the contributions of different stabilizing interactions to protein stability, including that of residual structure in the unfolded state, the small sweet protein monellin has been studied in both its two variant forms, the naturally occurring double-chain variant (dcMN) and the artificially created single-chain variant (scMN). Equilibrium guanidine hydrochloride-induced unfolding studies at pH 7 show that the standard free energy of unfolding, ΔG°(U), of dcMN to unfolded chains A and B and its dependence on guanidine hydrochloride (GdnHCl) concentration are both independent of protein concentration, while the midpoint of unfolding has an exponential dependence on protein concentration. Hence, the unfolding of dcMN like that of scMN can be described as two-state unfolding. The free energy of dissociation, ΔG°(d), of the two free chains, A and B, from dcMN, as measured by equilibrium binding studies, is significantly lower than ΔG°(U), apparently because of the presence of residual structure in free chain B. The value of ΔG°(U), at the standard concentration of 1 M, is found to be ∼5.5 kcal mol(-1) higher for dcMN than for scMN in the range from pH 4 to 9, over which unfolding appears to be two-state. Hence, dcMN appears to be more stable than scMN. It seems that unfolded scMN is stabilized by residual structure that is absent in unfolded dcMN and/or that native scMN is destabilized by strain that is relieved in native dcMN. The value of ΔG°(U) for both protein variants decreases with an increase in pH from 4 to 9, apparently because of the thermodynamic coupling of unfolding to the protonation of a buried carboxylate side chain whose pK(a) shifts from 4.5 in the unfolded state to 9 in the native state. Finally, it is shown that although the thermodynamic stabilities of dcMN and scMN are very different, their kinetic stabilities with respect to unfolding in GdnHCl are very similar.

Evidence for initial non-specific polypeptide chain collapse during the refolding of the SH3 domain of PI3 kinase

Refolding of the SH3 domain of PI3 kinase from the guanidine hydrochloride (GdnHCl)-unfolded state has been probed with millisecond (stopped flow) and sub-millisecond (continuous flow) measurements of the change in fluorescence, circular dichroism, ANS fluorescence and three-site fluorescence resonance energy transfer (FRET) efficiency. Fluorescence measurements are unable to detect structural changes preceding the rate-limiting step of folding, whereas measurements of changes in ANS fluorescence and FRET efficiency indicate that polypeptide chain collapse precedes the major structural transition. The initial chain collapse reaction is complete within 150 μs. The collapsed form at this time possesses hydrophobic clusters to which ANS binds. Each of the three measured intra-molecular distances has contracted to an extent predicted by the dependence of the FRET signal in completely unfolded protein on denaturant concentration, indicating that contraction is non-specific. The extent of contraction of each intra-molecular distance in the collapsed product of sub-millisecond folding increases continuously with a decrease in [GdnHCl]. The gradual contraction is continuous with the gradual contraction seen in completely unfolded protein, and its dependence on [GdnHCl] is not indicative of an all-or-none collapse reaction. The dependence of the extent of contraction on [GdnHCl] was similar for the three distances, indicating that chain collapse occurs in a synchronous manner across different segments of the polypeptide chain. The sub-millisecond measurements of folding in GdnHCl were unable to determine whether hydrophobic cluster formation, probed by ANS fluorescence measurement, precedes chain contraction probed by FRET. To determine whether hydrogen bonding plays a role in initial chain collapse, folding was initiated by dilution of the urea-unfolded state. The extent of contraction of at least one intra-molecular distance in the collapsed product of sub-millisecond folding in urea is similar to that seen in GdnHCl, and the initial contraction in urea too appears to be gradual.

Structure-based engineering of a monoclonal antibody for improved solubility

Protein aggregation is of great concern to pharmaceutical formulations and has been implicated in several diseases. We engineered an anti-IL-13 monoclonal antibody CNTO607 for improved solubility. Three structure-based engineering approaches were employed in this study: (i) modifying the isoelectric point (pI), (ii) decreasing the overall surface hydrophobicity and (iii) re-introducing an N-linked carbohydrate moiety within a complementarity-determining region (CDR) sequence. A mutant was identified with a modified pI that had a 2-fold improvement in solubility while retaining the binding affinity to IL-13. Several mutants with decreased overall surface hydrophobicity also showed moderately improved solubility while maintaining a similar antigen affinity. Structural studies combined with mutagenesis data identified an aggregation 'hot spot' in heavy-chain CDR3 (H-CDR3) that contains three residues ((99)FHW(100a)). The same residues, however, were found to be essential for high affinity binding to IL-13. On the basis of the spatial proximity and germline sequence, we reintroduced the consensus N-glycosylation site in H-CDR2 which was found in the original antibody, anticipating that the carbohydrate moiety would shield the aggregation 'hot spot' in H-CDR3 while not interfering with antigen binding. Peptide mapping and mass spectrometric analysis revealed that the N-glycosylation site was generally occupied. This variant showed greatly improved solubility and bound to IL-13 with affinity similar to CNTO607 without the N-linked carbohydrate. All three engineering approaches led to improved solubility and adding an N-linked carbohydrate to the CDR was the most effective route for enhancing the solubility of CNTO607.

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.

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.