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

Applications of the Newly Developed Force-Field Parameters Uncover a Dynamic Nature of Ω-Loop C in the Lys-Ligated Alkaline Form of Cytochrome c

Lys-ligated cytochromes make up an emerging family of heme proteins. Density functional theory calculations on the amine/imidazole-ligated c-type ferric heme were employed to develop force-field parameters for molecular dynamics (MD) simulations of structural and dynamic features of these proteins. The new force-field parameters were applied to the alkaline form of yeast iso-1 cytochrome c to rationalize discrepancies resulting from distinct experimental conditions in prior structural studies and to provide insights into the mechanisms of the alkaline transition. Our simulations have revealed the dynamic nature of Ω-loop C in the Lys-ligated protein and its unfolding in the Lys-ligated conformer having this loop in the same position as in the native Met-ligated protein. The proximity of Tyr67 or Tyr74 to the Lys ligand of ferric heme iron suggests a possible mechanism of the backward alkaline transition where a proton donor Tyr assists in Lys dissociation. The developed force-field parameters will be useful in structural and dynamic characterization of other native or engineered Lys-ligated heme proteins.

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.

Real-Time Observation of Conformational Changes and Translocation of Endogenous Cytochrome c within Intact Mitochondria

Cytochrome c (cyt c) is a multifunctional protein with varying conformations. However, the conformation of cyt c in its native environment, mitochondria, is still unclear. Here, we applied NMR spectroscopy to investigate the conformation and location of endogenous cyt c within intact mitochondria at natural isotopic abundance, mainly using widespread methyl groups as probes. By monitoring time-dependent chemical shift perturbations, we observed that most cyt c is located in the inner mitochondrial membrane and partially unfolded, which is distinct from its native conformation in solution. When suffering oxidative stress, cyt c underwent oxidative modifications due to increasing reactive oxygen species (ROS), weakening electrostatic interactions with the membrane, and gradually translocating into the inner membrane spaces of mitochondria. Meanwhile, the lethality of oxidatively modified cyt c to cells was reduced compared with normal cyt c. Our findings significantly improve the understanding of the molecular mechanisms underlying the regulation of ROS by cyt c in mitochondria. Moreover, it highlights the potential of NMR to monitor high-concentration molecules at a natural isotopic abundance within intact cells or organelles.

Effect on intrinsic peroxidase activity of substituting coevolved residues from Ω-loop C of human cytochrome c into yeast iso-1-cytochrome c

Naturally-occurring variants of human cytochrome c (Cytc) that induce thrombocytopenia IV occur within Ω-loop C (residues 40-57). These variants enhance the peroxidase activity of human Cytc apparently by facilitating access to the heme by destabilizing Ω-loops C and D (residues 70-85). Given the importance of peroxidase activity in the early stages of apoptosis, we identified three sites with the EVmutation algorithm in or near Ω-loop C that coevolve and differ between yeast iso-1-Cytc and human Cytc. We prepared iso-1-Cytc variants with all possible combinations of the S40T, V57I and N63T substitutions to determine if these residues decrease the peroxidase activity of iso-1-Cytc to that of human Cytc producing an effective off state for a peroxidase signaling switch. At pH 6 and above, all variants significantly decreased peroxidase activity. However, the correlation of peroxidase activity with local and global stability, expected if cooperative unfolding of Ω-loops C and D is required for peroxidase activity, was generally poor. The m-values derived from the guanidine hydrochloride dependence of the kinetics of imidazole binding to horse Cytc, which is well-characterized by native-state hydrogen exchange methods, and K72A/K73A/K79A iso-1-Cytc show that local structural fluctuations and not subglobal cooperative unfolding of Ω-loops C and D are sufficient to permit binding of a small molecule like peroxide to the heme. A 2.46 Å structure of N63T iso-1-Cytc identifies a change to a hydrogen bond network linking Ω-loops C and D that could modulate the local fluctuations needed for the intrinsic peroxidase activity of Cytc.

NMR Reveals the Conformational Changes of Cytochrome C upon Interaction with Cardiolipin

Conformational change of cytochrome c (cyt c) caused by interaction with cardiolipin (CL) is an important step during apoptosis, but the underlying mechanism is controversial. To comprehensively clarify the structural transformations of cyt c upon interaction with CL and avoid the unpredictable alias that might come from protein labeling or mutations, the conformation of purified yeast iso–1 cyt c with natural isotopic abundance in different contents of CL was measured by using NMR spectroscopy, in which the trimethylated group of the protein was used as a natural probe. The data demonstrate that cyt c has two partially unfolded conformations when interacted with CL: one with Fe–His33 coordination and the other with a penta–coordination heme. The Fe–His33 coordination conformation can be converted into a penta–coordination heme conformation in high content of CL. The structure of cyt c becomes partially unfolded with more exposed heme upon interaction with CL, suggesting that cyt c prefers a high peroxidase activity state in the mitochondria, which, in turn, makes CL easy to be oxidized, and causes the release of cyt c into the cytoplasm as a trigger in apoptosis.

Anion-Specific Effects on the Alkaline State of Cytochrome c

Specific effects of anions on the structure, thermal stability, and peroxidase activity of native (state III) and alkaline (state IV) cytochrome c (cyt c) have been studied by the UV-VIS absorbance spectroscopy, intrinsic tryptophan fluorescence, and circular dichroism. Thermal and isothermal denaturation monitored by the tryptophan fluorescence and circular dichroism, respectively, implied lower stability of cyt c state IV in comparison with the state III. The pK_a value of alkaline isomerization of cyt c depended on the present salts, i.e., kosmotropic anions increased and chaotropic anions decreased pK_a (Hofmeister effect on protein stability). The peroxidase activity of cyt c in the state III, measured by oxidation of guaiacol, showed clear dependence on the salt position in the Hofmeister series, while cyt c in the alkaline state lacked the peroxidase activity regardless of the type of anions present in the solution. The alkaline isomerization of cyt c in the presence of 8 M urea, measured by Trp59 fluorescence, implied an existence of a high-affinity non-native ligand for the heme iron even in a partially denatured protein conformation. The conformation of the cyt c alkaline state in 8 M urea was considerably modulated by the specific effect of anions. Based on the Trp59 fluorescence quenching upon titration to alkaline pH in 8 M urea and molecular dynamics simulation, we hypothesize that the Lys79 conformer is most likely the predominant alkaline conformer of cyt c. The high affinity of the sixth ligand for the heme iron is likely a reason of the lack of peroxidase activity of cyt c in the alkaline state.

Insights on the Conformational Ensemble of Cyt C Reveal a Compact State during Peroxidase Activity

Cytochrome c (cyt c) is known for its role in the electron transport chain but transitions to a peroxidase-active state upon exposure to oxidative species. The peroxidase activity ultimately results in the release of cyt c into the cytosol for the engagement of apoptosis. The accumulation of oxidative modifications that accompany the onset of the peroxidase function are well-characterized. However, the concurrent structural and conformational transitions of cyt c remain undercharacterized. Fast photochemical oxidation of proteins (FPOP) coupled with mass spectrometry is a protein footprinting technique used to structurally characterize proteins. FPOP coupled with native ion mobility separation shows that exposure to H₂O₂ results in the accumulation of a compact state of cyt c. Subsequent top-down fragmentation to localize FPOP modifications reveals changes in heme coordination between conformers. A time-resolved functional assay suggests that this compact conformer is peroxidase active. Altogether, combining FPOP, ion mobility separation, and top-down and bottom-up mass spectrometry allows us to discern individual conformations in solution and obtain a better understanding of the conformational ensemble and structural transitions of cyt c as it transitions from a respiratory role to a proapoptotic role.

A Heme Propionate Staples the Structure of Cytochrome c for Methionine Ligation to the Heme Iron

Ligand-switch reactions at the heme iron are common in biological systems, but their mechanisms and the features of the polypeptide fold that support dual ligation are not well understood. In cytochrome c (cyt c), two low-stability loops (Ω-loop C and Ω-loop D) are connected by the heme propionate HP6. At alkaline pH, the native Met80 ligand from Ω-loop D switches to a Lys residue from the same loop. Deprotonation of an as yet unknown group triggers the alkaline transition. We have created the two cyt c variants T49V/K79G and T78V/K79G with altered connections of these two loops to HP6. Electronic absorption, NMR, and EPR studies demonstrate that at pH 7.4 ferric forms of these variants are Lys-ligated, whereas ferrous forms maintain the native Met80 ligation. Measurements of protein stability, cyclic voltammetry, pH-jump and gated electron-transfer kinetics have revealed that these Thr to Val substitutions greatly affect the alkaline transition in both ferric and ferrous proteins. The substitutions modify the stability of the Met-ligated species and reduction potentials of the heme iron. The kinetics of ligand-switch processes are also altered, and analyses of these effects implicate redox-dependent differences in metal-ligand interactions and the role of the protein dynamics, including cross-talk between the two Ω-loops. With the two destabilized variants, it is possible to map energy levels for the Met- and Lys-ligated species in both ferric and ferrous proteins and assess the role of the protein scaffold in redox-dependent preferences for these two ligands. The estimated shift in the heme iron reduction potential upon deprotonation of the "trigger" group is consistent with those associated with deprotonation of an HP, suggesting that HP6, on its own or as a part of a hydrogen-bonded cluster, is a likely "trigger" for the Met to Lys ligand switch.

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.

Histidine-Lysine Axial Ligand Switching in a Hemoglobin: A Role for Heme Propionates

The hemoglobin of Synechococcus sp. PCC 7002, GlbN, is a monomeric group I truncated protein (TrHb1) that coordinates the heme iron with two histidine ligands at neutral pH. One of these is the distal histidine (His46), a residue that can be displaced by dioxygen and other small molecules. Here, we show with mutagenesis, electronic absorption spectroscopy, and nuclear magnetic resonance (NMR) spectroscopy that at high pH and exclusively in the ferrous state, Lys42 competes with His46 for the iron coordination site. When b heme is originally present, the population of the lysine-bound species remains too small for detailed characterization; however, the population can be increased significantly by using dimethyl-esterified heme. Electronic absorption and NMR spectroscopies showed that the reversible ligand switching process occurs with an apparent pK_a of 9.3 and a Lys-ligated population of ∼60% at the basic pH limit in the modified holoprotein. The switching rate, which is slow on the chemical shift time scale, was estimated to be 20-30 s^-1 by NMR exchange spectroscopy. Lys42-His46 competition and attendant conformational rearrangement appeared to be related to weakened bis-histidine ligation and enhanced backbone dynamics in the ferrous protein. The pH- and redox-dependent ligand exchange process observed in GlbN illustrates the structural plasticity allowed by the TrHb1 fold and demonstrates the importance of electrostatic interactions at the heme periphery for achieving axial ligand selection. An analogy is drawn to the alkaline transition of cytochrome c, in which Lys-Met competition is detected at alkaline pH, but, in contrast to GlbN, in the ferric state only.

Effect of a K72A Mutation on the Structure, Stability, Dynamics, and Peroxidase Activity of Human Cytochrome c

We test the hypothesis that Lys72 suppresses the intrinsic peroxidase activity of human cytochrome c, as observed previously for yeast iso-1-cytochrome c [McClelland, L. J., et al. (2014) Proc. Natl. Acad. Sci. U. S. A. 111, 6648-6653]. A 1.25 Å X-ray structure of K72A human cytochrome c shows that the mutation minimally affects structure. Guanidine hydrochloride denaturation demonstrates that the K72A mutation increases global stability by 0.5 kcal/mol. The K72A mutation also increases the apparent pKa of the alkaline transition, a measure of the stability of the heme crevice, by 0.5 unit. Consistent with the increase in the apparent pKa, the rate of formation of the dominant alkaline conformer decreases, and this conformer is no longer stabilized by proline isomerization. Peroxidase activity measurements show that the K72A mutation increases kcat by 1.6-4-fold at pH 7-10, an effect larger than that seen for the yeast protein. X-ray structures of wild type and K72A human cytochrome c indicate that direct interactions of Lys72 with the far side of Ω-loop D, which are seen in X-ray structures of horse and yeast cytochrome c and could suppress peroxidase activity, are lacking. Instead, we propose that the stronger effect of the K72A mutation on the peroxidase activity of human versus yeast cytochrome c results from relief of steric interactions between the side chains at positions 72 and 81 (Ile in human vs Ala in yeast), which suppress the dynamics of Ω-loop D necessary for the intrinsic peroxidase activity 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.

Dynamics of the His79-heme alkaline transition of yeast iso-1-cytochrome c probed by conformationally gated electron transfer with Co(II)bis(terpyridine)

Alkaline conformers of cytochrome c may be involved in both its electron transport and apoptotic functions. We use cobalt(II)bis(terpyridine), Co(terpy)2(2+), as a reagent for conformationally gated electron-transfer (gated ET) experiments to study the alkaline conformational transition of K79H variants of yeast iso-1-cytochrome c expressed in Escherichia coli , WT*K79H, with alanine at position 72 and Saccharomyces cerevisiae , yK79H, with trimethyllysine (Tml) at position 72. Co(terpy)2(2+) is well-suited to the 100 ms to 1 s time scale of the His79-mediated alkaline conformational transition of these variants. Reduction of the His79-heme alkaline conformer by Co(terpy)2(2+) occurs primarily by gated ET, which involves conversion to the native state followed by reduction, with a small fraction of the His79-heme alkaline conformer directly reduced by Co(terpy)2(2+). The gated ET experiments show that the mechanism of formation of the His79-heme alkaline conformer involves only two ionizable groups. In previous work, we showed that the mechanism of the His73-mediated alkaline conformational transition requires three ionizable groups. Thus, the mechanism of heme crevice opening depends upon the position of the ligand mediating the process. The microscopic rate constants provided by gated ET studies show that mutation of Tml72 (yK79H variant) in the heme crevice loop to Ala72 (WT*K79H variant) affects the dynamics of heme crevice opening through a small destabilization of both the native conformer and the transition state relative to the His79-heme alkaline conformer. Previous pH jump data had indicated that the Tml72→Ala mutation primarily stabilized the transition state for the His79-mediated alkaline conformational transition.

Structure and function of heme proteins in non-native states: a mini-review

Heme proteins perform various biological functions ranging from electron transfer, oxygen binding and transport, catalysis, to signaling. Although adopting proper native states is very important for these functions, progresses in representative heme proteins, including cytochrome c (cyt c), cytochrome b5 (cyt b5), myoglobin (Mb), neuroglobin (Ngb), cytochrome P450 (CYP) and heme-based sensor proteins such as CO sensor CooA, showed that various native functions, or new functions evolved, are also closely associated with non-native states. The structure and function relationship of heme proteins in non-native states is thus as important as that in native states for elucidating the precise roles of heme proteins in biological systems.

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.

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.

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.

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.

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.

Thermodynamics of the alkaline transition in phytocyanins

The thermodynamics of the alkaline transition which influences the spectral and redox properties of the type 1 copper center in phytocyanins has been determined spectroscopically. The proteins investigated include Rhus vernicifera stellacyanin, cucumber basic protein and its Met89Gln variant, and umecyanin, the stellacyanin from horseradish roots, along with its Gln95Met variant. The changes in reaction enthalpy and entropy within the protein series show partial compensatory behavior. Thus, the reaction free energy change (hence the pK (a) value) is rather variable. This indicates that species-dependent differences in reaction thermodynamics, although containing an important contribution from changes in the hydrogen-bonding network of water molecules in the hydration sphere of the protein (which feature enthalpy-entropy compensation), are to a large extent protein-based. The data for axial ligand variants are consistent with the hypothesis of a copper-binding His as the deprotonating residue responsible for this transition.

Conformational stability and dynamics of cytochrome c affect its alkaline isomerization

The alkaline isomerization of horse heart ferricytochrome c (cyt c) has been studied by electronic absorption spectroscopy in the presence of the Hofmeister series of anions: chloride, bromide, rhodanide and perchlorate. The anions significantly affect the apparent pK (a) value of the transition in a concentration-dependent manner according to their position in the Hofmeister series. The Soret region of the absorption spectra is not affected by the presence of the salts and shows no significant structural perturbation of the heme crevice. In the presence of perchlorate and rhodanide anions, the cyanide exchange rate between the bulk solvent and the binding site is increased. These results imply higher flexibility of the protein structure in the presence of chaotropic salts. The thermal and isothermal denaturations monitored by differential scanning calorimetry and circular dichroism, respectively, showed a decrease in the conformational stability of cyt c in the presence of the chaotropic salts. A positive correlation between the stability, DeltaG, of cyt c and the apparent pK (a) values that characterize the alkaline transition indicates the presence of a thermodynamic linkage between these conformational transitions. In addition, the rate constant of the cyanide binding and the partial molar entropies of anions negatively correlate with the pK (a) values. This indicates the important role of anion-induced solvent reorganization on the structural flexibility of cyt c in the alkaline transitions.

Structural Characterization of an Equilibrium Unfolding Intermediate in Cytochrome c

Although the denaturant-induced unfolding transition of cytochrome c was initially thought to be a cooperative process, recent spectroscopic studies have shown deviations from two-state behavior consistent with accumulation of an equilibrium intermediate. However, little is known about the structural and thermodynamic properties of this state, and whether it is stabilized by the presence of non-native heme ligands. We monitored the reversible denaturant-induced unfolding equilibrium of oxidized horse cytochrome c using various spectroscopic probes, including fluorescence, near and far-UV CD, heme absorbance bands in the Soret, visible and near-IR regions of the spectrum, as well as 2D NMR. Global fitting techniques were used for a quantitative interpretation of the results in terms of a three-state model, which enabled us to determine the intrinsic spectroscopic properties of the intermediate. A well-populated intermediate was observed in equilibrium experiments at pH 5 using either guanidine-HCl or urea as a denaturant, both for wild-type cytochrome c as well as an H33N mutant chosen to prevent formation of non-native His-heme ligation. For a more detailed structural characterization of the intermediate, we used 2D 1H-15N correlation spectroscopy to follow the changes in peak intensity for individual backbone amide groups. The equilibrium state observed in our optical and NMR studies contains many native-like structural features, including a well-structured alpha-helical sub-domain, a short Trp59-heme distance and solvent-shielded heme environment, but lacks the native Met80 sulfur-iron linkage and shows major perturbations in side-chain packing and other tertiary interactions. These structural properties are reminiscent of the A-state of cytochrome c, a compact denatured form found under acidic high-salt conditions, as well as a kinetic intermediate populated at a late stage of folding. The denaturant-induced intermediate also resembles alkaline forms of cytochrome c with altered heme ligation, suggesting that disruption of the native methionine ligand favors accumulation of structurally analogous states both in the presence and absence of non-native ligands.

Impact of pH on Bax α conformation, oligomerisation and mitochondrial integration

The change in the conformation of Bax at the onset of apoptosis is a determinant for the execution of this cell death programme. However, very few models can account for this modification and the factors involved in this process remain elusive. We have analysed the modifications in the conformation induced by a variation in pH using a cell-free assay. We show that a moderate basic or acidic pH can induce apoptotic-like changes in the conformation of Bax, such as the exposure of the N-terminal or the BH3 domain. These changes in the conformation are associated with the binding of Bax to mitochondria and an enhanced Bax homo- and oligomerisation. Our results suggest that variations in the pH, in a range consistent with that often observed during apoptosis, are sufficient to trigger Bax translocation to mitochondria and the subsequent release of apoptogenic factors from this organelle.

The 40s ?-loop plays a critical role in the stability and the alkaline conformational transition of cytochrome c

The structural and redox properties of a non-covalent complex reconstituted upon mixing two non-contiguous fragments of horse cytochrome c, the residues 1-38 heme-containing N-fragment with the residues 57-104 C-fragment, have been investigated. With respect to native cyt c, the complex lacks a segment of 18 residues, corresponding, in the native protein, to an omega (Omega)-loop region. The fragment complex shows compact structure, native-like alpha-helix content but a less rigid atomic packing and reduced stability with respect to the native protein. Structural heterogeneity is observed at pH 7.0, involving formation of an axially misligated low-spin species and consequent partial displacement of Met80 from the sixth coordination position of the heme-iron. Spectroscopic data suggest that a lysine (located in the Met80-containing loop, namely Lys72, Lys73, or Lys79) replaces the methionine residue. The residues 1-38/57-104 fragment complex shows an unusual biphasic alkaline titration characterized by a low (p K(a1)=6.72) and a high p K(a)-associated state transition (p K(a2)=8.56); this behavior differs from that of native cyt c, which shows a monophasic alkaline transition (p K(a)=8.9). The data indicate that the 40s Omega-loop plays an important role in the stability of cyt c and in ensuring a correct alkaline conformational transition of the protein.

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.

Protonation of two adjacent tyrosine residues influences the reduction of cytochrome c by diphenylacetaldehyde: a possible mechanism to select the reducer agent of heme iron

We have shown that diphenlacetaldehyde (DPAA) is able to promote mitochondrial DeltaPsi disruption accompanied by damage in mitochondrial DNA, lipids, and proteins [Almeida, A. M.; Bechara, E. J. H.; Vercesi, A. E.; Nantes, I. L. Free Radic. Biol. Med. 27:744-747; 1999]. In this work, DPAA was used as a model of carbonyl reagent for cytochrome c. The results suggest that DPAA is a redox cytochrome c modifier. Conversion of Fe(III) to Fe(II) cytochrome c promoted by DPAA is pH dependent. The second-order rate determined for heme iron reduction (k2) is 698 M(-1) s(-1) and this process occurs with an activation energy of 8.5 +/- 0.8 kcal/mol. Analysis of the pH profile suggests the presence of two ionizable cytochrome c groups (pKa1 = 8.9 and pKa2 = 11.4) related to the electron transfer from DPAA to heme iron. The heats of ionization of the two prototropic groups, pKa1 (DeltaH(ion) = 6.5 kcal/mol, DeltaS(ion) = -29.0 cal/mol.K), and pKa2 (DeltaH(ion) = 5.0 kcal/mol, DeltaS(ion) = -24.0 cal/mol.K), suggest involvement of two tyrosine residues, probably Y67 and Y74, related to DPAA-promoted heme iron reduction. The cytochrome c chemical modification by iodination of tyrosine groups significantly decreased the reduction rate promoted by DPAA, and shifted the pH(opt) value from 10.0 to 9.25. The cytochrome c-promoted DPAA electron abstraction quickly produces the expected enol-derived radical, as indicated by 3,5-dibromo-4-nitrosobenzenesulfonate (DBNBS) spin trapping EPR measurements. This radical reacts with molecular oxygen, producing a peroxyl intermediate radical that, via a putative dioxetane intermediate, promotes formation of benzophenone as the main final product of this reaction, detected by high-performance liquid chromatography coupled with tandem mass spectrometry.

Evaluation of cooperative interactions between substructures of iso-1-cytochrome c using double mutant cycles

A double mutant cycle has been used to evaluate interaction energies between the global stabilizer mutation asparagine 52 --> isoleucine (N52I) in iso-1-cytochrome c and mutations producing single surface histidines at positions 26, 33, 39, 54, 73, 89, and 100. These histidine mutation sites are distributed through the four cooperative folding units of cytochrome c. The double mutant cycle starts with the iso-1-cytochrome c variant AcTM, a variant with no surface histidines and with asparagine at position 52. Isoleucine is added singly at position 52, AcTMI52 variant, as are the surface histidines, AcHX variants, where X indicates the histidine sequence position. The double mutant variants, AcHXI52, provide the remaining corner of the double mutant cycle. The stabilities of all variants were determined by guanidine hydrochloride denaturation and interaction energies were calculated between position 52 and each histidine site. Six of the seven double mutants show additive (AcH33I52, AcH39I52, AcH54I52, AcH89I52, and AcH100I52) stability effects or weak interaction energies (AcH73I52) of the histidine mutations and the N52I mutation, consistent with cooperative effects on protein folding and stability being sparsely distributed through the protein structure. The AcH26I52 variant shows a strong favorable interaction energy, 2.0 +/- 0.5 kcal/mol, between the N52I mutation in one substructure and the addition of His 26 to an adjacent substructure. The data are consistent with an entropic stabilization of the intersubstructure hydrogen bond between His 26 and Glu 44 by the Ile 52 mutation.

Folding units govern the cytochrome c alkaline transition

The alkaline transition of cytochrome c is a model for protein structural switching in which the normal heme ligand is replaced by another group. Stopped flow data following a jump to high pH detect two slow kinetic phases, suggesting two rate-limiting structure changes. Results described here indicate that these events are controlled by the same structural unfolding reactions that account for the first two steps in the reversible unfolding pathway of cytochrome c. These and other results show that the cooperative folding-unfolding behavior of protein foldons can account for a variety of functional activities in addition to determining folding pathways.

NMR investigation of the alkaline-like conformational transition of horse heart cytochrome c in the presence of exogenous thiazole

The conformational transition of horse heart cyt c in the presence of exogenous thiazole is investigated by NMR spectroscopy. Surprisingly, besides the native form and the ligand-bound form, another species (species A) exists at neutral pH. Titration of thiazole at alkaline condition indicates that this form has Lys as the axial ligands and in fact is one of alkaline cyt c conformers. At a high concentration of thiazole, species A has a dominant population at a pH of approximately 7. It is the first time that a pure alkaline conformer of cyt c is obtained at neutral pH. The observed NOEs show that species A maintains the overall protein fold although large structural rearrangements are expected in the distal pocket. Furthermore, the different effect of imidazole, pyridine and thiazole on the alkaline transformation of cyt c is discussed.

NMR study of the conformational transition of cytochrome c upon the displacement of Met80 by exogenous ligand: structural and magnetic characterization of azidoferricytochrome c

As the exogenous ligand-cytochrome c complexes were purported to represent models for the unfolding intermediate of cytochrome c, NMR spectroscopy has been utilized to study the azide adduct of horse heart cytochrome c. The structure of azidoferricytochrome c was modeled by restrained energy minimization using paramagenetic pseudocontact shifts as constraints. The bound azide moiety was found to be tilted approximately 15 degrees from the heme normal. The displacement of Met80 by the exogenous azide molecule causes large structural rearrangement in the distal cavity. Furthermore, the conformation transition associated with the swing out of the loop containing Met80 and the shift of the 50s-helix increases the solvent accessibility of the heme group. To elucidate the heme electronic structure of the complex, the paramagnetic 13C shifts were analyzed in terms of a model based on the pi molecular orbitals of the heme under perturbed D(4) symmetry. It turned out that the His-Fe bonding provides the protein constraint that orients the in-plane anisotropy in the complex. The electronic properties are in accordance with the calculated magnetic susceptibility anisotropy and the structural information.

Conservation of the free energy change of the alkaline isomerization in mitochondrial and bacterial cytochromes c

The thermodynamic parameters of the alkaline transition for oxidized native yeast iso-1 cytochrome c and Rhodopseudomonas palustris cytochrome c(2) (cytc(2)) have been determined through direct electrochemistry experiments carried out at variable pH and temperature and compared to those for horse and beef heart cytochromes c. We have found that both transition enthalpy and entropy are remarkably species dependent, following the order R. palustris cytc(2) >> beef (horse) heart cytc>yeast iso-1 cytc. Considering the high homology at the heme-protein interface in the native species, this variability is likely to be mainly determined by differences in the structural and solvation properties and the relative abundance of the various alkaline conformers. Notably, changes in transition enthalpy and entropy among these cytochromes c are compensative and result in small variations in the free energy change of the process (which amounts approximately to +50 kJ mol(-1)) and consequently in the apparent pK(a) value. This compensation indicates that solvent reorganization effects play an important role in the thermodynamics of the transition. This mechanism is functional to ensure a relatively high pK(a) value for the alkaline transition, which is needed to preserve His,Met ligation to the heme iron in cytochrome c at physiological pH and temperature, hence the E(o) value required for the biological function.