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

Comparing Properties of Common Bioinorganic Ligands with Switchable Variants of Cytochrome c

Ligand substitution at the metal center is common in catalysis and signal transduction of metalloproteins. Understanding the effects of particular ligands, as well as the polypeptide surrounding, is critical for uncovering mechanisms of these biological processes and exploiting them in the design of bioinspired catalysts and molecular devices. A series of switchable K79G/M80X/F82C (X = Met, His, or Lys) variants of cytochrome (cyt) c was employed to directly compare the stability of differently ligated proteins and activation barriers for Met, His, and Lys replacement at the ferric heme iron. Studies of these variants and their nonswitchable counterparts K79G/M80X have revealed stability trends Met < Lys < His and Lys < His < Met for the protein Fe^III-X and Fe^II-X species, respectively. The differences in the hydrogen-bonding interactions in folded proteins and in solvation of unbound X in the unfolded proteins explain these trends. Calculations of free energy of ligand dissociation in small heme model complexes reveal that the ease of the Fe^III-X bond breaking increases in the series amine < imidazole < thioether, mirroring trends in hardness of these ligands. Experimental rate constants for X dissociation in differently ligated cyt c variants are consistent with this sequence, but the differences between Met and His dissociation rates are attenuated because the former process is limited by the heme crevice opening. Analyses of activation parameters and comparisons to those for the Lys-to-Met ligand switch in the alkaline transition suggest that ligand dissociation is entropically driven in all the variants and accompanied by Lys protonation at neutral pH. The described thiolate redox-linked switches have offered a wealth of new information about interactions of different protein-derived ligands with the heme iron in cyt c model proteins, and we anticipate that the strategy of employing these switches could benefit studies of other redox metalloproteins and model complexes.

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.

Mechanistic Studies of Proton-Coupled Electron Transfer in a Calorimetry Cell

Mechanistic studies of proton-coupled electron-transfer (PCET) reactions in proteins are complicated by the challenge of following proton transfer (PT) in these large molecules. Herein we describe the use of isothermal titration calorimetry (ITC) to establish proton involvement in protein redox reactions and the identity of PT sites. We validate this approach with three variants of a heme protein cytochrome c (cyt c) and show that the method yields a wealth of thermodynamic information that is important for characterizing PCET reactions, including reduction potentials, redox-dependent p K_a values, and reaction enthalpies for both electron-transfer (ET) and PT steps. We anticipate that this facile and label-free ITC approach will find widespread applications in studies of other redox proteins and enhance our knowledge of PCET reaction mechanisms.

Ligation and Reactivity of Methionine-Oxidized Cytochrome c

Met80, one of the heme iron ligands in cytochrome c (cyt c), is readily oxidized to Met sulfoxide (Met-SO) by several biologically relevant oxidants. The modification has been suggested to affect both the electron-transfer (ET) and apoptotic functions of this metalloprotein. The coordination of the heme iron in Met-oxidized cyt c (Met-SO cyt c) is critical for both of these functions but has remained poorly defined. We present electronic absorption, NMR, and EPR spectroscopic investigations as well as kinetic studies and mutational analyses to identify the heme iron ligands in yeast iso-1 Met-SO cyt c. Similar to the alkaline form of native cyt c, Lys73 and Lys79 ligate to the ferric heme iron in the Met80-oxidized protein, but this coordination takes place at much lower pH. The ferrous heme iron is ligated by Met-SO, implying the redox-linked ligand switch in the modified protein. Binding studies with the model peptide microperoxidase-8 provide a rationale for alterations in ligation and for the role of the polypeptide packing in native and Met-SO cyt c. Imidazole binding experiments have revealed that Lys dissociation from the ferric heme in K73A/K79G/M80K (M80K#) and Met-SO is more than 3 orders of magnitude slower than the opening of the heme pocket that limits Met80 replacement in native cyt c. The Lys-to-Met-SO ligand substitution gates ET of ferric Met-SO cyt c with Co(terpy)22+. Owing to the slow Lys dissociation step, ET reaction is slow but possible, which is not the case for nonswitchable M80A and M80K#. Acidic conditions cause Lys replacement by a water ligand in Met-SO cyt c (p Ka = 6.3 ± 0.1), increasing the intrinsic peroxidase activity of the protein. This pH-driven ligand switch may be a mechanism to boost peroxidase function of cyt c specifically in apoptotic cells.

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.

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.

Structure of a mitochondrial cytochrome c conformer competent for peroxidase activity

At the onset of apoptosis, the peroxidation of cardiolipin at the inner mitochondrial membrane by cytochrome c requires an open coordination site on the heme. We report a 1.45-Å resolution structure of yeast iso-1-cytochrome c with the Met80 heme ligand swung out of the heme crevice and replaced by a water molecule. This conformational change requires modest adjustments to the main chain of the heme crevice loop and is facilitated by a trimethyllysine 72-to-alanine mutation. This mutation also enhances the peroxidase activity of iso-1-cytochrome c. The structure shows a buried water channel capable of facilitating peroxide access to the active site and of moving protons produced during peroxidase activity to the protein surface. Alternate positions of the side chain of Arg38 appear to mediate opening and closing of the buried water channel. In addition, two buried water molecules can adopt alternate positions that change the network of hydrogen bonds in the buried water channel. Taken together, these observations suggest that low and high proton conductivity states may mediate peroxidase function. Comparison of yeast and mammalian cytochrome c sequences, in the context of the steric factors that permit opening of the heme crevice, suggests that higher organisms have evolved to inhibit peroxidase activity, providing a more stringent barrier to the onset of apoptosis.