Folding network of villin headpiece subdomain

Protein folding is a complex multidimensional process that is difficult to illustrate by the traditional analyses based on one- or two-dimensional profiles. Analyses based on transition networks have become an alternative approach that has the potential to reveal detailed features of protein folding dynamics. However, due to the lack of successful reversible folding of proteins from conventional molecular-dynamics simulations, this approach has rarely been utilized. Here, we analyzed the folding network from several 10 μs conventional molecular-dynamics reversible folding trajectories of villin headpiece subdomain (HP35). The folding network revealed more complexity than the traditional two-dimensional map and demonstrated a variety of conformations in the unfolded state, intermediate states, and the native state. Of note, deep enthalpic traps at the unfolded state were observed on the folding landscape. Furthermore, in contrast to the clear separation of the native state and the primary intermediate state shown on the two-dimensional map, the two states were mingled on the folding network, and prevalent interstate transitions were observed between these two states. A more complete picture of the folding mechanism of HP35 emerged when the traditional and network analyses were considered together.     

A Two-step Mechanism for the Folding of Actin by the Yeast Cytosolic Chaperonin

Actin requires the chaperonin containing TCP1 (CCT), a hexadecameric ATPase essential for cell viability in eukaryotes, to fold to its native state. Following binding of unfolded actin to CCT, the cavity of the chaperone closes and actin is folded and released in an ATP-dependent folding cycle. In yeast, CCT forms a ternary complex with the phosducin-like protein PLP2p to fold actin, and together they can return nascent or chemically denatured actin to its native state in a pure in vitro folding assay. The complexity of the CCT-actin system makes the study of the actin folding mechanism technically challenging. We have established a novel spectroscopic assay through selectively labeling the C terminus of yeast actin with acrylodan and observe significant changes in the acrylodan fluorescence emission spectrum as actin is chemically unfolded and then refolded by the chaperonin. The variation in the polarity of the environment surrounding the fluorescent probe during the unfolding/folding processes has allowed us to monitor actin as it folds on CCT. The rate of actin folding at a range of temperatures and ATP concentrations has been determined for both wild type CCT and a mutant CCT, CCT4anc2, defective in folding actin in vivo. Binding of the non-hydrolysable ATP analog adenosine 5′-(β,γ-imino)triphosphate to the ternary complex leads to 3-fold faster release of actin from CCT following addition of ATP, suggesting a two-step folding process with a conformational change occurring upon closure of the cavity and a subsequent final folding step involving packing of the C terminus to the native-like state.     

Protein Disulfide-Isomerase Interacts with a Substrate Protein at All Stages along Its Folding Pathway

In contrast to molecular chaperones that couple protein folding to ATP hydrolysis, protein disulfide-isomerase (PDI) catalyzes protein folding coupled to formation of disulfide bonds (oxidative folding). However, we do not know how PDI distinguishes folded, partly-folded and unfolded protein substrates. As a model intermediate in an oxidative folding pathway, we prepared a two-disulfide mutant of basic pancreatic trypsin inhibitor (BPTI) and showed by NMR that it is partly-folded and highly dynamic. NMR studies show that it binds to PDI at the same site that binds peptide ligands, with rapid binding and dissociation kinetics; surface plasmon resonance shows its interaction with PDI has a K_d of ca. 10⁻⁵ M. For comparison, we characterized the interactions of PDI with native BPTI and fully-unfolded BPTI. Interestingly, PDI does bind native BPTI, but binding is quantitatively weaker than with partly-folded and unfolded BPTI. Hence PDI recognizes and binds substrates via permanently or transiently unfolded regions. This is the first study of PDI''s interaction with a partly-folded protein, and the first to analyze this folding catalyst''s changing interactions with substrates along an oxidative folding pathway. We have identified key features that make PDI an effective catalyst of oxidative protein folding – differential affinity, rapid ligand exchange and conformational flexibility.     

Slow and Bimolecular Folding of a De Novo Designed Monomeric Protein DS119

De novo protein design offers a unique means to test and advance our understanding of how proteins fold. However, most current design methods are native structure eccentric and folding kinetics has rarely been considered in the design process. Here, we show that a de novo designed mini-protein DS119, which folds into a βαβ structure, exhibits unusually slow and concentration-dependent folding kinetics. For example, the folding time for 50 μM of DS119 was estimated to be ∼2 s. Stopped-flow fluorescence resonance energy transfer experiments further suggested that its folding was likely facilitated by a transient dimerization process. Taken together, these results highlight the need for consideration of the entire folding energy landscape in de novo protein design and provide evidence suggesting nonnative interactions can play a key role in protein folding.     

Tryptophan-tryptophan energy migration as a tool to follow apoflavodoxin folding

Submolecular details of Azotobacter vinelandii apoflavodoxin (apoFD) (un)folding are revealed by time-resolved fluorescence anisotropy using wild-type protein and variants lacking one or two of apoFD's three tryptophans. ApoFD equilibrium (un)folding by guanidine hydrochloride follows a three-state model: native ↔ unfolded ↔ intermediate. In native protein, W128 is a sink for Förster resonance energy transfer (FRET). Consequently, unidirectional FRET with a 50-ps transfer correlation time occurs from W167 to W128. FRET from W74 to W167 is much slower (6.9 ns). In the intermediate, W128 and W167 have native-like geometry because the 50-ps transfer time is observed. However, non-native structure exists between W74 and W167 because instead of 6.9 ns the transfer correlation time is 2.0 ns. In unfolded apoFD this 2.0-ns transfer correlation time is also detected. This decrease in transfer correlation time is a result of W74 and W167 becoming solvent accessible and randomly oriented toward one another. Apparently W74 and W167 are near-natively separated in the folding intermediate and in unfolded apoFD. Both tryptophans may actually be slightly closer in space than in the native state, even though apoFD's radius increases substantially upon unfolding. In unfolded apoFD the 50-ps transfer time observed for native and intermediate folding states becomes 200 ps as W128 and W167 are marginally further separated than in the native state. Apparently, apoFD's unfolded state is not a featureless statistical coil but contains well-defined substructures. The approach presented is a powerful tool to study protein folding.     

Stepwise Folding of an Autotransporter Passenger Domain Is Not Essential for Its Secretion

Autotransporters are a superfamily of virulence proteins produced by Gram-negative bacteria. They consist of an N-terminal β-helical domain (“passenger domain”) that is secreted into the extracellular space and a C-terminal β-barrel domain (“β-domain”) that anchors the protein to the outer membrane. Because the periplasm lacks ATP, vectorial folding of the passenger domain in a C-to-N-terminal direction has been proposed to drive the secretion reaction. Consistent with this hypothesis, mutations that disrupt the folding of the C terminus of the passenger domain of the Escherichia coli O157:H7 autotransporter EspP have been shown to cause strong secretion defects. Here, we show that point mutations introduced at specific locations near the middle or N terminus of the EspP β-helix that perturb folding also impair secretion, but to a lesser degree. Surprisingly, we found that even multiple mutations that potentially abolish the stability of several consecutive rungs of the β-helix only moderately reduce secretion efficiency. Although these results provide evidence that the free energy derived from passenger domain folding contributes to secretion efficiency, they also suggest that a significant fraction of the energy required for secretion is derived from another source.     

Protein Vivisection Reveals Elusive Intermediates in Folding

Although most folding intermediates escape detection, their characterization is crucial to the elucidation of folding mechanisms. Here, we outline a powerful strategy to populate partially unfolded intermediates: A buried aliphatic residue is substituted with a charged residue (e.g., Leu → Glu⁻) to destabilize and unfold a specific region of the protein. We applied this strategy to ubiquitin, reversibly trapping a folding intermediate in which the β5-strand is unfolded. The intermediate refolds to a native-like structure upon charge neutralization under mildly acidic conditions. Characterization of the trapped intermediate using NMR and hydrogen exchange methods identifies a second folding intermediate and reveals the order and free energies of the two major folding events on the native side of the rate-limiting step. This general strategy may be combined with other methods and have broad applications in the study of protein folding and other reactions that require trapping of high-energy states.     

The B domain of protein A retains residual structures in 6 M guanidinium chloride as revealed by hydrogen/deuterium‐exchange NMR spectroscopy

The characterization of residual structures persistent in unfolded proteins is an important issue in studies of protein folding, because the residual structures present, if any, may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the residual structures of the isolated B domain (BDPA) of staphylococcal protein A in 6 M guanidinium chloride. BDPA is a small three‐helix‐bundle protein, and until recently its folding/unfolding reaction has been treated as a simple two‐state process between the native and the fully unfolded states. We employed a dimethylsulfoxide (DMSO)‐quenched hydrogen/deuterium (H/D)‐exchange 2D NMR techniques with the use of spin desalting columns, which allowed us to investigate the H/D‐exchange behavior of individually identified peptide amide (NH) protons. We obtained H/D‐exchange protection factors of the 21 NH protons that form an α‐helical hydrogen bond in the native structure, and the majority of these NH protons were significantly protected with a protection factor of 2.0–5.2 in 6 M guanidinium chloride, strongly suggesting that these weakly protected NH protons form much stronger hydrogen bonds under native folding conditions. The results can be used to deduce the structure of an early folding intermediate, when such an intermediate is shown by other methods. Among three native helical regions, the third helix in the C‐terminal side was highly protected and stabilized by side‐chain salt bridges, probably acting as the folding initiation site of BDPA. The present results are discussed in relation to previous experimental and computational findings on the folding mechanisms of BDPA.     

Guanidine-HCl Dependent Structural Unfolding of M-Crystallin: Fluctuating Native State Like Topologies and Intermolecular Association

Numerous experimental techniques and computational studies, proposed in recent times, have revolutionized the understanding of protein-folding paradigm. The complete understanding of protein folding and intermediates are of medical relevance, as the aggregation of misfolding proteins underlies various diseases, including some neurodegenerative disorders. Here, we describe the unfolding of M-crystallin, a βγ-crystallin homologue protein from archaea, from its native state to its denatured state using multidimensional NMR and other biophysical techniques. The protein, which was earlier characterized to be a predominantly β-sheet protein in its native state, shows different structural propensities (α and β), under different denaturing conditions. In 2 M GdmCl, the protein starts showing two distinct sets of peaks, with one arising from a partially unfolded state and the other from a completely folded state. The native secondary structural elements start disappearing as the denaturant concentration approaches 4 M. Subsequently, the protein is completely unfolded when the denaturant concentration is 6 M. The ¹⁵N relaxation data (T₁/T₂), heteronuclear ¹H-¹⁵N Overhauser effects (nOes), NOESY data, and other biophysical data taken together indicate that the protein shows a consistent, gradual change in its structural and motional preferences with increasing GdmCl concentration.     

Single-Molecule Studies of the Im7 Folding Landscape

Under appropriate conditions, the four-helical Im7 (immunity protein 7) folds from an ensemble of unfolded conformers to a highly compact native state via an on-pathway intermediate. Here, we investigate the unfolded, intermediate, and native states populated during folding using diffusion single-pair fluorescence resonance energy transfer by measuring the efficiency of energy transfer (or proximity or P ratio) between pairs of fluorophores introduced into the side chains of cysteine residues placed in the center of helices 1 and 4, 1 and 3, or 2 and 4. We show that while the native states of each variant give rise to a single narrow distribution with high P values, the distributions of the intermediates trapped at equilibrium (denoted I^eqm) are fitted by two Gaussian distributions. Modulation of the folding conditions from those that stabilize the intermediate to those that destabilize the intermediate enabled the distribution of lower P value to be assigned to the population of the unfolded ensemble in equilibrium with the intermediate state. The reduced stability of the I^eqm variants allowed analysis of the effect of denaturant concentration on the compaction and breadth of the unfolded state ensemble to be quantified from 0 to 6 M urea. Significant compaction is observed as the concentration of urea is decreased in both the presence and absence of sodium sulfate, as previously reported for a variety of proteins. In the presence of Na₂SO₄ in 0 M urea, the P value of the unfolded state ensemble approaches that of the native state. Concurrent with compaction, the ensemble displays increased peak width of P values, possibly reflecting a reduction in the rate of conformational exchange among iso-energetic unfolded, but compact conformations. The results provide new insights into the initial stages of folding of Im7 and suggest that the unfolded state is highly conformationally constrained at the outset of folding.     

Assembly and Kinetic Folding Pathways of a Tetrameric β-Sheet Complex: Molecular Dynamics Simulations on Simplified Off-Lattice Protein Models

We have performed discontinuous molecular dynamic simulations of the assembly and folding kinetics of a tetrameric β-sheet complex that contains four identical four-stranded antiparallel β-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 a model protein's preference for its native state, and the intermolecular contact parameter η, which measures the ratio of intermolecular to intramolecular native attractions. The formation of the β-sheet complex and its equilibrium properties strongly depend on the size of the intermolecular contact parameter η. The ordered β-sheet complex in the folded state and nonaligned β-sheets or tangled chains in the misfolded state are distinguished by measuring the squared radius of gyration and the fraction of native contacts Q. The folding yield for the folded state is high at intermediate values of η, but is low at both small and large values of η. The folded state at small η is liquid-like, but is solid-like at both intermediate and large η. The misfolded state at small η contains nonaligned β-sheets and tangled chains with poor secondary structure at large η. Various folding pathways via dimeric and trimeric intermediates are observed, depending on η. Comparison with experimental results on protein aggregation indicates that intermediate η values are most appropriate for modeling fibril formation and small η values are most appropriate for modeling the formation of amorphous aggregates.     

Cunning Simplicity of a Hierarchical Folding

Abstract

A hierarchic scheme of protein folding does not solve the Levinthal paradox since it cannot provide a simultaneous explanation for major features observed for protein folding: (i) folding within non-astronomical time, (ii) independence of the native structure on large variations in the folding rates of given protein under different conditions, and (iii) co-existence, in a visible quantity, of only the native and the unfolded molecules during folding of moderate size (single-domain) proteins. On the contrary, a nucleation mechanism can account for all these major features simultaneously and resolves the Levinthal paradox.     

Role of Protein Stabilizers on the Conformation of the Unfolded State of Cytochrome c and Its Early Folding Kinetics: INVESTIGATION AT SINGLE MOLECULAR RESOLUTION*

An insight into the conformation and dynamics of unfolded and early intermediate states of a protein is essential to understand the mechanism of its aggregation and to design potent inhibitor molecules. Fluorescence correlation spectroscopy has been used to study the effects of several model protein stabilizers on the conformation of the unfolded state and early folding dynamics of tetramethyl rhodamine-labeled cytochrome c from Saccharomyces cerevisiae at single molecular resolution. Special attention has been given to arginine, which is a widely used stabilizer for improving refolding yield of different proteins. The value of the hydrodynamic radius (r_H) obtained by analyzing the intensity fluctuations of the diffusing molecules has been found to increase in a two-state manner as the protein is unfolded by urea. The results further show that the presence of arginine and other protein stabilizers favors a relatively structured conformation of the unfolded states (r_H of 29 Å) over an extended one (r_H of 40 Å), which forms in their absence. Also, the time constant of a kinetic component (τ_R) of about 30 μs has been observed by analyzing the correlation functions, which represents formation of a collapsed state. This time constant varies with urea concentration representing an inverted Chevron plot that shows a roll-over and behavior in the absence of arginine. To the best of our knowledge, this is one of the first applications of fluorescence correlation spectroscopy to study direct folding kinetics of a protein.     

Visualization of Protein Folding Funnels in Lattice Models

Protein folding occurs in a very high dimensional phase space with an exponentially large number of states, and according to the energy landscape theory it exhibits a topology resembling a funnel. In this statistical approach, the folding mechanism is unveiled by describing the local minima in an effective one-dimensional representation. Other approaches based on potential energy landscapes address the hierarchical structure of local energy minima through disconnectivity graphs. In this paper, we introduce a metric to describe the distance between any two conformations, which also allows us to go beyond the one-dimensional representation and visualize the folding funnel in 2D and 3D. In this way it is possible to assess the folding process in detail, e.g., by identifying the connectivity between conformations and establishing the paths to reach the native state, in addition to regions where trapping may occur. Unlike the disconnectivity maps method, which is based on the kinetic connections between states, our methodology is based on structural similarities inferred from the new metric. The method was developed in a 27-mer protein lattice model, folded into a 3×3×3 cube. Five sequences were studied and distinct funnels were generated in an analysis restricted to conformations from the transition-state to the native configuration. Consistent with the expected results from the energy landscape theory, folding routes can be visualized to probe different regions of the phase space, as well as determine the difficulty in folding of the distinct sequences. Changes in the landscape due to mutations were visualized, with the comparison between wild and mutated local minima in a single map, which serves to identify different trapping regions. The extension of this approach to more realistic models and its use in combination with other approaches are discussed.     

Folding of Active β-Lactamase in the Yeast Cytoplasm before Translocation into the Endoplasmic Reticulum

Polypeptides targeted to the yeast endoplasmic reticulum (ER) posttranslationally are thought to be kept in the cytoplasm in an unfolded state by Hsp70 chaperones before translocation. We show here that Escherichia coli β-lactamase associated with Hsp70, but adopted a native-like conformation before translocation in living Saccharomyces cerevisiae cells. β-Lactamase is a globular trypsin-resistant molecule in authentic form. For these studies, it was linked to the C terminus of a yeast polypeptide Hsp150Δ, which conferred posttranslational translocation and provided sites for O-glycosylation. We devised conditions to retard translocation of Hsp150Δ-β-lactamase. This enabled us to show by protease protection assays that an unglycosylated precursor was associated with the cytoplasmic surface of isolated microsomes, whereas a glycosylated form resided inside the vesicles. Both proteins were trypsin resistant and had similar β-lactamase activity and K_m values for nitrocefin. The enzymatically active cytoplasmic intermediate could be chased into the ER, followed by secretion of the activity to the medium. Productive folding in the cytoplasm occurred in the absence of disulfide formation, whereas in the ER lumen, proper folding required oxidation of the sulfhydryls. This suggests that the polypeptide was refolded in the ER and consequently, at least partially unfolded for translocation.     

Tubulin Subunits Exist in an Activated Conformational State Generated and Maintained by Protein Cofactors

The production of native α/β tubulin heterodimer in vitro depends on the action of cytosolic chaperonin and several protein cofactors. We previously showed that four such cofactors (termed A, C, D, and E) together with native tubulin act on β-tubulin folding intermediates generated by the chaperonin to produce polymerizable tubulin heterodimers. However, this set of cofactors generates native heterodimers only very inefficiently from α-tubulin folding intermediates produced by the same chaperonin. Here we describe the isolation, characterization, and genetic analysis of a novel tubulin folding cofactor (cofactor B) that greatly enhances the efficiency of α-tubulin folding in vitro. This enabled an integrated study of α- and β-tubulin folding: we find that the pathways leading to the formation of native α- and β-tubulin converge in that the folding of the α subunit requires the participation of cofactor complexes containing the β subunit and vice versa. We also show that sequestration of native α-or β-tubulins by complex formation with cofactors results in the destabilization and decay of the remaining free subunit. These data demonstrate that tubulin folding cofactors function by placing and/or maintaining α-and β-tubulin polypeptides in an activated conformational state required for the formation of native α/β heterodimers, and imply that each subunit provides information necessary for the proper folding of the other.     

N-Alpha-Acetylation of α-Synuclein Increases Its Helical Folding Propensity,GM1 Binding Specificity and Resistance to Aggregation

A switch in the conformational properties of α-synuclein (αS) is hypothesized to be a key step in the pathogenic mechanism of Parkinson’s disease (PD). Whereas the beta-sheet-rich state of αS has long been associated with its pathological aggregation in PD, a partially alpha-helical state was found to be related to physiological lipid binding; this suggests a potential role of the alpha-helical state in controlling synaptic vesicle cycling and resistance to β-sheet rich aggregation. N-terminal acetylation is the predominant post-translational modification of mammalian αS. Using circular dichroism, isothermal titration calorimetry, and fluorescence spectroscopy, we have analyzed the effects of N-terminal acetylation on the propensity of recombinant human αS to form the two conformational states in interaction with lipid membranes. Small unilamellar vesicles of negatively charged lipids served as model membranes. Consistent with previous NMR studies using phosphatidylserine, we found that membrane-induced α-helical folding was enhanced by N-terminal acetylation and that greater exothermic heat could be measured upon vesicle binding of the modified protein. Interestingly, the folding and lipid binding enhancements with phosphatidylserine in vitro were weak when compared to that of αS with GM1, a lipid enriched in presynaptic membranes. The resultant increase in helical folding propensity of N-acetylated αS enhanced its resistance to aggregation. Our findings demonstrate the significance of the extreme N-terminus for folding nucleation, for relative GM1 specificity of αS-membrane interaction, and for a protective function of N-terminal-acetylation against αS aggregation mediated by GM1.