Denatured state thermodynamics: residual structure, chain stiffness and scaling factors

Denatured State Conformational Biases in Three-Helix Bundles Containing Divergent Sequences Localize near Turns and Helix Capping Residues

Rhodopseudomonas palustris cytochrome c', a four-helix bundle, and the second ubiquitin-associated domain, UBA(2), a three-helix bundle from the human homologue of yeast Rad23, HHR23A, deviate from random coil behavior under denaturing conditions in a fold-specific manner. The random coil deviations in each of these folds occur near interhelical turns and loops in their tertiary structures. Here, we examine an additional three-helix bundle with an identical fold to UBA(2), but a highly divergent sequence, the first ubiquitin-associated domain, UBA(1), of HHR23A. We use histidine-heme loop formation methods, employing eight single histidine variants, to probe for denatured state conformational bias of a UBA(1) domain fused to the N-terminus of iso-1-cytochrome c (iso-1-Cytc). Guanidine hydrochloride (GuHCl) denaturation shows that the iso-1-Cytc domain unfolds first, followed by the UBA(1) domain. Denatured state (4 and 6 M GuHCl) histidine-heme loop formation studies show that as the size of the histidine-heme loop increases, loop stability decreases, as expected for the Jacobson-Stockmayer relationship. However, loops formed with His35, His31, and His15, of UBA(1), are 0.6-1.1 kcal/mol more stable than expected from the Jacobson-Stockmayer relationship, confirming the importance of deviations of the denatured state from random coil behavior near interhelical turns of helical domains for facilitating folding to the correct topology. For UBA(1) and UBA(2), hydrophobic clusters on either side of the turns partially explain deviations from random coil behavior; however, helix capping also appears to be important.

Helical Propensity Affects the Conformational Properties of the Denatured State of Cytochrome c'

Changing the helical propensity of a polypeptide sequence might be expected to affect the conformational properties of the denatured state of a protein. To test this hypothesis, alanines at positions 83 and 87 near the center of helix 3 of cytochrome c' from Rhodopseudomonas palustris were mutated to serine to decrease the stability of this helix. A set of 13 single histidine variants in the A83S/A87S background were prepared to permit assessment of the conformational properties of the denatured state using histidine-loop formation in 3 M guanidine hydrochloride. The data are compared with previous histidine-heme loop formation data for wild-type cytochrome c'. As expected, destabilization of helix 3 decreases the global stabilities of the histidine variants in the A83S/A87S background relative to the wild-type background. Loop stability versus loop size data yields a scaling exponent of 2.1 ± 0.2, similar to the value of 2.3 ± 0.2 obtained for wild-type cytochrome c'. However, the stabilities of all histidine-heme loops, which contain the helix 3 sequence segment, are increased in the A83S/A87S background compared to the wild-type background. Rate constants for histidine-heme loop breakage are similar for the wild-type and A83S/A87S variants. However, for histidine-heme loops that contain the helix 3 sequence segment, the rate constants for loop formation increase in the A83S/A87S background compared to the wild-type background. Thus, residual helical structure appears to stiffen the polypeptide chain slowing loop formation in the denatured state. The implications of these results for protein folding mechanisms are discussed.

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

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

Heightened Dynamics of the Oxidized Y48H Variant of Human Cytochrome c Increases Its Peroxidatic Activity

Proteins performing multiple biochemical functions are called "moonlighting proteins" or extreme multifunctional (EMF) proteins. Mitochondrial cytochrome c is an EMF protein that binds multiple partner proteins to act as a signaling molecule, transfers electrons in the respiratory chain, and acts as a peroxidase in apoptosis. Mutations in the cytochrome c gene lead to the disease thrombocytopenia, which is accompanied by enhanced apoptotic activity. The Y48H variant arises from one such mutation and is found in the 40-57 Ω-loop, the lowest-unfolding free energy substructure of the cytochrome c fold. A 1.36 Å resolution X-ray structure of the Y48H variant reveals minimal structural changes compared to the wild-type structure, with the axial Met80 ligand coordinated to the heme iron. Despite this, the intrinsic peroxidase activity is enhanced, implying that a pentacoordinate heme state is more prevalent in the Y48H variant, corroborated through determination of a Met80 "off rate" of >125 s^-1 compared to a rate of ∼6 s^-1 for the wild-type protein. Heteronuclear nuclear magnetic resonance measurements with the oxidized Y48H variant reveal heightened dynamics in the 40-57 Ω-loop and the Met80-containing 71-85 Ω-loop relative to the wild-type protein, illustrating communication between these substructures. Placed into context with the G41S cytochrome c variant, also implicated in thrombocytopenia, a dynamic picture associated with this disease relative to cytochrome c is emerging whereby increasing dynamics in substructures of the cytochrome c fold serve to facilitate an increased population of the peroxidatic pentacoordinate heme state in the following order: wild type < G41S < Y48H.

Disruption of a hydrogen bond network in human versus spider monkey cytochrome c affects heme crevice stability

The hypothesis that the recent rapid evolution of primate cytochromes c, which primarily involves residues in the least stable Ω-loop (Ω-loop C, residues 40-57), stabilizes the heme crevice of cytochrome c relative to other mammals, is tested. To accomplish this goal, we have compared the properties of human and spider monkey cytochrome c and a set of four variants produced in the process of converting human cytochrome c into spider monkey cytochrome c. The global stability of all variants has been measured by guanidine hydrochloride denaturation. The stability of the heme crevice has been assessed with the alkaline conformational transition. Structural insight into the effects of the five amino acid substitutions needed to convert human cytochrome c into spider monkey cytochrome c is provided by a 1.15Å resolution structure of spider monkey cytochrome c. The global stability for all variants is near 9.0kcal/mol at 25°C and pH7, which is higher than that observed for other mammalian cytochromes c. The heme crevice stability is more sensitive to the substitutions required to produce spider monkey cytochrome c with decreases of up to 0.5 units in the apparent pKa of the alkaline conformational transition relative to human cytochrome c. The structure of spider monkey cytochrome c indicates that the Y46F substitution destabilizes the heme crevice by disrupting an extensive hydrogen bond network that connects three surface loops including Ω-loop D (residues 70-85), which contains the Met80 heme ligand.

The response of Ω-loop D dynamics to truncation of trimethyllysine 72 of yeast iso-1-cytochrome c depends on the nature of loop deformation

Trimethyllysine 72 (tmK72) has been suggested to play a role in sterically constraining the heme crevice dynamics of yeast iso-1-cytochrome c mediated by the Ω-loop D cooperative substructure (residues 70-85). A tmK72A mutation causes a gain in peroxidase activity, a function of cytochrome c that is important early in apoptosis. More than one higher energy state is accessible for the Ω-loop D substructure via tier 0 dynamics. Two of these are alkaline conformers mediated by Lys73 and Lys79. In the current work, the effect of the tmK72A mutation on the thermodynamic and kinetic properties of wild-type iso-1-cytochrome c (yWT versus WT*) and on variants carrying a K73H mutation (yWT/K73H versus WT*/K73H) is studied. Whereas the tmK72A mutation confers increased peroxidase activity in wild-type yeast iso-1-cytochrome c and increased dynamics for formation of a previously studied His79-heme alkaline conformer, the tmK72A mutation speeds return of the His73-heme alkaline conformer to the native state through destabilization of the His73-heme alkaline conformer relative to the native conformer. These opposing behaviors demonstrate that the response of the dynamics of a protein substructure to mutation depends on the nature of the perturbation to the substructure. For a protein substructure which mediates more than one function of a protein through multiple non-native structures, a mutation could change the partitioning between these functions. The current results suggest that the tier 0 dynamics of Ω-loop D that mediates peroxidase activity has similarities to the tier 0 dynamics required to form the His79-heme alkaline conformer.

Effect of an Ala81His mutation on the Met80 loop dynamics of iso-1-cytochrome c

An A81H variant of yeast iso-1-cytochrome c is prepared to test the hypothesis that the steric size of the amino acid at sequence position 81 of cytochrome c, which has evolved from Ala in yeast to Ile in mammals, slows the dynamics of the opening of the heme crevice. The A81H mutation is used both to increase steric size and to provide a probe of the dynamics of the heme crevice through measurement of the thermodynamics and kinetics of the His81-mediated alkaline conformational transition of A81H iso-1-cytochrome c. Thermodynamic measurements show that the native conformer is more stable than the His81-heme alkaline conformer for A81H iso-1-cytochrome c. ΔGu°(H2O) is approximately 1.9 kcal/mol for formation of the His81-heme alkaline conformer. By contrast, for K79H iso-1-cytochrome c, the native conformer is less stable than the His79-heme alkaline conformer. ΔGu°(H2O) is approximately -0.34 kcal/mol for formation of the His79-heme alkaline conformer. pH jump and gated electron transfer kinetics demonstrate that this stabilization of the native conformer in A81H iso-1-cytochrome c arises primarily from a decrease in the rate constant for formation of the His81-heme alkaline conformer, kf,His81, relative to kf,His79 for formation of the His79-heme alkaline conformer, which forms by a mechanism similar to that observed for the His81-heme alkaline conformer. The result is discussed in terms of the effect of global protein stability on protein dynamics and in terms of optimization of the sequence of cytochrome c for its role as a peroxidase in the early stages of apoptosis in higher eukaryotes.

Mutation of trimethyllysine 72 to alanine enhances His79-heme-mediated dynamics of iso-1-cytochrome c

Trimethyllysine 72 (Tml72) of yeast iso-1-cytochrome c lies across the surface of the heme crevice loop (Ω-loop D, residues 70-85) like a brace. Lys72 is oriented similarly in horse cytochrome c (Cytc). To determine whether this residue affects the dynamics of opening the heme crevice loop, we have studied the effect of a Tml72 to Ala substitution on the formation of the His79-heme alkaline conformer near neutral pH using a variant of iso-1-Cytc including K72A and K79H mutations. Guanidine hydrochloride denaturation shows that the Tml72 to Ala substitution within error does not affect the global stability of the protein. The effect of the Tml72 to Ala substitution on the thermodynamics of the His79-heme alkaline transition is also small. However, pH-jump kinetic studies of the His79-heme alkaline transition show that both the forward and backward rates of conformational change are increased by the Tml72 to Ala substitution. The barrier for opening the heme crevice is reduced by 0.5 kcal/mol and for closing the heme crevice by 0.3 kcal/mol. The ability of Tml72 to modulate the heme crevice dynamics may indicate a crucial role in regulating function, such as in the peroxidase activity seen in the early stages of apoptosis.

Recombinant expression, biophysical characterization, and cardiolipin-induced changes of two Caenorhabditis elegans cytochrome c proteins

Cytochrome c (cyt c) is one of the most widely studied biomolecules, but not much is known about this protein from nematodes. Recombinant expression of Caenorhabditis elegans CYC-2.1 and CYC-2.2 allowed for detailed characterization of their structural features, redox properties, stabilities, and interactions with cardiolipin (CL)-containing liposomes. Using a variety of spectroscopic tools, we show that CYC-2.1 and CYC-2.2 adopt a globular α-helical fold with His/Met heme ligation. The longer CYC-2.2 has a lower thermodynamic stability than CYC-2.1 and lacks His residues to misligate to the heme in the protein's denatured state. Both C. elegans proteins bind to CL-containing liposomes, and these interactions promote the proteins' peroxidase activity but to a much greater degree for CYC-2.2. Dye-to-heme distance distributions from time-resolved fluorescence resonance energy transfer in bimane-labeled CYC-2.1 and CYC-2.2 revealed similar populations of extended and compact conformers for CL-bound proteins, suggesting that their distinct peroxidase activities in the presence of CL arise from differences in the local heme environments for the two polypeptide ensembles. Without inhibition from His misligation, a less stable and more prone to unfolding CYC-2.2 allows for better access of substrates to the heme and thus exhibits higher peroxidase activity. Similar features of the conformational ensembles of CYC-2.1 and CYC-2.2 to those of mammalian cyt c suggest that C. elegans proteins, particularly the former, could serve as useful models for examining the mechanism of cyt c-CL interactions in live organisms.

Conformational properties of polyglutamine sequences in guanidine hydrochloride solutions

Two sets of iso-1-cytochrome c variants have been prepared with N-terminal insertions of pure polyglutamine, i.e., PolyQ variants, or polyglutamine interrupted with lysine every sixth residue, i.e., Gln-rich variants. The polymer properties of these pure polyGln or Gln-rich sequences have been evaluated using equilibrium and kinetic His-heme loop formation methods for loop sizes ranging from 22 to 46 in 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). In 6.0 M GdnHCl, the scaling exponent, ν(3), for the pure polyGln sequences, is ~1.7--significantly less than ν(3) ≈ 2.15 for the Gln-rich sequences. The stability of the His-heme loops becomes progressively greater for the pure polyGln sequences relative to the Gln-rich sequences as GdnHCl concentration decreases from 6.0 to 1.5 M. Thus, the context of the sequence effects the polymer properties of Gln repeats even in denaturing concentrations of GdnHCl. Comparison of data for the Gln-rich variants with previous results for Gly-rich and Ala-rich variants shows that ν(3) ~ 2.2 for the Gln-rich, Gly-rich, and Ala-rich sequences in 6.0 M GdnHCl, whereas ν(3) remains unchanged at 3.0 M GdnHCl concentration for the Gln-rich and Ala-rich sequences but decreases to ~1.7 for the Gly-rich sequences. Thus, the polymer properties of Gln-rich and Ala-rich sequences are less sensitive to solvent quality in denaturing solutions of GdnHCl than Gly-rich sequences. Evaluation of Flory's characteristic ratio, C(n), for the Gln-rich and Ala-rich sequences relative to the Gly-rich sequences shows that Gln-rich sequences are stiffer than Ala-rich sequences at both 3.0 and 6.0 M GdnHCl.

Origin of the conformational heterogeneity of cardiolipin-bound cytochrome C

Interactions of cytochrome c (cyt c) with cardiolipin (CL) partially unfold the protein, activating its peroxidase function, a critical event in the execution of apoptosis. However, structural features of the altered protein species in the heterogeneous ensemble are difficult to probe with ensemble averaging. Analyses of the dye-to-heme distance distributions P(r) from time-resolved FRET (TR-FRET) have uncovered two distinct types of CL-bound cyt c conformations, extended and compact. We have combined TR-FRET, fluorescence correlation spectroscopy (FCS), and biolayer interferometry to develop a systematic understanding of the functional partitioning between the two conformations. The two subpopulations are in equilibrium with each other, with a submillisecond rate of conformational exchange reflecting the protein folding into a compact non-native state, as well as protein interactions with the lipid surface. Electrostatic interactions with the negatively charged lipid surface that correlate with physiologically relevant changes in CL concentrations strongly affect the kinetics of cyt c binding and conformational exchange. A predominantly peripheral binding mechanism, rather than deep protein insertion into the membrane, provides a rationale for the general denaturing effect of the CL surface and the large-scale protein unfolding. These findings closely relate to cyt c folding dynamics and suggest a general strategy for extending the time window in monitoring the kinetics of folding.

Scaling properties of glycine-rich sequences in guanidine hydrochloride solutions

The intrinsic polymer properties of glycine-rich sequences are evaluated with a set of iso-1-cytochrome c variants with N-terminal inserts of the sequence (GGGGGK)(n) for n = 1-5. The thermodynamics and kinetics of His-heme loop formation are measured as a function of guanidine hydrochloride (GdnHCl) concentration for loop sizes ranging from 22 to 46 residues. The scaling exponent for loop formation, ν(3), evaluated using the Jacobson-Stockmayer equation is near 1.8, at 1.5 and 3.0 M GdnHCl, but it increases to 2.2 in 6.0 M GdnHCl. Previous work on a set of iso-1-cytochrome c variants with (AAAAAK)(n) inserts gave ν(3) = 2.2 for alanine-rich sequences in both 3.0 and 6.0 M GdnHCl. Chain stiffness was evaluated from the relative magnitude of Flory's characteristic ratio, C(n), for alanine-rich versus glycine-rich sequences. In 3.0 M GdnHCl, C(n)(Ala)/C(n)(Gly) is 1.6, decreasing to 1.3 in 6.0 M GdnHCl. The data suggest that solvent-backbone interactions dominate polypeptide conformational properties under good solvent conditions whereas side-chain-dependent properties are more important under poor solvent conditions. The results provide a direct experimental assessment in terms of polymer properties of the distinct roles of Gly versus Ala in the folding code.

Tryptophan stabilizes His-heme loops in the denatured state only when it is near a loop end

We use a host-guest approach to evaluate the effect of Trp guest residues relative to Ala on the kinetics and thermodynamics of formation of His-heme loops in the denatured state of iso-1-cytochrome c at 1.5, 3.0, and 6.0 M guanidine hydrochloride (GdnHCl). Trp guest residues are inserted into an alanine-rich segment placed after a unique His near the N-terminus of iso-1-cytochrome c. Trp guest residues are either 4 or 10 residues from the His end of the 28-residue loop. We find the guest Trp stabilizes the His-heme loop at all GdnHCl concentrations when it is the 4th, but not the 10th, residue from the His end of the loop. Thus, residues near loop ends are most important in developing topological constraints in the denatured state that affect protein folding. In 1.5 M GdnHCl, the loop stabilization is ~0.7 kcal/mol, providing a thermodynamic rationale for the observation that Trp often mediates residual structure in the denatured state. Measurement of loop breakage rate constants, k(b,His), indicates that loop stabilization by the Trp guest residues occurs completely after the transition state for loop formation in 6.0 M GdnHCl. Under poorer solvent conditions, approximately half of the stabilization of the loop develops in the transition state, consistent with contacts in the denatured state being energetically downhill and providing evidence for funneling even near the rim of the folding funnel.

Diverse pathways of oxidative folding of disulfide proteins: underlying causes and folding models

The pathway of oxidative folding of disulfide proteins exhibits a high degree of diversity, which is manifested mainly by distinct structural heterogeneity and diverse rearrangement pathways of folding intermediates. During the past two decades, the scope of this diversity has widened through studies of more than 30 disulfide-rich proteins by various laboratories. A more comprehensive landscape of the mechanism of protein oxidative folding has emerged. This review will cover three themes. (1) Elaboration of the scope of diversity of disulfide folding pathways, including the two opposite extreme models, represented by bovine pancreatic trypsin inhibitor (BPTI) and hirudin. (2) Demonstration of experimental evidence accounting for the underlying mechanism of the folding diversity. (3) Discussion of the convergence between the extreme models of oxidative folding and models of conventional conformational folding (framework model, hydrophobic collapse model).

Manifestations of native topology in the denatured state ensemble of Rhodopseudomonas palustris cytochrome c'

To provide insight into the role of local sequence in the nonrandom coil behavior of the denatured state, we have extended our measurements of histidine-heme loop formation equilibria for cytochrome c' to 6 M guanidine hydrochloride. We observe that there is some reduction in the scatter about the best fit line of loop stability versus loop size data in 6 M versus 3 M guanidine hydrochloride, but the scatter is not eliminated. The scaling exponent, ν(3), of 2.5 ± 0.2 is also similar to that found previously in 3 M guanidine hydrochloride (2.6 ± 0.3). Rates of histidine-heme loop breakage in the denatured state of cytochrome c' show that some histidine-heme loops are significantly more persistent than others at both 3 and 6 M guanidine hydrochloride. Rates of histidine-heme loop formation more closely approximate random coil behavior. This observation indicates that heterogeneity in the denatured state ensemble results mainly from contact persistence. When mapped onto the structure of cytochrome c', the histidine-heme loops with slow breakage rates coincide with chain reversals between helices 1 and 2 and between helices 2 and 3. Molecular dynamics simulations of the unfolding of cytochrome c' at 498 K show that these reverse turns persist in the unfolded state. Thus, these portions of the primary structure of cytochrome c' set up the topology of cytochrome c' in the denatured state, predisposing the protein to fold efficiently to its native structure.

Propensities of aromatic amino acids versus leucine and proline to induce residual structure in the denatured-state ensemble of iso-1-cytochrome c

Histidine-heme loop formation in the denatured state of a protein is a sensitive means for probing residual structure under unfolding conditions. In this study, we use a host-guest approach to investigate the relative tendencies of different amino acids to promote residual structure under denaturing conditions. The host for this work is a 6-amino-acid insert of five alanines, followed by a lysine engineered immediately following a unique histidine near the N-terminus of yeast iso-1-cytochrome c. We substitute the fourth alanine in this sequence HAAAXAK (with X=Trp, Phe, Tyr, and Leu). The effects of proline are tested with substitutions at positions 1 and 5 in the insert (HPAAAAK and HAAAAPK, respectively). Thermodynamic studies on His-heme loop formation in 3 M guanidine hydrochloride reveal significant stabilization of residual structure by aromatic amino acids, particularly Trp and Phe, and minimal stabilization of residual structure by Leu. Prolines slightly disfavor His-heme loop formation, presumably due to enhanced chain stiffness. Kinetic studies reveal that much of the change in His-heme loop stability for the aromatic amino acids is caused by a slowdown in the rate of His-heme loop breakage, indicating that residual structure is preferentially stabilized in the closed-loop form of the denatured state.

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 A) over an extended one (r(H) of 40 A), which forms in their absence. Also, the time constant of a kinetic component (tau(R)) of about 30 micros 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.

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.

Probing early events in ferrous cytochrome c folding with time-resolved natural and magnetic circular dichroism spectroscopies

In a 1998 collaboration with Tony Fink, we coupled nanosecond circular dichroism methods (TRCD) with a CO-photolysis system for quickly triggering folding in cytochrome c (cyt c) in order to make the first time-resolved far-UV CD measurement of early secondary structure formation in a protein. The small signal observed in that initial study, approximately 10% of native helicity, became the seed for increasingly robust results from subsequent studies bringing additional natural and magnetic circular polarization dichroism and optical rotatory dispersion detection methods (e.g., TRORD, TRMCD, and TRMORD), coupled to fast photolysis and photoreduction triggers, to the study of early folding events. Nanosecond polarization methods are reviewed here in the context of the range of initiation methods and structure-sensitive probes currently available for fast folding studies. We also review the impact of experimental results from fast polarization studies on questions in folding dynamics such as the possibility of multiple folding pathways implied by energy landscape models, the sequence dependence of ultrafast helix formation, and the simultaneity of chain collapse and secondary structure formation implicit in molten globule models for kinetic folding intermediates.

Structural heterogeneity of 6 M GdmCl-denatured proteins: implications for the mechanism of protein folding

An in vitro experiment with protein folding is typically initiated with 6 M GdmCl-denatured proteins, which are generally considered fully unfolded. However, studies conducted by various laboratories have shown that many 6 M GdmCl-denatured proteins are structurally heterogeneous and still retain nativelike residual structures. The extent of conformational heterogeneity of the 6 M GdmCl-denatured protein has significant implications for the folding landscape as well as the interpretation of the observed early stage folding mechanism. Using the method of disulfide scrambling, we are able to gain rough insight into the diverse structural properties of 6 M GdmCl-denatured proteins. It demonstrates that most 6 M GdmCl-denatured proteins are approximately fully denatured, but partially unfolded. Most of them comprise diverse conformational isomers. We review here the cumulative evidence obtained from various laboratories and also provide experimental data obtained in our laboratory.

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.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

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.

Electrostatic effects on funneled landscapes and structural diversity in denatured protein ensembles

The denatured state of proteins is heterogeneous and susceptible to general hydrophobic and electrostatic forces, but to what extent does the funneled nature of protein energy landscapes play a role in the unfolded ensemble? We simulate the denatured ensemble of cytochrome c using a series of models. The models pinpoint the efficacy of incorporating energetic funnels toward the native state in contrast with models having no native structure-seeking tendency. These models also contain varying strengths of electrostatic effects and hydrophobic collapse. The simulations based on these models are compared with experimental distributions for the distances between a fluorescent donor and the heme acceptor that were extracted from time-resolved fluorescence energy transfer experiments on cytochrome c. Comparing simulations to detailed experimental data on several labeling sites allows us to quantify the dominant forces in denatured protein ensembles.

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.

Folding mechanism of reduced Cytochrome c: equilibrium and kinetic properties in the presence of carbon monoxide

Despite close structural similarity, the ferric and ferrous forms of cytochrome c differ greatly in terms of their ligand binding properties, stability, folding, and dynamics. The reduced heme iron binds diatomic ligands such as CO only under destabilizing conditions that promote weakening or disruption of native methionine-iron linkage. This makes CO a useful conformational probe for detecting partially structured states that cannot be observed in the absence of endogenous ligands. Heme absorbance, circular dichroism, and NMR were used to characterize the denaturant-induced unfolding equilibrium of ferrocytochrome c in the presence and in the absence of CO. In addition to the native state (N), which does not bind CO, and the unfolded CO complex (U-CO), a structurally distinct CO-bound form (M-CO) accumulates to high levels (approximately 75% of the population) at intermediate guanidine HCl concentrations. Comparison of the unfolding transitions for different conformational probes reveals that M-CO is a compact state containing a native-like helical core and regions of local disorder in the segment containing the native Met80 ligand and adjacent loops. Kinetic measurements of CO binding and dissociation under native, partially denaturing, and fully unfolded conditions indicate that a state M that is structurally analogous to M-CO is populated even in the absence of CO. The binding energy of the CO ligand lowers the free energy of this high-energy state to such an extent that it accumulates even under mildly denaturing equilibrium conditions. The thermodynamic and kinetic parameters obtained in this study provide a fully self-consistent description of the linked unfolding/CO binding equilibria of reduced cytochrome c.

Probing the bottom of a folding funnel using conformationally gated electron transfer reactions

The effect of global stability on the kinetics of interconversion between the native (N) and a compact, partially unfolded form (I) of iso-1-cytochrome c stabilized by His73-heme ligation is investigated using a novel conformationally gated ET method. For the K73H variant and the 2-fold less stable AcH73 variant, the N and I conformers are of nearly equal stability at pH 7.5. The pH jump kinetic data yield kobs = kNI + kIN of 35-40 s-1 at final pH values from 6 to 8 for the AcH73 variant, about 3-fold faster than for the more stable K73H variant. Gated ET measurements give kNI = 28 s-1 and kIN = 13 s-1 for the AcH73 variant, 10- and 2-fold greater than that for the more stable K73H variant. Thus, funneled landscapes have evolved such that loss of global stability lowers barriers at the bottom of a folding funnel, still allowing for efficient folding.

Non-native heme-histidine ligation promotes microsecond time scale secondary structure formation in reduced horse heart cytochrome c

Previous far-UV time-resolved optical rotatory dispersion (TRORD) studies of the sub-millisecond (burst) phase of secondary structure formation in horse and tuna cytochromes c after photoreduction in denaturant suggested that the non-native His18-Fe-His33 heme ligation dominant in the unfolded horse protein facilitated this fast folding better than did the His18-Fe-His26 coordination dominant in tuna [Chen, E., Goldbeck, R.A., and Kliger, D.S. (2003) J. Phys. Chem. A 107, 8149-8155; Chen, E., Goldbeck, R.A., and Kliger, D.S. (2004) J. Am. Chem. Soc. 126, 11175-11181]. Whether His18-Fe-His33 coordination actually facilitates fast secondary structure formation or just slows folding less than His18-Fe-His26 coordination is probed by examining the double histidine mutant H26QH33N of horse heart cytochrome c. The fast folding phase is absent in H26QH33N, indicating that His18-Fe-His33 misligation does promote fast secondary structure formation, as does His18-Fe-His26 to a lesser extent. His33 may be better able to facilitate folding because it is not as constrained by hydrogen bonding interactions in the denatured state as is His26. A greater flexibility, not only because of weakened or disrupted Van der Waals interactions in the presence of guanidine hydrochloride (GuHCl) but also because of its position relative to His18, may allow His33 to ligate to the heme group more easily than His26. These results are discussed along with the results of far-UV CD and Soret and visible region MCD measurements, which were performed to probe heme ligation in H26QH33N and to understand how GuHCl affects its folding stability and cooperativity.

Alkaline conformational transition and gated electron transfer with a Lys 79 --> his variant of iso-1-cytochrome c

To probe the mechanism of the alkaline conformational transition and its effect on the dynamics of gated electron transfer (ET) reactions, a Lys 79 --> His (K79H) variant of iso-1-cytochrome c has been prepared. Guanidine hydrochloride denaturation monitored by circular dichroism and absorbance at 695 nm indicates that this variant unfolds from a partially unfolded state. The conformation of the wild type (WT) and K79H proteins was monitored at 695 nm from pH 2 to 11. These data indicate that acid unfolding is multi-state for both K79H and WT proteins and that the His 79-heme alkaline conformer is more stable than a previously reported His 73-heme alkaline conformer. Fast and slow phases are observed in the kinetics of the alkaline transition of the K79H variant. The pH dependence of the fast phase kinetic data shows that ionizable groups with pKa values near 6.8 and 9 modulate the formation of the His 79-heme alkaline conformer. The slow phase kinetic data are consistent with a single ionizable group with a pKa near 9.5 promoting the Lys 73-heme alkaline transition. In the broader context of data on the alkaline transition, ionization of the ligand replacing Met 80 appears to play a primary role in promoting the formation of the alkaline conformer, with other ionizable groups acting as secondary modulators. Intermolecular ET with hexaammineruthenium(II) chloride shows conformational gating due to both His 79-heme and Lys 73-heme alkaline conformers. Both the position and the nature of the alkaline state ligand modulate the dynamics of ET gating.

Probing the cytochrome c' folding landscape

The folding kinetics of R. palustris cytochrome c' (cyt c') have been monitored by heme absorption and native Trp72 fluorescence at pH 5. The Trp72 fluorescence burst signal suggests early compaction of the polypeptide ensemble. Analysis of heme transient absorption spectra reveals deviations from two-state behavior, including a prominent slow phase that is accelerated by the prolyl isomerase cyclophilin. A nonnative proline configuration (Pro21) likely interferes with the formation of the helical bundle surrounding the heme.

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.

Role of the Charge–Charge Interactions in Defining Stability and Halophilicity of the CspB Proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.

Thermodynamics of protein denatured states

Recent work on the thermodynamics of protein denatured states is providing insight into the stability of residual structure and the conformational constraints that affect the disordered states of proteins. Current data from native state hydrogen exchange and the pH dependence of protein stability indicate that residual structure can modulate the stability of the denatured state by up to 4 kcal mol(-1). NMR structural data have emphasized the role of hydrophobic clusters in stabilizing denatured state residual structures, however recent results indicate that electrostatic interactions, both favorable and unfavorable, are also important modulators of the stability of the denatured state. Thermodynamics methods that take advantage of histidine-heme ligation chemistry have also been developed to probe the conformational constraints that act on denatured states. These methods have provided insights into the role of excluded volume, chain stiffness, and loop persistence in modulating the conformational preferences of highly disordered proteins. New insights into protein folding and novel methods to manipulate protein stability are emerging from this work.

Mutation of asparagine 52 to glycine promotes the alkaline form of iso-1-cytochrome c and causes loss of cooperativity in acid unfolding

The kinetics and thermodynamics of the alkaline and acid conformational transitions of a Lys 79 --> Ala/Asn 52 --> Gly (A79G52) variant of iso-1-cytochrome c are studied. The Lys 79 --> Ala mutation is designed to limit heme ligation in the alkaline conformer to Lys 73. The Asn 52 --> Gly mutation is intended to shift the population of the alkaline conformer to physiological pH based on the hierarchical nature of the cooperative substructures of this protein. The midpoint pH for formation of the alkaline conformer is approximately 7.45. The kinetics for the alkaline conformational transition of the A79G52 variant are consistent with the ionization constant, pK(H), for the trigger group controlling formation of the alkaline conformer being approximately 9.5. This pK(H) is low for alkaline conformers involving lysine-heme ligation but is consistent with the pK(a) of the highest of three ionizable groups which modulate formation of the histidine-heme alkaline conformer of a His 73 variant of iso-1-cytochrome c [Martinez, R. E., and Bowler, B. E. (2004) J. Am. Chem. Soc. 126, 6751-6758]. The acid transition of the A79G52 variant is split into two phases. Both the Lys 79 --> Ala and Asn 52 --> Gly mutations are expected to affect the buried hydrogen bond network of cytochrome c, suggesting that this network is an important modulator of the acid unfolding of cytochrome c.

Thermodynamics and kinetics of formation of the alkaline state of a Lys 79-->Ala/Lys 73-->His variant of iso-1-cytochrome c

The alkaline transition kinetics of a Lys 73-->His (H73) variant of iso-1-cytochrome c are triggered by three ionizable groups [Martinez, R. E., and Bowler, B. E. (2004) J. Am. Chem. Soc. 126, 6751-6758]. To eliminate ambiguities caused by overlapping phases due to formation of the Lys 79 alkaline conformer and proline isomerization associated with the His 73 alkaline conformer, we mutated Lys 79 to Ala in the H73 variant (A79H73). The stability and guanidineHCl m-values of the A79H73 and H73 variants at pH 7.5 are the same. The Ala 79 mutation causes formation of the alkaline conformer to depend on [NaCl]. The salt dependence saturates at 500 mM NaCl, and the thermodynamics of alkaline state formation for the A79H73 and H73 variants become identical. The salt dependence is consistent with loss of an electrostatic contact between Lys 79 and heme propionate D in the A79H73 variant. The kinetics of alkaline state formation for the A79H73 variant support the three trigger group model developed for the H73 variant, with the primary trigger, pK(HL), being ionization of His 73. The low pH ionization, pK(H1), is perturbed by the Ala 79 mutation indicating that this ionization is modulated by the buried hydrogen bond network involving heme propionate D. The A79H73 variant has a high spin heme above pH 9 suggesting that the high pH ionization, pK(H2), involves a high spin heme conformer. The proline isomerization phase is modulated by both pK(HL) and pK(H2) indicating that it is sensitive to protein conformation.

Kinetics of Loop Formation and Breakage in the Denatured State of Iso-1-cytochrome c

The earliest events in protein folding involve the formation of simple loops. Observing the rates of loop closure under denaturing conditions can provide direct insight into the relative probability and sequence determinants for formation of loops of different sizes. The persistence of these initial contacts is equally important for efficient folding, so measurement of rates of loop breakage under denaturing conditions is also essential. We have used stopped-flow and continuous-flow methods to measure the rates of histidine-heme loop formation and breakage in the denatured state of iso-1-cytochrome c (in the presence of 3 M guanidine HCl). The data indicate that the mechanism for forming loops is a two-step process, the first step being the deprotonation of the histidine, and the second step being the binding of the histidine to the heme. This mechanism makes it possible to extract both the rate constants of formation, k(f), and breakage, k(b), of loops from the pH dependence of the observed rate constant, k(obs). To determine the dependence of k(f) and k(b) on loop size, we have carried out kinetic measurements for seven single surface histidine variants of iso-1-cytochrome c. A scaling factor (the dependence of k(f) on log[loop size]) of approximately -1.8 is observed for loop formation, similar to that observed in other systems. The magnitude of k(b) varies from 30 s(-1) to 300 s(-1), indicating that the stability of different loops varies considerably. The implications of the kinetics of loop formation and breakage in the denatured state for the mechanism of protein folding are discussed.

Role of Hydrogen Bond Networks and Dynamics in Positive and Negative Cooperative Stabilization of a Protein

Cooperativity mediated through hydrogen bond networks in yeast iso-1-cytochrome c was studied using a thermodynamic triple mutant cycle. Three known stabilizing mutations, Asn 26 to His, Asn 52 to Ile, and Tyr 67 to Phe, were used to construct the triple mutant cycle. The side chain of His 26, a wild-type residue, forms two hydrogen bonds that bridge two substructures of the wild-type protein, and Tyr 67 and Asn 52 are part of an extensive buried hydrogen bond network. The stabilities of all variants in the triple mutant cycle were determined by guanidine hydrochloride denaturation methods and used to determine the pairwise, Delta(2)G(int), and triple interaction energies. His 26 and Ile 52 interact cooperatively (Delta(2)G(int) is 1-2 kcal/mol), whereas the two other pairs of mutations interact anticooperatively (Delta(2)G(int) is -0.5 to -1.5 kcal/mol). Previously reported structural data for iso-1-cytochrome c variants containing these mutations show that changes in the strength of the His 26 to Glu 44 hydrogen bond, apparently caused by changes in main chain dynamics, provide a mechanism for the long distance (His 26 to Phe 67 and His 26 to Ile 52) propagation of pairwise interaction energies. Opposing changes in the strength of the His 26 to Glu 44 hydrogen bond caused by the N52I and Y67F mutations generate a negative triple interaction energy (-0.9 +/-0.7 kcal/mol) that combined with cancellation of cooperative and anticooperative pairwise interactions produce apparent additivity for the stabilizing effects of the single mutations in the triple mutant variant.

Many faces of the unfolded state: conformational heterogeneity in denatured yeast cytochrome C

We have measured fluorescence energy-transfer (FET) kinetics from a dansyl fluorophore (Dns) introduced by derivatization of a Cys side-chain to the Fe(III) heme covalently attached to unfolded yeast iso-1 cytochrome c (cyt). To gain a global picture of the unfolded state, we examined variants with the fluorophore attached on three different helices (K4C, E66C, K99C) and in three different loops (H39C, D50C, L85C). Analysis of the FET kinetics data gave distributions of distances between the fluorescent donor and acceptor; these distributions demonstrate that the guanidine hydrochloride (GuHCl)-denatured polypeptide ensemble is not a simple random coil. Although misligation imposes some constraints, it is not the only source of structural complexity in the unfolded protein. Our FET kinetics data reveal a high degree of heterogeneity in the unfolded ensemble of cytochrome c. We detect relatively large populations of compact structures in unfolded Dns(C50)cyt, Dns(C39)cyt, and Dns(C66)cyt. These structures likely play a role in forming a hydrophobic core during the folding process.

Insertion of the cytochrome b5 heme-binding loop into an SH3 domain. Effects on structure and stability, and clues about the cytochrome's architecture

Under native conditions, apocytochrome b(5) exhibits a stable core and a disordered heme-binding region that refolds upon association with the cofactor. The termini of this flexible region are in close proximity, suggesting that loop closure may contribute to the thermodynamic properties of the apocytochrome. A chimeric protein containing 43 residues encompassing the cytochrome loop was constructed using the cyanobacterial photosystem I accessory protein E (PsaE) from Synechococcus sp. PCC 7002 as a structured scaffold. PsaE has the topology of an SH3 domain, and the insertion was engineered to replace its 14-residue CD loop. NMR and optical spectroscopies showed that the hybrid protein (named EbE1) was folded under native conditions and that it retained the characteristics of an SH3 domain. NMR spectroscopy revealed that structural and dynamic differences were confined near the site of loop insertion. Variable-temperature 1D NMR spectra of EbE1 confirmed the presence of a kinetic unfolding barrier. Thermal and chemical denaturations of PsaE and EbE1 demonstrated cooperative, two-state transitions; the stability of the PsaE scaffold was found only moderately compromised by the insertion, with a DeltaT(m) of 8.3 degrees C, a DeltaC(m) of 1.5 M urea, and a DeltaDeltaG degrees of 4.2 kJ/mole. The data implied that the penalty for constraining the ends of the inserted region was lower than the approximately 6.4 kJ/mole calculated for a self-avoiding chain. Extrapolation of these results to cytochrome b(5) suggested that the intrinsic stability of the folded portion of the apoprotein reflected only a small detrimental contribution from the large heme-binding domain.

The role of the unfolded state in hairpin stability

The effects of a T3S mutation on the stability of a 3:5 beta-hairpin forming peptide (YITNSNGTWT) are investigated. Molecular dynamics simulations in explicit water indicate that the wild-type peptide forms a stable hairpin whereas the T3S mutant does not, in agreement with the experimental data. Thermodynamic integration calculations for the mutation of Thr to Ser suggest that the free-energy changes in the folded state are small, but the corresponding changes in the unfolded state are large and favorable. One of the main reasons for the difference appears to be the formation of a stable cluster involving the Tyr1 and Ser3 hydroxyl groups and their interaction with the C-terminal carboxylate group, which was observed after unfolding of the T3S mutant. Further analysis of the side-chain preferences of Thr and Ser indicate that the corresponding cluster in the wild-type peptide is unstable due to the high preference of the Thr chi1 dihedral for g+ states, which appeared to be incompatible with formation of a stable cluster. The results suggest that one should consider the nature of the unfolded state when attempting to fully explain the effects of mutations on hairpin stability.