Formation of short-lived protein aggregates directly from the coil in two-state folding

Assessing the impact of choline chloride and benzyltrimethylammonium chloride-based deep eutectic solvents on the structure and conformational dynamics of bovine serum albumin: a combined steady-state, time-resolved fluorescence and fluorescence correlation spectroscopic study

Although deep eutectic solvents (DESs) are regarded as useful substitutes for both ionic liquids and common organic solvents for storage and applications of biomolecules, it is still unclear whether all DESs or only specific types of DESs will be suitable for the said purpose. In view of this, the current study aims to report on the structure and conformational dynamics of BSA in the presence of two DESs, namely ethaline (choline chloride:ethylene glycol) and BMEG (benzyltrimethyl ammonium chloride:ethylene glycol), having the same hydrogen bond donor but with a distinct hydrogen bond acceptor, so that how small changes in one constituent of a DES alter the protein-DES interaction at the molecular level can be understood. The protein-DES interaction is investigated by exploiting both ensemble-averaged measurements like steady-state and time-resolved fluorescence spectroscopy, circular dichroism (CD) spectroscopy, and single-molecule sensitive techniques based on fluorescence correlation spectroscopy (FCS). Interestingly, the results obtained from these studies have demonstrated that while a very small quantity of BMEG completely unfolds the native structure of the protein, it remains in a partially unfolded state even at very high ethaline content. More interestingly, it has been found that at very high concentrations of BMEG, the unfolded protein undergoes enhanced protein-protein interaction resulting in the aggregation of BSA. All of the results obtained from these investigations have essentially suggested that both protein-DES interaction and interspecies interaction among the constituent of DESs play a crucial role in governing the overall stability and conformational dynamics of the protein in DESs.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

Quantum Dots as Promising Theranostic Tools Against Amyloidosis: A Review

Amyloids are highly ordered beta sheet rich stable protein aggregates, which have been found to play a significant role in the onset of several degenerative diseases such as Alzheimer's disease, Huntington's disease, Parkinson's disease, Type II diabetes mellitus and so on. Aggregation of proteins leading to amyloid fibril formation via intermediate(s), is thought to be a nucleated condensation polymerization process associated with many pathological conditions. There has been extensive research to identify inhibitors of these disease oriented aggregation processes. In recent times, quantum dots, with their unique physico-chemical properties have grabbed the attention of scientific community due to its applications in medical sciences. Quantum dots are nano-particles usually made of semiconductor materials which emit fluorescence upon radiation. The wavelength of fluorescence emission varies with changes in size of quantum dots. Several studies have reported significant inhibitory effects of these quantum dots towards amyloidogenesis, thereby presenting themselves as promising candidates against amyloidosis. Further, studies have also revealed amyloid detection capacity of quantum dots with sensitivity and specificity better than conventional probes. In the current review, we will discuss the various effects of quantum dots on protein aggregation pathways, their mechanism of actions and their potentials as effective therapeutics against amyloidosis.

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.

Anesthetic 2,2,2-trifluoroethanol induces amyloidogenesis and cytotoxicity in human serum albumin

Trifluoroethanol (TFE) mimics the membrane environments as it simulates the hydrophobic environment and better stabilizes the secondary structures in peptides owing to its hydrophobicity and hydrogen bond-forming properties. Its dielectric constant approximates that of the interior of proteins and is one-third of that of water. Human serum albumin (HSA) is a biological transporter. The effect of TFE on HSA gives the clue about the conformational changes taking place in HSA on transport of ligands across the biological membranes. At 25% (v/v) and 60% (v/v) TFE, HSA exhibits non-native β-sheet, altered tryptophan fluorescence, exposed hydrophobic clusters, increased thioflavin T fluorescence and prominent red shifted Congo red absorbance, and large hydrodynamic radii suggesting the aggregate formation. Isothermal titration calorimetric results indicate weak binding of TFE and HSA. This suggests that solvent-mediated effects dominate the interaction of TFE and HSA. TEM confirmed prefibrillar at 25% (v/v) and fibrillar aggregates at 60% (v/v) TFE. Comet assay of prefibrillar aggregates showed DNA damage causing cell necrosis hence confirming cytotoxic nature. On increasing concentration of TFE to 80% (v/v), HSA showed retention of native-like secondary structure, increased Trp and ANS fluorescence, a transition from β-sheet to α-helix. Thus, TFE at high concentration possess anti- aggregating potency.

Structural transformation of bovine serum albumin induced by dimethyl sulfoxide and probed by fluorescence correlation spectroscopy and additional methods

Determining the structure of a protein and its transformation under different conditions is key to understanding its activity. The structural stability and activity of proteins in aqueous-organic solvent mixtures, which is an intriguing topic of research in biochemistry, is dependent on the nature of the protein and the properties of the medium. Herein, the effect of a commonly used cosolvent, dimethyl sulfoxide (DMSO), on the structure and conformational dynamics of bovine serum albumin (BSA) protein is studied by fluorescence correlation spectroscopy (FCS) measurements on fluorescein isothiocyanate (FITC)-labeled BSA. The FCS study reveals a change of the hydrodynamic radius of BSA from 3.7 nm in the native state to 7.0 nm in the presence of 40% DMSO, which suggests complete unfolding of the protein under these conditions. Fluorescence self-quenching of FITC has been exploited to understand the conformational dynamics of BSA. The time constant of the conformational dynamics of BSA is found to change from 35 μs in its native state to 50 μs as the protein unfolds with increasing DMSO concentration. The FCS results are corroborated by the near-UV circular dichroism spectra of the protein, which suggest a loss of its tertiary structure with increasing concentration of DMSO. The intrinsic fluorescence of BSA and the fluorescence response of 1-anilinonaphthalene-8-sulfonic acid, used as a probe molecule, provide information that is consistent with the FCS measurements, except that aggregation of BSA is observed in the presence of 40% DMSO in the ensemble measurements.

Early aggregated States in the folding of interleukin-1β

Kinetic data measured from folding of the protein interleukin-1β fits best to three exponential phases when studied with tryptophan fluorescence but only two exponential phases when measured using other methods. The technique of ANS fluorescence was used to determine whether the additional phase observed in tryptophan fluorescence was also detected with ANS dye binding. Unlike trytophan fluorescence, the ANS fluorescence was highly dependent on the concentration of protein present during the folding experiment. Experimental controls provide evidence that ANS binds to protein aggregates, present at higher concentrations and absent at lower concentrations. Protein concentration-dependent folding studies demonstrate that, at lower interleukin-1β concentrations, tryptophan fluorescence kinetics can be fit adequately with a two exponential fit. This study indicates that (1) measured interleukin-1β folding kinetics fit to a 2 phase model and (2) at higher protein concentrations, transient association of IL-1β may result in a kinetic fit of 3 phases.

Frustration in the energy landscapes of multidomain protein misfolding

Frustration from strong interdomain interactions can make misfolding a more severe problem in multidomain proteins than in single-domain proteins. On the basis of bioinformatic surveys, it has been suggested that lowering the sequence identity between neighboring domains is one of nature's solutions to the multidomain misfolding problem. We investigate folding of multidomain proteins using the associative-memory, water-mediated, structure and energy model (AWSEM), a predictive coarse-grained protein force field. We find that reducing sequence identity not only decreases the formation of domain-swapped contacts but also decreases the formation of strong self-recognition contacts between β-strands with high hydrophobic content. The ensembles of misfolded structures that result from forming these amyloid-like interactions are energetically disfavored compared with the native state, but entropically favored. Therefore, these ensembles are more stable than the native ensemble under denaturing conditions, such as high temperature. Domain-swapped contacts compete with self-recognition contacts in forming various trapped states, and point mutations can shift the balance between the two types of interaction. We predict that multidomain proteins that lack these specific strong interdomain interactions should fold reliably.

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.

Complexes of native ubiquitin and dodecyl sulfate illustrate the nature of hydrophobic and electrostatic interactions in the binding of proteins and surfactants

A previous study, using capillary electrophoresis (CE) [J. Am. Chem. Soc. 2008, 130, 17384-17393], reported that six discrete complexes of ubiquitin (UBI) and sodium dodecyl sulfate (SDS) form at different concentrations of SDS along the pathway to unfolding of UBI in solutions of SDS. One complex (which formed between 0.8 and 1.8 mM SDS) consisted of native UBI associated with approximately 11 molecules of SDS. The current study used CE and (15)N/(13)C-(1)H heteronuclear single quantum coherence (HSQC) NMR spectroscopy to identify residues in folded UBI that associate specifically with SDS at 0.8-1.8 mM SDS, and to correlate these associations with established biophysical and structural properties of this well-characterized protein. The ability of the surface charge and hydrophobicity of folded UBI to affect the association with SDS (at concentrations below the CMC) was studied, using CE, by converting lys-ε-NH(3)(+) to lys-ε-NHCOCH(3) groups. According to CE, the acetylation of lysine residues inhibited the binding of 11 SDS ([SDS] < 2 mM) and decreased the number of complexes of composition UBI-(NHAc)(8)·SDS(n) that formed on the pathway of unfolding of UBI-(NHAc)(8) in SDS. A comparison of (15)N-(1)H HSQC spectra at 0 mM and 1 mM SDS with calculated electrostatic surface potentials of folded UBI (e.g., solutions to the nonlinear Poisson-Boltzmann (PB) equation) suggested, however, that SDS binds preferentially to native UBI at hydrophobic residues that are formally neutral (i.e., Leu and Ile), but that have positive electrostatic surface potential (as predicted from solutions to nonlinear PB equations); SDS did not uniformly interact with residues that have formal positive charge (e.g., Lys or Arg). Cationic functional groups, therefore, promote the binding of SDS to folded UBI because these groups exert long-range effects on the positive electrostatic surface potential (which extend beyond their own van der Waals radii, as predicted from PB theory), and not because cationic groups are necessarily the site of ionic interactions with sulfate groups. Moreover, SDS associated with residues in native UBI without regard to their location in α-helix or β-sheet structure (although residues in hydrogen-bonded loops did not bind SDS). No correlation was observed between the association of an amino acid with SDS and the solvent accessibility of the residue or its rate of amide H/D exchange. This study establishes a few (of perhaps several) factors that control the simultaneous molecular recognition of multiple anionic amphiphiles by a folded cytosolic protein.

Nucleation process in the folding of a domain-swapped dimer

Nucleation processes are important for the understanding in protein dynamics. To evaluate the effect of nucleation mechanism in dimerization process, a domain-swapped dimer (Esp8) is simulated with the symmetrized Gō model and the classical Gō model. The pathways of the dimerization are analyzed with computational phi -analysis method. It is found out that some nuclei are observed in the kinetic steps of the dimeric association though the whole pathway is a process with multiple intermediate states. The key residues in the nuclei are rather similar to those observed in the monomeric folding. The differences with the monomeric cases are also discussed. These differences illustrate the effects of dimeric feature on the nucleation process. Besides, manual mutations are carried out to illustrate the importance of the interactions related to the nuclei. It is observed that the mutations in the nuclei-related interactions apparently change the dynamics while other mutations have little effect on the kinetics. All of these results outline a picture that the nucleation processes act as the fundamental steps of high-order organization of protein systems.

Competition between reversible aggregation and loop formation in denatured iso-1-cytochrome c

The competition between intramolecular histidine-heme loop formation and ligand-mediated oligomer formation in the denatured state is investigated for two yeast iso-1-cytochrome c variants, AcH26I52 and AcA25H26I52. Besides the native His 18 heme ligand, both variants contain a single His at position 26. The AcA25H26I52 variant has Pro 25 mutated to Ala. The concentration dependence of the apparent pK(a) for His 26-heme binding in 3 M guanidine hydrochloride indicates that the P25A mutation disfavors oligomerization mediated by intermolecular heme ligation by 10-fold. Single- and double-pH-jump stopped-flow experiments with the AcH26I52 variant show that fast phases for His-heme bond formation and breakage are due to intramolecular loop formation and slow phases for His-heme bond formation and breakage are due to intermolecular aggregation. The presence of two closely spaced slow phases in the kinetics of loop formation for both variants suggests that intermolecular His 26-heme ligation results in both dimers and higher-order aggregates. The P25A mutation slows formation and speeds breakdown of an initial dimer, demonstrating a strong effect of local sequence on aggregation. Analysis of the kinetic data yields equilibrium constants for intramolecular loop formation and intermolecular dimerization at pH 7.1 and indicates that the rate constant for intermolecular aggregation is very fast at this pH (10(7)-10(8) M(-1) s(-1)). In light of the very fast rates of aggregation in the denatured state, comparison of models involving reversible or irreversible oligomerization steps suggests that equilibrium control of the partitioning between folding and aggregation is advantageous for productive protein folding in vivo.

SOD1-associated ALS: a promising system for elucidating the origin of protein-misfolding disease

Amyotropic lateral sclerosis (ALS) is a neurodegenerative disease linked to misfolding and aggregation of the homodimeric enzyme superoxide dismutase (SOD1). In contrast to the precursors of other neurodegenerative diseases, SOD1 is a soluble and simple-to-study protein with immunoglobulin-like structure. Also, there are more than 120 ALS-provoking SOD1 mutations at the disposal for detailed elucidation of the disease-triggering factors at molecular level. In this article, we review recent progress in the characterization of the folding and assembly pathway of the SOD1 dimer and how this is affected by ALS-provoking mutations. Despite the diverse nature of these mutations, the results offer so far a surprising simplicity. The ALS-provoking mutations decrease either protein stability or net repulsive charge: the classical hallmarks for a disease mechanism triggered by association of non-native protein. In addition, the mutant data identifies immature SOD1 monomers as the species from which the cytotoxic pathway emerges, and point at compromised folding cooperativity as a key disease determinant. The relative ease by which these data can be obtained makes SOD1 a promising model for elucidating also the origin of other neurodegenerative diseases where the precursor proteins are structurally more elusive.

The hyperthermophilic nature of the metallo-oxidase from Aquifex aeolicus

The stability of the Aquifex aeolicus multicopper oxidase (McoA) was studied by spectroscopy, calorimetry and chromatography to understand its thermophilic nature. The enzyme is hyperthermostable as deconvolution of the differential scanning calorimetry trace shows that thermal unfolding is characterized by temperature values at the mid-point of 105, 110 and 114 degrees C. Chemical denaturation revealed however a very low stability at room temperature (2.8 kcal/mol) because copper bleaching/depletion occur before the unfolding of the tertiary structure and McoA is highly prone to aggregate. Indeed, unfolding kinetics measured with the stopped-flow technique quantified the stabilizing effect of copper on McoA (1.5 kcal/mol) and revealed quite an uncommon observation further confirmed by light scattering and gel filtration chromatography: McoA aggregates in the presence of guanidinium hydrochloride, i.e., under unfolding conditions. The aggregation process results from the accumulation of a quasi-native state of McoA that binds to ANS and is the main determinant of the stability curve of McoA. Kinetic partitioning between aggregation and unfolding leads to a very low heat capacity change and determines a flat dependence of stability on temperature.

Contrasting disease and nondisease protein aggregation by molecular simulation

[Figurre: see text]. Protein aggregation can be defined as the sacrifice of stabilizing intrachain contacts of the functional state that are replaced with interchain contacts to form non-functional states. The resulting aggregate morphologies range from amorphous structures without long-range order typical of nondisease proteins involved in inclusion bodies to highly structured fibril assemblies typical of amyloid disease proteins. In this Account, we describe the development and application of computational models for the investigation of nondisease and disease protein aggregation as illustrated for the proteins L and G and the Alzheimer's Abeta systems. In each case, we validate the models against relevant experimental observables and then expand on the experimental window to better elucidate the link between molecular properties and aggregation outcomes. Our studies show that each class of protein exhibits distinct aggregation mechanisms that are dependent on protein sequence, protein concentration, and solution conditions. Nondisease proteins can have native structural elements in the denatured state ensemble or rapidly form early folding intermediates, which offers avenues of protection against aggregation even at relatively high concentrations. The possibility that early folding intermediates may be evolutionarily selected for their protective role against unwanted aggregation could be a useful strategy for reengineering sequences to slow aggregation and increase folding yield in industrial protein production. The observed oligomeric aggregates that we see for nondisease proteins L and G may represent the nuclei for larger aggregates, not just for large amorphous inclusion bodies, but potentially as the seeds of ordered fibrillar aggregates, since most nondisease proteins can form amyloid fibrils under conditions that destabilize the native state. By contrast, amyloidogenic protein sequences such as Abeta 1-40,42 and the familial Alzheimer's disease (FAD) mutants favor aggregation into ordered fibrils once the free-energy barrier for forming a critical nucleus is crossed. However, the structural characteristics and oligomer size of the soluble nucleation species have yet to be determined experimentally for any disease peptide sequence, and the molecular mechanism of polymerization that eventually delineates a mature fibril is unknown. This is in part due to the limited experimental access to very low peptide concentrations that are required to characterize these early aggregation events, providing an opportunity for theoretical studies to bridge the gap between the monomer and fibril end points and to develop testable hypotheses. Our model shows that Abeta 1-40 requires as few as 6-10 monomer chains (depending on sequence) to begin manifesting the cross-beta order that is a signature of formation of amyloid filaments or fibrils assessed in dye-binding kinetic assays. The richness of the oligomeric structures and viable filament and fibril polymorphs that we observe may offer structural clues to disease virulence variations that are seen for the WT and hereditary mutants.

Simulation-based fitting of protein-protein interaction potentials to SAXS experiments

We present a new method for computing interaction potentials of solvated proteins directly from small-angle x-ray scattering data. An ensemble of proteins is modeled by Monte Carlo or molecular dynamics simulation. The global x-ray scattering of the whole model ensemble is then computed at each snapshot of the simulation, and averaged to obtain the x-ray scattering intensity. Finally, the interaction potential parameters are adjusted by an optimization algorithm, and the procedure is iterated until the best agreement between simulation and experiment is obtained. This new approach obviates the need for approximations that must be made in simplified analytical models. We apply the method to lambda repressor fragment 6-85 and fyn-SH3. With the increased availability of fast computer clusters, Monte Carlo and molecular dynamics analysis using residue-level or even atomistic potentials may soon become feasible.

N-terminal domains of native multidomain proteins have the potential to assist de novo folding of their downstream domains in vivo by acting as solubility enhancers

The fusion of soluble partner to the N terminus of aggregation-prone polypeptide has been popularly used to overcome the formation of inclusion bodies in the E. coli cytosol. The chaperone-like functions of the upstream fusion partner in the artificial multidomain proteins could occur in de novo folding of native multidomain proteins. Here, we show that the N-terminal domains of three E. coli multidomain proteins such as lysyl-tRNA synthetase, threonyl-tRNA synthetase, and aconitase are potent solubility enhancers for various C-terminal heterologous proteins. The results suggest that the N-terminal domains could act as solubility enhancers for the folding of their authentic C-terminal domains in vivo. Tandem repeat of N-terminal domain or insertion of aspartic residues at the C terminus of the N-terminal domain also increased the solubility of fusion proteins, suggesting that the solubilizing ability correlates with the size and charge of N-terminal domains. The solubilizing ability of N-terminal domains would contribute to the autonomous folding of multidomain proteins in vivo, and based on these results, we propose a model of how N-terminal domains solubilize their downstream domains.

The Highly Cooperative Folding of Small Naturally Occurring Proteins Is Likely the Result of Natural Selection

To illuminate the evolutionary pressure acting on the folding free energy landscapes of naturally occurring proteins, we have systematically characterized the folding free energy landscape of Top7, a computationally designed protein lacking an evolutionary history. Stopped-flow kinetics, circular dichroism, and NMR experiments reveal that there are at least three distinct phases in the folding of Top7, that a nonnative conformation is stable at equilibrium, and that multiple fragments of Top7 are stable in isolation. These results indicate that the folding of Top7 is significantly less cooperative than the folding of similarly sized naturally occurring proteins, suggesting that the cooperative folding and smooth free energy landscapes observed for small naturally occurring proteins are not general properties of polypeptide chains that fold to unique stable structures but are instead a product of natural selection.

Detection of an intermediate during the unfolding process of the dimeric ketosteroid isomerase

Failure to detect the intermediate in spite of its existence often leads to the conclusion that two-state transition in the unfolding process of the protein can be justified. In contrast to the previous equilibrium unfolding experiment fitted to a two-state model by circular dichroism and fluorescence spectroscopies, an equilibrium unfolding intermediate of a dimeric ketosteroid isomerase (KSI) could be detected by small angle X-ray scattering (SAXS) and analytical ultracentrifugation. The sizes of KSI were determined to be 18.7A in 0M urea, 17.3A in 5.2M urea, and 25.1A in 7M urea by SAXS. The size of KSI in 5.2M urea was significantly decreased compared with those in 0M and 7M urea, suggesting the existence of a compact intermediate. Sedimentation velocity as obtained by ultracentrifugation confirmed that KSI in 5.2M urea is distinctly different from native and fully-unfolded forms. The sizes measured by pulse field gradient nuclear magnetic resonance (NMR) spectroscopy were consistent with those obtained by SAXS. Discrepancy of equilibrium unfolding studies between size measurement methods and optical spectroscopies might be due to the failure in detecting the intermediate by optical spectroscopic methods. Further characterization of the intermediate using (1)H NMR spectroscopy and Kratky plot supported the existence of a partially-folded form of KSI which is distinct from those of native and fully-unfolded KSIs. Taken together, our results suggest that the formation of a compact intermediate should precede the association of monomers prior to the dimerization process during the folding of KSI.

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

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

Binary and ternary aggregation within tethered protein constructs

The free energy per monomer of a protein aggregate varies with the number of participating monomers n. The change of this free energy with aggregate size, DeltaDeltaG(n), is difficult to determine by sedimentation or concentration studies. We introduce a kinetic approach to quantitate the free energy of aggregates in the presence of tethers. By linking the protein U1A into dimers and trimers, a high effective concentration of the monomers is achieved, together with exact size control of the aggregates. We found that the free energy of the aggregate relative to the native monomer reached a maximum for n = 2, and decreased by DeltaDeltaG(2) = -3.1 kT between dimer and trimer.

The importance of sequence diversity in the aggregation and evolution of proteins

Incorrect folding of proteins, leading to aggregation and amyloid formation, is associated with a group of highly debilitating medical conditions including Alzheimer's disease and late-onset diabetes. The issue of how unwanted protein association is normally avoided in a living system is particularly significant in the context of the evolution of multidomain proteins, which account for over 70% of all eukaryotic proteins, where the effective local protein concentration in the vicinity of each domain is very high. Here we describe the aggregation kinetics of multidomain protein constructs of immunoglobulin domains and the ability of different homologous domains to aggregate together. We show that aggregation of these proteins is a specific process and that the efficiency of coaggregation between different domains decreases markedly with decreasing sequence identity. Thus, whereas immunoglobulin domains with more than about 70% identity are highly prone to coaggregation, those with less than 30-40% sequence identity do not detectably interact. A bioinformatics analysis of consecutive homologous domains in large multidomain proteins shows that such domains almost exclusively have sequence identities of less than 40%, in other words below the level at which coaggregation is likely to be efficient. We propose that such low sequence identities could have a crucial and general role in safeguarding proteins against misfolding and aggregation.

Influence of denatured and intermediate states of folding on protein aggregation

We simulate the aggregation thermodynamics and kinetics of proteins L and G, each of which self-assembles to the same alpha/beta [corrected] topology through distinct folding mechanisms. We find that the aggregation kinetics of both proteins at an experimentally relevant concentration exhibit both fast and slow aggregation pathways, although a greater proportion of protein G aggregation events are slow relative to those of found for protein L. These kinetic differences are correlated with the amount and distribution of intrachain contacts formed in the denatured state ensemble (DSE), or an intermediate state ensemble (ISE) if it exists, as well as the folding timescales of the two proteins. Protein G aggregates more slowly than protein L due to its rapidly formed folding intermediate, which exhibits native intrachain contacts spread across the protein, suggesting that certain early folding intermediates may be selected for by evolution due to their protective role against unwanted aggregation. Protein L shows only localized native structure in the DSE with timescales of folding that are commensurate with the aggregation timescale, leaving it vulnerable to domain swapping or nonnative interactions with other chains that increase the aggregation rate. Folding experiments that characterize the structural signatures of the DSE, ISE, or the transition state ensemble (TSE) under nonaggregating conditions should be able to predict regions where interchain contacts will be made in the aggregate, and to predict slower aggregation rates for proteins with contacts that are dispersed across the fold. Since proteins L and G can both form amyloid fibrils, this work also provides mechanistic and structural insight into the formation of prefibrillar species.

Time-dependent changes in the denatured state(s) influence the folding mechanism of an all beta-sheet protein

Newt fibroblast growth factor (nFGF-1) is an approximately 15-kDa all beta-sheet protein devoid of disulfide bonds. Urea-induced equilibrium unfolding of nFGF-1, monitored by steady state fluorescence and far-UV circular dichroism spectroscopy, is cooperative with no detectable intermediate(s). Urea-induced unfolding of nFGF-1 is reversible, but the percentage of the protein recovered in the native state depends on the time of incubation of the protein in the denaturant. The yield of the protein in the native state decreases with the increase in time of incubation in the denaturant. The failure of the protein to refold to its native state is not due to trivial chemical reactions that could possibly occur upon prolonged incubation in the denaturant. (1)H-(15)N heteronuclear single quantum coherence (HSQC) spectra, limited proteolytic digestion, and fluorescence data suggest that the misfolded state(s) of nFGF-1 has structural features resembling that of the denatured state(s). GroEL, in the presence of ATP, is observed to rescue the protein from being trapped in the misfolded state(s). (1)H-(15)N HSQC data of nFGF-1, acquired in the denatured state(s) (in 8 m urea), suggest that the protein undergoes subtle time-dependent structural changes in the denaturant. To our knowledge, this report for the first time demonstrates that the commitment to adapt unproductive pathways leading to protein misfolding/aggregation occurs in the denatured state ensemble.

Pathway complexity of Alzheimer's beta-amyloid Abeta16-22 peptide assembly

Recent studies suggest that both soluble oligomers and insoluble fibrils have toxic effects in cell cultures, raising the interest in determining the first steps of the assembly process. We have determined the aggregation mechanisms of Abeta(16-22) dimer using the activation-relaxation technique and an approximate free energy model. Consistent with the NMR solid-state analysis, the dimer is predicted to prefer an antiparallel beta sheet structure with the expected registry of intermolecular hydrogen bonds. The simulations, however, locate three other antiparallel minima with nonnative beta sheet registries and one parallel beta sheet structure, slightly destabilized with respect to the ground state. This result is significant because it can explain the observed dependency of beta sheet registry on pH conditions. We also find that assembly of Abeta(16-22) into dimers follows multiple routes, but alpha-helical intermediates are not obligatory. This indicates that destabilization of alpha-helical intermediates is unlikely to abolish oligomerization of Abeta peptides.

Detection and characterization of partially unfolded oligomers of the SH3 domain of alpha-spectrin

For the purpose of equilibrium and kinetic folding-unfolding studies, the SH3 domain of alpha-spectrin (spc-SH3) has long been considered a classic two-state folding protein. In this work we have indeed observed that the thermal unfolding curves of spc-SH3 measured at pH 3.0 by differential scanning calorimetry, circular dichroism, and NMR follow apparently the two-state model when each unfolding profile is considered individually. Nevertheless, we have found that protein concentration has a marked effect upon the thermal unfolding profiles. This effect cannot be properly explained in terms of the two-state unfolding model and can only be interpreted in terms of the accumulation of intermediate associated states in equilibrium with the monomeric native and unfolded states. By chemical cross-linking and pulsed-field gradient NMR diffusion experiments we have been able to confirm the existence of associated states formed during spc-SH3 unfolding. A three-state model, in which a dimeric intermediate state is assumed to be significantly populated, provides the simplest interpretation of the whole set of thermal unfolding data and affords a satisfactory explanation for the concentration effects observed. Whereas at low concentrations the population of the associated intermediate state is negligible and the unfolding process consequently takes place in a two-state fashion, at concentrations above approximately 0.5 mM the population of the intermediate state becomes significant at temperatures between 45 degrees C and 80 degrees C and reaches up to 50% at the largest concentration investigated. The thermodynamic properties of the intermediate state implied by this analysis fall in between those of the unfolded state and the native ones, indicating a considerably disordered conformation, which appears to be stabilized by oligomerization.

Trapping the on-pathway folding intermediate of Im7 at equilibrium

The four-helical protein Im7 folds via a rapidly formed on-pathway intermediate (k(UI)=3000 s(-1) at pH 7.0, 10 degrees C) that contains three (helices I, II and IV) of the four native alpha-helices. The relatively slow (k(IN)=300 s(-1)) conversion of this intermediate into the native structure is driven by the folding and docking of the six residue helix III onto the developing hydrophobic core. Here, we describe the structural properties of four Im7* variants designed to trap the protein in the intermediate state by disrupting the stabilising interactions formed between helix III and the rest of the protein structure. In two of these variants (I54A and L53AI54A), hydrophobic residues within helix III have been mutated to alanine, whilst in the other two mutants the sequence encompassing the native helix III was replaced by a glycine linker, three (H3G3) or six (H3G6) residues in length. All four variants were shown to be monomeric, as judged by analytical ultracentrifugation, and highly helical as measured by far-UV CD. In addition, all the variants denature co-operatively and have a stability (DeltaG(UF)) and buried hydrophobic surface area (M(UF)) similar to those of the on-pathway kinetic intermediate. Structural characterisation of these variants using 1-anilino-8-napthalene sulphonic acid (ANS) binding, near-UV CD and 1D (1)H NMR demonstrate further that the trapped intermediate ensemble is highly structured with little exposed hydrophobic surface area. Interestingly, however, the structural properties of the variants I54A and L53AI54A differ in detail from those of H3G3 and H3G6. In particular, the single tryptophan residue, located near the end of helix IV, and distant from helix III, is in a distinct environment in the two sets of mutants as judged by fluorescence, near-UV CD and the sensitivity of tryptophan fluorescence to iodide quenching. Overall, the results confirm previous kinetic analysis that demonstrated the hierarchical folding of Im7 via an on-pathway intermediate, and show that this species is a highly helical ensemble with a well-formed hydrophobic core. By contrast with the native state, however, the intermediate ensemble is flexible enough to change in response to mutation, its structural properties being tailored by residues in the sequence encompassing the native helix III.

The "two-state folder" MerP forms partially unfolded structures that show temperature dependent hydrogen exchange

We have analysed the folding energy landscape of the 72 amino acid protein MerP by monitoring native state hydrogen exchange as a function of temperature in the range of 7-55 degrees C. The temperature dependence of the hydrogen exchange has allowed us to determine DeltaG, DeltaH and DeltaC(p) values for the conformational processes that permit hydrogen exchange. When studied with the traditional probes, fluorescence and CD, MerP appears to behave as a typical two-state protein, but the results from the hydrogen exchange analysis reveal a much more complex energy landscape. Analysis at the individual amino acid level show that exchange is allowed from an ensemble of partially unfolded structures (i.e. intermediates) in which the stabilities at the amino acid level form a broad distribution throughout the protein. The formation of partially unfolded structures might contribute to the unusually slow folding of MerP.

Folding lambda-repressor at its speed limit

We show that the five-helix bundle lambda(6-85) can be engineered and solvent-tuned to make the transition from activated two-state folding to downhill folding. The transition manifests itself as the appearance of additional dynamics faster than the activated kinetics, followed by the disappearance of the activated kinetics when the bias toward the native state is increased. Our fastest value of 1 micros for the "speed" limit of lambda(6-85) is measured at low concentrations of a denaturant that smooths the free-energy surface. Complete disappearance of the activated phase is obtained in stabilizing glucose buffer. Langevin dynamics on a rough free-energy surface with variable bias toward the native state provides a robust and quantitative description of the transition from activated to downhill folding. Based on our simulation, we estimate the residual energetic frustration of lambda(6-85) to be delta(2) G approximately 0.64 k(2)T(2). We show that lambda(6-86), as well as very fast folding proteins or folding intermediates estimated to lie near the speed limit, provide a better rate-topology correlation than proteins with larger energetic frustration. A limit of beta > or = 0.7 on any stretching of lambda(6-85) barrier-free dynamics suggests that a low-dimensional free-energy surface is sufficient to describe folding.

The regions of the sequence most exposed to the solvent within the amyloidogenic state of a protein initiate the aggregation process

Formation of misfolded aggregates is an essential part of what proteins can do. The process of protein aggregation is central to many human diseases and any aggregating event needs to be prevented within a cell and in protein design. In order to aggregate, a protein needs to unfold its native state, at least partially. The conformational state that is prone to aggregate is difficult to study, due to its aggregating potential and heterogeneous nature. Here, we use a systematic approach of limited proteolysis, in combination with electrospray ionisation mass spectrometry, to investigate the regions that are most flexible and solvent-exposed within the native, ligand-bound and amyloidogenic states of muscle acylphosphatase (AcP), a protein previously shown to form amyloid fibrils in the presence of trifluoroethanol. Seven proteases with different degrees of specificity have been used for this purpose. Following exposure to the aggregating conditions, a number of sites along the sequence of AcP become susceptible to proteolytic digestion. The pattern of proteolytic cleavages obtained under these conditions is considerably different from that of the native and ligand-bound conformations and includes a portion within the N-terminal tail of the protein (residues 6-7), the region of the sequence 18-23 and the position 94 near the C terminus. There is a significant overlap between the regions of the sequence found to be solvent-exposed from the present study and those previously identified to be critical in the rate-determining steps of aggregation from protein engineering approaches. This indicates that a considerable degree of solvent exposure is a feature of the portions of a protein that initiate the process of aggregation.

Intermediates Control Domain Swapping during Folding of p13

The 13-kDa protein p13(suc1) has two folded states, a monomer and a structurally similar domain-swapped dimer formed by exchange of a beta-strand. The refolding reaction of p13(suc1) is multiphasic, and in this paper we analyze the kinetics as a function of denaturant and protein concentration and compare the behavior of wild type and a set of mutants previously designed with dimerization propensities that span 9 orders of magnitude. We show that the folding reactions of wild type and all mutants produce the monomer predominantly despite their very different equilibrium behavior. However, the addition of low concentrations of denaturant in the refolding buffer leads to thermodynamic control of the folding reaction with products that correspond to the wild type and mutant equilibrium dimerization propensities. We present evidence that the kinetic control in the absence of urea arises because of the population of the folding intermediates. Intermediates are usually considered to be detrimental to folding because they slow down the reaction; however, our work shows that intermediates buffer the monomeric folding pathway against the effect of mutations that favor the nonfunctional, dimeric state at equilibrium.

Cystatin forms a Tetramer through Structural Rearrangement of Domain-swapped Dimers prior to Amyloidogenesis

The cystatins were the first amyloidogenic proteins to be shown to oligomerize through a 3D domain swapping mechanism. Here we show that, under conditions leading to the formation of amyloid deposits, the domain-swapped dimer of chicken cystatin further oligomerizes to a tetramer, prior to fibrillization. The tetramer has a very similar circular dichroism and fluorescence signature to the folded monomer and dimer structures, but exhibits some loss of dispersion in the 1H-NMR spectrum. 8-Anilino-1-naphthalene sulfonate fluorescence enhancement indicates an increase in the degree of disorder. While the dimerization reaction is bimolecular and most likely limited by the availability of a predominantly unfolded form of the monomer, the tetramerization reaction is first-order. The tetramer is formed slowly (t(1/2)=six days at 85 degrees C), dimeric cystatin is the precursor to tetramer formation, and thus the rate is limited by structural rearrangement within the dimer. Some higher-order oligomerization events parallel tetramer formation while others follow from the tetrameric form. Thus, the tetramer is a transient intermediate within the pathway of large-scale oligomerization.

Thermodynamics and stability of a beta-sheet complex: molecular dynamics simulations on simplified off-lattice protein models

We have performed discontinuous molecular dynamics simulations of the thermodynamics and stability of a tetrameric beta-sheet complex that contains four identical four-stranded antiparallel beta-sheet peptides. The potential used in the simulation is a hybrid Go-type potential characterized by the bias gap parameter g, an artificial measure of the preference of a model protein for its native state, and the intermolecular contact parameter eta, which measures the ratio of intermolecular to intramolecular native attractions. Despite the simplicity of the model, a complex set of thermodynamic transitions for the beta-sheet complex is revealed that shows there are three distinct oligomer (partially ordered, ordered, and highly ordered beta-sheet complex) states and four noninteracting monomers phases. The thermodynamic properties of the three oligomer states strongly depend on both the size of the intermolecular contact parameter eta and the temperature. The partially ordered beta-sheet complex is made up of four ordered globules and is observed at intermediate to large eta at high temperatures. The ordered beta-sheet complex contains four native beta-sheets and is located at small to intermediate eta at low temperatures in the phase diagram. The highly ordered beta-sheet complex has fully-stiff beta-sheet strands, the same as the global energy minimum structure, and is observed for all eta at low temperatures.

The unfolding story of three-dimensional domain swapping

Three-dimensional domain swapping is the event by which a monomer exchanges part of its structure with identical monomers to form an oligomer where each subunit has a similar structure to the monomer. The accumulating number of observations of this phenomenon in crystal structures has prompted speculation as to its biological relevance. Domain swapping was originally proposed to be a mechanism for the emergence of oligomeric proteins and as a means for functional regulation, but also to be a potentially harmful process leading to misfolding and aggregation. We highlight experimental studies carried out within the last few years that have led to a much greater understanding of the mechanism of domain swapping and of the residue- and structure-specific features that facilitate the process. We discuss the potential biological implications of domain swapping in light of these findings.

High concentrations of viscogens decrease the protein folding rate constant by prematurely collapsing the coil

In several studies, viscogenic osmolytes have been suggested to decrease the folding rate constant of polypeptides by slowing their motion through the solvent. Here, we show that osmolytes may slow protein folding by prematurely collapsing the coil. At low or moderate concentrations of osmolytes (<30%), folding of the two-state protein CI2 becomes faster with increasing osmolyte concentrations, suggesting that the kinetics are governed by protein stability. However, at higher concentrations of osmolyte, the coil collapses in the dead-time of the refolding experiment, causing a dramatic drop in the folding rate. The collapsed state is non-native and appears to be different for different osmolytes.

Protein folding in the post-genomic era

Protein folding is a topic of fundamental interest since it concerns the mechanisms by which the genetic message is translated into the three-dimensional and functional structure of proteins. In these post-genomic times, the knowledge of the fundamental principles are required in the exploitation of the information contained in the increasing number of sequenced genomes. Protein folding also has practical applications in the understanding of different pathologies and the development of novel therapeutics to prevent diseases associated with protein misfolding and aggregation. Significant advances have been made ranging from the Anfinsen postulate to the "new view" which describes the folding process in terms of an energy landscape. These new insights arise from both theoretical and experimental studies. The problem of folding in the cellular environment is briefly discussed. The modern view of misfolding and aggregation processes that are involved in several pathologies such as prion and Alzheimer diseases. Several approaches of structure prediction, which is a very active field of research, are described.

Residue-specific real-time NMR diffusion experiments define the association states of proteins during folding

Characterizing the association states of proteins during folding is critical for understanding the nature of protein-folding intermediates and protein-folding pathways, protein aggregation, and disease-related aggregation. To study the association states of unfolded, folded, and intermediate species during protein folding, we have introduced a novel residue-specific real-time NMR diffusion experiment. This experiment, a combination of NMR real-time folding experiments and 3D heteronuclear pulsed field gradient NMR diffusion experiments (LED-HSQC), measures hydrodynamic properties, or molecular sizes, of kinetic species directly during the folding process. Application of the residue-specific real-time NMR diffusion experiments to characterize the folding of the collagen triple helix motif shows that this experiment can be used to determine the association states of unfolded, folded, and kinetic intermediates with transient lifetimes simultaneously. The ratio of the apparent translational diffusion coefficients of the unfolded to the folded form of the triple helix is 0.59, which correlates very well with a theoretical ratio for monomer to linear trimer. The apparent diffusion coefficients of the kinetic intermediates formed during triple helix folding indicate the formation of trimer-like associates which is consistent with previously published kinetic and relaxation data. The residue-specific time dependence of apparent diffusion coefficients of monomer and trimer peaks also illustrates the ability to use diffusion data to probe the directionality of triple helix formation. NMR diffusion experiments provide a new strategy for the investigation of protein-folding mechanisms, both to understand the role of kinetic intermediates and to determine the time-dependent aggregation processes in human diseases.

The structure of the transition state for folding of domain-swapped dimeric p13suc1

suc1 has two native states, a monomer and a domain-swapped dimer, in which one molecule exchanges a beta strand with an identical partner. Thus, monomer and dimer have the same structures but are topologically distinct. Importantly, residues that exchange are part of the folding nucleus of the monomer and therefore forming these interactions in the dimer would be expected to incur a large entropic cost. Here we present the transition state for folding/unfolding of domain-swapped dimeric suc1 and compare it with its monomeric counterpart. The same overall structure is observed in the two transition states but the phi values are consistently higher for the domain-swapped dimer. Thus, a greater entropic penalty for bringing together the key interactions in the dimer is overcome by mobilizing more contacts in the transition state, thereby achieving a greater enthalpic gain.

Exploring protein aggregation and self-propagation using lattice models: phase diagram and kinetics

Many seemingly unrelated neurodegenerative disorders, such as amyloid and prion diseases, are associated with propagating fibrils whose structures are dramatically different from the native states of the corresponding monomers. This observation, along with the experimental demonstration that any protein can aggregate to form either fibrils or amorphous structures (inclusion bodies) under appropriate external conditions, suggest that there must be general principles that govern aggregation mechanisms. To probe generic aspects of prion-like behavior we use the model of Harrison, Chan, Prusiner, and Cohen. In this model, aggregation of a structure, that is conformationally distinct from the native state of the monomer, occurs by three parallel routes. Kinetic partitioning, which leads to parallel assembly pathways, occurs early in the aggregation process. In all pathways transient unfolding precedes polymerization and self-propagation. Chain polymerization is consistent with templated assembly, with the dimer being the minimal nucleus. The kinetic effciency of R(n-1) + G --> R(n) (R is the aggregation prone state and G is either U, the unfolded state, or N, the native state of the monomer) is increased when polymerization occurs in the presence of a "seed" (a dimer). These results support the seeded nucleated-polymerization model of fibril formation in amyloid peptides. To probe generic aspects of aggregation in two-state proteins, we use lattice models with side chains. The phase diagram in the (T,C) plane (T is the temperature and C is the polypeptide concentration) reveals a bewildering array of "phases" or structures. Explicit computations for dimers show that there are at least six phases including ordered structures and amorphous aggregates. In the ordered region of the phase diagram there are three distinct structures. We find ordered dimers (OD) in which each monomer is in the folded state and the interaction between the monomers occurs via a well-defined interface. In the domain-swapped structures a certain fraction of intrachain contacts are replaced by interchain contacts. In the parallel dimers the interface is stabilized by favorable intermolecular hydrophobic interactions. The kinetics of folding to OD shows that aggregation proceeds directly from U in a dynamically cooperative manner without populating partially structured intermediates. These results support the experimental observation that ordered aggregation in the two-state folders U1A and CI2 takes place from U. The contrasting aggregation processes in the two models suggest that there are several distinct mechanisms for polymerization that depend not only on the polypeptide sequence but also on external conditions (such as C, T, pH, and salt concentration).

Conserved and nonconserved features of the folding pathway of hisactophilin, a beta-trefoil protein

Based on previous studies of interleukin-1beta (IL-1beta) and both acidic and basic fibroblast growth factors (FGFs), it has been suggested that the folding of beta-trefoil proteins is intrinsically slow and may occur via the formation of essential intermediates. Using optical and NMR-detected quenched-flow hydrogen/deuterium exchange methods, we have measured the folding kinetics of hisactophilin, another beta-trefoil protein that has < 10% sequence identity and unrelated function to IL-1beta and FGFs. We find that hisactophilin can fold rapidly and with apparently two-state kinetics, except under the most stabilizing conditions investigated where there is evidence for formation of a folding intermediate. The hisactophilin intermediate has significant structural similarities to the IL-1beta intermediate that has been observed experimentally and predicted theoretically using a simple, topology-based folding model; however, it appears to be different from the folding intermediate observed experimentally for acidic FGF. For hisactophilin and acidic FGF, intermediates are much less prominent during folding than for IL-1beta. Considering the structures of the different beta-trefoil proteins, it appears that differences in nonconserved loops and hydrophobic interactions may play an important role in differential stabilization of the intermediates for these proteins.

Ultrafast folding of WW domains without structured aromatic clusters in the denatured state

Ultrafast-folding proteins are important for combining experiment and simulation to give complete descriptions of folding pathways. The WW domain family comprises small proteins with a three-stranded antiparallel beta-sheet topology. Previous studies on the 57-residue YAP 65 WW domain indicate the presence of residual structure in the chemically denatured state. Here we analyze three minimal core WW domains of 38-44 residues. There was little spectroscopic or thermodynamic evidence for residual structure in either their chemically or thermally denatured states. Folding and unfolding kinetics, studied by using rapid temperature-jump and continuous-flow techniques, show that each domain folds and unfolds very rapidly in a two-state transition through a highly compact transition state. Folding half-times were as short as 17 micros at 25 degrees C, within an order of magnitude of the predicted maximal rate of loop formation. The small size and topological simplicity of these domains, in conjunction with their very rapid two-state folding, may allow us to reduce the difference in time scale between experiment and theoretical simulation.

Folding of horse cytochrome c in the reduced state

Equilibrium and kinetic folding studies of horse cytochrome c in the reduced state have been carried out under strictly anaerobic conditions at neutral pH, 10 degrees C, in the entire range of aqueous solubility of guanidinium hydrochloride (GdnHCl). Equilibrium unfolding transitions observed by Soret heme absorbance, excitation energy transfer from the lone tryptophan residue to the ferrous heme, and far-UV circular dichroism (CD) are all biphasic and superimposable, implying no accumulation of structural intermediates. The thermodynamic parameters obtained by two-state analysis of these transitions yielded DeltaG(H2O)=18.8(+/-1.45) kcal mol(-1), and C(m)=5.1(+/-0.15) M GdnHCl, indicating unusual stability of reduced cytochrome c. These results have been used in conjunction with the redox potential of native cytochrome c and the known stability of oxidized cytochrome c to estimate a value of -164 mV as the redox potential of the unfolded protein. Stopped-flow kinetics of folding and unfolding have been recorded by Soret heme absorbance, and tryptophan fluorescence as observables. The refolding kinetics are monophasic in the transition region, but become biphasic as moderate to strongly native-like conditions are approached. There also is a burst folding reaction unobservable in the stopped-flow time window. Analyses of the two observable rates and their amplitudes indicate that the faster of the two rates corresponds to apparent two-state folding (U<-->N) of 80-90 % of unfolded molecules with a time constant in the range 190-550 micros estimated by linear extrapolation and model calculations. The remaining 10-20 % of the population folds to an off-pathway intermediate, I, which is required to unfold first to the initial unfolded state, U, in order to refold correctly to the native state, N (I<-->U<-->N). The slower of the two observable rates, which has a positive slope in the linear functional dependence on the denaturant concentration indicating that an unfolding process under native-like conditions indeed exists, originates from the unfolding of I to U, which rate-limits the overall folding of these 10-20 % of molecules. Both fast and slow rates are independent of protein concentration and pH of the refolding milieu, suggesting that the off-pathway intermediate is not a protein aggregate or trapped by heme misligation. The nature or type of unfolded-state heme ligation does not interfere with refolding. Equilibrium pH titration of the unfolded state yielded coupled ionization of the two non-native histidine ligands, H26 and H33, with a pK(a) value of 5.85. A substantial fraction of the unfolded population persists as the six-coordinate form even at low pH, suggesting ligation of the two methionine residues, M65 and M80. These results have been used along with the known ligand-binding properties of unfolded cytochrome c to propose a model for heme ligation dynamics. In contrast to refolding kinetics, the unfolding kinetics of reduced cytochrome c recorded by observation of Soret absorbance and tryptophan fluorescence are all slow, simple, and single-exponential. In the presence of 6.8 M GdnHCl, the unfolding time constant is approximately 300(+/-125) ms. There is no burst unfolding reaction. Simulations of the observed folding-unfolding kinetics by numerical solutions of the rate equations corresponding to the three-state I<-->U<-->N scheme have yielded the microscopic rate constants.

Bergerac-SH3: "frustation" induced by stabilizing the folding nucleus

The influence of an inserted exogenous independent folding element on the thermodynamics and folding properties of SH3 domain from alpha-spectrin has been investigated by creating a fused form between this small all-beta domain and a stable beta-hairpin (BH19). NMR analysis of synthetic peptides shows that insertion of BH19 nucleates formation of the original natural beta-hairpin (distal loop) that is part of the SH3 folding nucleus. The resulting protein (Bergerac-SHH) is more stable, folds faster and contains an elongated hairpin protruding from the globular domain as determined by 2D-NMR. "Protein engineering" analysis of the inserted region shows that it is folded in the transition state. Interestingly, stabilisation by insertion of the distal loop region results in the appearance of a compact intermediate revealed by a curved chevron plot at low denaturant concentration. This effect is eliminated at low salt concentrations by a single mutation of a hydrophobic residue within BH19 sequence, which is most probably involved in non-native interactions. Local stabilisation by enlargement and reinforcement of the folding nucleus, global compaction by the addition of salt and non-native interactions are shown to contribute to the observed deviation from the two-state behaviour.