High populations of non-native structures in the denatured state are compatible with the formation of the native folded state

Implications of the unfolded state in the folding energetics of heterogeneous-backbone protein mimetics

Sequence-encoded folding is the foundation of protein structure and is also possible in synthetic chains of artificial chemical composition. In natural proteins, the characteristics of the unfolded state are as important as those of the folded state in determining folding energetics. While much is known about folded structures adopted by artificial protein-like chains, corresponding information about the unfolded states of these molecules is lacking. Here, we report the consequences of altered backbone composition on the structure, stability, and dynamics of the folded and unfolded states of a compact helix-rich protein. Characterization through a combination of biophysical experiments and atomistic simulation reveals effects of backbone modification that depend on both the type of artificial monomers employed and where they are applied in sequence. In general, introducing artificial connectivity in a way that reinforces characteristics of the unfolded state ensemble of the prototype natural protein minimizes the impact of chemical changes on folded stability. These findings have implications in the design of protein mimetics and provide an atomically detailed picture of the unfolded state of a natural protein and artificial analogues under non-denaturing conditions.

Biophysical experiments and atomistic simulation reveal impacts of protein backbone alteration on the ensemble that defines the unfolded state. These effects have implications on folded stability of protein mimetics.

Thermodynamic characterization for the denatured state of bovine prion protein and the BSE Associated variant E211K

Propagation of transmissible spongiform encephalopathies involves the conversion of cellular prion protein, PrP^C, into a misfolded oligomeric form, PrP^Sc. The most common hereditary prion disease is a genetic form of Creutzfeldt-Jakob disease in humans, in which a mutation in the prion gene results in a glutamic acid to lysine substitution at position 200 (E200K) in PrP. In cattle, the analogous amino acid substitution is found at residue 211 (E211K) and has been associated with a case of bovine spongiform encephalopathy. Here, we have compared the secondary structure of E211K to that of wild type using circular dichroism and completed a thermodynamic analysis of the folding of recombinant wild type and E211K variants of the bovine prion protein. The secondary structure of the E211K variant was essentially indistinguishable from that of wild type. The thermodynamic stability of E211K substitution showed a slight destabilization relative to the wild type consistent with results reported for recombinant human prion protein and its mutant E200K. In addition, the E211K variant exhibits a similarly compact denatured state to that of wild type based upon similar m-value and change in heat capacity of unfolding for the proteins. Together these results indicate that residual structure in the denatured state of bPrP is present in both the wild type protein and BSE associated variant E211K. Given this observation, as well as folding similarities reported for other disease associated variants of PrP it is worth consideration that functional aspects of PrP conformation may play a role in the misfolding process.

Effect of an Imposed Contact on Secondary Structure in the Denatured State of Yeast Iso-1-cytochrome c

There is considerable evidence that long-range interactions stabilize residual protein structure under denaturing conditions. However, evaluation of the effect of a specific contact on structure in the denatured state has been difficult. Iso-1-cytochrome c variants with a Lys54 → His mutation form a particularly stable His-heme loop in the denatured state, suggestive of loop-induced residual structure. We have used multidimensional nuclear magnetic resonance methods to assign ¹H and ¹⁵N backbone amide and ¹³C backbone and side chain chemical shifts in the denatured state of iso-1-cytochrome c carrying the Lys54 → His mutation in 3 and 6 M guanidine hydrochloride and at both pH 6.4, where the His54-heme loop is formed, and pH 3.6, where the His54-heme loop is broken. Using the secondary structure propensity score, with the 6 M guanidine hydrochloride chemical shift data as a random coil reference state for data collected in 3 M guanidine hydrochloride, we found residual helical structure in the denatured state for the 60s helix and the C-terminal helix, but not in the N-terminal helix in the presence or absence of the His54-heme loop. Non-native helical structure is observed in two regions that form Ω-loops in the native state. There is more residual helical structure in the C-terminal helix at pH 6.4 when the loop is formed. Loop formation also appears to stabilize helical structure near His54, consistent with induction of helical structure observed when His-heme bonds form in heme-peptide model systems. The results are discussed in the context of the folding mechanism of cytochrome c.

Insights into Unfolded Proteins from the Intrinsic ϕ/ψ Propensities of the AAXAA Host-Guest Series

Various host-guest peptide series are used by experimentalists as reference conformational states. One such use is as a baseline for random-coil NMR chemical shifts. Comparison to this random-coil baseline, through secondary chemical shifts, is used to infer protein secondary structure. The use of these random-coil data sets rests on the perception that the reference chemical shifts arise from states where there is little or no conformational bias. However, there is growing evidence that the conformational composition of natively and nonnatively unfolded proteins fail to approach anything that can be construed as random coil. Here, we use molecular dynamics simulations of an alanine-based host-guest peptide series (AAXAA) as a model of unfolded and denatured states to examine the intrinsic propensities of the amino acids. We produced ensembles that are in good agreement with the experimental NMR chemical shifts and confirm that the sampling of the 20 natural amino acids in this peptide series is be far from random. Preferences toward certain regions of conformational space were both present and dependent upon the environment when compared under conditions typically used to denature proteins, i.e., thermal and chemical denaturation. Moreover, the simulations allowed us to examine the conformational makeup of the underlying ensembles giving rise to the ensemble-averaged chemical shifts. We present these data as an intrinsic backbone propensity library that forms part of our Structural Library of Intrinsic Residue Propensities to inform model building, to aid in interpretation of experiment, and for structure prediction of natively and nonnatively unfolded states.

Mapping the Native Conformational Ensemble of Proteins from a Combination of Simulations and Experiments: New Insight into the src-SH3 Domain

The biological function of a protein is strongly tied to the ensemble of three-dimensional conformations populated at physiological temperature, and dynamically transforming into each other. Experimental techniques such as nuclear magnetic resonance spectroscopy (NMR) provide a wealth of structural and dynamical information, which, in combination with an accurate atomic-level computational modeling, can disclose the details of protein behavior. We here propose a fast and efficient protocol employing molecular dynamics (MD) simulations and NMR chemical shifts, which allows one to reconstruct the detailed conformational ensemble of small globular proteins. In the case of the well-studied src-SH3 domain, we are able to obtain new important insight including the existence of a helical state in the RT loop and a pathway for single-file water diffusion in and out of the core.

Improvement of structure-based potentials for protein folding by native and nonnative hydrogen bonds

Pure Gō models (where every native interaction equally stabilizes the folded state) have widely proved their convenience in the computational investigation of protein folding. However, a chemistry-based description of the real interactions also provides a desirable tune in the analysis of the folding process, and thus some hybrid Gō potentials that combine both aspects have been proposed. Among all the noncovalent interactions that contribute to protein folding, hydrogen bonds are the only ones with a partial covalent character. This feature makes them directional and, thus, more difficult to model as part of the coarse-grained descriptions that are typically employed in Gō models. Thanks to a simplified but rigorous representation of backbone hydrogen bonds that we have recently proposed, we present in this article a combined potential (Gō + backbone hydrogen bond) to study the thermodynamics of protein folding in the frame of very simple simulation models. We show that the explicit inclusion of hydrogen bonds leads to a systematic improvement in the description of protein folding. We discuss a representative set of examples (from two-state folders to downhill proteins, with different types of native structures) that reveal a relevant agreement with experimental data.

Cooperative formation of native-like tertiary contacts in the ensemble of unfolded states of a four-helix protein

In studies of the ensembles of unfolded structures of a four-helix bundle protein, we have detected the presence of potential precursors of native tertiary structures. These observations were based on the perturbation of NMR chemical shifts of the protein backbone atoms by single site mutations. Some mutations change the chemical shifts of residues remote from the site of mutation indicating the presence of an interaction between the mutated and the remote residues, suggesting that the formation of helix segments and helix-helix interactions is cooperative. We can begin to track down the folding mechanism of this protein using only experimental data by combining the information available for the rate limiting structure formation during the folding process with measurements of the site specific hydrogen bond formation in the burst phase, and with the existence prior to the folding reaction of tertiary structures in the ensemble of otherwise unfolded structures observed in the present study.

Residue-specific description of non-native transient structures in the ensemble of acid-denatured structures of the all-beta protein c-src SH3

Secondary chemical shift analysis has been used to characterize the unfolded state of acid-denatured c-src SH3. Even though native c-src SH3 adopts an all-beta fold, we found evidence of transient helicity in regions corresponding to native loops. In particular, residues 40-46, connecting the n-src loop to the third beta-strand, exhibited an apparent helicity of nearly 45%. Furthermore, the RT loop and the diverging turn appeared to adopt non-native-like helical conformations. Interestingly, none of the residues found in transient helical conformations exhibited significant varphi-values [Riddle, D. S., et al. (1999) Nat. Struct. Biol. 6, 1016-1024]. This indicated that the transient helicity has no influence or only a weak influence on the actual protein folding reaction. The residual structural propensities were compared to those of other SH3 domains, revealing heterogeneity in the unfolded ensemble that clearly contrasts with the conserved character of the topology of native state and transition state ensembles typical for SH3 domains.

Characterization of intrinsically disordered proteins with electrospray ionization mass spectrometry: conformational heterogeneity of alpha-synuclein

Conformational heterogeneity of alpha-synuclein was studied with electrospray ionization mass spectrometry by analyzing protein ion charge state distributions, where the extent of multiple charging reflects compactness of the protein conformations in solution. Although alpha-synuclein lacks a single well-defined structure under physiological conditions, it was found to sample four distinct conformational states, ranging from a highly structured one to a random coil. The compact highly structured state of alpha-synuclein is present across the entire range of conditions tested (pH ranging from 2.5 to 10, alcohol content from 0% to 60%), but is particularly abundant in acidic solutions. The only other protein state populated in acidic solutions is a partially folded intermediate state lacking stable tertiary structure. Another, more compact intermediate state is induced by significant amounts of ethanol used as a co-solvent and appears to represent a partially folded conformation with high beta-sheet content. Protein dimerization is observed throughout the entire range of conditions tested, although only acidic solutions favor formation of highly structured dimers of alpha-synuclein. These dimers are likely to present the earliest stages in protein aggregation leading to globular oligomers and, subsequently, protofibrils.

Identification of the folding inhibitors of hen-egg lysozyme: gathering the right tools

The unfolded state of proteins displays a surprisingly rich amount of local native structure, which appears to be critical for driving the protein to its native state. Peptides with the same sequence of the corresponding structured segments can be used to interfere with the correct folding of the protein. Using model simulations, we investigate the folding of hen-egg lysozyme, identifying its key segments. Activity assays, NMR and circular dichroism experiments are used to screen the peptides which are able to inhibit the folding of lysozyme. Few peptides, corresponding to the segments of the protein which are structured in the unfolded state, are identified to have significant inhibitory effects.

Fluorescence characterization of denatured proteins

Characterization of unfolded states, while critical to a complete understanding of protein folding, is inherently difficult due to structural heterogeneity and dynamic interchange between states. The growing body of work focusing on single molecule fluorescence techniques for the study of protein folding, also highlights their potential for studies of unfolded proteins. These methods can obtain conformational information about individual subpopulations of molecules in an ensemble, and measure dynamics without the need for synchronization. The studies highlighted here demonstrate the promise of these techniques for obtaining novel information about unfolded states in vitro and in more physiologically relevant milieu.

NMR evidence for forming highly populated helical conformations in the partially folded hNck2 SH3 domain

Recent studies of several proteins implied that the folding of beta-proteins may follow a nonhierarchical mechanism in which two major transitions are essential, i.e., the collapse of a random coil to form a nonnative helical intermediate, followed by a transformation into the native beta-structure. We report that the first hNck2 SH3 domain, assuming an all-beta barrel in the native form, can be reversibly transformed into a stable and nonnative helical state by acid-unfolding. We also conducted extensive NMR and mutagenesis studies that led to two striking findings: 1), NMR analysis reveals that in the helical state formed at pH 2.0, the first and last beta-strands in the native form become unstructured, whereas the rest is surprisingly converted into two highly populated helices with a significantly limited backbone motion; and 2), a conserved four-residue sequence is identified on the second beta-strand, a mutation of which suddenly renders the SH3 domain into a helical state even at pH 6.5, with NMR conformational and dynamic properties highly similar to those of the wild-type at pH 2.0. This observation implies that the region might contribute key interactions to disrupt the helical state, and to facilitate a further transformation into the native SH3 fold in the second transition.

Characterization of protein-folding pathways by reduced-space modeling

Ab initio simulations of the folding pathways are currently limited to very small proteins. For larger proteins, some approximations or simplifications in protein models need to be introduced. Protein folding and unfolding are among the basic processes in the cell and are very difficult to characterize in detail by experiment or simulation. Chymotrypsin inhibitor 2 (CI2) and barnase are probably the best characterized experimentally in this respect. For these model systems, initial folding stages were simulated by using CA-CB-side chain (CABS), a reduced-space protein-modeling tool. CABS employs knowledge-based potentials that proved to be very successful in protein structure prediction. With the use of isothermal Monte Carlo (MC) dynamics, initiation sites with a residual structure and weak tertiary interactions were identified. Such structures are essential for the initiation of the folding process through a sequential reduction of the protein conformational space, overcoming the Levinthal paradox in this manner. Furthermore, nucleation sites that initiate a tertiary interactions network were located. The MC simulations correspond perfectly to the results of experimental and theoretical research and bring insights into CI2 folding mechanism: unambiguous sequence of folding events was reported as well as cooperative substructures compatible with those obtained in recent molecular dynamics unfolding studies. The correspondence between the simulation and experiment shows that knowledge-based potentials are not only useful in protein structure predictions but are also capable of reproducing the folding pathways. Thus, the results of this work significantly extend the applicability range of reduced models in the theoretical study of proteins.

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.

Statistical properties and kinetics of end-end contact formation of unfolded polypeptides: a systematic molecular dynamics study

The authors have systematically examined the statistical properties of the unfolded states of series of polypeptides and the kinetics of their end-to-end contact (ring closure) formation by molecular dynamics simulations. The formation of an end-to-end contact follows a single-exponential decay as measured by the first-passage time. It is shown that the shifted Gaussian chain model can be applied to describe the dimensions of glycine-rich polypeptides at high temperature. However, notable deviation from the ideal Gaussian chain model was observed at lower temperatures particularly for those polypeptides without glycines, due to the tendency to form local structures.

Spectrin organization and dynamics: new insights

Spectrin is the major constituent protein of the erythrocyte cytoskeleton which forms a filamentous network on the cytoplasmic face of the membrane by providing a scaffold for a variety of proteins. In this review, several aspects of spectrin organization are highlighted, particularly with respect to its ability to bind hydrophobic ligands and its interaction with membrane surfaces. The characteristic binding of the fluorescent hydrophobic probes Prodan and pyrene to spectrin, which allows an estimation of the polarity of the hydrophobic probe binding site, is illustrated. In addition, the contribution of uniquely localized and conserved tryptophan residues in the 'spectrin repeats' in these processes is discussed. A functional implication of the presence of hydrophobic binding sites in spectrin is its recently discovered chaperone-like activity. Interestingly, spectrin exhibits residual structural integrity even after denaturation which could be considered as a hallmark of cytoskeletal proteins. Future research could provide useful information about the possible role played by spectrin in cellular physiology in healthy and diseased states.

Characterization of the unfolded state of bovine alpha-lactalbumin and comparison with unfolded states of homologous proteins

The unfolded states of three homologous proteins with a very similar fold have been investigated by heteronuclear NMR spectroscopy. Secondary structure propensities as derived from interpretation of chemical shifts and motional restrictions as evidenced by heteronuclear (15)N relaxation rates have been analyzed in the reduced unfolded states of hen lysozyme and the calcium-binding proteins bovine alpha-lactalbumin and human alpha-lactalbumin. For all three proteins, significant deviations from random-coil predictions can be identified; in addition, the unfolded states also differ from each other, despite the fact that they possess very similar structures in their native states. Deviations from random-coil motional properties are observed in the alpha- and the beta-domain in bovine alpha-lactalbumin and lysozyme, while only regions within the alpha-domain deviate in human alpha-lactalbumin. The motional restrictions and residual secondary structure are determined both by the amino acid sequence of the protein and by residual long-range interactions. Even a conservative single point mutation from I to L in a highly conserved region between the two alpha-lactalbumins results in considerable differences in the motional properties. Given the differences in oxidative folding between hen lysozyme and alpha-lactalbumin, the results obtained on the unfolded states suggest that residual long-range interactions, i.e., those between the alpha- and the beta-domain of lysozyme, may act as nucleation sites for protein folding, while this property of residual structure is replaced by the calcium-binding site between the domains in alpha-lactalbumin.

NMR analysis of Caenorhabditis elegans FLP-18 neuropeptides: implications for NPR-1 activation

Phe-Met-Arg-Phe-NH2 (FMRFamide)-like peptides (FLPs) are the largest neuropeptide family in animals, particularly invertebrates. FLPs are characterized by a C-N-terminal gradient of decreasing amino acid conservation. Neuropeptide receptor 1 (NPR-1) is a G-protein coupled receptor (GPCR), which has been shown to be a strong regulator of foraging behavior and aggregation responses in Caenorhabditis elegans. Recently, ligands for NPR-1 were identified as neuropeptides coded by the precursor genes flp-18 and flp-21 in C. elegans. The flp-18 gene encodes eight FLPs including DFDGAMPGVLRF-NH2 and EMPGVLRF-NH2. These peptides exhibit considerably different activities on NPR-1, with the longer one showing a lower potency. We have used nuclear magnetic resonance and biological activity to investigate structural features that may explain these activity differences. Our data demonstrate that long-range electrostatic interactions exist between N-terminal aspartates and the C-terminal penultimate arginine as well as N-terminal hydrogen-bonding interactions that form transient loops within DFDGAMPGVLRF-NH2. We hypothesize that these loops, along with peptide charge, diminish the activity of this peptide on NPR-1 relative to that of EMPGVLRF-NH2. These results provide some insight into the large amino acid diversity in FLPs.

Characterization of the hydrodynamic properties of the folding transition state of an SH3 domain by magnetization transfer NMR spectroscopy

Protein folding kinetic data have been obtained for the marginally stable N-terminal Src homology 3 domain of the Drosophila protein drk (drkN SH3) in an investigation of the hydrodynamic properties of its folding transition state. Due to the presence of NMR resonances of both folded and unfolded states at equilibrium, kinetic data can be derived from NMR magnetization transfer techniques under equilibrium conditions. Kinetic analysis as a function of urea (less than approximately 1 M) and glycerol enables determination of alpha values, measures of the energetic sensitivity of the transition state to the perturbation relative to the end states of the protein folding reaction (the folded and unfolded states). Both end states have previously been studied experimentally by NMR spectroscopic and other biophysical methods in great detail and under nondenaturing conditions. Combining these results with the kinetic folding data obtained here, we can characterize the folding transition state without requiring empirical models for the unfolded state structure. We are thus able to give a reliable measure of the solvent-accessible surface area of the transition state of the drkN SH3 domain (4730 +/- 360 A(2)) based on urea titration data. Glycerol titration data give similar results and additionally demonstrate that folding of this SH3 domain is dependent on solvent viscosity, which is indicative of at least partial hydration of the transition state. Because SH3 domains appear to fold by a common folding mechanism, the data presented here provide valuable insight into the transition states of the drkN and other SH3 domains.

Motional properties of unfolded ubiquitin: a model for a random coil protein

The characterization of unfolded states of proteins has recently attracted considerable interest, as the residual structure present in these states may play a crucial role in determining their folding and misfolding behavior. Here, we investigated the dynamics in the denatured state of ubiquitin in 8 M urea at pH2. Under these conditions, ubiquitin does not have any detectable local residual structure, and uniform 15N relaxation rates along the sequence indicate the absence of motional restrictions caused by residual secondary structure and/or long-range interactions. A comparison of different models to predict relaxation data in unfolded proteins suggests that the subnanosecond dynamics in unfolded states depend on segmental motions only and do not show a dependence on the residue type but for proline and glycine residues.

Extending the folding nucleus of ubiquitin with an independently folding beta-hairpin finger: hurdles to rapid folding arising from the stabilisation of local interactions

The N-terminal beta-hairpin sequence of ubiquitin has been implicated as a folding nucleation site. To extend and stabilise the ubiquitin folding nucleus, we have inserted an autonomously folding 14-residue peptide sequence beta4 which in isolation forms a highly populated beta-hairpin (>70%) stabilised by local interactions. NMR structural analysis of the ubiquitin mutant (Ubeta4) shows that the hairpin finger is fully structured and stabilises ubiquitin by approximately 8kJmol(-1). Protein engineering and kinetic (phi(F)-value) analysis of a series of Ubeta4 mutants shows that the hairpin extension of Ubeta4 is also significantly populated in the transition state (phi(F)-values >0.7) and has the effect of templating the formation of native contacts in the folding nucleus of ubiquitin. However, at low denaturant concentrations the chevron plot of Ubeta4 shows a small deviation from linearity (roll-over effect), indicative of the population of a compact collapsed state, which appears to arise from over-stabilisation of local interactions. Destabilising mutations within the native hairpin sequence and within the engineered hairpin extension, but not elsewhere, eliminate this non-linearity and restore apparent two-state behaviour. The pitfall to stabilising local interactions is to present hurdles to the rapid and efficient folding of small proteins down a smooth folding funnel by trapping partially folded or misfolded states that must unfold or rearrange before refolding.

Determination of an ensemble of structures representing the denatured state of the bovine acyl-coenzyme a binding protein

The denatured state of a protein contains important information about the determinants of the folding process. By combining site-directed spin-labeling NMR experiments and restrained computer simulations, we have determined ensembles of conformations that represent the denatured state of the bovine acyl-coenzyme A binding protein (ACBP) at three different concentrations of guanidine hydrochloride. As the experimentally determined distance information corresponds to weighted averages over a broad ensemble of structures, we applied the experimental restraints to a system of noninteracting replicas of the protein by using a Monte Carlo sampling scheme. This procedure permits us to sample ensembles of conformations that are compatible with the experimental data and thus to obtain information regarding the distribution of structures in the denatured state. Our results show that the denatured state of ACBP is highly heterogeneous. The high sensitivity of the computational method that we present, however, enabled us to identify long-range interactions between two regions, located near the N- and C-termini, that include both native and non-native elements. The preferential formation of these contacts suggests that the sequence-dependent patterns of helical propensity and hydrophobicity are important determinants of the structure in the denatured state of ACBP.

Direct Demonstration of Structural Similarity between Native and Denatured Eglin C

To characterize the long-range structure that persists in the unfolded form of the 70-residue protein eglin C, residual dipolar couplings (RDCs) for HN-N and HA-CA bond vectors were measured by NMR spectroscopy for both its low pH, urea denatured state and its native state. When the data sets for the two different structural states were compared, a statistically significant correlation was found, with both sets of dipolar couplings yielding a correlation coefficient of r = 0.47 to 0.51. This finding directly demonstrates that the denatured state of eglin C has a nativelike global structure, a conclusion reached indirectly for staphylococcal nuclease by combining two different types of NMR data. A simple computer simulation showed that the degree of variation in phi and psi angles that yields the RDC correlation of r = 0.5 was inversely dependent on the statistical segment length, ranging from +/-6 to +/-30 degrees at the upper limit. Stable nativelike topologies that persist on unfolding would explain the rapid refolding kinetics displayed by many proteins and might provide a natural barrier against amyloid fibril formation.

The kinetics of side chain stabilization during protein folding

The rate of stabilization of side chains during protein folding has never been carefully studied. Recent developments in labeling proteins with (19)F-labeled amino acids coupled with real-time NMR measurements have allowed such measurements to be made. This paper describes the application of this method to the study of several proteins using 6-(19)F-tryptophan as the reporting group. It is found that these side chains adopt their final stable state at the last stages of the folding process and that the stabilization of side chains into their final conformation is a highly cooperative process. It is also possible to show the presence of intermediates in which the side chains are not correctly packed. The technique should be applicable to many systems.

Folding a protein in a computer: an atomic description of the folding/unfolding of protein A

We study the folding mechanism of a three-helix bundle protein at atomic resolution, including effects of explicit water. Using replica exchange molecular dynamics we perform enough sampling over a wide range of temperatures to obtain the free energy, entropy, and enthalpy surfaces as a function of structural reaction coordinates. Simulations were started from different configurations covering the folded and unfolded states. Because many transitions between all minima at the free energy surface are observed, a quantitative determination of the free energy barriers and the ensemble of configurations associated with them is now possible. The kinetic bottlenecks for folding can be determined from the thermal ensembles of structures on the free energy barriers, provided the kinetically determined transition-state ensembles are similar to those determined from free energy barriers. A mechanism incorporating the interplay among backbone ordering, sidechain packing, and desolvation arises from these calculations. Large Phi values arise not only from native contacts, which mostly form at the transition state, but also from contacts already present in the unfolded state that are partially destroyed at the transition.

The present view of the mechanism of protein folding

We can track the positions and movements of all the atoms in small proteins as they fold and unfold by combining experimental studies with atomic-resolution molecular dynamics simulations. General principles as to how such complex architectures form so rapidly are now emerging from in-depth studies of a few proteins.

The complete folding pathway of a protein from nanoseconds to microseconds

Combining experimental and simulation data to describe all of the structures and the pathways involved in folding a protein is problematical. Transition states can be mapped experimentally by phi values, but the denatured state is very difficult to analyse under conditions that favour folding. Also computer simulation at atomic resolution is currently limited to about a microsecond or less. Ultrafast-folding proteins fold and unfold on timescales accessible by both approaches, so here we study the folding pathway of the three-helix bundle protein Engrailed homeodomain. Experimentally, the protein collapses in a microsecond to give an intermediate with much native alpha-helical secondary structure, which is the major component of the denatured state under conditions that favour folding. A mutant protein shows this state to be compact and contain dynamic, native-like helices with unstructured side chains. In the transition state between this and the native state, the structure of the helices is nearly fully formed and their docking is in progress, approximating to a classical diffusion-collision model. Molecular dynamics simulations give rate constants and structural details highly consistent with experiment, thereby completing the description of folding at atomic resolution.

beta-Hairpins, alpha-helices, and the intermediates among the secondary structures in the energy landscape of a peptide from a distal beta-hairpin of SH3 domain

Energy landscape of a peptide, extracted from a distal beta-hairpin of src SH3 domain, in explicit water was obtained with the multicanonical molecular dynamics. A variety of beta-hairpins with various strand-strand hydrogen bonds were found in the energy landscape at 300 K. There was no energy barrier between random-coil and hairpins. Thus, the peptide conformation can easily change from the random-coil to the hairpins in the thermal fluctuations at 300 K. The landscape also included two clusters of alpha-helices, among which an energy barrier existed, and besides, these helix clusters were separated from the other conformations. Thus, the free-energy barrier exists among the helices and the other conformations. Intermediate clusters were found between the helix and the hairpin clusters. The current study showed that the isolated state of this peptide in water fluctuates among random-coil, beta-hairpin, and alpha-helix. In SH3 domain, which has a topology of mainly beta-protein, the whole-protein folding may proceed when the segment is folded in the beta-hairpin and the other parts of the protein are coupled with the beta-hairpin in an energetically or kinetically favorite way.

An NMR view of the folding process of a CheY mutant at the residue level

The folding of CheY mutant F14N/V83T was studied at 75 residues by NMR. Fluorescence, NMR, and sedimentation equilibrium studies at different urea and protein concentrations reveal that the urea-induced unfolding of this CheY mutant includes an on-pathway molten globule-like intermediate that can associate off-pathway. The populations of native and denatured forms have been quantified from a series of 15N-1H HSQC spectra recorded under increasing concentrations of urea. A thermodynamic analysis of these data provides a detailed picture of the mutant's unfolding at the residue level: (1) the transition from the native state to the molten globule-like intermediate is highly cooperative, and (2) the unfolding of this state is sequential and yields another intermediate showing a collapsed N-terminal domain and an unfolded C-terminal tail. This state presents a striking similarity to the kinetic transition state of the CheY folding pathway.

Structural dynamics of the membrane translocation domain of colicin E9 and its interaction with TolB

In order for the 61 kDa colicin E9 protein toxin to enter the cytoplasm of susceptible cells and kill them by hydrolysing their DNA, the colicin must interact with the outer membrane BtuB receptor and Tol translocation pathway of target cells. The translocation function is located in the N-terminal domain of the colicin molecule. (1)H, (1)H-(1)H-(15)N and (1)H-(13)C-(15)N NMR studies of intact colicin E9, its DNase domain, minimal receptor-binding domain and two N-terminal constructs containing the translocation domain showed that the region of the translocation domain that governs the interaction of colicin E9 with TolB is largely unstructured and highly flexible. Of the expected 80 backbone NH resonances of the first 83 residues of intact colicin E9, 61 were identified, with 43 of them being assigned specifically. The absence of secondary structure for these was shown through chemical shift analyses and the lack of long-range NOEs in (1)H-(1)H-(15)N NOESY spectra (tau(m)=200 ms). The enhanced flexibility of the region of the translocation domain containing the TolB box compared to the overall tumbling rate of the protein was identified from the relatively large values of backbone and tryptophan indole (15)N spin-spin relaxation times, and from the negative (1)H-(15)N NOEs of the backbone NH resonances. Variable flexibility of the N-terminal region was revealed by the (15)N T(1)/T(2) ratios, which showed that the C-terminal end of the TolB box and the region immediately following it was motionally constrained compared to other parts of the N terminus. This, together with the observation of inter-residue NOEs involving Ile54, indicated that there was some structural ordering, resulting most probably from the interactions of side-chains. Conformational heterogeneity of parts of the translocation domain was evident from a multiplicity of signals for some of the residues. Im9 binding to colicin E9 had no effect on the chemical shifts or other NMR characteristics of the region of colicin E9 containing the TolB recognition sequence, though the interaction of TolB with intact colicin E9 bound to Im9 did affect resonances from this region. The flexibility of the translocation domain of colicin E9 may be connected with its need to recognise protein partners that assist it in crossing the outer membrane and in the translocation event itself.

Unspecific hydrophobic stabilization of folding transition states

Here we present a method for determining the inference of non-native conformations in the folding of a small domain, alpha-spectrin Src homology 3 domain. This method relies on the preservation of all native interactions after Tyr/Phe exchanges in solvent-exposed, contact-free positions. Minor changes in solvent exposure and free energy of the denatured ensemble are in agreement with the reverse hydrophobic effect, as the Tyr/Phe mutations slightly change the polypeptide hydrophilic/hydrophobic balance. Interestingly, more important Gibbs energy variations are observed in the transition state ensemble (TSE). Considering the small changes induced by the H/OH replacements, the observed energy variations in the TSE are rather notable, but of a magnitude that would remain undetected under regular mutations that alter the folded structure free energy. Hydrophobic residues outside of the folding nucleus contribute to the stability of the TSE in an unspecific nonlinear manner, producing a significant acceleration of both unfolding and refolding rates, with little effect on stability. These results suggest that sectors of the protein transiently reside in non-native areas of the landscape during folding, with implications in the reading of phi values from protein engineering experiments. Contrary to previous proposals, the principle that emerges is that non-native contacts, or conformations, could be beneficial in evolution and design of some fast folding proteins.

The chicken-egg scenario of protein folding revisited

What is the first step in protein folding - hydrophobic collapse (compaction) or secondary structure formation? It is still not clear if the major driving force in protein folding is hydrogen bonding or hydrophobic interactions or both. We analyzed data on the conformational characteristics of 41 globular proteins in native and partially folded conformational states. Our analysis shows that a good correlation exists between relative decrease in hydrodynamic volume and increase in secondary structure content. No compact equilibrium intermediates lacking secondary structure, or highly ordered non-compact species, were found. This correlation provides experimental support for the hypothesis that hydrophobic collapse occurs simultaneously with formation of secondary structure in the early stages of the protein folding.

Long-range interactions within a nonnative protein

Protein folding and unfolding are coupled to a range of biological phenomena, from the regulation of cellular activity to the onset of neurodegenerative diseases. Defining the nature of the conformations sampled in nonnative proteins is crucial for understanding the origins of such phenomena. We have used a combination of nuclear magnetic resonance (NMR) spectroscopy and site-directed mutagenesis to study unfolded states of the protein lysozyme. Extensive clusters of hydrophobic structure exist within the wild-type protein even under strongly denaturing conditions. These clusters involve distinct regions of the sequence but are all disrupted by a single point mutation that replaced residue Trp62 with Gly located at the interface of the two major structural domains in the native state. Thus, nativelike structure in the denatured protein is stabilized by the involvement of Trp62 in nonnative and long-range interactions.

NMR identification of local structural preferences in HIV-1 protease tethered heterodimer in 6 M guanidine hydrochloride

Understanding protein folding requires complete characterization of all the states of the protein present along the folding pathways. For this purpose nuclear magnetic resonance (NMR) has proved to be a very powerful technique because of the great detail it can unravel regarding the structure and dynamics of protein molecules. We report here NMR identification of local structural preferences in human immunodeficiency virus-1 protease in the 'unfolded state'. Analyses of the chemical shifts revealed the presence of local structural preferences many of which are native-like, and there are also some non-native structural elements. Three-bond H(N)-H(alpha) coupling constants that could be measured for some of the N-terminal and C-terminal residues are consistent with the native-like beta-structure. Unusually shifted 15N and amide proton chemical shifts of residues adjacent to some prolines and tryptophans also indicate the presence of some structural elements. These conclusions are supported by amide proton temperature coefficients and nuclear Overhauser enhancement data. The locations of the residues exhibiting preferred structural propensities on the crystal structure of the protein, give useful insights into the folding mechanism of this protein.

Slow dynamics in folded and unfolded states of an SH3 domain

(15)N relaxation dispersion experiments were applied to the isolated N-terminal SH3 domain of the Drosophila protein drk (drkN SH3) to study microsecond to second time scale exchange processes. The drkN SH3 domain exists in equilibrium between folded (F(exch)) and unfolded (U(exch)) states under nondenaturing conditions in a ratio of 2:1 at 20 degrees C, with an average exchange rate constant, k(ex), of 2.2 s(-1) (slow exchange on the NMR chemical shift time scale). Consequently a discrete set of resonances is observed for each state in NMR spectra. Within the U(exch) ensemble there is a contiguous stretch of residues undergoing conformational exchange on a micros/ms time scale, likely due to local, non-native hydrophobic collapse. For these residues both the F(exch) <--> U(exch) conformational exchange process and the micros/ms exchange event within the U(exch) state contribute to the (15)N line width and can be analyzed using CPMG-based (15)N relaxation dispersion measurements. The contribution of both processes to the apparent relaxation rate can be deconvoluted numerically by combining the experimental (15)N relaxation dispersion data with results from an (15)N longitudinal relaxation experiment that accurately quantifies exchange rates in slow exchanging systems (Farrow, N. A.; Zhang, O.; Forman-Kay, J. D.; Kay, L. E. J. Biomol. NMR 1994, 4, 727-734). A simple, generally applicable analytical expression for the dependence of the effective transverse relaxation rate constant on the pulse spacing in CPMG experiments has been derived for a two-state exchange process in the slow exchange limit, which can be used to fit the experimental data on the global folding/unfolding transition. The results illustrate that relaxation dispersion experiments provide an extremely sensitive tool to probe conformational exchange processes in unfolded states and to obtain information on the free energy landscape of such systems.

Energy landscape of a peptide consisting of alpha-helix, 3(10)-helix, beta-turn, beta-hairpin, and other disordered conformations

The energy landscape of a peptide [Ace-Lys-Gln-Cys-Arg-Glu-Arg-Ala-Nme] in explicit water was studied with a multicanonical molecular dynamics simulation, and the AMBER parm96 force field was used for the energy calculation. The peptide was taken from the recognition helix of the DNA-binding protein, c-MYB: A rugged energy landscape was obtained, in which the random-coil conformations were dominant at room temperature. The CD spectra of the synthesized peptide revealed that it is in the random state at room temperature. However, the 300 K canonical ensemble, Q(300K), contained alpha-helix, 3(10)-helix, beta-turn, and beta-hairpin structures with small but notable probabilities of existence. The complete alpha-helix, imperfect alpha-helix, and random-coil conformations were separated from one another in the conformational space. This means that the peptide must overcome energy barriers to form the alpha-helix. The overcoming process may correspond to the hydrogen-bond rearrangements from peptide-water to peptide-peptide interactions. The beta-turn, imperfect 3(10)-helix, and beta-hairpin structures, among which there are no energy barriers at 300 K, were embedded in the ensemble of the random-coil conformations. Two types of beta-hairpin with different beta-turn regions were observed in Q(300K). The two beta-hairpin structures may have different mechanisms for the beta-hairpin formation. The current study proposes a scheme that the random state of this peptide consists of both ordered and disordered conformations. In contrast, the energy landscape obtained from the parm94 force field was funnel like, in which the peptide formed the helical conformation at room temperature and random coil at high temperature.

NMR and SAXS characterization of the denatured state of the chemotactic protein CheY: implications for protein folding initiation

The denatured state of a double mutant of the chemotactic protein CheY (F14N/V83T) has been analyzed in the presence of 5 M urea, using small angle X-ray scattering (SAXS) and heteronuclear magnetic resonance. SAXS studies show that the denatured protein follows a wormlike chain model. Its backbone can be described as a chain composed of rigid elements connected by flexible links. A comparison of the contour length obtained for the chain at 5 M urea with the one expected for a fully expanded chain suggests that approximately 25% of the residues are involved in residual structures. Conformational shifts of the alpha-protons, heteronuclear (15)N-[(1)H] NOEs and (15)N relaxation properties have been used to identify some regions in the protein that deviate from a random coil behavior. According to these NMR data, the protein can be divided into two subdomains, which largely coincide with the two folding subunits identified in a previous kinetic study of the folding of the protein. The first of these subdomains, spanning residues 1-70, is shown here to exhibit a restricted mobility as compared to the rest of the protein. Two regions, one in each subdomain, were identified as deviating from the random coil chemical shifts. Peptides corresponding to these sequences were characterized by NMR and their backbone (1)H chemical shifts were compared to those in the intact protein under identical denaturing conditions. For the region located in the first subdomain, this comparison shows that the observed deviation from random coil parameters is caused by interactions with the rest of the molecule. The restricted flexibility of the first subdomain and the transient collapse detected in that subunit are consistent with the conclusions obtained by applying the protein engineering method to the characterization of the folding reaction transition state.

Dramatic stabilization of an SH3 domain by a single substitution: roles of the folded and unfolded states

The N-terminal SH3 domain of the Drosophila drk protein (drkN SH3) exists in equilibrium between folded and unfolded states under non-denaturing buffer conditions. In order to examine the origins of this instability, we have made mutations in the domain and characterized the thermodynamics and kinetics of folding. Results of substitutions of negatively charged residues to neutral amino acid residues suggest that the large electrostatic potential of the domain does not play a dominant role in the instability of the domain. Sequence alignment of a large number of SH3 domains reveals that the drkN SH3 domain has a threonine (T22) at a position corresponding to an otherwise highly conserved glycine residue in the diverging beta-turn connecting the beta3 and beta4 strands. Mutation of T22 to glycine results in significant stabilization of the drkN SH3 domain by 2.5 kcal/mole. To further characterize the basis for the stabilization of the T22 mutant relative to wild-type, we made additional mutant proteins with substitutions of residue T22. A strong correlation is seen between protein stability or folding rate and propensity for native beta-turn structure at this position. Correlation of folding rates with AGADIR predictions of non-native helical structure in the diverging turn region, along with our previous NMR evidence for non-native structure in this region of the unfolded state of the drkN SH3 domain, suggests that the free energy of the unfolded state also plays a role in stability. This result highlights the importance of both folded and unfolded states for understanding protein stability.

Polymer principles of protein calorimetric two-state cooperativity

The experimental calorimetric two-state criterion requires the van't Hoff enthalpy DeltaH(vH) around the folding/unfolding transition midpoint to be equal or very close to the calorimetric enthalpy DeltaH(cal) of the entire transition. We use an analytical model with experimental parameters from chymotrypsin inhibitor 2 to elucidate the relationship among several different van't Hoff enthalpies used in calorimetric analyses. Under reasonable assumptions, the implications of these DeltaH(vH)'s being approximately equal to DeltaH(cal) are equivalent: Enthalpic variations among denatured conformations in real proteins are much narrower than some previous lattice-model estimates, suggesting that the energy landscape theory "folding to glass transition temperature ratio" T(f) /T(g) may exceed 6.0 for real calorimetrically two-state proteins. Several popular three-dimensional lattice protein models, with different numbers of residue types in their alphabets, are found to fall short of the high experimental standard for being calorimetrically two-state. Some models postulate a multiple-conformation native state with substantial pre-denaturational energetic fluctuations well below the unfolding transition temperature, or predict a significant post-denaturational continuous conformational expansion of the denatured ensemble at temperatures well above the transition point, or both. These scenarios either disagree with experiments on protein size and dynamics, or are inconsistent with conventional interpretation of calorimetric data. However, when empirical linear baseline subtractions are employed, the resulting DeltaH(vH)/DeltaH(cal)'s for some models can be increased to values closer to unity, and baseline subtractions are found to correspond roughly to an operational definition of native-state conformational diversity. These results necessitate a re-assessment of theoretical models and experimental interpretations.

Partially unfolded species populated during equilibrium denaturation of the beta-sheet protein Y74W apo-pseudoazurin

Apo-pseudoazurin is a single domain cupredoxin. We have engineered a mutant in which a unique tryptophan replaces the tyrosine residue found in the tyrosine corner of this Greek key protein, a region that has been proposed to have an important role in folding. Equilibrium denaturation of Y74W apo-pseudoazurin demonstrated multistate unfolding in urea (pH 7.0, 0.5 M Na(2)SO(4) at 15 degrees C), in which one or more partially folded species are populated in 4. 3 M urea. Using a variety of biophysical techniques, we show that these species, on average, have lost a substantial portion of the native secondary structure, lack fixed tertiary packing involving tryptophan and tyrosine residues, are less compact than the native state as determined by fluorescence lifetimes and time-resolved anisotropy, but retain significant residual structure involving the trytophan residue. Peptides ranging in length from 11 to 30 residues encompassing this region, however, did not contain detectable nonrandom structure, suggesting that long-range interactions are important for stabilizing the equilibrium partially unfolded species in the intact protein. On the basis of these results, we suggest that the equilibrium denaturation of Y74W apo-pseudoazurin generates one or more partially unfolded species that are globally collapsed and retain elements of the native structure involving the newly introduced tryptophan residue. We speculate on the role of such intermediates in the generation of the complex Greek key fold.

Similarities between the spectrin SH3 domain denatured state and its folding transition state

We have expanded our description of the energy landscape for folding of the SH3 domain of chicken alpha-spectrin by a detailed structural characterization of its denatured state ensemble (DSE). This DSE is significantly populated under mildly acidic conditions in equilibrium with the folded state. Evidence from heteronuclear nuclear magnetic resonance (NMR) experiments on (2)H, (15)N-labeled protein suggests the presence of conformers whose residual structure bears some resemblence to the structure of the folding transition state of this protein. NMR analysis in a mutant with an engineered, non-native alpha-helical tendency shows a significant amount of local non-native structure in the mutant, while the overall characteristics of the DSE are unchanged. Comparison with recent theoretical predictions of SH3 domain folding reactions reveals an interesting correlation with the predicted early events. Based on these results and recent data from other systems, we propose that the DSE of a protein will resemble the intermediate or transition state of its nearest rate-limiting step, as a consequence of simple energetic and kinetic principles.

Protein folding mechanisms: new methods and emerging ideas

During the past year, advances in our understanding of folding mechanisms have been made through detailed experimental and theoretical studies of a number of proteins. The development of new methods has allowed the earliest events in folding to be probed and the measurement of folding at the level of individual molecules is now possible, opening the door to exciting new experiments.

A tale of two secondary structure elements: when a beta-hairpin becomes an alpha-helix

In this work, we have analyzed the relative importance of secondary versus tertiary interactions in stabilizing and guiding protein folding. For this purpose, we have designed four different mutants to replace the alpha-helix of the GB1 domain by a sequence with strong beta-hairpin propensity in isolation. In particular, we have chosen the sequence of the second beta-hairpin of the GB1 domain, which populates the native conformation in aqueous solution to a significant extent. The resulting protein has roughly 30 % of its sequence duplicated and maintains the 3D-structure of the wild-type protein, but with lower stability (up to -5 kcal/mol). The loss of intrinsic helix stability accounts for about 80 % of the decrease in free energy, illustrating the importance of local interactions in protein stability. Interestingly enough, all the mutant proteins, included the one with the duplicated beta-hairpin sequence, fold with similar rates as the GB1 domain. Essentially, it is the nature of the rate-limiting step in the folding reaction that determines whether a particular interaction will speed up, or not, the folding rates. While local contacts are important in determining protein stability, residues involved in tertiary contacts in combination with the topology of the native fold, seem to be responsible for the specificity of protein structures. Proteins with non-native secondary structure tendencies can adopt stable folds and be as efficient in folding as those proteins with native-like propensities.