Stability and folding properties of a model beta-sheet protein, Escherichia coli CspA

From Alpha to Beta - a co-translational way to fold?

Protein folding in the cell is largely a co-translational process occurring during protein synthesis on the ribosome. It has become evident that co-translational folding is characteristic to almost every protein in the cell of pro- and eukaryotic origin that are single and multidomain, single and multisubunit, cytosolic, secretory and membrane. Co-translational protein folding begins very early during the process of polypeptide chain synthesis on the ribosome, with some secondary structure elements forming inside the ribosomal tunnel and some tertiary structures forming inside the vestibule (lower/wider) region of the ribosomal exit tunnel. However, many details of co-translational folding remains incompletely understood. New data show that folding of a β-barrel protein begins with formation of an α-helix inside the ribosome that rearranges into a β-hairpin structure as the growing peptide reaches the wider/vestibule region of the exit tunnel. While it was previously suggested that such scenario can take place on the ribosome, the new data provide the first experimental evidence in support of this notion.

A switch from α-helical to β-strand conformation during co-translational protein folding

Cellular proteins begin to fold as they emerge from the ribosome. The folding landscape of nascent chains is not only shaped by their amino acid sequence but also by the interactions with the ribosome. Here, we combine biophysical methods with cryo‐EM structure determination to show that folding of a β‐barrel protein begins with formation of a dynamic α‐helix inside the ribosome. As the growing peptide reaches the end of the tunnel, the N‐terminal part of the nascent chain refolds to a β‐hairpin structure that remains dynamic until its release from the ribosome. Contacts with the ribosome and structure of the peptidyl transferase center depend on nascent chain conformation. These results indicate that proteins may start out as α‐helices inside the tunnel and switch into their native folds only as they emerge from the ribosome. Moreover, the correlation of nascent chain conformations with reorientation of key residues of the ribosomal peptidyl‐transferase center suggest that protein folding could modulate ribosome activity.

Domain topology, stability, and translation speed determine mechanical force generation on the ribosome

The concomitant folding of a nascent protein domain with its synthesis can generate mechanical forces that act on the ribosome and alter translation speed. Such changes in speed can affect the structure and function of the newly synthesized protein as well as cellular phenotype. The domain properties that govern force generation have yet to be identified and understood, and the influence of translation speed is unknown because all reported measurements have been carried out on arrested ribosomes. Here, using coarse-grained molecular simulations and statistical mechanical modeling of protein synthesis, we demonstrate that force generation is determined by a domain's stability and topology, as well as translation speed. The statistical mechanical models we create predict how force profiles depend on these properties. These results indicate that force measurements on arrested ribosomes will not always accurately reflect what happens in a cell, especially for slow-folding domains, and suggest the possibility that certain domain properties may be enriched or depleted across the structural proteome of organisms through evolutionary selection pressures to modulate protein synthesis speed and posttranslational protein behavior.

Estimating the mean first passage time of protein misfolding

Most theoretical and experimental studies confirm that proteins fold in the time scale of microseconds to milliseconds, but the kinetics of the protein misfolding remains largely unexplored. The kinetics of unfolding-folding-misfolding equilibrium in proteins is formulated in the analytical framework of the Master equation. The folded, unfolded and the misfolded state are characterized in terms of their respective contacts. The Mean First Passage Time (MFPT) to acquire the misfolded conformation from the native or folded state is derived from this equation with different boundary conditions. The MFPT is found to be practically independent of the length of the protein, the number of native contacts and the rate constant for the misfolded to the folded state. The results obtained from the survival probability are directly correlated to the age of onset and appearance of misfolding diseases in humans.

The Role of Electrostatic Interactions in Folding of β-Proteins

Atomic-level molecular dynamic simulations are capable of fully folding structurally diverse proteins; however, they are limited in their ability to accurately represent electrostatic interactions. Here we have experimentally tested the role of charged residues on stability and folding kinetics of one of the most widely simulated β-proteins, the WW domain. The folding of wild type Pin1 WW domain, which has two positively charged residues in the first turn, was compared to the fast folding mutant FiP35 Pin1, which introduces a negative charge into the first turn. A combination of FTIR spectroscopy and laser-induced temperature-jump coupled with infrared spectroscopy was used to probe changes in the amide I region. The relaxation dynamics of the peptide backbone, β-sheets and β-turns, and negatively charged aspartic acid side chain of FiP35 were measured independently by probing the corresponding bands assigned in the amide I region. Folding is initiated in the turns and the β-sheets form last. While the global folding mechanism is in good agreement with simulation predictions, we observe changes in the protonation state of aspartic acid during folding that have not been captured by simulation methods. The protonation state of aspartic acid is coupled to protein folding; the apparent pKa of aspartic acid in the folded protein is 6.4. The dynamics of the aspartic acid follow the dynamics of the intermediate phase, supporting assignment of this phase to formation of the first hairpin. These results demonstrate the importance of electrostatic interactions in turn stability and formation of extended β-sheet structures.

Folding simulations for proteins with diverse topologies are accessible in days with a physics-based force field and implicit solvent

The millisecond time scale needed for molecular dynamics simulations to approach the quantitative study of protein folding is not yet routine. One approach to extend the simulation time scale is to perform long simulations on specialized and expensive supercomputers such as Anton. Ideally, however, folding simulations would be more economical while retaining reasonable accuracy, and provide feedback on structure, stability and function rapidly enough if partnered directly with experiment. Approaches to this problem typically involve varied compromises between accuracy, precision, and cost; the goal here is to address whether simple implicit solvent models have become sufficiently accurate for their weaknesses to be offset by their ability to rapidly provide much more precise conformational data as compared to explicit solvent. We demonstrate that our recently developed physics-based model performs well on this challenge, enabling accurate all-atom simulated folding for 16 of 17 proteins with a variety of sizes, secondary structure, and topologies. The simulations were carried out using the Amber software on inexpensive GPUs, providing ∼1 μs/day per GPU, and >2.5 ms data presented here. We also show that native conformations are preferred over misfolded structures for 14 of the 17 proteins. For the other 3, misfolded structures are thermodynamically preferred, suggesting opportunities for further improvement.

WW domain folding complexity revealed by infrared spectroscopy

Although the intrinsic tryptophan fluorescence of proteins offers a convenient probe of protein folding, interpretation of the fluorescence spectrum is often difficult because it is sensitive to both global and local changes. Infrared (IR) spectroscopy offers a complementary measure of structural changes involved in protein folding, because it probes changes in the secondary structure of the protein backbone. Here we demonstrate the advantages of using multiple probes, infrared and fluorescence spectroscopy, to study the folding of the FBP28 WW domain. Laser-induced temperature jumps coupled with fluorescence or infrared spectroscopy have been used to probe changes in the peptide backbone on the submillisecond time scale. The relaxation dynamics of the β-sheets and β-turn were measured independently by probing the corresponding IR bands assigned in the amide I region. Using these wavelength-dependent measurements, we observe three kinetics phases, with the fastest process corresponding to the relaxation kinetics of the turns. In contrast, fluorescence measurements of the wild-type WW domain and tryptophan mutants exhibit single-exponential kinetics with a lifetime that corresponds to the slowest phase observed by infrared spectroscopy. Mutant sequences provide evidence of an intermediate dry molten globule state. The slowest step in the folding of this WW domain is the tight packing of the side chains in the transition from the dry molten globule intermediate to the native structure. This study demonstrates that using multiple complementary probes enhances the interpretation of protein folding dynamics.

Dynamics of an ultrafast folding subdomain in the context of a larger protein fold

Small fast folding subdomains with low contact order have been postulated to facilitate the folding of larger proteins. We have tested this idea by determining how the fastest folding linear β-hairpin, CLN025, which folds on the nanosecond time scale, folds within the context of a two-hairpin WW domain system, which folds on the microsecond time scale. The folding of the wild type FBP28 WW domain was compared to constructs in which each of the loops was replaced by CLN025. A combination of FTIR spectroscopy and laser-induced temperature-jump coupled with infrared spectroscopy was used to probe changes in the peptide backbone. The relaxation dynamics of the β-sheets and β-turn were measured independently by probing the corresponding bands assigned in the amide I region. The folding rate of the CLN025 β-hairpin is unchanged within the larger protein. Insertion of the β-hairpin into the second loop results in an overall stabilization of the WW domain and a relaxation lifetime five times faster than the parent WW domain. In both mutants, folding is initiated in the turns and the β-sheets form last. These results demonstrate that fast folding subdomains can be used to speed the folding of more complex proteins, and that the folding dynamics of the subdomain is unchanged within the context of the larger protein.

The entropic nature of protein thermal stabilization

We performed thermodynamic analysis of temperature-induced unfolding of mesophilic and thermophilic proteins. It was shown that the variability in protein thermostability associated with pH-dependent unfolding or linked to the substitution of amino acid residues on the protein surface is evidence of the governing role of the entropy factor. Numerical values of conformational components in enthalpy, entropy and free energy which characterize protein unfolding in the "gas phase" were obtained. Based on the calculated absolute values of entropy and free energy, a model of protein unfolding is proposed in which the driving force is the conformational entropy of native protein, as an energy of the heat motion (T·S(NC)) increasing with temperature and acting as an factor devaluating the energy of intramolecular weak bonds in the transition state.

Early turn formation and chain collapse drive fast folding of the major cold shock protein CspA of Escherichia coli

The folding mechanism of the β-sheet protein CspA, the major cold shock protein of Escherichia coli, was previously reported to be a concerted, two-state process. We have reexamined the folding of CspA using multiple spectroscopic probes of the equilibrium transition and laser-induced temperature jump (T-jump) to achieve better time resolution of the kinetics. Equilibrium temperature-dependent Fourier transform infrared (1634 cm(-1)) and tryptophan fluorescence measurements reveal probe-dependent thermal transitions with midpoints (T(m)) of 66 ± 1 and 61 ± 1 °C, respectively. Singular-value decomposition analysis with global fitting of the temperature-dependent infrared (IR) difference spectra reveals two spectral components with distinct melting transitions with different midpoints. T-jump relaxation measurements of CspA probed by IR and fluorescence spectroscopy show probe-dependent multiexponential kinetics characteristic of non-two-state folding. The frequency-dependent IR transients all show biphasic relaxation with average time constants of 50 ± 7 and 225 ± 25 μs at a T(f) of 77 °C and almost equal amplitudes. Similar biphasic kinetics are observed using Trp fluorescence of the wild-type protein and the Y42W and T68W mutants, with comparable lifetimes. All of these observations support a model for the folding of CspA through a compact intermediate state. The transient IR and fluorescence spectra are consistent with a diffuse intermediate having β-turns and substantial β-sheet structure. The loop β3-β4 structure is likely not folded in the intermediate state, allowing substantial solvent penetration into the barrel structure.

Is there an en route folding intermediate for Cold shock proteins?

Cold shock proteins (Csps) play an important role in cold shock response of a diverse number of organisms ranging from bacteria to humans. Numerous studies of the Csp from various species showed that a two-state folding mechanism is conserved and the transition state (TS) appears to be very compact. However, the atomic details of the folding mechanism of Csp remain unclear. This study presents the folding mechanism of Csp in atomic detail using an all-atom Go model-based simulations. Our simulations predict that there may exist an en route intermediate, in which β strands 1-2-3 are well ordered and the contacts between β1 and β4 are almost developed. Such an intermediate might be too unstable to be detected in the previous fluorescence energy transfer experiments. The transition state ensemble has been determined from the P(fold) analysis and the TS appears even more compact than the intermediate state.

The effect of context on the folding of β-hairpins

Small β-hairpin peptides have been widely used as models for the folding of β-sheets. But how applicable is the folding of such models to β-structure in larger proteins with conventional hydrophobic cores? Here we present multiple unfolding simulations of three such proteins that contain the WW domain double hairpin β-sheet motif: cold shock protein A (CspA), cold shock protein B (CspB) and glucose permease IIA domain. We compare the behavior of the free motif in solution and in the context of proteins of different size and architecture. Both Csp proteins lost contacts between the double-hairpin motif and the protein core as the first step of unfolding and proceeded to unfold with loss of the third β-strand, similar to the isolated WW domain. The glucose permease IIA domain is a larger protein and the contacts between the motif and the core were not lost as quickly. Instead the unfolding pathway of glucose permease IIA followed a different pathway with β1 pulling away from the sheet first. Interestingly, when the double hairpin motif was excised from the glucose permease IIA domain and simulated in isolation in water it unfolded by the same pathway as the WW domain, indicating that it is tertiary interactions with the protein that alter the motif's unfolding not a sequence dependent effect on its intrinsic unfolding behavior. With respect to the unfolding of the hairpins, there was no consistent order to the loss of hydrogen bonds between the β-strands in the hairpins in any of the systems. Our results show that while the folding behavior of the isolated WW domain is generally consistent with the double hairpin motif's behavior in the cold shock proteins, it is not the case for the glucose permease IIA domain. So, one must be cautious in extrapolating findings from model systems to larger more complicated proteins where tertiary interactions can overwhelm intrinsic behavior.

Dynameomics: a consensus view of the protein unfolding/folding transition state ensemble across a diverse set of protein folds

The Dynameomics project aims to simulate a representative sample of all globular protein metafolds under both native and unfolding conditions. We have identified protein unfolding transition state (TS) ensembles from multiple molecular dynamics simulations of high-temperature unfolding in 183 structurally distinct proteins. These data can be used to study individual proteins and individual protein metafolds and to mine for TS structural features common across all proteins. Separating the TS structures into four different fold classes (all proteins, all-alpha, all-beta, and mixed alpha/beta and alpha +beta) resulted in no significant difference in the overall protein properties. The residues with the most contacts in the native state lost the most contacts in the TS ensemble. On average, residues beginning in an alpha-helix maintained more structure in the TS ensemble than did residues starting in beta-strands or any other conformation. The metafolds studied here represent 67% of all known protein structures, and this is, to our knowledge, the largest, most comprehensive study of the protein folding/unfolding TS ensemble to date. One might have expected broad distributions in the average global properties of the TS relative to the native state, indicating variability in the amount of structure present in the TS. Instead, the average global properties converged with low standard deviations across metafolds, suggesting that there are general rules governing the structure and properties of the TS.

Urea impedes the hydrophobic collapse of partially unfolded proteins

Proteins are denatured in aqueous urea solution. The nature of the molecular driving forces has received substantial attention in the past, whereas the question how urea acts at different phases of unfolding is not yet well understood at the atomic level. In particular, it is unclear whether urea actively attacks folded proteins or instead stabilizes unfolded conformations. Here we investigated the effect of urea at different phases of unfolding by molecular dynamics simulations, and the behavior of partially unfolded states in both aqueous urea solution and in pure water was compared. Whereas the partially unfolded protein in water exhibited hydrophobic collapses as primary refolding events, it remained stable or even underwent further unfolding steps in aqueous urea solution. Further, initial unfolding steps of the folded protein were found not to be triggered by urea, but instead, stabilized. The underlying mechanism of this stabilization is a favorable interaction of urea with transiently exposed, less-polar residues and the protein backbone, thereby impeding back-reactions. Taken together, these results suggest that, quite generally, urea-induced protein unfolding proceeds primarily not by active attack. Rather, thermal fluctuations toward the unfolded state are stabilized and the hydrophobic collapse of partially unfolded proteins toward the native state is impeded. As a result, the equilibrium is shifted toward the unfolded state.

Integrase-directed recovery of functional genes from genomic libraries

Large population sizes, rapid growth and 3.8 billion years of evolution firmly establish microorganisms as a major source of the planet's biological and genetic diversity. However, up to 99% of the microorganisms in a given environment cannot be cultured. Culture-independent methods that directly access the genetic potential of an environmental sample can unveil new proteins with diverse functions, but the sequencing of random DNA can generate enormous amounts of extraneous data. Integrons are recombination systems that accumulate open reading frames (gene cassettes), many of which code for functional proteins with enormous adaptive potential. Some integrons harbor hundreds of gene cassettes and evidence suggests that the gene cassette pool may be limitless in size. Accessing this genetic pool has been hampered since sequence-based techniques, such as hybridization or PCR, often recover only partial genes or a small subset of those present in the sample. Here, a three-plasmid genetic strategy for the sequence-independent recovery of gene cassettes from genomic libraries is described and its use by retrieving functional gene cassettes from the chromosomal integron of Vibrio vulnificus ATCC 27562 is demonstrated. By manipulating the natural activity of integrons, we can gain access to the caches of functional genes amassed by these structures.

Confocal spectrofluorimetric evidence for the hetero-aggregation of sequence-scrambled forms of two model all-beta sheet proteins

Can two unrelated proteins with deliberately compromised folding abilities, marked propensities to aggregate upon increase of protein concentration, and proclivities towards beta sheet formation, be caused to hetero-aggregate? We address this question here using the 'designer' backbone-reversed forms of two model all-beta sheet proteins, E. coli CspA and C. elegans HSP12.6, both earlier created and characterized by our group. These were covalently labeled with fluorescent dyes of well-resolved spectral characteristics [retro-CspA with FITC, and retro-HSP12.6 with TRITC] and then allowed to aggregate within the same reaction vessel. The resultant aggregates are shown by spectrofluorimetry-coupled confocal laser scanning microscopy to constitute uniform mixtures of both proteins, existing within every cylindrical volume element of approximately 200nm diameter, and comparable height, in all sections of the co-aggregated material suggesting that the two proteins do not selectively associate with copies of themselves during aggregation. Thus, it would appear that aggregation can occur without reference to protein molecular identity.

Structure and function of cold shock proteins in archaea

Archaea are abundant and drive critical microbial processes in the Earth's cold biosphere. Despite this, not enough is known about the molecular mechanisms of cold adaptation and no biochemical studies have been performed on stenopsychrophilic archaea (e.g., Methanogenium frigidum). This study examined the structural and functional properties of cold shock proteins (Csps) from archaea, including biochemical analysis of the Csp from M. frigidum. csp genes are present in most bacteria and some eucarya but absent from most archaeal genome sequences, most notably, those of all archaeal thermophiles and hyperthermophiles. In bacteria, Csps are small, nucleic acid binding proteins involved in a variety of cellular processes, such as transcription. In this study, archaeal Csp function was assessed by examining the ability of csp genes from psychrophilic and mesophilic Euryarchaeota and Crenarchaeota to complement a cold-sensitive growth defect in Escherichia coli. In addition, an archaeal gene with a cold shock domain (CSD) fold but little sequence identity to Csps was also examined. Genes encoding Csps or a CSD structural analog from three psychrophilic archaea rescued the E. coli growth defect. The three proteins were predicted to have a higher content of solvent-exposed basic residues than the noncomplementing proteins, and the basic residues were located on the nucleic acid binding surface, similar to their arrangement in E. coli CspA. The M. frigidum Csp was purified and found to be a single-domain protein that folds by a reversible two-state mechanism and to exhibit a low conformational stability typical of cold-adapted proteins. Moreover, M. frigidum Csp was characterized as binding E. coli single-stranded RNA, consistent with its ability to complement function in E. coli. The studies show that some Csp and CSD fold proteins have retained sufficient similarity throughout evolution in the Archaea to be able to function effectively in the Bacteria and that the function of the archaeal proteins relates to cold adaptation. The initial biochemical analysis of M. frigidum Csp has developed a platform for further characterization and demonstrates the potential for expanding molecular studies of proteins from this important archaeal stenopsychrophile.

BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

We exploit the availability of recent experimental data on a variety of proteins to develop a Web-based prediction algorithm (BPPred) to calculate several biophysical parameters commonly used to describe the folding process. These parameters include the equilibrium m-values, the length of proteins, and the changes upon unfolding in the solvent-accessible surface area, in the heat capacity, and in the radius of gyration. We also show that the knowledge of any one of these quantities allows an estimate of the others to be obtained, and describe the confidence limits with which these estimations can be made. Furthermore, we discuss how the kinetic m-values, or the Beta Tanford values, may provide an estimate of the solvent-accessible surface area and the radius of gyration of the transition state for protein folding. Taken together, these results suggest that BPPred should represent a valuable tool for interpreting experimental measurements, as well as the results of molecular dynamics simulations.

Folding rate prediction using n-order contact distance for proteins with two- and three-state folding kinetics

It is a challenging task to understand the relationship between sequences and folding rates of proteins. Previous studies are found that one of contact order (CO), long-range order (LRO), total contact distance (TCD), chain topology parameter (CTP), and effective length (Leff) has a significant correlation with folding rate of proteins. In this paper, we introduce a new parameter called n-order contact distance (nOCD) and use it to predict folding rate of proteins with two- and three-state folding kinetics. A good linear correlation between the folding rate logarithm lnkf and nOCD with n=1.2, alpha=0.6 is found for two-state folders (correlation coefficient is -0.809, P-value<0.0001) and n=2.8, alpha=1.5 for three-state folders (correlation coefficient is -0.816, P-value<0.0001). However, this correlation is completely absent for three-state folders with n=1.2, alpha=0.6 (correlation coefficient is 0.0943, P-value=0.661) and for two-state folders with n=2.8, alpha=1.5 (correlation coefficient is -0.235, P-value=0.2116). We also find that the average number of contacts per residue Pm in the interval of m for two-state folders is smaller than that for three-state folders. The probability distribution P(gamma) of residue having gamma pairs of contacts fits a Gaussian distribution for both two- and three-state folders. We observe that the correlations between square radius of gyration S2 and number of residues for two- and three-state folders are both good, and the correlation coefficient is 0.908 and 0.901, and the slope of the fitting line is 1.202 and 0.795, respectively. Maybe three-state folders are more compact than two-state folders. Comparisons with nTCD and nCTP are also made, and it is found that nOCD is the best one in folding rate prediction.

Entropy capacity determines protein folding

Search and study of the general principles that govern kinetics and thermodynamics of protein folding generate a new insight into the factors controlling this process. Here, based on the known experimental data and using theoretical modeling of protein folding, we demonstrate that there exists an optimal relationship between the average conformational entropy and the average energy of contacts per residue-that is, an entropy capacity-for fast protein folding. Statistical analysis of conformational entropy and number of contacts per residue for 5829 protein structures from four general structural classes (all-alpha, all-beta, alpha/beta, alpha+beta) demonstrates that each class of proteins has its own class-specific average number of contacts (class alpha/beta has the largest number of contacts) and average conformational entropy per residue (class all-alpha has the largest number of rotatable angles phi, psi, and chi per residue). These class-specific features determine the folding rates: alpha proteins are the fastest folding proteins, then follow beta and alpha+beta proteins, and finally alpha/beta proteins are the slowest ones. Our result is in agreement with the experimental folding rates for 60 proteins. This suggests that structural and sequence properties are important determinants of protein folding rates.

Heat stable ssDNA/RNA-binding activity of a wheat cold shock domain protein

The cold-induced wheat WCSP1 protein belongs to the cold shock domain (CSD) protein family. In prokaryotes and eukaryotes, the CSD functions as a nucleic acid-binding domain. Here, we demonstrated that purified recombinant WCSP1 is boiling soluble and binds ss/dsDNA and mRNA. Furthermore, boiled-WCSP1 retained its characteristic nucleic acid-binding activity. A WCSP1 deletion mutant, containing only a CSD, lost ssDNA/RNA-binding activity; while a mutant containing the CSD and the first glycine-rich region (GR) displayed the activity. These data indicated that the first GR of WCSP1 is necessary for the binding activity but is not for the heat stability of the protein.

Folding and Stability of a Primitive Protein

We have previously attempted to simulate domain creation in early protein evolution by recombining polypeptide segments from non-homologous proteins, and we have described the structure of one such de novo protein, 1b11, a segment-swapped tetramer with novel architecture. Here, we have analyzed the thermodynamic stability and folding kinetics of the 1b11 tetramer and its monomeric and dimeric intermediates, and of 1b11 mutants with changes at the domain interface. Denatured 1b11 polypeptides fold into transient, folded monomers with marginal stability (DeltaG<1kcalmol(-1)) which convert rapidly ( approximately 6x10(4)M(-1)s(-1)) into dimers (DeltaG=9.8kcal/mol) and then more slowly ( approximately 3M(-1)s(-1)) into tetramers (DeltaG=28kcalmol(-1)). Segment swapping takes place during dimerization, as suggested by mass spectroscopic analysis of covalently linked peptides derived from proteolysis of a disulfide-linked dimer. Our results confirm that segment swapping and associated oligomerization are both powerful ways of stabilizing proteins, and we suggest that this may have been a feature of early protein evolution.

Desolvation is a likely origin of robust enthalpic barriers to protein folding

Experimental data from global analyses of temperature (T) and denaturant dependence of the folding rates of small proteins led to an intrinsic enthalpic folding barrier hypothesis: to a good approximation, the T-dependence of folding rate under constant native stability conditions is Arrhenius. Furthermore, for a given protein, the slope of isostability folding rate versus 1/T is essentially independent of native stability. This hypothesis implies a simple relationship between chevron and Eyring plots of folding that is easily discernible when both sets of rates are expressed as functions of native stability. Using experimental data in the literature, we verify the predicted chevron-Eyring relationship for 14 proteins and determine their intrinsic enthalpic folding barriers, which vary approximately from 15 kcal/mol to 40 kcal/mol for different proteins. These enthalpic barriers do not appear to correlate with folding rates, but they exhibit correlation with equilibrium unfolding enthalpy at room temperature. Intrinsic enthalpic barriers with similarly high magnitudes apply as well to at least two cases of peptide-peptide and peptide-protein association, suggesting that these barriers are a hallmark of certain general and fundamental kinetic processes during folding and binding. Using a class of explicit-chain C(alpha) protein models with constant elementary enthalpic desolvation barriers between C(alpha) positions, we show that small microscopic pairwise desolvation barriers, which are a direct consequence of the particulate nature of water, can act cooperatively to give rise to a significant overall enthalpic barrier to folding. This theoretical finding provides a physical rationalization for the high intrinsic enthalpic barriers in protein folding energetics. Ramifications of entropy-enthalpy compensation in hydrophobic association for the height of enthalpic desolvation barrier are discussed.

Correspondence between anomalous m- and DeltaCp-values in protein folding

Proteins folding according to a classical two-state system characteristically show V-shaped chevron plots. We have previously interpreted the symmetrically curved chevron plot of the protein U1A as denaturant-dependent movements in the position of the transition state ensemble (TSE). S6, a structural analog of U1A, shows a classical V-shaped chevron plot indicative of straightforward two-state kinetics, but the mutant LA30 has a curved unfolding limb, which is most consistent with TSE mobility. The kinetic m-values (derivatives of the rate constants with respect to denaturant concentration) in themselves depend on denaturant concentration. To obtain complementary information about putative mobile TSEs, we have carried out a thermodynamic analysis of the three proteins, based on data for refolding and unfolding over the range 10 degrees C to 70 degrees C. The data at all temperatures can be fitted to two-state model systems. Importantly, for all three proteins the activation heat capacities are, within error, identical to the heat capacities measured in independent experiments under equilibrium conditions. Although the equilibrium heat capacities are essentially invariant with regard to denaturant concentration, the activation heat capacities, similar to the structurally equivalent kinetic m-values, show marked denaturant dependence. Furthermore, the values of beta++ at different denaturant concentrations measured by m-values and by heat capacity values are very similar. These observations are consistent with significant transition state movements within the framework of two-state folding. The basis for TSE movement appears to be enthalpic rather than entropic, suggesting that the binding energy of denaturant-protein interactions is a major determinant of the response of energy landscape contours to changing environments.

A Statistical Model for Predicting Protein Folding Rates from Amino Acid Sequence with Structural Class Information

Prediction of protein folding rates from amino acid sequences is one of the most important challenges in molecular biology. In this work, I have related the protein folding rates with physical-chemical, energetic and conformational properties of amino acid residues. I found that the classification of proteins into different structural classes shows an excellent correlation between amino acid properties and folding rates of two- and three-state proteins, indicating the importance of native state topology in determining the protein folding rates. I have formulated a simple linear regression model for predicting the protein folding rates from amino acid sequences along with structural class information and obtained an excellent agreement between predicted and experimentally observed folding rates of proteins; the correlation coefficients are 0.99, 0.96 and 0.95, respectively, for all-alpha, all-beta and mixed class proteins. This is the first available method, which is capable of predicting the protein folding rates just from the amino acid sequence with the aid of generic amino acid properties and structural class information.

Clitocypin, a new cysteine proteinase inhibitor, is monomeric: impact on the mechanism of folding

The molecular mass of clitocypin, a new type of cysteine proteinase inhibitor from the mushroom Clitocybe nebularis, has been determined by analytical ultracentrifugation and gel exclusion chromatography. The result is in agreement with the formula mass of 16.8 kDa, demonstrating that the inhibitor is a monomer in aqueous solution. This enables the kinetics of unfolding and refolding to be interpreted in terms of folding in a kinetically two state, highly cooperative transition from the thermally unfolded state.

Mechanism of prion propagation: amyloid growth occurs by monomer addition

Abundant nonfibrillar oligomeric intermediates are a common feature of amyloid formation, and these oligomers, rather than the final fibers, have been suggested to be the toxic species in some amyloid diseases. Whether such oligomers are critical intermediates for fiber assembly or form in an alternate, potentially separable pathway, however, remains unclear. Here we study the polymerization of the amyloidogenic yeast prion protein Sup35. Rapid polymerization occurs in the absence of observable intermediates, and both targeted kinetic and direct single-molecule fluorescence measurements indicate that fibers grow by monomer addition. A three-step model (nucleation, monomer addition, and fiber fragmentation) accurately accounts for the distinctive kinetic features of amyloid formation, including weak concentration dependence, acceleration by agitation, and sigmoidal shape of the polymerization time course. Thus, amyloid growth can occur by monomer addition in a reaction distinct from and competitive with formation of potentially toxic oligomeric intermediates.

Identification and transcriptional control of Caulobacter crescentus genes encoding proteins containing a cold shock domain

The cold shock proteins are small peptides that share a conserved domain, called the cold shock domain (CSD), that is important for nucleic acid binding. The Caulobacter crescentus genome has four csp genes that encode proteins containing CSDs. Three of these (cspA, cspB, and cspC) encode peptides of about 7 kDa and are very similar to the cold shock proteins of other bacteria. Analysis by reverse transcription-PCR of the fourth gene (cspD), which was previously annotated as encoding a 7-kDa protein, revealed that the mRNA is larger and probably encodes a putative 21-kDa protein, containing two CSDs. A search in protein sequences databases revealed that this new domain arrangement has thus far only been found among deduced peptides of alpha-proteobacteria. Expression of each Caulobacter csp gene was studied both in response to cold shock and to growth phase, and we have found that only cspA and cspB are induced by cold shock, whereas cspC and cspD are induced at stationary phase, with different induction rates. The transcription start sites were determined for each gene, and a deletion mapping of the cspD promoter region defined a sequence required for maximal levels of expression, indicating that regulation of this gene occurs at the transcriptional level. Deletion of cspA, but not cspD, caused a reduction in viability when cells were incubated at 10 degrees C for prolonged times, suggesting that cspA is important for adaptation to a low temperature.

Unification of the Folding Mechanisms of Non-two-state and Two-state Proteins

We have collected the kinetic folding data for non-two-state and two-state globular proteins reported in the literature, and investigated the relationships between the folding kinetics and the native three-dimensional structure of these proteins. The rate constants of formation of both the intermediate and the native state of non-two-state folders were found to be significantly correlated with protein chain length and native backbone topology, which is represented by the absolute contact order and sequence-distant native pairs. The folding rate of two-state folders, which is known to be correlated with the native backbone topology, apparently does not correlate significantly with protein chain length. On the basis of a comparison of the folding rates of the non-two-state and two-state folders, it was found that they are similarly dependent on the parameters that reflect the native backbone topology. This suggests that the mechanisms behind non-two-state and two-state folding are essentially identical. The present results lead us to propose a unified mechanism of protein folding, in which folding occurs in a hierarchical manner, reflecting the hierarchy of the native three-dimensional structure, as embodied in the case of non-two-state folding with an accumulation of the intermediate. Apparently, two-state folding is merely a simplified version of hierarchical folding caused either by an alteration in the rate-limiting step of folding or by destabilization of the intermediate.

Rapid collapse precedes the fast two-state folding of the cold shock protein

The cold shock protein Bc-Csp folds very rapidly in a reaction that is well described by a kinetic two-state mechanism without intermediates. We measured the shortening of six intra-protein distances during folding by Förster resonance energy transfer (FRET) in combination with stopped-flow experiments. Single tryptophan residues were engineered into the protein as the donors, and single 5-(((acetylamino)ethyl)amino)naphthalene-1-sulfonate (AEDANS) residues were placed as the acceptors at solvent-exposed sites of Bc-Csp. Their R0 value of about 22 A was well suited for following distance changes during the folding of this protein with a high sensitivity. The mutagenesis and the labeling did not alter the refolding kinetics. The changes in energy transfer during folding were monitored by both donor and acceptor emission and reciprocal effects were found. In two cases the donor-acceptor distances were similar in the unfolded and the folded state and, as a consequence, the kinetic changes in energy transfer upon folding were very small. For four donor/acceptor pairs we found that > or =50% of the increase in energy transfer upon folding occurred prior to the rate-limiting step of folding. This reveals that about half of the shortening of the intra-molecular distances upon folding has occurred already before the rate-limiting step and suggests that the fast two-state folding reaction of Bc-Csp is preceded by a very rapid collapse.

The importance of explicit chain representation in protein folding models: an examination of Ising-like models

A class of models that represents a protein chain as a sequence of "folded" and "unfolded" residues has recently been used to correlate rates and mechanisms of protein folding with the protein native structure. In order to better understand the conditions under which these "Ising-like" models apply, we compare results from this model to those obtained from an off-lattice model which uses the same potential function. We find that Ising-like models by construction impose folding via a highly sequential nucleation-condensation mechanism, which in turn leads to more rugged energy landscapes, fewer "pathways" to the native state, and in the specific case examined here, the cold shock protein A from Escherichia coli, a qualitative difference in the most likely order of events in folding.

Contact order revisited: influence of protein size on the folding rate

Guided by the recent success of empirical model predicting the folding rates of small two-state folding proteins from the relative contact order (CO) of their native structures, by a theoretical model of protein folding that predicts that logarithm of the folding rate decreases with the protein chain length L as L(2/3), and by the finding that the folding rates of multistate folding proteins strongly correlate with their sizes and have very bad correlation with CO, we reexamined the dependence of folding rate on CO and L in attempt to find a structural parameter that determines folding rates for the totality of proteins. We show that the Abs_CO = CO x L, is able to predict rather accurately folding rates for both two-state and multistate folding proteins, as well as short peptides, and that this Abs_CO scales with the protein chain length as L(0.70 +/- 0.07) for the totality of studied single-domain proteins and peptides.

Chain length is the main determinant of the folding rate for proteins with three-state folding kinetics

We demonstrate that chain length is the main determinant of the folding rate for proteins with the three-state folding kinetics. The logarithm of their folding rate in water (k(f)) strongly anticorrelates with their chain length L (the correlation coefficient being -0.80). At the same time, the chain length has no correlation with the folding rate for two-state folding proteins (the correlation coefficient is -0.07). Another significant difference of these two groups of proteins is a strong anticorrelation between the folding rate and Baker's "relative contact order" for the two-state folders and the complete absence of such correlation for the three-state folders.

Prediction of folding rates and transition-state placement from native-state geometry

A variety of experimental and theoretical studies have established that the folding process of monomeric proteins is strongly influenced by the topology of the native state. In particular, folding times have been shown to correlate well with the contact order, a measure of contact locality. Our investigation focuses on identifying additional topologic properties that correlate with experimentally measurable quantities, such as folding rates and transition-state placement, for both two- and three-state folders. The validation against data from 40 experiments shows that a particular topological property that measures the interdependence of contacts, termed cliquishness or clustering coefficient, can account with statistically significant accuracy both for the transition state placement and especially for folding rates. The observed correlations can be further improved by optimally combining the distinct topological information captured by cliquishness and contact order.

Structural characterization of thyroglobulin type-1 domains of equistatin

Equistatin is a protein composed of three thyroglobulin type-1 domains. It inhibits papain-like cysteine proteinases and the aspartic proteinase, cathepsin D. To determine the structural basis for this inhibition we cloned and expressed the separated domains (eq d-1, eq d-2, eq d-3) in Pichia pastoris. Kinetic constants for the interaction of eq d-1 with papain and that of eq d-2 with cathepsin D are of similar order (subnanomolar) and are comparable to the constants obtained for full-length equistatin. The target proteinase for the third domain remains unknown. Thus, we demonstrate here that thyroglobulin type-1 motifs per se are able to support specific structural features that enable them to inhibit proteases from different classes. The overall conformation of three domains in equistatin is such that the interaction of domains 1 or 2 with their respective target enzymes is not hindered sterically by either domain. In addition, we show that the interaction of eq d-2 with cathepsin D results in conformational changes, which is not the case for the eq d-1/papain interaction.

Bacterial cold shock responses

As a measure for molecular motion, temperature is one of the most important environmental factors for life as it directly influences structural and hence functional properties of cellular components. After a sudden increase in ambient temperature, which is termed heat shock, bacteria respond by expressing a specific set of genes whose protein products are designed to mainly cope with heat-induced alterations of protein conformation. This heat shock response comprises the expression of protein chaperones and proteases, and is under central control of an alternative sigma factor (sigma 32) which acts as a master regulator that specifically directs RNA polymerase to transcribe from the heat shock promotors. In a similar manner, bacteria express a well-defined set of proteins after a rapid decrease in temperature, which is termed cold shock. This protein set, however, is different from that expressed under heat shock conditions and predominantly comprises proteins such as helicases, nucleases, and ribosome-associated components that directly or indirectly interact with the biological information molecules DNA and RNA. Interestingly, in contrast to the heat shock response, to date no cold-specific sigma factor has been identified. Rather, it appears that the cold shock response is organized as a complex stimulon in which post-transcriptional events play an important role. In this review, we present a summary of research results that have been acquired in recent years by examinations of bacterial cold shock responses. Important processes such as cold signal perception, membrane adaptation, and the modification of the translation apparatus are discussed together with many other cold-relevant aspects of bacterial physiology and first attempts are made to dissect the cold shock stimulon into less complex regulatory subunits. Special emphasis is placed on findings concerning the nucleic acid-binding cold shock proteins which play a fundamental role not only during cold shock adaptation but also under optimal growth conditions.

Sequential vs. parallel protein-folding mechanisms: experimental tests for complex folding reactions

The recent emphasis on rough energy landscapes for protein folding reactions by theoreticians, and the many observations of complex folding kinetics by experimentalists provide a rationale for a brief literature survey of various empirical approaches for validating the underlying mechanisms. The determination of the folding mechanism is a key step in defining the energy surface on which the folding reactions occurs and in interpreting the effects of amino acid replacements on this reaction. Case studies that illustrate methods for differentiating between sequential and parallel channel folding mechanisms are presented. The ultimate goal of such efforts is to understand how the one-dimensional information contained in the amino acid sequence is rapidly and efficiently translated into three-dimensional structure.

Folding kinetics of the lipoic acid-bearing domain of human mitochondrial branched chain alpha-ketoacid dehydrogenase complex

A reversible two-step (native state<-->denatured state) folding mechanism based on equilibrium and stopped flow experiments is proposed for urea denaturation of the lipoyl-bearing domain (hbLBD) of human mitochondrial branched chain alpha-ketoacid dehydrogenase (BCKD) complex. The results from this circular dichroism (CD) and fluorescence study have ruled out populated kinetic or equilibrium intermediates on folding pathway of this beta-barrel domain under experimental conditions. Both studies suggested mono-exponential kinetics without any burst phases. Moreover the thermodynamic parameters DeltaG(H(2)O) and m obtained from the kinetic analysis are consistent with the equilibrium measurements.

Probing the free-energy surface for protein folding with single-molecule fluorescence spectroscopy

Protein folding is inherently a heterogeneous process because of the very large number of microscopic pathways that connect the myriad unfolded conformations to the unique conformation of the native structure. In a first step towards the long-range goal of describing the distribution of pathways experimentally, Förster resonance energy transfer (FRET) has been measured on single, freely diffusing molecules. Here we use this method to determine properties of the free-energy surface for folding that have not been obtained from ensemble experiments. We show that single-molecule FRET measurements of a small cold-shock protein expose equilibrium collapse of the unfolded polypeptide and allow us to calculate limits on the polypeptide reconfiguration time. From these results, limits on the height of the free-energy barrier to folding are obtained that are consistent with a simple statistical mechanical model, but not with the barriers derived from simulations using molecular dynamics. Unlike the activation energy, the free-energy barrier includes the activation entropy and thus has been elusive to experimental determination for any kinetic process in solution.

Folding, stability, and secondary structure of a new dimeric cysteine proteinase inhibitor

Clitocypin, a new type of cysteine proteinase inhibitor from the mushroom Clitocybe nebularis, is a 34-kDa homodimer lacking disulphide bonds, reported to have unusual stability properties. Sequence similarity is limited solely to certain proteins from mushrooms. Infrared spectroscopy shows that clitocypin is a high beta-structure protein which was lost at high temperatures. The far UV circular dichroism spectrum is not that of classical beta-structure, but similar to those of a group of small beta-strand proteins, with a peak at 189nm and a trough at 202nm. An aromatic peak at 232nm and infrared bands at 1633 and 1515cm(-1) associated with the peptide backbone and the tyrosine microenvironment, respectively, were used to characterize the thermal unfolding. The reversible transition has a midpoint at 67 degrees C, with DeltaG=34kJ/mol and DeltaH=300kJ/mol, and is, unusually, independent of protein concentration. The kinetics of thermal unfolding and refolding are slow, with activation energies of 167 and 44kJ/mol, respectively. A model for folding and assembly is discussed.

Single-domain antibody fragments with high conformational stability

A variety of techniques, including high-pressure unfolding monitored by Fourier transform infrared spectroscopy, fluorescence, circular dichroism, and surface plasmon resonance spectroscopy, have been used to investigate the equilibrium folding properties of six single-domain antigen binders derived from camelid heavy-chain antibodies with specificities for lysozymes, beta-lactamases, and a dye (RR6). Various denaturing conditions (guanidinium chloride, urea, temperature, and pressure) provided complementary and independent methods for characterizing the stability and unfolding properties of the antibody fragments. With all binders, complete recovery of the biological activity after renaturation demonstrates that chemical-induced unfolding is fully reversible. Furthermore, denaturation experiments followed by optical spectroscopic methods and affinity measurements indicate that the antibody fragments are unfolded cooperatively in a single transition. Thus, unfolding/refolding equilibrium proceeds via a simple two-state mechanism (N <--> U), where only the native and the denatured states are significantly populated. Thermally-induced denaturation, however, is not completely reversible, and the partial loss of binding capacity might be due, at least in part, to incorrect refolding of the long loops (CDRs), which are responsible for antigen recognition. Most interestingly, all the fragments are rather resistant to heat-induced denaturation (apparent T(m) = 60-80 degrees C), and display high conformational stabilities (DeltaG(H(2)O) = 30-60 kJ mole(-1)). Such high thermodynamic stability has never been reported for any functional conventional antibody fragment, even when engineered antigen binders are considered. Hence, the reduced size, improved solubility, and higher stability of the camelid heavy-chain antibody fragments are of special interest for biotechnological and medical applications.

Examination of the folding of E. coli CspA through tryptophan substitutions

Escherichia coli cold shock protein, CspA, folds very rapidly (time constant, tau = 4 msec) by an apparent two-state mechanism. However, recent time-resolved infrared (IR) temperature-jump experiments indicate that the folding trajectory of CspA may be more complicated. The sole tryptophan of wild-type CspA (Trp11), which is used to monitor the folding process by fluorescence spectroscopy, is located in an unusual aromatic cluster on the surface of CspA within the nucleic acid binding site. To gain a more global picture of the folding kinetics of CspA and to determine if there are any previously undetected intermediates, we have introduced a second tryptophan at three different surface locations in the protein. The three mutations did not significantly alter the tertiary structure of CspA, although two of the substitutions were found to be slightly stabilizing. Two-state folding, as detected by stopped-flow fluorescence spectroscopy, is preserved in all three mutants. These results indicate that the fast folding of CspA is driven by a concerted mechanism.

On the relationship between protein stability and folding kinetics: a comparative study of the N-terminal domains of RNase HI, E. coli and Bacillus stearothermophilus L9

There is currently a great deal of interest in proteins that fold in a single highly cooperative step. Particular attention has been focused on elucidating the factors that govern folding rates of simple proteins. Recently, the topology of the native state has been proposed to be the most important determinant of their folding rates. Here we report a comparative study of the folding of three topologically equivalent proteins that adapt a particularly simple alpha/beta fold. The folding kinetics of the N-terminal domain of RNase HI and the N-terminal domain of the ribosomal protein L9 from Escherichia coli (eNTL9) were compared to the previously characterized N-terminal domain of L9 from Bacillus stearothermophilus (bNTL9). This 6.2 kDa protein, which is one of simplest examples of the ABCalphaD motif, folds via a two-state mechanism on the millisecond to submillisecond time scale. The RNase HI domain and bNTL9 have very similar tertiary structures but there is little similarity in primary sequence. bNTL9 and eNTL9 share the same biological function and a similar primary sequence but differ significantly in stability. Fluorescence-detected stopped-flow experiments showed that the three proteins fold in a two-state fashion. The folding rates in the absence of denaturant were found to be very different, ranging form 21 s(-1) to 790 s(-1) at 10 degrees C. The diverse folding rates appear to reflect large differences in the stability of the proteins. When compared at an isostability point, the folding rates converged to a similar value and there is a strong linear correlation between the log of the folding rate and stability for this set of proteins. These observations are consistent with the idea that stability can play an important role in dictating relative folding rates among topologically equivalent proteins.

Solution NMR structure of the cold-shock protein from the hyperthermophilic bacterium Thermotoga maritima

Cold-shock proteins (Csps) are a subgroup of the cold-induced proteins preferentially expressed in bacteria and other organisms on reduction of the growth temperature below the physiological temperature. They are related to the cold-shock domain found in eukaryotes and are some of the most conserved proteins known. Their exact function is still not known, but translational regulation, possibly via RNA chaperoning, has been discussed. Here we present the structure of a hyperthermophilic member of the Csp family. The NMR solution structure of TmCsp from Thermotoga maritima, the hyperthermophilic member of this class of proteins, was solved on the basis of 1015 conformational constraints. It contains five beta strands combined in two antiparallel beta sheets making up a beta barrel structure, in which beta strands 1-4 are arranged in a Greek-key topology. The side chain of R2, which is exclusively found in thermophilic members of the Csp family, probably participates in a peripheral ion cluster involving residues D20, R2, E47 and K63, suggesting that the thermostability of TmCsp is based on the peripheral ion cluster around the side chain of R2.