Kinetic studies on isomerization of ferricytochrome c in alkaline and acid pH ranges by the circular dichroism stopped-flow method

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.

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.

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.

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.

Characterization of alkaline transitions in ferricytochrome c using carbon-deuterium infrared probes

The alkaline-induced structural transitions of ferricytochrome c have been studied intensively as a model for how changes in metal ligation contribute to protein function and folding. Previous studies have demonstrated that multiple non-native species accumulate with increasing pH. Here, we used a combination of experiments and simulations to provide a high-resolution view of the changes associated with increasing alkaline conditions. Alkaline-induced transitions were characterized under equilibrium conditions by following changes in the IR absorptions of carbon-deuterium chromophores incorporated at Leu68, Lys72, Lys73, Lys79, and Met80. The data suggest that at least four intermediates are formed as the pH is increased prior to complete unfolding of the protein. The first alkaline transition observed appears to be driven by a single deprotonation and occurs with a midpoint of pH 8.8, but surprisingly, the intermediate formed does not appear to be one of the well-characterized lysine misligates. At higher pH, second and third deprotonations, with a combined apparent midpoint pH of 10.2, induce transitions to Lys73- or Lys79-misligated species. Interestingly, the lysine misligates appear to undergo iron reduction by the coordinated amine. A transition from the lysine misligates to another intermediate, likely a hydroxide-misligated species, is associated with a fourth deprotonation and a midpoint of pH 10.7. Finally, the protein loses tertiary structure with a fifth deprotonation that occurs with a midpoint of pH 12.7. Native topology-based models with enforced misligation are employed to help understand the structures of the observed intermediates.

A millisecond infrared stopped-flow apparatus

In this paper, we have described a stopped-flow apparatus that is capable of measuring infrared kinetics in the amide I' region of a protein's vibrational spectrum. The dead time of this setup, determined by the reducing reaction of 2,6-Dichlorophenolindophenol by L-ascorbic acid, is between 6 to 15 ms, depending on the flow rate. Therefore, this stopped-flow IR method provides a means of measuring infrared kinetics in a time window that is not easily accessible to other mixing-based IR techniques. Using this apparatus, we have studied the alkaline transition of cytrochrome c and have found that this conformational event proceeds in a biphasic manner. The characteristic time constants of these two phases were determined to be 68 +/- 20 ms and 624 +/- 37 ms, respectively.

Proton-mediated dynamics of the alkaline conformational transition of yeast iso-1-cytochrome c

The kinetics of the alkaline conformational transition of a Lys 73-->His variant of iso-1-cytochrome c have been investigated using pH jump stopped-flow methods to probe the nature of the ionizable "trigger" group for this conformational change. This mutation moves the pK(a) of the ligand replacing Met 80 from about 10.5 to approximately 6.6 and has unmasked two other ionizable groups, besides the ligand replacing Met 80, that modulate the kinetics of this process. The results are discussed in terms of the impact of ionization equilibria on protein folding mechanisms.

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.

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.

Overexpression of Desulfovibrio vulgaris Hildenborough cytochrome c553 in Desulfovibrio desulfuricans G200. Evidence of conformational heterogeneity in the oxidized protein by NMR

Plasmid pRC41, containing the cyf gene encoding cytochrome c533 from Desulfovibrio vulgaris Hildenborough, was transferred by conjugation from Escherichia coli to Desulfovibrio desulfuricans G200. The structural properties of the purified protein were studied by one-dimensional and two-dimensional NMR. A heterogeneity in the folding of the cytochrome isolated from D. vulgaris Hildenborough and from D. desulfuricans G200 was observed for the oxidized from. Temperature, pH and salt-dependence studies indicated that the heterogeneity does not result from an intermediate in the protein unfolding process, but derives from two conformations which are not in dynamic equilibrium.

NMR study of the alkaline isomerization of ferricytochrome c

The pH-induced isomerization of horse heart cytochrome c has been studied by 1H NMR. We find that the transition occurring in D2O with a pKa measured as 9.5 +/- 0.1 is from the native species to a mixture of two basic forms which have very similar NMR spectra. The heme methyl peaks of these two forms have been assigned by 2D exchange NMR. The forward rate constant (native to alkaline cytochrome c) has a value of 4.0 +/- 0.6 s-1 at 27 degrees C and is independent of pH; the reverse rate constant is pH-dependent. The activation parameters are delta H not equal to = 12.8 +/- 0.8 kcal.mol1, delta S not equal to = -12.9 +/- 2.0 e.u. for the forward reaction and delta H not equal to = 6.0 +/- 0.3 kcal.mol-1, delta S not equal to = -35.1 +/- 1.3 e.u. for the reverse reaction (pH* = 9.28). delta H degree and delta S degree for the isomerization are 6.7 +/- 0.6 kcal.mol-1 and 21.9 +/- 1.0 e.u., respectively.

Kinetic aspects of the interaction of horse heart cytochrome c with sodium dodecyl sulfate

The sodium dodecyl sulfate (SDS) concentration dependence of spectral changes in circular dichroism (CD) and in absorbance of cytochrome c were examined in the far-ultraviolet region, aromatic region, and the Soret band. The Soret peak obtained in 0.60 mM SDS was nine times greater than that of the native state. (The critical micelle concentration, CMC, of SDS was 2.2 mM in the phosphate buffer used.) The results indicated that the drastic change at the Soret band did not accompany the corresponding large-scale change in secondary structure of the protein. In the stopped-flow measurements, two and three processes were followed at 406 nm below and above the CMC, respectively. At 289 nm only one process was observed, and this corresponded to the second process at 406 nm. Therefore, the second process at 406 nm was considered to be a change in tertiary structure around the heme group. The first process and the third process seemed to reflect a change in the heme environment; the former appeared to be due to a solvent effect and the latter due to a binding effect of a large number of dodecyl sulfate ions.

A spin label study of conformational changes in cytochrome c

Spin-labeled pig heart cytochromes c singly modified at Met-65, Tyr-74 and at one of the lysine residues, Lys-72 or Lys-73, were investigated by the ESR method under conditions of different ligand and redox states of the heme and at various pH values. Replacement of Met-80 by the external ligand, cyanide, was shown to produce a sharp increase in the mobility of all the three bound labels while reduction of the spin-labeled ferricytochromes c did not cause any marked changes in their ESR spectra. In the pH range 6-13, two conformational transitions in ferricytochrome c were observed which preceded its alkaline denaturation: the first with pK 9.3 registered by the spin label at the Met-65 position, and the second with pK 11.1 registered by the labels bound to Tyr-74 and Lys-72(73). The conformational changes in the 'left-hand part' of ferricytochrome c are most probably induced in both cases by the exchange of internal protein ligands at the sixth coordination site of the heme.