Evidence for the principle of minimal frustration in the evolution of protein folding landscapes.

Theoretical and experimental studies have firmly established that protein folding can be described by a funneled energy landscape. This funneled energy landscape is the result of foldable protein sequences evolving following the principle of minimal frustration, which allows proteins to rapidly fold to their native biologically functional conformations. For a protein family with a given functional fold, the principle of minimal frustration suggests that, independent of sequence, all proteins within this family should fold with similar rates. However, depending on the optimal living temperature of the organism, proteins also need to modulate their thermodynamic stability. Consequently, the difference in thermodynamic stability should be primarily caused by differences in the unfolding rates. To test this hypothesis experimentally, we performed comprehensive thermodynamic and kinetic analyses of 15 different proteins from the thioredoxin family. Eight of these thioredoxins were extant proteins from psychrophilic, mesophilic, or thermophilic organisms. The other seven protein sequences were obtained using ancestral sequence reconstruction and can be dated back over 4 billion years. We found that all studied proteins fold with very similar rates but unfold with rates that differ up to three orders of magnitude. The unfolding rates correlate well with the thermodynamic stability of the proteins. Moreover, proteins that unfold slower are more resistant to proteolysis. These results provide direct experimental support to the principle of minimal frustration hypothesis.

Intramolecular diffusion controls aggregation of the PAPf39 peptide.

The 39-residue fragment of human prostatic acidic phosphatase (PAP) is found in high concentrations in semen and easily form fibrils. Previous work has shown that fibrillization is accelerated with a deletion of the first 8, mostly charged residues and it was hypothesized that fibrillization depended on the dynamics of these peptides. To test this hypothesis we have measured the intramolecular diffusion of the full length and 8-residue deletion peptides at two different pHs and found a correlation with fibrillization lag time. These results can be explained by a simple kinetic model of the early stages of aggregation in which oligomerization is controlled by the rate of peptide reconfiguration.

Modulation of folding energy landscape by charge-charge interactions: linking experiments with computational modeling.

The kinetics of folding-unfolding of a structurally diverse set of four proteins optimized for thermodynamic stability by rational redesign of surface charge-charge interactions is characterized experimentally. The folding rates are faster for designed variants compared with their wild-type proteins, whereas the unfolding rates are largely unaffected. A simple structure-based computational model, which incorporates the Debye-Hückel formalism for the electrostatics, was used and found to qualitatively recapitulate the experimental results. Analysis of the energy landscapes of the designed versus wild-type proteins indicates the differences in refolding rates may be correlated with the degree of frustration of their respective energy landscapes. Our simulations indicate that naturally occurring wild-type proteins have frustrated folding landscapes due to the surface electrostatics. Optimization of the surface electrostatics seems to remove some of that frustration, leading to enhanced formation of native-like contacts in the transition-state ensembles (TSE) and providing a less frustrated energy landscape between the unfolded and TS ensembles. Macroscopically, this results in faster folding rates. Furthermore, analyses of pairwise distances and radii of gyration suggest that the less frustrated energy landscapes for optimized variants are a result of more compact unfolded and TS ensembles. These findings from our modeling demonstrates that this simple model may be used to: (i) gain a detailed understanding of charge-charge interactions and their effects on modulating the energy landscape of protein folding and (ii) qualitatively predict the kinetic behavior of protein surface electrostatic interactions.

Structural and thermodynamic characterization of the recognition of the S100-binding peptides TRTK12 and p53 by calmodulin.

Calmodulin (CaM) is a multifunctional messenger protein that activates a wide variety of signaling pathways in eukaryotic cells in a calcium-dependent manner. CaM has been proposed to be functionally distinct from the S100 proteins, a related family of eukaryotic calcium-binding proteins. Previously, it was demonstrated that peptides derived from the actin-capping protein, TRTK12, and the tumor-suppressor protein, p53, interact with multiple members of the S100 proteins. To test the specificity of these peptides, they were screened using isothermal titration calorimetry against 16 members of the human S100 protein family, as well as CaM, which served as a negative control. Interestingly, both the TRTK12 and p53 peptides were found to interact with CaM. These interactions were further confirmed by both fluorescence and nuclear magnetic resonance spectroscopies. These peptides have distinct sequences from the known CaM target sequences. The TRTK12 peptide was found to independently interact with both CaM domains and bind with a stoichiometry of 2:1 and dissociations constants Kd,C-term = 2 ± 1 µM and Kd,N-term = 14 ± 1 µM. In contrast, the p53 peptide was found to interact only with the C-terminal domain of CaM, Kd,C-term = 2 ± 1 µM, 25°C. Using NMR spectroscopy, the locations of the peptide binding sites were mapped onto the structure of CaM. The binding sites for both peptides were found to overlap with the binding interface for previously identified targets on both domains of CaM. This study demonstrates the plasticity of CaM in target binding and may suggest a possible overlap in target specificity between CaM and the S100 proteins.

Novel interactions of the TRTK12 peptide with S100 protein family members: specificity and thermodynamic characterization.

The S100 protein family consists of small, dimeric proteins that exert their biological functions in response to changing calcium concentrations. S100B is the best-studied member and has been shown to interact with more than 20 binding partners in a calcium-dependent manner. The TRTK12 peptide, derived from the consensus binding sequence for S100B, has previously been found to interact with S100A1 and has been proposed to be a general binding partner of the S100 family. To test this hypothesis and gain a better understanding of the specificity of binding for the S100 proteins, 16 members of the human S100 family were screened against this peptide and its alanine variants. Novel interactions were found with only two family members, S100P and S100A2, indicating that TRTK12 selectively interacts with a small subset of the S100 proteins. Substantial promiscuity was observed in the binding site of S100B thereby accommodating variations in the peptide sequence, while S100A1, S100A2, and S100P exhibited larger differences in the binding constants for the TRTK12 alanine variants. This suggests that single-point substitutions can be used to selectively modulate the affinity of TRTK12 peptides for individual S100 proteins. This study has important implications for the rational drug design of inhibitors for the S100 proteins, which are involved in a variety of cancers and neurodegenerative diseases.

Denatured states of low-complexity polypeptide sequences differ dramatically from those of foldable sequences.

How the primary sequence of a protein encodes conformational preferences that operate early in folding to promote efficient formation of the correct native topology is still poorly understood. To address this issue, we have prepared a set of yeast iso-1-cytochrome c variants that contain polyalanine inserts ranging from 6 to 30 residues in length near the N terminus of the protein. We study the thermodynamics and kinetics of His-heme loop formation in the denatured state at 3 and 6 M guanidine-HCl concentration. We find that polyalanine closely approximates a random coil with excluded volume giving scaling exponents, nu(3), for equilibrium loop formation of 2.26 +/- 0.13 and 1.97 +/- 0.04 in 3 and 6 M guanidine-HCl, respectively. The rate of loop breakage initially decreases and then becomes independent of loop size as would be expected for a random coil. Comparison with previously reported data for denatured state His-heme loop formation for iso-1-cytochrome c and Rhodopseudomonas palustris cytochrome c', shows that foldable sequences deviate significantly from random coil behavior and that the deviation is fold-dependent.

Thermodynamics of loop formation in the denatured state of rhodopseudomonas palustris cytochrome c': scaling exponents and the reconciliation problem.

The observation that denatured proteins yield scaling exponents, nu, consistent with random-coil behavior and yet can also have pockets of residual or nonrandom structure has been termed the "reconciliation problem". To provide greater insight into the denatured state of a foldable sequence, we have measured histidine-heme loop formation equilibria in the denatured state of a class II c-type cytochrome, cytochrome c' from Rhodopseudomonas palustris. We have prepared a series of variants that provide His-heme loop stabilities, pK(loop)(His), for loop sizes ranging from 10 to 111 residues at intervals of 7 to 11 residues along the sequence of the protein. We observe a scaling exponent for loop formation, nu(3), of 2.5+/-0.3. Theoretical values for nu(3) range from 1.8 to 2.4; thus, the observed nu(3) is consistent with random-coil behavior. However, in contrast to data for loop formation as a function of loop size obtained with peptides of homogeneous sequence, we observe considerable scatter about the linear dependence of loop stability on loop size. Thus, foldable sequences behave very differently from homogeneous peptide sequences. The observed scatter suggests that there is considerable variation in the conformational properties along the backbone of a foldable sequence, consistent with alternating compact and extended regions. With regard to the reconciliation problem, it is evident that a scaling exponent consistent with a random coil is necessary but not sufficient to demonstrate random-coil behavior.

Importance of contact persistence in denatured state loop formation: kinetic insights into sequence effects on nucleation early in folding.

Protein folding is dependent on the formation and persistence of simple loops early in folding. Ease of loop formation and persistence is believed to be dependent on the steric interactions of the residues involved in loop formation. We have previously investigated this factor in the denatured state of iso-1-cytochrome c using a five-amino-acid insert in front of a unique histidine in the N-terminal region of the protein. Previously, we reported that the apparent pK(a) values of loop formation for the most flexible (all Gly) and least flexible (all Ala) insert were, within error, the same. We evaluate whether this observation is due to differences in the persistence of loop contacts or due to effects of local sequence sterics and main-chain hydration on the persistence length of the chain. We also test whether sequence order affects loop formation. Here, we report kinetic results coupled to further mutagenesis of the insert to discern between these possibilities. We find that the amino acid-glycine versus alanine-next to the loop forming histidine has a dominant effect on loop kinetics and equilibria. A glycine in this position speeds loop breakage relative to alanine, resulting in less stable loops. At high percentage of Gly in the insert, rates of loop formation and breakage exactly compensate, leading to a leveling out in loop stability. Loop formation rates also increase with glycine content, inconsistent with poly-Gly segments being more extended than previously suspected due to main-chain hydration or local sterics. Unlike loop breakage rates, loop formation rates are insensitive to local sequence. Together, these observations suggest that contact persistence plays a more important role in defining the "folding code" than rates of loop formation.

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

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

Sequence composition effects on denatured state loop formation in iso-1-cytochrome c variants: polyalanine versus polyglycine inserts.

Protein folding is dependent on the formation and persistence of simple loops during the earliest events of the folding process. Ease of loop formation and persistence is believed to be dependent on the steric properties of the residues involved in loop formation. We have investigated this conformational factor in the denatured state of iso-1-cytchrome c using a five alanine insert in front of a unique histidine in the N-terminal region of the protein. The alanine residues have then been progressively substituted with sterically less-constrained glycine residues. Guanidine-HCl unfolding shows that all variants have a free energy of unfolding of approximately 2 kcal/mol. The low stability of these variants is well accounted for by stabilization of the denatured state by histidine-heme loop formation. The stability of the 22 residue histidine-heme loop has been measured in 3 M guanidine hydrochloride for all variants. Surprisingly, relative to alanine, glycine has only a very modest effect on equilibrium loop stability. Thus, the greater flexibility that glycine confers on the main-chain provides no advantage in terms of the persistence of simple loops early in folding. The underlying basis for the similar behavior of loops with polyalanine versus polyglycine inserts is discussed in terms of the current knowledge of the structure and loop formation kinetics of glycine versus alanine-rich peptides.