Apparent two-state tendamistat folding is a sequential process along a defined route

The assembly platform FimD is required to obtain the most stable quaternary structure of type 1 pili

Type 1 pili are important virulence factors of uropathogenic Escherichia coli that mediate bacterial attachment to epithelial cells in the urinary tract. The pilus rod is comprised of thousands of copies of the main structural subunit FimA and is assembled in vivo by the assembly platform FimD. Although type 1 pilus rods can self-assemble from FimA in vitro, this reaction is slower and produces structures with lower kinetic stability against denaturants compared to in vivo-assembled rods. Our study reveals that FimD-catalysed in vitro-assembled type 1 pilus rods attain a similar stability as pilus rods assembled in vivo. Employing structural, biophysical and biochemical analyses, we show that in vitro assembly reactions lacking FimD produce pilus rods with structural defects, reducing their stability against dissociation. Overall, our results indicate that FimD is not only required for the catalysis of pilus assembly, but also to control the assembly of the most stable quaternary structure.

Determining folding and binding properties of the C-terminal SH2 domain of SHP2

SH2 domains are a class of protein–protein interaction modules with the function to recognize and bind sequences characterized by the presence of a phosphorylated tyrosine. SHP2 is a protein phosphatase involved in the Ras‐ERK1/2 signaling pathway that possess two SH2 domains, namely, N‐SH2 and C‐SH2, that mediate the interaction of SHP2 with various partners and determine the regulation of its catalytic activity. One of the main interactors of the SH2 domains of SHP2 is Gab2, a scaffolding protein with critical role in determining cell differentiation. Despite their key biological role and the importance of a correct native fold to ensure it, the mechanism of binding of SH2 domains with their ligands and the determinants of their stability have been poorly characterized. In this article, we present a comprehensive kinetic study of the folding of the C‐SH2 domain and the binding mechanism with a peptide mimicking a region of Gab2. Our data, obtained at different pH and ionic strength conditions and supported by site‐directed mutagenesis, highlight the role of electrostatic interactions in the early events of recognition. Interestingly, our results suggest a key role of a highly conserved histidine residue among SH2 family in the interaction with negative charges carried by the phosphotyrosine of Gab2. Moreover, the analysis of the equilibrium and kinetic folding data of C‐SH2 describes a complex mechanism implying a change in rate‐limiting step at high denaturant concentrations. Our data are discussed under the light of previous works on N‐SH2 domain of SHP2 and other SH2 domains.

Folding and Stability of Ankyrin Repeats Control Biological Protein Function

Ankyrin repeat proteins are found in all three kingdoms of life. Fundamentally, these proteins are involved in protein-protein interaction in order to activate or suppress biological processes. The basic architecture of these proteins comprises repeating modules forming elongated structures. Due to the lack of long-range interactions, a graded stability among the repeats is the generic properties of this protein family determining both protein folding and biological function. Protein folding intermediates were frequently found to be key for the biological functions of repeat proteins. In this review, we discuss most recent findings addressing this close relation for ankyrin repeat proteins including DARPins, Notch receptor ankyrin repeat domain, IκBα inhibitor of NFκB, and CDK inhibitor p19^INK4d. The role of local folding and unfolding and gradual stability of individual repeats will be discussed during protein folding, protein-protein interactions, and post-translational modifications. The conformational changes of these repeats function as molecular switches for biological regulation, a versatile property for modern drug discovery.

Investigating the Effect of Chain Connectivity on the Folding of a Beta-Sheet Protein On and Off the Ribosome

Determining the relationship between protein folding pathways on and off the ribosome remains an important area of investigation in biology. Studies on isolated domains have shown that alteration of the separation of residues in a polypeptide chain, while maintaining their spatial contacts, may affect protein stability and folding pathway. Due to the vectorial emergence of the polypeptide chain from the ribosome, chain connectivity may have an important influence upon cotranslational folding. Using MATH, an all β-sandwich domain, we investigate whether the connectivity of residues and secondary structure elements is a key determinant of when cotranslational folding can occur on the ribosome. From Φ-value analysis, we show that the most structured region of the transition state for folding in MATH includes the N and C terminal strands, which are located adjacent to each other in the structure. However, arrest peptide force-profile assays show that wild-type MATH is able to fold cotranslationally, while some C-terminal residues remain sequestered in the ribosome, even when destabilized by 2–3 kcal mol⁻¹. We show that, while this pattern of Φ-values is retained in two circular permutants in our studies of the isolated domains, one of these permutants can fold only when fully emerged from the ribosome. We propose that in the case of MATH, onset of cotranslational folding is determined by the ability to form a sufficiently stable folding nucleus involving both β-sheets, rather than by the location of the terminal strands in the ribosome tunnel.

•

Adjacent N and C terminal strands are most structured region in the transition state.

•

Two circular permutants retain the same folding pathway as wild-type MATH.

•

On the ribosome, early emergence of terminal strands does not promote earlier folding.

•

Formation of both β-sheets is energetically critical for folding on the ribosome.

•

Folding pathway minimizes formation of partly structured states prone to mis-folding.

The Cryoelectron Microscopy Structure of the Type 1 Chaperone-Usher Pilus Rod

Adhesive chaperone-usher pili are long, supramolecular protein fibers displayed on the surface of many bacterial pathogens. The type 1 and P pili of uropathogenic Escherichia coli (UPEC) play important roles during urinary tract colonization, mediating attachment to the bladder and kidney, respectively. The biomechanical properties of the helical pilus rods allow them to reversibly uncoil in response to flow-induced forces, allowing UPEC to retain a foothold in the unique and hostile environment of the urinary tract. Here we provide the 4.2-Å resolution cryo-EM structure of the type 1 pilus rod, which together with the previous P pilus rod structure rationalizes the remarkable "spring-like" properties of chaperone-usher pili. The cryo-EM structure of the type 1 pilus rod differs in its helical parameters from the structure determined previously by a hybrid approach. We provide evidence that these structural differences originate from different quaternary structures of pili assembled in vivo and in vitro.

Three-State Fluorescence of a 2-Functionalized Pyrene-Based RNA Label

The pyrene-based RNA-fluorescence label 2-(2-pyrenylethynyl) adenosine (2PyA) shows triexponential fluorescence, which depends strongly on the excitation wavelength. Most strikingly, a structured, long-lived fluorescence is observed in solution at room temperature after excitation into the S₂ state, which is shifted hypsochromically by 30 nm compared to excitation into the S₁ state. This very unusual behavior is investigated in detail with steady-state and time-resolved emission spectroscopy, ultrafast transient absorption spectroscopy, and quantum chemical calculations with both wave functions (CC2-level) and density-functional theory (DFT). 2PyA is found to emit simultaneously from two different intramolecular charge transfer states (mesomeric and twisted, MICT and TICT) which are populated most efficiently via the S₁ state and a pyrene-like locally excited (LE) state. Rotational momentum derived from excess excitation energy is required to populate twisted LE configurations. Therefore, the LE state is most efficiently accessible via excitation to the S₂. The stabilization of the different substates is related to two distinct reaction coordinates: the adenine-pyrene distance and the adenine-pyrene tilt angle, respectively.

Perplexing cooperative folding and stability of a low-sequence complexity, polyproline 2 protein lacking a hydrophobic core

The burial of hydrophobic side chains in a protein core generally is thought to be the major ingredient for stable, cooperative folding. Here, we show that, for the snow flea antifreeze protein (sfAFP), stability and cooperativity can occur without a hydrophobic core, and without α-helices or β-sheets. sfAFP has low sequence complexity with 46% glycine and an interior filled only with backbone H-bonds between six polyproline 2 (PP2) helices. However, the protein folds in a kinetically two-state manner and is moderately stable at room temperature. We believe that a major part of the stability arises from the unusual match between residue-level PP2 dihedral angle bias in the unfolded state and PP2 helical structure in the native state. Additional stabilizing factors that compensate for the dearth of hydrophobic burial include shorter and stronger H-bonds, and increased entropy in the folded state. These results extend our understanding of the origins of cooperativity and stability in protein folding, including the balance between solvent and polypeptide chain entropies.

Disorder drives cooperative folding in a multidomain protein

Many human proteins contain intrinsically disordered regions, and disorder in these proteins can be fundamental to their function-for example, facilitating transient but specific binding, promoting allostery, or allowing efficient posttranslational modification. SasG, a multidomain protein implicated in host colonization and biofilm formation in Staphylococcus aureus, provides another example of how disorder can play an important role. Approximately one-half of the domains in the extracellular repetitive region of SasG are intrinsically unfolded in isolation, but these E domains fold in the context of their neighboring folded G5 domains. We have previously shown that the intrinsic disorder of the E domains mediates long-range cooperativity between nonneighboring G5 domains, allowing SasG to form a long, rod-like, mechanically strong structure. Here, we show that the disorder of the E domains coupled with the remarkable stability of the interdomain interface result in cooperative folding kinetics across long distances. Formation of a small structural nucleus at one end of the molecule results in rapid structure formation over a distance of 10 nm, which is likely to be important for the maintenance of the structural integrity of SasG. Moreover, if this normal folding nucleus is disrupted by mutation, the interdomain interface is sufficiently stable to drive the folding of adjacent E and G5 domains along a parallel folding pathway, thus maintaining cooperative folding.

Human Frataxin Folds Via an Intermediate State. Role of the C-Terminal Region

The aim of this study is to investigate the folding reaction of human frataxin, whose deficiency causes the neurodegenerative disease Friedreich's Ataxia (FRDA). The characterization of different conformational states would provide knowledge about how frataxin can be stabilized without altering its functionality. Wild-type human frataxin and a set of mutants, including two highly destabilized FRDA-associated variants were studied by urea-induced folding/unfolding in a rapid mixing device and followed by circular dichroism. The analysis clearly indicates the existence of an intermediate state (I) in the folding route with significant secondary structure content but relatively low compactness, compared with the native ensemble. However, at high NaCl concentrations I-state gains substantial compaction, and the unfolding barrier is strongly affected, revealing the importance of electrostatics in the folding mechanism. The role of the C-terminal region (CTR), the key determinant of frataxin stability, was also studied. Simulations consistently with experiments revealed that this stretch is essentially unstructured, in the most compact transition state ensemble (TSE2). The complete truncation of the CTR drastically destabilizes the native state without altering TSE2. Results presented here shed light on the folding mechanism of frataxin, opening the possibility of mutating it to generate hyperstable variants without altering their folding kinetics.

Fast Folding Dynamics of an Intermediate State in RNase H Measured by Single-Molecule FRET

We have studied the folding kinetics of the core intermediate (I) state of RNase H by using a combination of single-molecule FRET (smFRET) and hidden Markov model analysis. To measure fast dynamics in thermal equilibrium as a function of the concentration of the denaturant GdmCl, a special FRET labeled variant, RNase H 60-113, which is sensitive to folding of the protein core, was immobilized on PEGylated surfaces. Conformational transitions between the unfolded (U) state and the I state could be described by a two-state model within our experimental time resolution, with millisecond mean residence times. The I state population was always a minority species in the entire accessible range of denaturant concentrations. By introducing the measured free energy differences between the U and I states as constraints in global fits of the GdmCl dependence of FRET histograms of a differently labeled RNase H variant (RNase H 3-135), we were able to reveal the free energy differences and, thus, population ratios of all three macroscopic state ensembles, U, I and F (folded state) as a function of denaturant concentration.

Acceleration of protein folding by four orders of magnitude through a single amino acid substitution

Cis prolyl peptide bonds are conserved structural elements in numerous protein families, although their formation is energetically unfavorable, intrinsically slow and often rate-limiting for folding. Here we investigate the reasons underlying the conservation of the cis proline that is diagnostic for the fold of thioredoxin-like thiol-disulfide oxidoreductases. We show that replacement of the conserved cis proline in thioredoxin by alanine can accelerate spontaneous folding to the native, thermodynamically most stable state by more than four orders of magnitude. However, the resulting trans alanine bond leads to small structural rearrangements around the active site that impair the function of thioredoxin as catalyst of electron transfer reactions by more than 100-fold. Our data provide evidence for the absence of a strong evolutionary pressure to achieve intrinsically fast folding rates, which is most likely a consequence of proline isomerases and molecular chaperones that guarantee high in vivo folding rates and yields.

Redefining the dry molten globule state of proteins

Dynamics and function of proteins are governed by the structural and energetic properties of the different states they adopt and the barriers separating them. In earlier work, native-state triplet-triplet energy transfer (TTET) on the villin headpiece subdomain (HP35) revealed an equilibrium between a locked native state and an unlocked native state, which are structurally similar but have different dynamic properties. The locked state is restricted to low amplitude motions, whereas the unlocked state shows increased conformational flexibility and undergoes local unfolding reactions. This classified the unlocked state as a dry molten globule (DMG), which was proposed to represent an expanded native state with loosened side-chain interactions and a solvent-shielded core. To test whether the unlocked state of HP35 is actually expanded compared to the locked state, we performed high-pressure TTET measurements. Increasing pressure shifts the equilibrium from the locked toward the unlocked state, with a small negative reaction volume for unlocking (ΔV(0)=-1.6±0.5cm(3)/mol). Therefore, rather than being expanded, the unlocked state represents an alternatively packed, compact state, demonstrating that native proteins can exist in several compact folded states, an observation with implications for protein function. The transition state for unlocking/locking, in contrast, has a largely increased volume relative to the locked and unlocked state, with respective activation volumes of 7.1±0.4cm(3)/mol and 8.7±0.9cm(3)/mol, indicating an expansion of the protein during the locking/unlocking transition. The presented results demonstrate the existence of both compact, low-energy and expanded, high-energy DMGs, prompting a broader definition of this state.

Analyses of the folding properties of ferredoxin-like fold proteins by means of a coarse-grained Gō model: relationship between the free energy profiles and folding cores

The folding mechanisms of proteins with multi-state transitions, the role of the intermediate states, and the precise mechanism how each transition occurs are significant on-going research issues. In this study, we investigate ferredoxin-like fold proteins which have a simple topology and multi-state transitions. We analyze the folding processes by means of a coarse-grained Gō model. We are able to reproduce the differences in the folding mechanisms between U1A, which has a high-free-energy intermediate state, and ADA2h and S6, which fold into the native structure through two-state transitions. The folding pathways of U1A, ADA2h, S6, and the S6 circular permutant, S6_p54-55, are reproduced and compared with experimental observations. We show that the ferredoxin-like fold contains two common regions consisting folding cores as predicted in other studies and that U1A produces an intermediate state due to the distinct cooperative folding of each core. However, because one of the cores of S6 loses its cooperativity and the two cores of ADA2h are tightly coupled, these proteins fold into the native structure through a two-state mechanism.

Protein folding: adding a nucleus to guide helix docking reduces landscape roughness

The elongated three-helix‐bundle spectrin domains R16 and R17 fold and unfold unusually slowly over a rough energy landscape, in contrast to the homologue R15, which folds fast over a much smoother, more typical landscape. R15 folds via a nucleation–condensation mechanism that guides the docking of the A and C-helices. However, in R16 and R17, the secondary structure forms first and the two helices must then dock in the correct register. Here, we use variants of R16 and R17 to demonstrate that substitution of just five key residues is sufficient to alter the folding mechanism and reduce the landscape roughness. We suggest that, by providing access to an alternative, faster, folding route over their landscape, R16 and R17 can circumvent their slow, frustrated wild-type folding mechanism.

Separating the effects of internal friction and transition state energy to explain the slow, frustrated folding of spectrin domains

The elongated three-helix bundle domains spectrin R16 and R17 fold some two to three orders of magnitude more slowly than their homologue R15. We have shown that this slow folding is due, at least in part, to roughness in the free-energy landscape of R16 and R17. We have proposed that this roughness is due to a frustrated search for the correct docking of partly preformed helices. However, this accounts for only a small part of the slowing of folding and unfolding. Five residues on the A helix of R15, when inserted together into R16 or R17, increase the folding rate constants, reduce landscape roughness, and alter the folding mechanism to one resembling R15. The effect of each of these mutations individually is investigated here. No one mutation causes the behavior seen for the five in combination. However, two mutations, E18F and K25V, significantly increase the folding and unfolding rates of both R16 and R17 but without a concomitant loss in landscape roughness. E18F has the greatest effect on the kinetics, and a Φ-value analysis of the C helix reveals that the folding mechanism is unchanged. For both E18F and K25V the removal of the charge and resultant transition state stabilization is the main origin of the faster folding. Consequently, the major cause of the unusually slow folding of R16 and R17 is the non-native burial of the two charged residues in the transition state. The slowing due to landscape roughness is only about fivefold.

Sequential four-state folding/unfolding of goat α-lactalbumin and its N-terminal variants

Equilibria and kinetics of folding/unfolding of α-lactalbumin and its two N-terminal variants were studied by circular dichroism spectroscopy. The two variants were wild-type recombinant and Glu1-deletion (E1M) variants expressed in Escherichia coli. The presence of an extra methionine at the N terminus in recombinant α-lactalbumin destabilized the protein by 2 kcal/mol, while the stability was recovered in the E1M variant in which Glu1 was replaced by Met1. Kinetic folding/unfolding reactions of the proteins, induced by stopped-flow concentration jumps of guanidine hydrochloride, indicated the presence of a burst-phase in refolding, and gave chevron plots with significant curvatures in both the folding and unfolding limbs. The folding-limb curvature was interpreted in terms of accumulation of the burst-phase intermediate. However, there was no burst phase observed in the unfolding kinetics to interpret the unfolding-limb curvature. We thus assumed a sequential four-state mechanism, in which the folding from the burst-phase intermediate takes place via two transition states separated by a high-energy intermediate. We estimated changes in the free energies of the burst-phase intermediate and two transition states, caused by the N-terminal variations and also by the presence of stabilizing calcium ions. The Φ values at the N terminus and at the Ca(2+)-binding site thus obtained increased successively during folding, demonstrating the validity of the sequential mechanism. The stability and the folding behavior of the E1M variant were essentially identical to those of the authentic protein, allowing us to use this variant as a pseudo-wild-type α-lactalbumin in future studies.

Folding mechanism of an extremely thermostable (βα)(8)-barrel enzyme: a high kinetic barrier protects the protein from denaturation

HisF, the cyclase subunit of imidazole glycerol phosphate synthase (ImGPS) from Thermotoga maritima, is an extremely thermostable (βα)(8)-barrel protein. We elucidated the unfolding and refolding mechanism of HisF. Its unfolding transition is reversible and adequately described by the two-state model, but 6 weeks is necessary to reach equilibrium (at 25 °C). During refolding, initially a burst-phase off-pathway intermediate is formed. The subsequent productive folding occurs in two kinetic phases with time constants of ~3 and ~20 s. They reflect a sequential process via an on-pathway intermediate, as revealed by stopped-flow double-mixing experiments. The final step leads to native HisF, which associates with the glutaminase subunit HisH to form the functional ImGPS complex. The conversion of the on-pathway intermediate to the native protein results in a 10(6)-fold increase of the time constant for unfolding from 89 ms to 35 h (at 4.0 M GdmCl) and thus establishes a high energy barrier to denaturation. We conclude that the extra stability of HisF is used for kinetic protection against unfolding. In its refolding mechanism, HisF resembles other (βα)(8)-barrel proteins.

Overview of protein folding mechanisms: experimental and theoretical approaches to probing energy landscapes

We present an overview of the current experimental and theoretical approaches to studying protein folding mechanisms, set against current models of the folding energy landscape. We describe how stability and folding kinetics can be determined experimentally and how this data can be interpreted in terms of the characteristic features of various models from the simplest two-state pathway to a multi-state mechanism. We summarize the pros and cons of a range of spectroscopic methods for measuring folding rates and present a theoretical framework, coupled with protein engineering approaches, for elucidating folding mechanisms and structural features of folding transition states. A series of case studies are used to show how experimental kinetic data can be interpreted in the context of non-native interactions, populated intermediates, parallel folding pathways, and sequential transition states. We also show how computational methods now allow transient species of high energy, such as folding transition states, to be modeled on the basis of experimental Φ-value analysis derived from the effects of point mutations on folding kinetics.

GB1 is not a two-state folder: identification and characterization of an on-pathway intermediate

The folding pathway of the small α/β protein GB1 has been extensively studied during the past two decades using both theoretical and experimental approaches. These studies provided a consensus view that the protein folds in a two-state manner. Here, we reassessed the folding of GB1, both by experiments and simulations, and detected the presence of an on-pathway intermediate. This intermediate has eluded earlier experimental characterization and is distinct from the collapsed state previously identified using ultrarapid mixing. Failure to identify the presence of an intermediate affects some of the conclusions that have been drawn for GB1, a popular model for protein folding studies.

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.

The stability of intramolecular DNA G-quadruplexes compared with other macromolecules

DNA quadruplexes are often conceived as very stable structures. However, most of the free energy of stabilization derives from specific ion binding via inner sphere coordination of the GO6 of the guanine residues comprising the basic quartet. When compared with other nucleic acid structures such as DNA or RNA duplexes and hairpins, or proteins of the same number of atoms, metal-coordinated intramolecular quadruplexes are found to be of comparable or lower thermodynamic stability under similar solution conditions. Furthermore, intramolecular quadruplexes are actually less stable kinetically, than DNA duplexes or hairpins of the same size. Although the literature is incomplete, it is clear that polyelectrolyte ion effects, the influence of solvation and steric crowding on stability are qualitatively different between intramolecular quadruplexes and DNA duplexes. For example, decreasing water activity destabilizes DNA duplexes, whereas quadruplexes are stabilized. The variety of folded conformations accessible to a single sequence further implies strong sensitivity of the conformational ensemble to the solution conditions, compared with DNA duplexes or small single domain proteins. These considerations may have relevance to the conditions prevailing inside cell nuclei and therefore the structures that potentially might form in vivo.

The effect of glycosylation on the folding kinetics of erythropoietin

Glycosylation is a common posttranslational modification that generally increases protein solubility and thermodynamic stability. Less is known about how this modification influences protein folding, particularly folding processes involving intermediate species. In the present report, folding comparisons of a nonglycosylated erythropoietin (EPO) mutant are made with the fully glycosylated EPO, which was recently shown to fold by a three-state on-pathway mechanism. The absence of glycosylation did not alter the folding mechanism of EPO but did greatly decrease the stability of the intermediate species, change the rate-limiting step of the folding reaction, and accelerate the folding kinetics to both the intermediate state and the native state. Surprisingly, glycosylation stabilized the intermediate species to a greater extent than it increased the EPO equilibrium stability. These results suggest that glycosylation impedes the latter EPO folding steps rather than accelerating them by biasing particular folding pathways, as previously proposed for folding reactions initiated from unfolded ensembles with minimal residual structure. Due to the specific biological processes modulated by EPO glycosylation, however, there may be little evolutionary pressure to fold on a faster, more direct pathway at the expense of biological function, particularly given the protective role glycosylation has at preventing EPO aggregation. Lastly, evidence that is consistent with glycosylation destabilizing the unfolded state to some degree and contributing to the greater equilibrium stability of the glycosylated EPO is presented.

Protein folding mechanisms studied by pulsed oxidative labeling and mass spectrometry

Deciphering the mechanisms of protein folding remains a considerable challenge. In this review we discuss the application of pulsed oxidative labeling for tracking protein structural changes in a time-resolved fashion. Exposure to a microsecond OH pulse at selected time points during folding induces the oxidation of solvent-accessible side chains, whereas buried residues are protected. Oxidative modifications can be detected by mass spectrometry. Folding is associated with dramatic accessibility changes, and therefore this method can provide detailed mechanistic insights. Solvent accessibility patterns are complementary to H/D exchange investigations, which report on the extent of hydrogen bonding. This review highlights the application of pulsed OH labeling to soluble proteins as well as membrane proteins.

Folding kinetics of the cooperatively folded subdomain of the IκBα ankyrin repeat domain

The ankyrin repeat (AR) domain of IκBα consists of a cooperative folding unit of roughly four ARs (AR1-AR4) and of two weakly folded repeats (AR5 and AR6). The kinetic folding mechanism of the cooperative subdomain, IκBα(67-206), was analyzed using rapid mixing techniques. Despite its apparent architectural simplicity, IκBα(67-206) displays complex folding kinetics, with two sequential on-pathway high-energy intermediates. The effect of mutations to or away from the consensus sequences of ARs on folding behavior was analyzed, particularly the GXTPLHLA motif, which have not been examined in detail previously. Mutations toward the consensus generally resulted in an increase in folding stability, whereas mutations away from the consensus resulted in decreased overall stability. We determined the free energy change upon mutation for three sequential transition state ensembles along the folding route for 16 mutants. We show that folding initiates with the formation of the interface of the outer helices of AR3 and AR4, and then proceeds to consolidate structure in these repeats. Subsequently, AR1 and AR2 fold in a concerted way in a single kinetic step. We show that this mechanism is robust to the presence of AR5 and AR6 as they do not strongly affect the folding kinetics. Overall, the protein appears to fold on a rather smooth energy landscape, where the folding mechanism conforms a one-dimensional approximation. However, we note that the AR does not necessarily act as a single folding element.

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.

Fundamental link between folding states and functional states of proteins

Folding and function of proteins are two aspects of proteins which are usually considered as basically unrelated phenomena that are optimized by evolution independently. From the funnel model of folding/unfolding and the associated energy landscape, we infer the paradigm that the minimum number of folding intermediates is determined by the number of all functional states of a protein ("essential" folding intermediates). Here, we demonstrate the supposed fundamental link using the Ras protein complexed with the GTP analogue GppNHp that occurs in two structural states coexisting in solution. State 2 was shown earlier to represent the effector interacting state, and the function of state 1 was hitherto unknown. By (31)P NMR spectroscopy, we demonstrate that state 1 represents the conformation interacting with guanine nucleotide exchange factors (GEFs). Denaturation experiments of the protein with a chaotropic reagent show that both functional states coexist during folding and unfolding. Application of high pressure represents another perturbation of the energy landscape, leading to an increased population of the state 1 as observed by NMR spectroscopy. The specific volume difference between the two states DeltaV(12) is 17.2 +/- 0.5 mL mol(-1), indicating that state 1 represents a more open conformation of the protein. The free energies of stabilization for state 1 and state 2 at 278 K can be determined as 8.3 and 9.8 kJ mol(-1), respectively.

Experimental evidence for a frustrated energy landscape in a three-helix-bundle protein family

Energy landscape theory is a powerful tool for understanding the structure and dynamics of complex molecular systems, in particular biological macromolecules. The primary sequence of a protein defines its free-energy landscape and thus determines the folding pathway and the rate constants of folding and unfolding, as well as the protein's native structure. Theory has shown that roughness in the energy landscape will lead to slower folding, but derivation of detailed experimental descriptions of this landscape is challenging. Simple folding models show that folding is significantly influenced by chain entropy; proteins in which the contacts are local fold quickly, owing to the low entropy cost of forming stabilizing, native contacts during folding. For some protein families, stability is also a determinant of folding rate constants. Where these simple metrics fail to predict folding behaviour, it is probable that there are features in the energy landscape that are unusual. Such general observations cannot explain the folding behaviour of the R15, R16 and R17 domains of alpha-spectrin. R15 folds approximately 3,000 times faster than its homologues, although they have similar structures, stabilities and, as far as can be determined, transition-state stabilities. Here we show that landscape roughness (internal friction) is responsible for the slower folding and unfolding of R16 and R17. We use chimaeric domains to demonstrate that this internal friction is a property of the core, and suggest that frustration in the landscape of the slow-folding spectrin domains may be due to misdocking of the long helices during folding. Theoretical studies have suggested that rugged landscapes will result in slower folding; here we show experimentally that such a phenomenon directly influences the folding kinetics of a 'normal' protein, that is, one with a significant energy barrier that folds on a relatively slow, millisecond-second, timescale.

Distinguishing between smooth and rough free energy barriers in protein folding

Analysis of curved chevron plots is a powerful tool in investigating protein folding pathways, as the curvatures can be used to gain information about both early and late folding events. When and if accumulation of low-energy intermediates can be ruled out, two different models have classically been applied to describe curved chevron plots, namely , (i) Hammond effects along smooth barrier profiles and (ii) changes in the rate-limiting step between two discrete transition states. The two models lead to very similar numerical solutions, which are generally indistinguishable. This is not surprising, since the smooth barrier assumption approximates barrier profiles with a more complex topology involving multiple local maxima that are too close, or too broad, to yield clear-cut kinks in the chevron data. In this work, we have reconstructed the transition state shifts as a function of protein stability over a wide stability range for three small globular proteins, to screen for fingerprints more sensitive for different barrier profiles. We show that such an analysis represents a valuable test for the discrimination between the two different scenarios.

Multiple Unfolding Intermediates Obtained by Molecular Dynamic Simulations under Stretching for Immunoglobulin-Binding Domain of Protein G

We have studied the mechanical properties of the immunoglobulin-binding domain of protein G at the atomic level under stretching at constant velocity using molecular dynamics simulations. We have found that the unfolding process can occur either in a single step or through intermediate states. Analysis of the trajectories from the molecular dynamic simulations showed that the mechanical unfolding of the immunoglobulin-binding domain of protein G is triggered by the separation of the terminal beta-strands and the order in which the secondary-structure elements break is practically the same in two- and multi-state events and at the different extension velocities studied. It is seen from our analysis of 24 trajectories that the theoretical pathway of mechanical unfolding for the immunoglobulin-binding domain of protein G does not coincide with that proposed in denaturant studies in the absence of force.

Folding mechanism of reduced Cytochrome c: equilibrium and kinetic properties in the presence of carbon monoxide

Despite close structural similarity, the ferric and ferrous forms of cytochrome c differ greatly in terms of their ligand binding properties, stability, folding, and dynamics. The reduced heme iron binds diatomic ligands such as CO only under destabilizing conditions that promote weakening or disruption of native methionine-iron linkage. This makes CO a useful conformational probe for detecting partially structured states that cannot be observed in the absence of endogenous ligands. Heme absorbance, circular dichroism, and NMR were used to characterize the denaturant-induced unfolding equilibrium of ferrocytochrome c in the presence and in the absence of CO. In addition to the native state (N), which does not bind CO, and the unfolded CO complex (U-CO), a structurally distinct CO-bound form (M-CO) accumulates to high levels (approximately 75% of the population) at intermediate guanidine HCl concentrations. Comparison of the unfolding transitions for different conformational probes reveals that M-CO is a compact state containing a native-like helical core and regions of local disorder in the segment containing the native Met80 ligand and adjacent loops. Kinetic measurements of CO binding and dissociation under native, partially denaturing, and fully unfolded conditions indicate that a state M that is structurally analogous to M-CO is populated even in the absence of CO. The binding energy of the CO ligand lowers the free energy of this high-energy state to such an extent that it accumulates even under mildly denaturing equilibrium conditions. The thermodynamic and kinetic parameters obtained in this study provide a fully self-consistent description of the linked unfolding/CO binding equilibria of reduced cytochrome c.

Shifting transition states in the unfolding of a large ankyrin repeat protein

The 33-amino-acid ankyrin motif comprises a beta-turn followed by two anti-parallel alpha-helices and a loop and tandem arrays of the motif pack in a linear fashion to produce elongated structures characterized by short-range interactions. In this article we use site-directed mutagenesis to investigate the kinetic unfolding mechanism of D34, a 426-residue, 12-ankyrin repeat fragment of the protein ankyrinR. The data are consistent with a model in which the N-terminal half of the protein unfolds first by unraveling progressively from the start of the polypeptide chain to form an intermediate; in the next step, the C-terminal half of the protein unfolds via two pathways whose transition states have either the early or the late C-terminal ankyrin repeats folded. We conclude that the two halves of the protein unfold by different mechanisms because the N-terminal moiety folds and unfolds in the context of a folded C-terminal moiety, which therefore acts as a "seed" and confers a unique directionality on the process, whereas the C-terminal moiety folds and unfolds in the context of an unfolded N-terminal moiety and therefore behaves like a single-domain ankyrin repeat protein, having a high degree of symmetry and consequently more than one unfolding pathway accessible to it.

Conformational stability and folding mechanisms of dimeric proteins

The folding of multisubunit proteins is of tremendous biological significance since the large majority of proteins exist as protein-protein complexes. Extensive experimental and computational studies have provided fundamental insights into the principles of folding of small monomeric proteins. Recently, important advances have been made in extending folding studies to multisubunit proteins, in particular homodimeric proteins. This review summarizes the equilibrium and kinetic theory and models underlying the quantitative analysis of dimeric protein folding using chemical denaturation, as well as the experimental results that have been obtained. Although various principles identified for monomer folding also apply to the folding of dimeric proteins, the effects of subunit association can manifest in complex ways, and are frequently overlooked. Changes in molecularity typically give rise to very different overall folding behaviour than is observed for monomeric proteins. The results obtained for dimers have provided key insights pertinent to understanding biological assembly and regulation of multisubunit proteins. These advances have set the stage for future advances in folding involving protein-protein interactions for natural multisubunit proteins and unnatural assemblies involved in disease.

The influence of proline isomerization and off-pathway intermediates on the folding mechanism of eukaryotic UMP/CMP Kinase

The globular 22-kDa protein UMP/CMP from Dictyostelium discoideum (UmpK) belongs to the family of nucleoside monophosphate (NMP) kinases. These enzymes not only show high sequence and structure similarities but also share the alpha/beta-fold, a very common protein topology. We investigated the protein folding mechanism of UmpK as a representative for this ubiquitous enzyme class. Equilibrium stability towards urea and the unfolding and refolding kinetics were studied by means of fluorescence and far-UV CD spectroscopy. Although the unfolding can be described by a two-state process, folding kinetics are rather complex with four refolding phases that can be resolved and an additional burst phase. Moreover, two of these phases exhibit a pronounced rollover in the refolding limb that cannot be explained by aggregation. Whilst secondary structure formation is not observed in the burst phase reaction, folding to the native structure is strongly influenced by the slowest phase, since 30% of the alpha-helical CD signal is restored therein. This process can be assigned to proline isomerization and is strongly accelerated by the Escherichia coli peptidyl-prolyl isomerase trigger factor. The analysis of our single-mixing and double-mixing experiments suggests the occurrence of an off-pathway intermediate and an unproductive collapsed structure, which appear to be rate limiting for the folding of UmpK.

The folding pathway of a single domain in a multidomain protein is not affected by its neighbouring domain

Domains are the structural, functional, and evolutionary components of proteins. Most folding studies to date have concentrated on the folding of single domains, but more than 70% of human proteins contain more than one domain, and interdomain interactions can affect both the stability and the folding kinetics. Whether the folding pathway is altered by interdomain interactions is not yet known. Here we investigated the effect of a folded neighbouring domain on the folding pathway of spectrin R16 (the 16th α-helical repeat from chicken brain α-spectrin) by using the two-domain construct R1516. The R16 folds faster and unfolds more slowly in the presence of its folded neighbour R15 (the 15th α-helical repeat from chicken brain α-spectrin). An extensive Φ-value analysis of the R16 domain in R1516 was completed to compare the transition state of the R16 domain alone with that of the R16 domain in a multidomain construct. The results indicate that the folding pathways are the same. This result validates the current approach of breaking up larger proteins into domains for the study of protein folding pathways.

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.