Kinetic study of isomerization of ferricytochrome c at alkaline pH

Naturally Occurring A51V Variant of Human Cytochrome c Destabilizes the Native State and Enhances Peroxidase Activity

The A51V variant of human cytochrome c is linked to thrombocytopenia 4 (THC4), a condition that causes decreased blood platelet counts. A 1.82 Å structure of the A51V variant shows only minor changes in tertiary structure relative to the wild-type (WT) protein. Guanidine hydrochloride denaturation demonstrates that the global stability of the A51V variant is 1.3 kcal/mol less than that of the WT protein. The midpoint pH, pH_1/2, of the alkaline transition of the A51V variant is 1 unit less than that of the WT protein. Stopped-flow pH jump experiments show that the A51V substitution affects the triggering ionization for one of two kinetically distinguishable alkaline conformers and enhances the accessibility of a high-spin heme transient. The pH_1/2 for acid unfolding of the A51V variant is 0.7 units higher than for that of the WT protein. Consistent with the greater accessibility of non-native conformers for the A51V variant, the k_cat values for its peroxidase activity increase by 6- to 15-fold in the pH range of 5-8 versus those of the WT protein. These data along with previously reported data for the other THC4-linked variants, G41S and Y48H, underscore the role of Ω-loop C (residues 40-57) in modulating the peroxidase activity of cytochrome c early in apoptosis.

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.

Binding of imidazole to the heme of cytochrome c1 and inhibition of the bc1 complex from Rhodobacter sphaeroides: II. Kinetics and mechanism of binding

The kinetics of imidazole (Im) and N-methylimidazole (MeIm) binding to oxidized cytochrome (cyt) c(1) of detergent-solubilized bc(1) complex from Rhodobacter sphaeroides are described. The rate of formation of the cyt c(1)-Im complex exhibited three separated regions of dependence on the concentration of imidazole: (i) below 8 mM Im, the rate increased with concentration in a parabolic manner; (ii) above 20 mM, the rate leveled off, indicating a rate-limiting conformational step with lifetime approximately 1 s; and (iii) at Im concentrations above 100 mM, the rate substantially increased again, also parabolically. In contrast, binding of MeIm followed a simple hyperbolic concentration dependence. The temperature dependences of the binding and release kinetics of Im and MeIm were also measured and revealed very large activation parameters for all reactions. The complex concentration dependence of the Im binding rate is not consistent with the popular model for soluble c-type cytochromes in which exogenous ligand binding is preceded by spontaneous opening of the heme cleft, which becomes rate-limiting at high ligand concentrations. Instead, binding of ligand to the heme is explained by a model in which an initial and superficial binding facilitates access to the heme by disruption of hydrogen-bonded structures in the heme domain. For imidazole, two separate pathways of heme access are indicated by the distinct kinetics at low and high concentration. The structural basis for ligand entry to the heme cleft is discussed.

Chemical frustration in the protein folding landscape: grand canonical ensemble simulations of cytochrome c

A grand canonical formalism is developed to combine discrete simulations for chemically distinct species in equilibrium. Each simulation is based on a perturbed funneled landscape. The formalism is illustrated using the alkaline-induced transitions of cytochrome c as observed by FTIR spectroscopy and with various other experimental approaches. The grand canonical simulation method accounts for the acid/base chemistry of deprotonation, the inorganic chemistry of heme ligation and misligation, and the minimally frustrated folding energy landscape, thus elucidating the physics of protein folding involved with an acid/base titration of a protein. The formalism combines simulations for each of the relevant chemical species, varying by protonation and ligation states. In contrast to models based on perfectly funneled energy landscapes that contain only contacts found in the native structure, this study introduces "chemical frustration" from deprotonation and misligation that gives rise to many intermediates at alkaline pH. While the nature of these intermediates cannot be easily inferred from available experimental data, this study provides specific structural details of these intermediates, thus extending our understanding of how cytochrome c changes with an increase in pH. The results demonstrate the importance of chemical frustration for understanding biomolecular energy landscapes.

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.

Protein folding and misfolding: mechanism and principles

Two fundamentally different views of how proteins fold are now being debated. Do proteins fold through multiple unpredictable routes directed only by the energetically downhill nature of the folding landscape or do they fold through specific intermediates in a defined pathway that systematically puts predetermined pieces of the target native protein into place? It has now become possible to determine the structure of protein folding intermediates, evaluate their equilibrium and kinetic parameters, and establish their pathway relationships. Results obtained for many proteins have serendipitously revealed a new dimension of protein structure. Cooperative structural units of the native protein, called foldons, unfold and refold repeatedly even under native conditions. Much evidence obtained by hydrogen exchange and other methods now indicates that cooperative foldon units and not individual amino acids account for the unit steps in protein folding pathways. The formation of foldons and their ordered pathway assembly systematically puts native-like foldon building blocks into place, guided by a sequential stabilization mechanism in which prior native-like structure templates the formation of incoming foldons with complementary structure. Thus the same propensities and interactions that specify the final native state, encoded in the amino-acid sequence of every protein, determine the pathway for getting there. Experimental observations that have been interpreted differently, in terms of multiple independent pathways, appear to be due to chance misfolding errors that cause different population fractions to block at different pathway points, populate different pathway intermediates, and fold at different rates. This paper summarizes the experimental basis for these three determining principles and their consequences. Cooperative native-like foldon units and the sequential stabilization process together generate predetermined stepwise pathways. Optional misfolding errors are responsible for 3-state and heterogeneous kinetic folding.

Functional role of a protein foldon-An Ω-loop foldon controls the alkaline transition in ferricytochrome c

Hydrogen exchange results for cytochrome c and several other proteins show that they are composed of a number of foldon units which continually unfold and refold and account for some functional properties. Previous work showed that one Omega-loop foldon controls the rate of the structural switching and ligand exchange behavior of cytochrome c known as the alkaline transition. The present work tests the role of foldons in the alkaline transition equilibrium. We measured the effects of denaturant and 14 destabilizing mutations. The results show that the ligand exchange equilibrium is controlled by the stability of the same foldon unit implicated before. In addition, the results obtained confirm the epsilon-amino group of Lys79 and Lys73 as the alkaline replacement ligands and bear on the search for a triggering group.

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.

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.

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.

Identification of the ligand-exchange process in the alkaline transition of horse heart cytochrome c

Magnetic-circular-dichroism (m.c.d.) spectra over the wavelength range 300-2000 nm at room temperature and at 4.2K of horse heart cytochrome c are reported at a series of pH values between 7.8 and 11.0, encompassing the alkaline transition. The effect of glassing agents on the e.p.r. spectrum at various pH values is also reported. Comparison of these results with spectra obtained for the n-butylamine adduct of soybean leghaemoglobin support the hypothesis that lysine is the sixth ligand in the alkaline form of horse heart cytochrome c. The m.c.d. and e.p.r. spectra of horse heart cytochrome c in the presence of 1-methylimidazole have also been examined. These studies strongly suggest that histidine-18, the proximal ligand of the haem, is the ionizing group that triggers the alkaline transition. Low-temperature m.c.d. and e.p.r. spectra are also reported for Pseudomonas aeruginosa cytochrome c551. It is shown that no ligand exchange takes place at the haem in this species over the pH range 6.0-11.3.

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.

Kinetic and equilibrium studies of alkaline isomerization of vertebrate cytochromes c

Equilibria and kinetics of alkaline isomerization of seven ferricytochromes c from vertebrates were studied by pH-titration and pH-jump methods in the pH region of 7-12. In the equilibrium behavior, no significant difference was detected among the cytochromes c, whereas marked differences in the kinetic behavior were observed. According to the kinetic behavior of the isomerization, the cytochromes c examined fall into three classes: Group I (horse, sheep, dog and pigeon cytochromes c), Group II (tuna and bonito cytochromes c) and Group III (rhesus monkey cytochrome c). The kinetic results are interpreted in terms of the sequential scheme: Neutral form in equilibrium with fast Transient form in equilibrium with slow Alkaline form where the neutral and alkaline forms are the species stable at neutral and alkaline pH, respectively, and the transient form is a kinetic intermediate. From comparison of the primary sequences of the seven cytochromes c and the classification of these cytochromes c, it is concluded that the amino acid substitution Phe/Tyr at the 46-th position has a major influence on the kinetic behavior. In Group II and III cytochromes c, the ionization of Tyr-46 is suggested to bring about loosening of the heme crevice and thus facilitate the ligand replacement involved in the isomerization.

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

The isomerization of horse-heart ferricytochrome c caused by varying pH was kinetically studied by using circular dichroism (CD) and optical absorption stopped-flow techniques. In the pH range of 7--13, the existence of the three different forms of ferricytochrome c (pH less than 10, pH 10--12, and pH greater than 12) was indicated from the statistical difference CD spectra. On the basis of analyses of the stopped-flow traces in the near-ultraviolet and Soret wavelength regions, the isomerization of ferricytochrome c from neutral form to the above three alkaline forms was interpreted as follows (1) below pH 10, the replacement of the intrinsic ligand of methionine residue by lysine residue occurs; (2) between pH 10 and 12, the uncoupling of the polypeptide chain from close proximity of the heme group occurs first, followed by the interconversion of the intrinsic ligands; and (3) above pH 12, hydroxide form of ferricytochrome c is formed, though the replacement of the intrinsic ligand by extrinsic ligands may occur via different routes from those below pH 12. The CD changes at 288 nm and in the Soret region caused by the pH-jump (down) from pH 6.0 to 1.6 were compared with the appearance of the 620-nm absorption band ascribed to the formation of the high-spin form of ferricytochrome c. Both CD and absorption changes indicated that the isomerization at pH 1.6 consisted of two processes: one proceeded within the dead-time (about 2 ms) of the stopped-flow apparatus and the other proceeded at a determinable rate with the apparatus. On the basis of these results, the isomerization of ferricytochrome c at pH 1.6 was explained as follows: (1) the transition from the low-spin form to the high-spin forms occurs within about 2 ms, the dead-time of the stopped-flow apparatus; and (2) the polypeptide chain is unfolded after the formation of the high-spin form.

Equilibrium and kinetic measurements of the redox potentials of cytochromes c2 in vitro and in vivo

The equilibrium oxidation-reduction mipoint potential (Em) of isolated Rhodopseudomonas sphaeroides cytochrome c2 exhibits a pH-dependent behavior which can be ascribed to a pK on the oxidized form at pH 8.0 (Pettigrew et al. (1975) Biochim. Biophys. Acta 430, 197-208). However, as with mammalian cytochrome c (Brandt, K.G. Parks, P.C., Czerlinski, G.H. and Hess, G.P. (1966) J. Biol. Chem. 241, 4180-4185) this pK can more properly be attributed to the combination of a pK beyond pH 11, and a slow conformational change of the ferricytochrome. This has been demonstrated by resolving the Em of cytochrome c2 before and after the conformational change. The Em of the unaltered form is essentially pH independent between pH 7 and 11.5, and the lower equilibrium Em is due solely to the conformational change. In vivo the conformational change is prevented by the binding of the cytochrome c2 to the photochemical reaction center, and the cytochrome exhibits an essentially pH-independent Em from pH 5 to 11. The alkaline transition thus has little physiological significance, and it is unlikely that the redox reactions of cytochrome c2 in vivo involve protons.

Alkaline isomerization of thermoresistant cytochrome c-552 and horse heart cytochrome c studied by absorption and resonance Raman spectroscopy

The structure of the thermoresistant cytochrome c (552, Thermus thermophilus) has been investigated at neutral and alkaline pH by absorption and resonance Raman spectroscopy and compared with that of horse heart cytochrome c. The ligands of the ferricytochrome c-552 at neutral pH are considered to be histidine and methionine, whereas the ligands of ferrocytochrome c-552 are histidine and another nitrogen base, histidine or lysine. Ferric cytochrome c-552 undergoes an alkaline isomerization with a pK of 12.3 (25 degrees C), accompanied by a ligand exchange. Horse heart cytochrome c has at least three isomerization states at alkaline pH (pK 9.3, 12.9 and greater than 13.5 at 25 degrees C). The replacement of the sixth ligand may not be involved in the second isomerization. The thermodynamic parameters for the isomerization were also estimated. The entropy change upon isomerization of cytochrome c-552 is negative, whereas for that of horse heart cytochrome c the entropy change is positive.

Kinetic studies on redox reactions of hemoproteins. I. Reduction of thermoresistant cytochrome c-552 and horse heart cytochrome c by ferrocyanide

The oxidation-reduction reaction of horse heart cytochrome c and cytochrome c (552, Thermus thermophilus), which is highly thermoresistant, was studied by temperature-jump method. Ferrohexacyanide was used as reductant. (Formula: see text.) Thermodynamic and activation parameters of the reaction obtained for both cytochromes were compared with each other. The results of this showed that (1) the redox potential of cytochrome c-552, + 0.19 V, is markedly less than that of horse heart cytochrome c. (2) deltaHox of cytochrome c-552 is considerably lower than that of horse heart cytochrome c. (3) deltaSox and deltaSred of cytochrome c-552 are more negative than those of horse heart cytochrome c. (4) kred of cytochrome c-552 is much lower than that of horse heart cytochrome c at room temperature.