Measuring denatured state energetics: deviations from random coil behavior and implications for the folding of iso-1-cytochrome c

Effect of proline content and histidine ligation on the dynamics of Ω-loop D and the peroxidase activity of iso-1-cytochrome c

To study how proline residues affect the dynamics of Ω-loop D (residues 70 to 85) of cytochrome c, we prepared G83P and G83A variants of yeast iso-1-cytochrome c (iso-1-Cytc) in the presence and absence of a K73H mutation. Ω-loop D is important in controlling both the electron transfer function of Cytc and the peroxidase activity of Cytc used in apoptosis because it provides the Met80 heme ligand. The G83P and G83A mutations have no effect on the global stability of iso-1-Cytc in presence or absence of the K73H mutation. However, both mutations destabilize the His73-mediated alkaline conformer relative to the native state. pH jump stopped-flow experiments show that the dynamics of the His73-mediated alkaline transition are significantly enhanced by the G83P mutation. Gated electron transfer studies show that the enhanced dynamics result from an increased rate of return to the native state, whereas the rate of loss of Met80 ligation is unchanged by the G83P mutation. Thus, the G83P substitution does not stiffen the conformation of the native state. Because bis-His heme ligation occurs when Cytc binds to cardiolipin-containing membranes, we studied the effect of His73 ligation on the peroxidase activity of Cytc, which acts as an early signal in apoptosis by causing oxygenation of cardiolipin. We find that the His73 alkaline conformer suppresses the peroxidase activity of Cytc. Thus, the bis-His ligated state of Cytc formed upon binding to cardiolipin is a negative effector for the peroxidase activity of Cytc early in apoptosis.

Influence of the Interdomain Interface on Structural and Redox Properties of Multiheme Proteins

Multiheme proteins are important in energy conversion and biogeochemical cycles of nitrogen and sulfur. A diheme cytochrome c4 (c4) was used as a model to elucidate roles of the interdomain interface on properties of iron centers in its hemes A and B. Isolated monoheme domains c4-A and c4-B, together with the full-length diheme c4 and its Met-to-His ligand variants, were characterized by a variety of spectroscopic and stability measurements. In both isolated domains, the heme iron is Met/His-ligated at pH 5.0, as in the full-length c4, but becomes His/His-ligated in c4-B at higher pH. Intradomain contacts in c4-A are minimally affected by the separation of c4-A and c4-B domains, and isolated c4-A is folded. In contrast, the isolated c4-B is partially unfolded, and the interface with c4-A guides folding of this domain. The c4-A and c4-B domains have the propensity to interact even without the polypeptide linker. Thermodynamic cycles have revealed properties of monomeric folded isolated domains, suggesting that ferrous (FeII), but not ferric (FeIII) c4-A and c4-B, is stabilized by the interface. This study illustrates the effects of the interface on tuning structural and redox properties of multiheme proteins and enriches our understanding of redox-dependent complexation.

Humanlike substitutions to Ω-loop D of yeast iso-1-cytochrome c only modestly affect dynamics and peroxidase activity

Structural studies of yeast iso-1-cytochrome c (L.J. McClelland, T.-C. Mou, M.E. Jeakins-Cooley, S.R. Sprang, B.E. Bowler, Proc. Natl. Acad. Sci. U.S.A. 111 (2014) 6648-6653) show that modest movement of Ω-loop D (residues 70-85, average RMSD versus the native structure: 0.81 Å) permits loss of Met80-heme ligation creating an available coordination site to catalyze the peroxidase activity mediated by cytochrome c early in apoptosis. However, Ala81 and Gly83 move significantly (RMSDs of 2.18 and 1.26 Å, respectively). Ala81 and Gly83 evolve to Ile and Val, respectively, in human cytochrome c and peroxidase activity decreases 25-fold relative to the yeast protein at pH 7. To test the hypothesis that these residues evolved to restrict the peroxidase activity of cytochrome c, A81I and G83V variants of yeast iso-1-cytochrome c were prepared. For both variants, the apparent pK_a of the alkaline transition increases by 0.2 to 0.3 relative to the wild type (WT) protein and the rate of opening the heme crevice is slowed. The cooperativity of acid unfolding is decreased for the G83V variant. At pH 7 and 8, the catalytic rate constant, k_cat, for the peroxidase activity of both variants decreases relative to WT, consistent with the effects on alkaline isomerization. Below pH 7, the loss in the cooperativity of acid unfolding causes k_cat for peroxidase activity to increase for the G83V variant relative to WT. Neither variant decreases k_cat to the level of the human protein, indicating that other residues also contribute to the low peroxidase activity of human cytochrome c.

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.

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.

Redox-dependent stability, protonation, and reactivity of cysteine-bound heme proteins

Cysteine-bound hemes are key components of many enzymes and biological sensors. Protonation (deprotonation) of the Cys ligand often accompanies redox transformations of these centers. To characterize these phenomena, we have engineered a series of Thr78Cys/Lys79Gly/Met80X mutants of yeast cytochrome c (cyt c) in which Cys78 becomes one of the axial ligands to the heme. At neutral pH, the protonation state of the coordinated Cys differs for the ferric and ferrous heme species, with Cys binding as a thiolate and a thiol, respectively. Analysis of redox-dependent stability and alkaline transitions of these model proteins, as well as comparisons to Cys binding studies with the minimalist heme peptide microperoxidase-8, demonstrate that the protein scaffold and solvent interactions play important roles in stabilizing a particular Cys-heme coordination. The increased stability of ferric thiolate compared with ferrous thiol arises mainly from entropic factors. This robust cyt c model system provides access to all four forms of Cys-bound heme, including the ferric thiol. Protein motions control the rates of heme redox reactions, and these effects are amplified at low pH, where the proteins are less stable. Thermodynamic signatures and redox reactivity of the model Cys-bound hemes highlight the critical role of the protein scaffold and its dynamics in modulating redox-linked transitions between thiols and thiolates.

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.

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.

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.

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.

pH-dependent surface-enhanced resonance Raman scattering of yeast iso-1-cytochrome c adsorbed on silver nanoparticle surfaces under denaturing conditions at pH<3

We measured the pH-induced spectral changes of yeast iso-1-cytochrome c on silver nanoparticle surfaces using surface-enhanced resonance Raman scattering (SERRS) at 457.9 nm. At a pH of 3, the Met80 ligand in yeast iso-1-cytochrome c is assumed to dissociate, leading to a marked conformational change as evidenced by the vibrational spectral shifts. The Soret band at 410 nm in the UV-Vis spectrum shifted to 396 nm at pH 3, indicating a transition from a low spin state to a high spin state from a weak interaction with a water molecule. Thus, SERRS spectroscopy can measure the pH-induced denaturalization of cyt c adsorbed on metal nanoparticle surfaces at a lower concentration with a better sensitivity than ordinary resonance Raman spectroscopy.

Stability of the heme Fe-N-terminal amino group coordination bond in denatured cytochrome c

In the denatured states of Hydrogenobacter thermophilus cytochrome c(552) (HT) and Pseudomonas aeruginosa cytochrome c(551) (PA), and their mutants, the N-terminal amino group of the polypeptide chain is coordinated to heme Fe in place of the axial Met, the His-N(term) form being formed. The coordination of the N-terminal amino group to heme Fe leads to loop formation by the N-terminal stretch preceding the first Cys residue bound to the heme, and the N-terminal stretches of HT and PA are different from each other in terms of both the sequence and the number of constituent amino acid residues. The His-N(term) form was shown to be rather stable, and hence it can influence the stability of the denatured state. We have investigated the heme Fe coordination structures and stabilities of the His-N(term) forms emerging upon guanidine hydrochloric acid-induced unfolding of the oxidized forms of the proteins. The Fe-N(term) coordination bond in the His-N(term) form with a 9-residue N-terminal stretch of HT proteins was found to be tilted to some extent away from the heme normal, as reflected by the great heme methyl proton shift spread. On the other hand, the small heme methyl proton shift spread of the His-N(term) form with an 11-residue stretch of PA proteins indicated that its Fe-N(term) bond is nearly parallel with the heme normal. The stability of the His-N(term) form was found to be affected by the structural properties of the N-terminal stretch, such as its length and the N-terminal residue. With a given N-terminal residue, the stability of the His-N(term) form is higher for a 9-residue N-terminal stretch than an 11-residue one. In addition, with a given length of the N-terminal stretch, the His-N(term) form with an N-terminal Glu is stabilized by a few kJ mol(-1) relative to that with an N-terminal Asn. These results provide a novel insight into the stabilizing interactions in the denatured cyts c that will facilitate elucidation of the folding/unfolding mechanisms of the proteins.

Zinc porphyrin: a fluorescent acceptor in studies of Zn-cytochrome c unfolding by fluorescence resonance energy transfer

FRET between the zinc porphyrin (ZnP) chromophore in zinc-substituted cytochrome c (Zn-cyt c) and an Alexa Fluor dye attached to specific surface sites was used to characterize Zn-cyt c unfolding. The use of ZnP as a fluorescent acceptor eliminates the need to doubly label the protein with exogenous dyes to perform FRET experiments in which both donor and acceptor fluorescence is monitored. The requirement for attachment of only one dye also minimizes perturbation to the protein. This sensitive technique allowed for the determination of distances between the label placed at six different sites and ZnP through a range of denaturant concentrations. Fitting of the data to a three-state model provides distances in the unfolding intermediate. The use of ZnP as a fluorescent acceptor of energy in FRET has a significant potential for application to a range of other systems including heme-binding proteins and proteins to which a covalently attached heme tag may be added.

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.

Characterization of N-terminal amino group–heme ligation emerging upon guanidine hydrochloric acid induced unfolding of Hydrogenobacter thermophilus ferricytochrome c 552

Nonnative heme coordination structures emerging upon guanidine hydrochloric acid (GdnHCl) induced unfolding of Hydrogenobacter thermophilus ferricytochrome c552 were characterized by means of paramagnetic NMR. The heme coordination structure possessing the N-terminal amino group of the peptide chain in place of axial Met (His-Nterm form) was determined in the presence of GdnHCl concentrations in excess of 1.5 M at neutral pH. The stability of the His-Nterm form at pH 7.0 was found to be comparable with that of the bis-His form which has been recognized as a major nonnative heme coordination structure in cytochrome c folding/unfolding. Consequently, in addition to the bis-His form, the His-Nterm form is a substantial intermediate which affects the pathway and kinetics of the folding/unfolding of cytochromes c, of which the N-terminal amino groups are not acetylated.

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.

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.

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.

Site-specific collapse dynamics guide the formation of the cytochrome c' four-helix bundle

The evolution of tryptophan-to-heme (W/heme) distance distributions extracted from analysis of fluorescence energy transfer kinetics during the refolding of Rhodopseudomonas palustris cytochrome c' reveals dramatic differences between two variants [W32 (Q1A/F32W/W72F) and W72 (Q1A)]. Both W32/heme and W72/heme distance distributions measured at the earliest time point attainable with a continuous-flow mixer (150 mus) confirm that the polypeptide ensemble is not uniformly collapsed and that native structure is not formed. Time-resolved fluorescence spectra indicate that W32 is sequestered from the aqueous solution during the first 700 mus of folding, whereas W72 remains exposed to solvent. The first moment of the W32/heme distance distribution evolves to its native value faster than that of W72, suggesting that the approach of W32 to the heme precedes that of W72.

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.

Gold Nanoparticles as a Colorimetric Sensor for Protein Conformational Changes

Spherical gold nanoparticles and flat gold films are prepared in which yeast iso-1-cytochrome c (Cyt c) is covalently bound to the gold surface by a thiol group in the cystein 102 residue. Upon exposure to solutions of different pH, bound Cyt c unfolds at low pH and refolds at high pH. This conformational change causes measurable shifts in the color of the coated nanoparticle solutions detected by UV-VIS absorption spectroscopy and in the refractive index (RI) of the flat gold films detected by surface plasmon resonance (SPR) spectroscopy. Both experiments demonstrate the same trend with pH, suggesting the use of protein-covered gold nanoparticles as a simple colorimetric sensor for conformational change.

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.

Communication of Stabilizing Energy between Substructures of a Protein

Thermodynamic communication between protein substructures has been investigated by determining the stabilizing effect of mutations at position 52 in the least stable, N-yellow, substructure of cytochrome c on the second least stable, Red, and most stable, Blue, substructures of the protein. A Lys 73 --> His (H73) variant of iso-1-cytochrome c, containing these mutations was used to measure the stability of the Red substructure of cytochrome c through the pH and guanidine hydrochloride (gdnHCl) dependence of the His 73-mediated alkaline conformational transition. The stability of the Blue substructure was measured by global unfolding with gdnHCl and increased by 1 to 3.5 kcal/mol versus the H73 variant. The data demonstrate that the increase in stability of the Red substructure is similar to the increase in global stability, consistent with upward propagation of stabilizing energy from less (N-yellow) to more stable (Red and Blue) protein substructures. The result also supports sequential rather than independent unfolding of the N-yellow and Red substructures of cytochrome c. The data indicate that a leucine at position 52 alters the nature of partial unfolding of the Red substructure, a surprising effect for a single-site mutation. For all variants, the thermodynamics of formation of the Lys 79 alkaline state, which does not unfold the entire Red substructure, shows less stabilization of the portion of the protein unfolded relative to the stabilization of the Blue substructure, indicating that propagation of energy between substructures is somewhat disrupted when unfolding does not correspond to a natural substructure.

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.

Conformational Properties of the Iso-1-Cytochrome c Denatured State: Dependence on Guanidine Hydrochloride Concentration

Production of seven single surface histidine variants of yeast iso-1-cytochrome c allowed measurement of the apparent pK(a), pK(a)(obs), for histidine-heme loop formation for loops of nine to 83 amino acid residues under varying denaturing conditions (2 M to 6 M guanidine hydrochloride, gdnHCl). A linear correlation between pK(a)(obs) and the log of the loop size is expected for a random coil, pK(a)(obs) proportional to k log(n), where k is a scaling factor and n is the number of monomers in the loop. For small loops of nine, 16, and 22 monomers, no dependence of pK(a)(obs) on loop size was observed at any denaturant concentration indicating effects from chain stiffness. For larger loops of 37, 56, 72, and 83 monomers, the dependence of pK(a)(obs) on log(n) was linear and the slope of that dependence decreased with increasing concentration of denaturant. The scaling factor obtained at 5 M and 6 M gdnHCl for the larger loop sizes was approximately -2.0, close to the value of -2.2 expected for a random coil with excluded volume. However, scaling factors obtained under less harsh denaturing conditions (2 M to 4.5 M gdnHCl) deviated strongly from that expected for a random coil, being in the range -3 to -4. The gdnHCl dependence of pK(a)(obs) at each loop size was also evaluated to obtain denaturant m-values. Short loops where chain stiffness dominates had similar m-values of approximately 0.25 kcal/mol M. For larger loops m-values decrease with increasing loop size indicating that less hydrophobic area is sequestered when larger loops form. It is known that the earliest events in protein folding involve the formation of simple loops. The data from these studies provide direct insight into the relative probability with which loops of different sizes will form, as well as the factors which affect loop formation.

The protein-folding speed limit: intrachain diffusion times set by electron-transfer rates in denatured Ru(NH3)5(His-33)-Zn-cytochrome c

The kinetics of electron transfer from the triplet-excited Zn-porphyrin to a Ru(NH(3))(5)(His-33)(3+) complex have been measured in Zn-substituted ruthenium-modified cytochrome c under denaturing conditions. In the folded protein, the electron-tunneling rate constant is 7.5 x 10(5) s(-1). As the protein is denatured with guanidine hydrochloride, a faster adiabatic electron-transfer reaction appears (4.0 x 10(6) s(-1), [guanidine hydrochloride] = 5.4 M) that is limited by the rate of intrachain diffusion to bring the Zn-porphyrin and Ru complex into contact. The 250-ns contact time for formation of a 15-residue loop in denatured cytochrome c is in accord with a statistical model developed by Camacho and Thirumalai [Camacho, C. J. & Thirumalai, D. (1995) Proc. Natl. Acad. Sci. USA 92, 1277-1281] that predicts that the most probable transient loops formed in denatured proteins are comprised of 10 amino acids. Extrapolation of the cytochrome c contact time to a 10-residue loop sets the folding speed limit at approximately 10(7) s(-1).

Effect of pH on the iso-1-cytochrome c denatured state: changing constraints due to heme ligation

The effect of pH on the denatured state (3 M guanidine hydrochloride) was evaluated with fluorescence spectroscopy for four variants of iso-1-cytochrome c, AcTM (no surface histidines), AcH26 (surface histidine at position 26), AcH54 (surface histidine at position 54), and AcH54I52 (stabilizing I52 mutation added to AcH54). Changes in the compactness and the heme ligation of the denatured state, as a function of pH, were monitored through changes in Trp 59-heme fluorescence quenching. With the AcTM and AcH26 variants, no change in the fluorescence intensity occurs from pH 4 to 10. However, for the AcH54 and AcH54I52 variants the fluorescence intensity drops significantly between pH 4 and 6, consistent with His 54 binding to the heme of cytochrome c. Between pH 8 and 10 fluorescence intensity increases again, indicating that the His 54 is displaced from the heme. The data are consistent with lysines 4 and 5 being the primary heme ligands at alkaline pH, under denaturing conditions. This conclusion was confirmed by site-directed mutagenesis. Thermodynamic analysis indicates that heme-ligand affinity in the denatured state is controlled primarily by sequence position (loop size) and that when histidines are present they inhibit lysine ligation until approximately pH 8.5-9.0 as compared to pH 7.5 with the AcTM variant. Thus, at physiological pH, histidine ligands provide the primary constraint on the denatured state of cytochrome c. The heme-Trp 59 distance in the denatured state of iso-1-cytochrome c, derived from analysis by Förster energy transfer theory, is approximately 26 A at pH 4 and 10, much shorter than the random coil prediction of 56 A. Surprisingly, the heme-Trp 59 distance in the His 54 bound conformation only drops to approximately 21 A, consistent with an extended conformation for the short polypeptide segment separating heme and Trp 59.

Peroxidase activity as a tool for studying the folding of c-type cytochromes

The peroxidase activity of c-type cytochromes increases substantially by unfolding. This phenomenon was used to study the equilibrium unfolding of ferricytochrome c. The peroxidase activity is already enhanced at low denaturant concentrations. The lowest free energy folding intermediate is easily detected by this method, while it is invisible using fluorescence or optical spectroscopy. The free energy difference between this folding intermediate and the native state depends on the strength of the sixth ligand of the heme-iron and the increase in peroxidase activity upon unfolding is shown to be a sensitive indicator of the strength of this ligand. Under fully denaturing conditions, the peroxidase activity is inhibited by protein-based ligands. It is shown that at least three different ligand groups can be responsible for this inhibition, and that at neutral or alkaline pH, the predominant ligand is not histidine. The use of peroxidase activity assays as a method to study the unfolding of cytochrome c is evaluated.

Effects of topology and excluded volume on protein denatured state conformational properties

The conformational constraints on protein denatured states are of prime importance in modulating early events in protein folding. Although structural studies have demonstrated residual structure in protein denatured states, much remains poorly understood with regard to the conformational properties of this state. Here, we investigate topological effects on loop formation probabilities in denatured iso-1-cytochrome c by comparing histidine-heme binding affinities for histidines on the N- versus the C-terminal side of the heme. For histidines N-terminal to the heme (preceding cysteine 14), the polypeptide emerges from the edge of the heme and must simply fold over to bind to the heme. For histidines C-terminal to the heme (following histidine 18), the polypeptide emerges from the back side of the heme and must wrap around the heme for the histidine to bind to the heme. Thus, the steric constraints on this wrap-around topology are expected to be much more demanding than for the heme-edge topology of the N-terminal histidines. Evaluation of loop formation probabilities in 3 M guanidine hydrochloride, conditions that fully denature the variants studied, demonstrates that N-terminal histidine-heme loop formation is 10-25-fold more favorable than C-terminal histidine-heme loop formation, for similar loop sizes. A two-dimensional square lattice model indicates that excluded volume is important in this topological preference. These data provide direct evidence that denatured state topology affects contact probability, and thus probable folding pathways, in a disordered protein.

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

A set of nine variants of yeast iso-1-cytochrome c with zero or one surface histidine have been engineered such that the N-terminal amino group is acetylated in vivo. N-terminal acetylation has been confirmed by mass spectral analysis of intact and proteolytically digested protein. The histidine-heme loop-forming equilibrium, under denaturing conditions (3 M guanidine hydrochloride), has been measured by pH titration providing an observed pK(a), pK(a)(obs), for each variant. N-terminal acetylation prevents the N-terminal amino group-heme binding equilibrium from interfering with measurements of histidine-heme affinity. Significant deviation is observed from the linear dependence of pK(a)(obs) on the log of the number of monomers in the loop formed, expected for a random coil denatured state. The maximum histidine-heme affinity occurs for a loop size of 37 monomers. For loop sizes of 37-83 monomers, histidine-heme pK(a)(obs) values are consistent with a scaling factor of -4.2+/-0.3. This value is much larger than the scaling factor of -1.5 for a freely jointed random coil, which is commonly used to represent the conformational properties of protein denatured states. For loop sizes of nine to 22 monomers, chain stiffness is likely responsible for the decreases in histidine-heme affinity relative to a loop size of 37. The results are discussed in terms of residual structure and sequence composition effects on the conformational properties of the denatured states of proteins.

Rate of intrachain contact formation in an unfolded protein: temperature and denaturant effects

We have measured the effect of temperature and denaturant concentration on the rate of intrachain diffusion in an unfolded protein. After photodissociating a ligand from the heme iron of unfolded horse cytochrome c, we use transient optical absorption spectroscopy to measure the time scale of the diffusive motions that bring the heme, located at His18, into contact with its native ligand, Met80. Measuring the rate at which this 62 residue intrachain loop forms under both folding and unfolding conditions, we find a significant effect of denaturant on the chain dynamics. The diffusion of the chain accelerates as denaturant concentration decreases, with the contact formation rate approaching a value near approximately 6x10(5) s(-1) in the absence of denaturant. This result agrees well with an extrapolation from recent loop formation measurements in short synthetic peptides. The temperature dependence of the rate of contact formation indicates an Arrhenius activation barrier, Ea approximately 20 kJ/mol, at high denaturant concentrations, comparable to what is expected from solvent viscosity effects alone. Although Ea increases by several kBT as denaturant concentration decreases, the overall rate of diffusion nevertheless increases. These results indicate that inter-residue energetic interactions do not control conformational diffusion in unfolded states, even under folding conditions.

pH dependence of formation of a partially unfolded state of a Lys 73 --> His variant of iso-1-cytochrome c: implications for the alkaline conformational transition of cytochrome c

The alkaline conformational transition of a lysine 73 --> histidine variant of iso-1-cytochrome c has been studied. The transition has been monitored at 695 nm, a band sensitive to the presence of the heme-methionine 80 bond, at the heme Soret band which is sensitive to the nature of the heme ligand, and by NMR methods. The guanidine hydrochloride dependence of the alkaline conformational transition has also been monitored. The histidine 73 protein has an unusual biphasic alkaline conformational transition at both 695 nm and the heme Soret band, consistent with a three-state process. The conformational transition is fully reversible. An equilibrium model has been developed to account for this behavior. With this model, it has been possible to obtain the acid constant for the trigger group, pK(H), of the low-pH phase from the equilibrium data. A pK(H) value of 6.6 +/- 0.1 in H(2)O was obtained, consistent with a histidine acting as the trigger group. The NMR data for the low-pH phase of the alkaline conformational transition are consistent with an imidazole ligand replacing Met 80. For the high-pH phase of the biphasic alkaline transition, the NMR data are consistent with lysine 79 being the heme ligand. Guanidine hydrochloride m values of 1.67 +/- 0.08 and 1.1 +/- 0.2 kcal mol(-1) M(-1) were obtained for the low- and high-pH phases of the biphasic alkaline transition of the histidine 73 protein, respectively, consistent with a greater structural disruption for the low-pH phase of the transition.

Characterization of an alkaline transition intermediate stabilized in the Phe82Trp variant of yeast iso-1-cytochrome c

In general, mutation of the phylogenetically conserved residue Phe82 in yeast iso-1-cytochrome c destabilizes the native conformation of the protein by facilitating the ligand exchange reactions that are associated with the alkaline conformational transitions of the ferricytochrome. Of the Phe82 variants surveyed thus far, Phe82Trp is unique in that it adopts a thermodynamically stable, high-spin conformation at mildly alkaline pH. This species exhibits spectroscopic features that can only be detected transiently in other ferricytochromes c within the first 100 ms immediately after a pH-jump from neutrality to pH >10. Spectroscopic characterization of this high-spin reaction intermediate suggests that in addition to an obligatory pentacoordinate heme iron, a group within the heme pocket coordinates the heme iron but is then replaced either by Met80, to revert to the native conformation, or by Lys73 or Lys79, to yield one of the conventional alkaline conformers. Evidence is presented to suggest that this group is either a hydroxide ion or Tyr67 rather than a loosely bound Met80.

Coupled kinetic traps in cytochrome c folding: His-heme misligation and proline isomerization

The effect of His-heme misligation on folding has been investigated for a triple mutant of yeast iso-2 cytochrome c (N26H,H33N,H39K iso-2). The variant contains a single misligating His residue at position 26, a location at which His residues are found in several cytochrome c homologues, including horse, tuna, and yeast iso-1. The amplitude for fast phase folding exhibits a strong initial pH dependence. For GdnHCl unfolded protein at an initial pH<5, the observed refolding at final pH 6 is dominated by a fast phase (tau(2f)=20 ms, alpha(2f)=90 %) that represents folding in the absence of misligation. For unfolded protein at initial pH 6, folding at final pH 6 occurs in a fast phase of reduced amplitude (alpha(2f) approximately 20 %) but the same rate (tau(2f)=20 ms), and in two slower phases (tau(m)=6-8 seconds, alpha(m) approximately 45 %; and tau(1b)=16-20 seconds, alpha(1b) approximately 35 %). Double jump experiments show that the initial pH dependence of the folding amplitudes results from a slow pH-dependent equilibrium between fast and slow folding species present in the unfolded protein. The slow equilibrium arises from coupling of the His protonation equilibrium to His-heme misligation and proline isomerization. Specifically, Pro25 is predominantly in trans in the unligated low-pH unfolded protein, but is constrained in a non-native cis isomerization state by His26-heme misligation near neutral pH. Refolding from the misligated unfolded form proceeds slowly due to the large energetic barrier required for proline isomerization and displacement of the misligated His26-heme ligand.