Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Dynamics of unfolded polypeptide chains as model for the earliest steps in protein folding

The rate of formation of intramolecular interactions in unfolded proteins determines how fast conformational space can be explored during folding. Characterization of the dynamics of unfolded proteins is therefore essential for the understanding of the earliest steps in protein folding. We used triplet-triplet energy transfer to measure formation of intrachain contacts in different unfolded polypeptide chains. The time constants (1/k) for contact formation over short distances are almost independent of chain length, with a maximum value of about 5 ns for flexible glycine-rich chains and of 12 ns for stiffer chains. The rates of contact formation over longer distances decrease with increasing chain length, indicating different rate-limiting steps for motions over short and long chain segments. The effect of the amino acid sequence on local chain dynamics was probed by using a series of host-guest peptides. Formation of local contacts is only sixfold slower around the stiffest amino acid (proline) compared to the most flexible amino acid (glycine). Good solvents for polypeptide chains like EtOH, GdmCl and urea were found to slow intrachain diffusion and to decrease chain stiffness. These data allow us to determine the time constants for formation of the earliest intrachain contacts during protein folding.     

Dynamics of unfolded polypeptide chains in crowded environment studied by fluorescence correlation spectroscopy

Proteins have evolved to fold and function within a cellular environment that is characterized by high macromolecular content. The earliest step of protein folding represents intrachain contact formation of amino acid residues within an unfolded polypeptide chain. It has been proposed that macromolecular crowding can have significant effects on rates and equilibria of biomolecular processes. However, the kinetic consequences on intrachain diffusion of polypeptides have not been tested experimentally, yet. Here, we demonstrate that selective fluorescence quenching of the oxazine fluorophore MR121 by the amino acid tryptophan (Trp) in combination with fast fluorescence correlation spectroscopy (FCS) can be used to monitor end-to-end contact formation rates of unfolded polypeptide chains. MR121 and Trp were incorporated at the terminal ends of polypeptides consisting of repetitive units of glycine (G) and serine (S) residues. End-to-end contact formation and dissociation result in "off" and "on" switching of MR121 fluorescence and underlying kinetics can be revealed in FCS experiments with nanosecond time resolution. We revisit previous experimental studies concerning the dependence of end-to-end contact formation rates on polypeptide chain length, showing that kinetics can be described by Gaussian chain theory. We further investigate effects of solvent viscosity and temperature on contact formation rates demonstrating that intrachain diffusion represents a purely diffusive, entropy-controlled process. Finally, we study the influence of macromolecular crowding on polypeptide chain dynamics. The data presented demonstrate that intrachain diffusion is fast in spite of hindered diffusion caused by repulsive interactions with macromolecules. Findings can be explained by effects of excluded volume reducing chain entropy and therefore accelerating the loop search process. Our results suggest that within a cellular environment the early formation of structural elements in unfolded proteins can still proceed quite efficiently in spite of hindered diffusion caused by high macromolecular content.     

A comparison of the folding characteristics of free and ribosome-tethered polypeptide chains using limited proteolysis and mass spectrometry

The kinetics and thermodynamics of protein folding are commonly studied in vitro by denaturing/renaturing intact protein sequences. How these folding mechanisms relate to de novo folding that occurs as the nascent polypeptide emerges from the ribosome is much less well understood. Here, we have employed limited proteolysis followed by mass spectrometry analyses to compare directly free and ribosome-tethered polypeptide chains of the Src-homology 3 (SH3) domain and its unfolded variant, SH3-m10. The disordered variant was found to undergo faster proteolysis than SH3. Furthermore, the trypsin cleavage patterns observed show minor, but significant, differences for the free and ribosome-bound nascent chains, with significantly fewer tryptic peptides detected in the presence of ribosome. The results highlight the utility of limited proteolysis coupled with mass spectrometry for the structural analysis of these complex systems, and pave the way for detailed future analyses by combining this technique with chemical labeling methods (for example, hydrogen-deuterium exchange, photochemical oxidation) to analyze protein folding in real time, including in the presence of additional ribosome-associated factors.     

Nascent β-Hairpin Formation of a Natively Unfolded Peptide Reveals the Role of Hydrophobic Contacts

Despite the important role of the unfolded states in protein stability, folding, and aggregation, they remain poorly understood due to the lack of residue-specific experimental data. Here, we explore features of the unfolded state of the NTL9 protein by applying all-atom replica-exchange simulations to the two fragment peptides NTL9(1–22) and NTL9(6–17). We found that while NTL9(6–17) is unstructured, NTL9(1–22) transiently folds as various β-hairpins, a fraction of which contain a native β-sheet. Interestingly, despite a large number of charged residues, the formation of backbone hydrogen bonds is concomitant with hydrophobic but not electrostatic contacts. Although the fragment peptides lack a proposed specific contact between Asp⁸ and Lys¹², the individually weak, nonspecific interactions with lysines together stabilize the charged Asp⁸, leading to a pK_a shift of nearly 0.5 units, in agreement with the NMR data. Taken together, our data suggest that the unfolded state of NTL9 likely contains a β-hairpin in segment 1–22 with sequence-distant hydrophobic contacts, thus lending support to a long-standing hypothesis that the unfolded states of proteins exhibit native-like topology with hydrophobic clusters.     

Conformational properties of a peptide model for unfolded alpha-helices

Models of protein folding often hypothesize that the first step is local secondary structure formation. The assumption is that unfolded polypeptide chains possess an intrinsic propensity to form these local secondary structures. On the basis of this idea, it is tempting to model the local conformational properties of unfolded proteins using well-established residue secondary structure propensities, in particular, alpha-helix forming propensities. We have used spectroscopic methods to investigate the conformational behavior of a host-guest series of peptides designed to model unfolded alpha-helices. A suitable peptide model for unfolded alpha-helices was determined from studies of the length dependence of the conformational properties of alanine-based peptides. The chosen host peptide possessed a small, detectable, alpha-helix content. Substituting various representative guest residues into the central position of the host peptide at times changed the conformational behavior dramatically, and often in ways that could not be predicted from known alpha-helix forming propensities. The data presented can be used to rationalize some of these propensities. However, it is clear that secondary structure propensities cannot be used to predict the local conformational properties of unfolded proteins.     

Reversible thermal denaturation of a 60-kDa genetically engineered beta-sheet polypeptide

A de novo 687-amino-acid residue polypeptide with a regular 32-amino-acid repeat sequence, (GA)₃GY(GA)₃GE(GA)₃GH(GA)₃GK, forms large β-sheet assemblages that exhibit remarkable folding properties and, as well, form fibrillar structures. This construct is an excellent tool to explore the details of β-sheet formation yielding intimate folding information that is otherwise difficult to obtain and may inform folding studies of naturally occurring materials. The polypeptide assumes a fully folded antiparallel β-sheet/turn structure at room temperature, and yet is completely and reversibly denatured at 125°C, adopting a predominant polyproline II conformation. Deep ultraviolet Raman spectroscopy indicated that melting/refolding occurred without any spectroscopically distinct intermediates, yet the relaxation kinetics depend on the initial polypeptide state, as would be indicative of a non-two-state process. Thermal denaturation and refolding on cooling appeared to be monoexponential with characteristic times of ~1 and ~60 min, respectively, indicating no detectable formation of hairpin-type nuclei in the millisecond timescale that could be attributed to nonlocal “nonnative” interactions. The polypeptide folding dynamics agree with a general property of β-sheet proteins, i.e., initial collapse precedes secondary structure formation. The observed folding is much faster than expected for a protein of this size and could be attributed to a less frustrated free-energy landscape funnel for folding. The polypeptide sequence suggests an important balance between the absence of strong nonnative contacts (salt bridges or hydrophobic collapse) and limited repulsion of charged side chains.     

Reversible denaturation of the gene V protein of bacteriophage f1

The guanidine hydrochloride (GuHCl)-induced denaturation of the gene V protein of bacteriophage f1 has been studied, using the chemical reactivity of a cysteine residue that is buried in the folded protein and the circular dichroism (CD) at 211 and 229 nm as measures of the fraction of polypeptide chains in the folded form. It is found that this dimeric protein unfolds in a single cooperative transition from a folded dimer to two unfolded monomers. A folded, monomeric form of the gene V protein was not detected at equilibrium. The kinetics of unfolding of the gene V protein in 3 M GuHCl and the refolding in 2 M GuHCl are also consistent with a transition between a folded dimer and two unfolded monomers. The GuHCl concentration dependence of the rates of folding and unfolding suggests that the transition state for folding is near the folded conformation.     

Kinetics of internal-loop formation in polypeptide chains: a simulation study

The speed of simple diffusional motions, such as the formation of loops in the polypeptide chain, places one physical limit on the speed of protein folding. Many experimental studies have explored the kinetics of formation of end-to-end loops in polypeptide chains; however, protein folding more often requires the formation of contacts between interior points on the chain. One expects that, for loops of fixed contour length, interior loops will form more slowly than end-to-end loops, owing to the additional excluded volume associated with the "tails". We estimate the magnitude of this effect by generating ensembles of randomly coiled, freely jointed chains, and then using the theory of Szabo, Schulten, and Schulten to calculate the corresponding contact formation rates for these ensembles. Adding just a few residues, to convert an end-to-end loop to an internal loop, sharply decreases the contact rate. Surprisingly, the relative change in rate increases for a longer loop; sufficiently long tails, however, actually reverse the effect and accelerate loop formation slightly. Our results show that excluded volume effects in real, full-length polypeptides may cause the rates of loop formation during folding to depart significantly from the values derived from recent loop-formation experiments on short peptides.     

High temperature unfolding simulations of the TRPZ1 peptide

We report high temperature molecular dynamics simulations of the unfolding of the TRPZ1 peptide using an explicit model for the solvent. The system has been simulated for a total of 6 μs with 100-ns minimal continuous stretches of trajectory. The populated states along the simulations are identified by monitoring multiple observables, probing both the structure and the flexibility of the conformations. Several unfolding and refolding transition pathways are sampled and analyzed. The unfolding process of the peptide occurs in two steps because of the accumulation of a metastable on-pathway intermediate state stabilized by two native backbone hydrogen bonds assisted by nonnative hydrophobic interactions between the tryptophan side chains. Analysis of the un/folding kinetics and classical commitment probability calculations on the conformations extracted from the transition pathways show that the rate-limiting step for unfolding is the disruption of the ordered native hydrophobic packing (Trp-zip motif) leading from the native to the intermediate state. But, the speed of the folding process is mainly determined by the transition from the completely unfolded state to the intermediate and specifically by the closure of the hairpin loop driven by formation of two native backbone hydrogen bonds and hydrophobic contacts between tryptophan residues. The temperature dependence of the unfolding time provides an estimate of the unfolding activation enthalpy that is in agreement with experiments. The unfolding time extrapolated to room temperature is in agreement with the experimental data as well, thus providing a further validation to the analysis reported here.     

Measuring “Unmeasurable” Folding Kinetics of Proteins by Single-Molecule Force Spectroscopy

Folding and unfolding are fundamental biological processes in cell and are important for the biological functions of proteins. Characterizing the folding and unfolding kinetics of proteins is important for understanding the energetic landscape leading to the active native conformations of these molecules. However, the thermal or chemical-induced unfolding of many proteins is irreversible in vitro, precluding characterization of the folding kinetics of such proteins, just as it is impossible to “un-boil” an egg. Irreversible unfolding often manifests as irreversible aggregation of unfolded polypeptide chains, which typically occurs between denatured protein molecules in response to the exposure of hydrophobic residues to solvent. An example of such a protein where thermal denaturation results in irreversible aggregation is the β-1,4 endoxylanase from Bacillus circulans (BCX). Here, we report the use of single-molecule atomic force microscopy to directly measure the folding kinetics of BCX in vitro. By mechanically unfolding BCX, we essentially allowed only one unfolded molecule to exist in solution at a given time, effectively eliminating the possibility for aggregation. We found that BCX can readily refold back to the native state, allowing us to measure its folding kinetics for the first time. Our results demonstrate that single-molecule force-spectroscopy-based methods can adequately tackle the challenge of “un-boiling eggs”, providing a general methodology to characterize the folding kinetics of many proteins that suffer from irreversible denaturation and thus cannot be characterized using traditional equilibrium methodologies.     

Helix nucleation kinetics from molecular simulations in explicit solvent

We study the reversible folding/unfolding of short Ala and Gly-based peptides by molecular dynamics simulations of all-atom models in explicit water solvent. A kinetic analysis shows that the formation of a first alpha-helical turn occurs within 0.1-1 ns, in agreement with the analyses of laser temperature jump experiments. The unfolding times exhibit Arrhenius temperature dependence. For a rapidly nucleating all-Ala peptide, the helix nucleation time depends only weakly on temperature. For a peptide with enthalpically competing turn-like structures, helix nucleation exhibits an Arrhenius temperature dependence, corresponding to the unfolding of enthalpic traps in the coil ensemble. An analysis of structures in a "transition-state ensemble" shows that helix-to-coil transitions occur predominantly through breaking of hydrogen bonds at the helix ends, particularly at the C-terminus. The temperature dependence of the transition-state ensemble and the corresponding folding/unfolding pathways illustrate that folding mechanisms can change with temperature, possibly complicating the interpretation of high-temperature unfolding simulations. The timescale of helix formation is an essential factor in molecular models of protein folding. The rapid helix nucleation observed here suggests that transient helices form early in the folding event.     

Hydrogen-bond energetics drive helix formation in membrane interfaces

The free energy cost ΔG of partitioning many unfolded peptides into membrane interfaces is unfavorable due to the cost of partitioning backbone peptide bonds. The partitioning cost is dramatically reduced if the peptide bonds participate in hydrogen bonds. The reduced cost underlies secondary structure formation by amphiphilic peptides partitioned into membrane interfaces through a process referred to as partitioning-folding coupling. This coupling is characterized by the free energy reduction per residue, ?G(res) that drives folding. There is some debate about the correct value of ?G(res) and its dependence on the hydrophobic moment (μ(H)) of amphiphilic α-helical peptides. We show how to compute ?G(res) correctly. Using published data for two families of peptides with different hydrophobic moments and charges, we find that ?G(res) does not depend upon μ(H). The best estimate of ?G(res) is -0.37 ± 0.02 kcal mol(-1). This article is part of a Special Issue entitled: Membrane protein structure and function.     

Kinetics control preferential heterodimer formation of platelet-derived growth factor from unfolded A- and B-chains

The folding and assembly of platelet-derived growth factor (PDGF), a potent mitogen involved in wound-healing processes and member of the cystine knot growth factor family, was studied. The kinetics of the formation of disulfide-bonded dimers were investigated under redox reshuffling conditions starting either from unfolded and reduced PDGF-A- or B-chains or an equimolar mixture of both chains. It is shown that in all cases the formation of disulfide-bonded dimers is a very slow process occurring in the time scale of hours with a first-order rate-determining step. The formation of disulfide-bonded PDGF-AA or PDGF-BB homodimers displayed identical kinetics, indicating that both monomeric forms as well as the dimerized homodimer have similar folding and assembly pathways. In contrast, the formation of the heterodimer occurred three times more rapidly compared with the formation of the homodimers. As both monomeric forms revealed similar renaturation kinetics, it can be concluded that the first-order rate-determining folding step does not occur during monomer folding but must be attributed to conformational rearrangements of the dimerized, not yet disulfide-bonded protein. These structural rearrangements allow a more rapid formation of intermolecular disulfide bonds between the two different monomers of a heterodimer compared with the formation of the disulfide bonds between two identical monomers. The preferential formation of disulfide-bonded heterodimers from an equimolar mixture of unfolded A- and B-chains is thus a kinetically controlled process. Moreover, similar activation enthalpies for the formation of all different isoforms suggest that faster heterodimerization is controlled by entropic factors.     

Effects of urea and guanidine · HCl on the folding and unfolding of pancreatic trypsin inhibitor

Concentrated solutions of urea and of guanidine · HCl produced a random spectrum of single-disulphide forms of the polypeptide chain of the pancreatic trypsin inhibitor. Guanidine · HCl also unfolded completely, with accompanying interchange of disulphide bonds, the two-disulphide form of this protein in the native-like conformation; urea produced an equilibrium mixture in which one-quarter of the molecules had the native-like conformation and disulphide bonds. The unfolded forms of the protein in the denaturants were very flexible polypeptide chains. The observations suggest that urea and guanidine · HCl are denaturants because they produce essentially equally favourable solvation of all portions of a polypeptide.The energetics of the conformational transitions involved in folding and unfolding of the inhibitor were determined in urea and compared with those observed in its absence. The denaturant lowers the stability of the native, folded inhibitor relative to that of the reduced, unfolded state by 6.5 kilocalories per mole; the greatest part of this apparent free-energy difference was expressed at the two-disulphide stage of folding. The results are consistent with other indications that most of the favourable interactions stabilizing the native conformation of this protein are not encountered until the final stage of folding, when all may occur simultaneously.The unfolded one- and two-disulphide species produced in guanidine · HCl were trapped, and their rearrangement to the normal intermediates followed after removal of the denaturant. The random single-disulphide species, with one exception, reverted very rapidly to the non-random spectrum of intermediates normally observed during folding; this confirms that these species are normally rapidly interconverted and that normal refolding of the reduced protein is not dependent kinetically upon residual stable conformation in the reduced protein. The unfolded two-disulphide species refolded to the native-like conformation more slowly, but appeared to pass through the same intermediates normally observed during refolding from the fully reduced state.     

Residual Structure of Unfolded Ubiquitin as Revealed by Hydrogen/Deuterium-Exchange 2D NMR

The characterization of residual structures persistent in unfolded proteins in concentrated denaturant solution is currently an important issue in studies of protein folding because the residual structure present, if any, in the unfolded state may form a folding initiation site and guide the subsequent folding reactions. Here, we studied the hydrogen/deuterium (H/D)-exchange behavior of unfolded human ubiquitin in 6 M guanidinium chloride. We employed a dimethylsulfoxide (DMSO)-quenched H/D-exchange NMR technique with the use of spin desalting columns, which allowed us to perform a quick medium exchange from 6 M guanidinium chloride to a quenching DMSO solution. Based on the backbone resonance assignment of ubiquitin in the DMSO solution, we successfully investigated the H/D-exchange kinetics of 60 identified peptide amide groups in the ubiquitin sequence. Although a majority of these amide groups were not protected, certain amide groups involved in a middle helix (residues 23–34) and an N-terminal β-hairpin (residues 2–16) were significantly protected with a protection factor of 2.1–4.2, indicating that there were residual structures in unfolded ubiquitin and that these amide groups were more than 52% hydrogen bonded in the residual structures. We show that the hydrogen-bonded residual structures in the α-helix and the β-hairpin are formed even in 6 M guanidinium chloride, suggesting that these residual structures may function as a folding initiation site to guide the subsequent folding reactions of ubiquitin.     

Relationship between the native-state hydrogen exchange and folding pathways of a four-helix bundle protein

The hydrogen exchange behavior of a four-helix bundle protein in low concentrations of denaturant reveals some partially unfolded forms that are significantly more stable than the fully unfolded state. Kinetic folding of the protein, however, is apparently two-state in the absence of the accumulation of early folding intermediates. The partially unfolded forms are either as folded as or more folded than the rate-limiting transition state and appear to represent the major intermediates in a folding and unfolding reaction. These results are consistent with the suggestion that partially unfolded intermediates may form after the rate-limiting step for small proteins with apparent two-state folding kinetics.     

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.     

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.     

NMR and CD spectroscopy show that imino acid restriction of the unfolded state leads to efficient folding

Protein folding is determined by molecular features in the unfolded state, as well as the native folded structure. In the unfolded state, imino acids both restrict conformational space and present cis-trans isomerization barriers to folding. Because of its high proline and hydroxyproline content, the collagen triple-helix offers an opportunity to characterize the impact of imino acids on the unfolded state and folding kinetics. Here, NMR and CD spectroscopy are used to characterize the role of imino acids in a triple-helical peptide, T1-892, which contains an 18-residue sequence from type I collagen and a C-terminal (Gly-Pro-Hyp)(4) domain. The replacement of Pro or Hyp by an Ala in the (Gly-Pro-Hyp)(4) region significantly decreases the folding rate at low but not high concentrations, consistent with less efficient nucleation. To understand the molecular basis of the decreased folding rate, changes in the unfolded as well as the folded states of the peptides were characterized. While the trimer states of the peptides are all similar, NMR dynamics studies show monomers with all trans (Gly-Pro-Hyp)(4) are less flexible than monomers containing Pro --> Ala or Hyp --> Ala substitutions. Nucleation requires all trans bonds in the (Gly-Pro-Hyp)(4) domain and the constrained monomer state of the all trans nucleation domain in T1-892 increases its competency to initiate triple-helix formation and illustrates the impact of the unfolded state on folding kinetics.