Conformational change and human cytochrome c function: mutation of residue 41 modulates caspase activation and destabilizes Met-80 coordination

Altered conformational dynamics contribute to species-specific effects of cytochrome c mutations on caspase activation

Variants in the gene encoding human cytochrome c (CYCS) cause mild autosomal dominant thrombocytopenia. Despite high sequence conservation between mouse and human cytochrome c, this phenotype is not recapitulated in mice for the sole mutant (G41S) that has been investigated. The effect of the G41S mutation on the in vitro activities of cytochrome c is also not conserved between human and mouse. Peroxidase activity is increased in both mouse and human G41S variants, whereas apoptosome activation is increased for human G41S cytochrome c but decreased for mouse G41S cytochrome c. These apoptotic activities of cytochrome c are regulated at least in part by conformational dynamics of the main chain. Here we use computational and in vitro approaches to understand why the impact of the G41S mutation differs between mouse and human cytochromes c. The G41S mutation increases the inherent entropy and main chain mobility of human but not mouse cytochrome c. Exclusively in human G41S cytochrome c this is accompanied by a decrease in occupancy of H-bonds between protein and heme during simulations. These data demonstrate that binding of cytochrome c to Apaf-1 to trigger apoptosome formation, but not the peroxidase activity of cytochrome c, is enhanced by increased mobility of the native protein conformation.

Effects of naturally occurring S47F/A mutations on the structure and function of human cytochrome c

The sequence and structure of human cytochrome c (hCyt c) exhibit evolutionary conservations, with only a limited number of naturally occurring mutations in humans. Herein, we investigated the effects of the naturally occurring S47F/A mutations on the structure and function of hCyt c in the oxidized form. Although the naturally occurring S47F/A mutations did not largely alter the protein structure, the S47F and S47A variants exhibited a small fraction of high-spin species. Kinetic studies showed that the peroxidase activity of the variants was enhanced by ∼2.5-fold under neutral pH conditions, as well as for the rate in reaction with H₂O₂, when compared to those of wild-type hCyt c. In addition, we evaluated the interaction between hCyt c and human neuroglobin (hNgb) by isothermal titration calorimetry (ITC) studies, which revealed that the binding constant was reduced by ∼8-fold as result of the mutation of the hydrophilic Ser to the hydrophobic Phe/Ala. These findings provide valuable insights into the role of Ser47 in Ω-loop C in sustaining the structure and function of hCyt c.

Molecular insights on the conformational dynamics of a P76C mutant of human cytochrome c and the enhancement on its peroxidase activity

In apoptotic pathway, the interaction of Cytochrome c (Cytc) with cardiolipin in vivo is a key process to induce peroxidase activity of Cytc and trigger the release of Cytc in the inner mitochondria into cytosol. The peroxidase active form of Cytc occurs due to local conformational changes that support the opening of the heme crevice and the loss of an axial ligand between Met80 and heme Fe. Structural adjustments at the Ω-loop segments of Cytc are required for such process. To study the role of the distal Ω-loop segments comprising residues 71-85 in human Cytc (hCytc), we investigated a cysteine mutation at Pro76, one of the highly conserved residues in this loop. The effect of P76C mutant was explored by the combination of experimental characterizations and molecular dynamics (MD) simulations. The peroxidase activity of the P76C mutant was found to be significantly increased by ∼13 folds relative to the wild type. Experimental data on global denaturation, alkaline transition, heme bleaching, and spin-labeling Electron Spin Resonance were in good agreement with the enhancement of peroxidase activity. The MD results of hCytc in the hexacoordinate form suggest the important changes in P76C mutant occurred due to the unfolding at the central Ω-loop (residues 40-57), and the weakening of H-bond between Tyr67 and Met80. Whereas the experimental data implied that the P76C mutant tend to be in equilibrium between the pentacoordinate and hexacoordinate forms, the MD and experimental information are complementary and were used to support the mechanisms of peroxidase active form of hCytc.

Correlated electric field modulation of electron transfer parameters and the access to alternative conformations of multifunctional cytochrome c

Cytochrome c (Cytc) is a multifunctional protein that, in its native conformation, shuttles electrons in the mitochondrial respiratory chain. Conformational transitions that involve replacement of the heme distal ligand lead to the gain of alternative peroxidase activity, which is crucial for membrane permeabilization during apoptosis. Using a time-resolved SERR spectroelectrochemical approach, we found that the key physicochemical parameters that characterize the electron transfer (ET) canonic function and those that determine the transition to alternative conformations are strongly correlated and are modulated by local electric fields (LEF) of biologically meaningful magnitude. The electron shuttling function is optimized at moderate LEFs of around 1 V nm^-1. A decrease of the LEF is detrimental for ET as it rises the reorganization energy. Moreover, LEF values below and above the optimal for ET favor alternative conformations with peroxidase activity and downshifted reduction potentials. The underlying proposed mechanism is the LEF modulation of the flexibility of crucial protein segments, which produces a differential effect on the kinetic ET and conformational parameters of Cytc. These findings might be related to variations in the mitochondrial membrane potential during apoptosis, as the basis for the switch between canonic and alternative functions of Cytc. Moreover, they highlight the possible role of variable LEFs in determining the function of other moonlighting proteins through modulation of the protein dynamics.

Activation of Cytochrome C Peroxidase Function Through Coordinated Foldon Loop Dynamics upon Interaction with Anionic Lipids

Cardiolipin (CL) is a mitochondrial anionic lipid that plays important roles in the regulation and signaling of mitochondrial apoptosis. CL peroxidation catalyzed by the assembly of CL-cytochrome c (cyt c) complexes at the inner mitochondrial membrane is a critical checkpoint. The structural changes in the protein, associated with peroxidase activation by CL and different anionic lipids, are not known at a molecular level. To better understand these peripheral protein-lipid interactions, we compare how phosphatidylglycerol (PG) and CL lipids trigger cyt c peroxidase activation, and correlate functional differences to structural and motional changes in membrane-associated cyt c. Structural and motional studies of the bound protein are enabled by magic angle spinning solid state NMR spectroscopy, while lipid peroxidase activity is assayed by mass spectrometry. PG binding results in a surface-bound state that preserves a nativelike fold, which nonetheless allows for significant peroxidase activity, though at a lower level than binding its native substrate CL. Lipid-specific differences in peroxidase activation are found to correlate to corresponding differences in lipid-induced protein mobility, affecting specific protein segments. The dynamics of omega loops C and D are upregulated by CL binding, in a way that is remarkably controlled by the protein:lipid stoichiometry. In contrast to complete chemical denaturation, membrane-induced protein destabilization reflects a destabilization of select cyt c foldons, while the energetically most stable helices are preserved. Our studies illuminate the interplay of protein and lipid dynamics in the creation of lipid peroxidase-active proteolipid complexes implicated in early stages of mitochondrial 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.

Altered structure and dynamics of pathogenic cytochrome c variants correlate with increased apoptotic activity

Mutation of cytochrome c in humans causes mild autosomal dominant thrombocytopenia. The role of cytochrome c in platelet formation, and the molecular mechanism underlying the association of cytochrome c mutations with thrombocytopenia remains unknown, although a gain-of-function is most likely. Cytochrome c contributes to several cellular processes, with an exchange between conformational states proposed to regulate changes in function. Here, we use experimental and computational approaches to determine whether pathogenic variants share changes in structure and function, and to understand how these changes might occur. Three pathogenic variants (G41S, Y48H, A51V) cause an increase in apoptosome activation and peroxidase activity. Molecular dynamics simulations of these variants, and two non-naturally occurring variants (G41A, G41T), indicate that increased apoptosome activation correlates with the increased overall flexibility of cytochrome c, particularly movement of the Ω loops. Crystal structures of Y48H and G41T complement these studies which overall suggest that the binding of cytochrome c to apoptotic protease activating factor-1 (Apaf-1) may involve an 'induced fit' mechanism which is enhanced in the more conformationally mobile variants. In contrast, peroxidase activity did not significantly correlate with protein dynamics. Thus, the mechanism by which the variants increase peroxidase activity is not related to the conformational dynamics of the native hexacoordinate state of cytochrome c. Recent molecular dynamics data proposing conformational mobility of specific cytochrome c regions underpins changes in reduction potential and alkaline transition pK was not fully supported. These data highlight that conformational dynamics of cytochrome c drive some but not all of its properties and activities.

The importance of Asn52 in the structure-function relationship of human cytochrome c

The function of the highly conserved residue Asn52 in human cytochrome c (H-Cyt c) is not fully understood. Herein, we show that the naturally occurring variant N52S H-Cyt c has a perturbed secondary structure, with a small fraction of high-spin species. Remarkably, it exhibits an enhanced peroxidase activity by 3-8-fold at neutral pH, as well as self-oxidation in reaction with H₂O₂. This study suggests that the H-bond network mediated by Asn52 is essential to suppress the apoptotic activity of H-Cyt c under physiological conditions.

The Human Cytochrome c Domain-Swapped Dimer Facilitates Tight Regulation of Intrinsic Apoptosis

Oxidation of cardiolipin (CL) by cytochrome c (cytc) has been proposed to initiate the intrinsic pathway of apoptosis. Domain-swapped dimer (DSD) conformations of cytc have been reported both by our laboratory and by others. The DSD is an alternate conformer of cytc that could oxygenate CL early in apoptosis. We demonstrate here that the cytc DSD has a set of properties that would provide tighter regulation of the intrinsic pathway. We show that the human DSD is kinetically more stable than horse and yeast DSDs. Circular dichroism data indicate that the DSD has a less asymmetric heme environment, similar to that seen when the monomeric protein binds to CL vesicles at high lipid:protein ratios. The dimer undergoes the alkaline conformational transition near pH 7.0, 2.5 pH units lower than that of the monomer. Data from fluorescence correlation spectroscopy and fluorescence anisotropy suggest that the alkaline transition of the DSD may act as a switch from a high affinity for CL nanodiscs at pH 7.4 to a much lower affinity at pH 8.0. Additionally, the peroxidase activity of the human DSD increases 7-fold compared to that of the monomer at pH 7 and 8, but by 14-fold at pH 6 when mixed Met80/H2O ligation replaces the lysine ligation of the alkaline state. We also present data that indicate that cytc binding shows a cooperative effect as the concentration of cytc is increased. The DSD appears to have evolved into a pH-inducible switch that provides a means to control activation of apoptosis near pH 7.0.

Distal Unfolding of Ferricytochrome c Induced by the F82K Mutation

It is well known that axial coordination of heme iron in mitochondrial cytochrome c has redox-dependent stability. The Met80 heme iron axial ligand in the ferric form of the protein is relatively labile and can be easily replaced by alternative amino acid side chains under non-native conditions induced by alkaline pH, high temperature, or denaturing agents. Here, we showed a redox-dependent destabilization induced in human cytochrome c by substituting Phe82-conserved amino acid and a key actor in cytochrome c intermolecular interactions-with a Lys residue. Introducing a positive charge at position 82 did not significantly affect the structure of ferrous cytochrome c but caused localized unfolding of the distal site in the ferric state. As revealed by ¹H NMR fingerprint, the ferric form of the F82K variant had axial coordination resembling the renowned alkaline species, where the detachment of the native Met80 ligand favored the formation of multiple conformations involving distal Lys residues binding to iron, but with more limited overall structural destabilization.

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.

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.

The proportion of Met80-sulfoxide dictates peroxidase activity of human cytochrome c

The peroxidase activity of cytochrome c is proposed to contribute to apoptosis by peroxidation of cardiolipin in the mitochondrial inner membrane. However, cytochrome c heme is hexa-coordinate with a methionine (Met80) on the distal side, stopping it from acting as an efficient peroxidase. The first naturally occurring variant of cytochrome c discovered, G41S, has higher peroxidase activity than wild-type. To understand the basis for this increase and gain insight into the peroxidase activity of wild-type, we have studied wild-type, G41S and the unnatural variant G41T. Through a combined kinetic and mass spectrometric analysis, we have shown that hydrogen peroxide specifically oxidizes Met80 to the sulfoxide. In the absence of substrate this can be further oxidized to the sulfone, leading to a decrease in peroxidase activity. Peroxidase activity can be correlated with the proportion of sulfoxide present and if fully in that form, all variants have the same activity without a lag phase caused by activation of the protein.

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.

X-ray structure of bovine heart cytochrome c at high ionic strength

Bovine heart cytochrome c (bCyt c) is an extensively studied hemoprotein of only 104 residues. Due to the existence of isoforms generated by non-enzymatic deaminidation, crystallization of bCyt c is difficult and involves extensive purification and the use of microseeding or the presence of an electric field. Taking advantage of the capacity of cytochrome c (cyt c) to bind anions on its protein surface, the commercially available bCyt c was crystallized without extra purifications, using ammonium sulfate as precipitant and nitrate ions as additives. The structure of the ferric bCyt c in a new crystal form is described and compared with that previously solved at low ionic strength and with those of human and horse cyt c. The overall structure of bCyt c is conserved, while the side chains of several residues that play a role in the interactions of cyt c with its partners have different rotamers in the two structures. The effect of the presence of nitrate ions on the structure of the protein is then evaluated and compared with that observed in the case of ferrous and ferric horse heart cyt c.

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.

Investigation of structural dynamics of Thrombocytopenia Cargeeg mutants of human apoptotic cytochrome c: A molecular dynamics simulation approach

Naturally occurring mutations to cytochrome c (cyt-c) have been identified recently in patients with mild autosomal dominant thrombocytopenia (low platelet levels), which yield cyt-c mutants with enhanced apoptotic activity. However, the molecular mechanism underlying this low platelet production and enhanced apoptosis remain unclear. Therefore, an attempt is made herein for the first time to investigate the effects of mutations of glycine 41 by serine (G41S) and tyrosine 48 by histidine (Y48H) on the conformational and dynamic changes of apoptotic (Fe³⁺) cyt-c using all atom molecular dynamics (MD) simulations in explicit water solvent. Our 30ns MD simulations demonstrate considerable structural differences in G41S and Y48H compared to wild type (WT) cyt-c, such as increasing distances between the critical electron transfer residues results in open conformation at the heme active site, large fluctuations in β-turns and α-helices. Additionally, although the β-sheets remain mostly unaffected in all the three cyt-c simulations, the α-helices undergo conformational switch to β-turns in both the mutant simulations. Importantly, this conformational switch of α-helix to β-turn around heme active site should attributes to the loss of intraprotein H-bonds in the mutant simulations especially between NE2 (His26) and O (Pro44) in agreement with the experimental report. Further, essential dynamics analysis reveals that overall motions of WT cyt-c is mainly involved only in the first eigenvector, but in G41S and Y48H the overall motions are mainly in three and two eigenvectors respectively. Overall, the detailed atomistic level information provide a unifying description for the molecular mechanism of structural destabilization, disregulation of platelet formation and enhanced peroxidase activity of the mutant cyt-c's in the pathology of intrinsic apoptosis.

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.

Remote Perturbations in Tertiary Contacts Trigger Ligation of Lysine to the Heme Iron in Cytochrome c

Perturbations in protein structure define the mechanism of allosteric regulation and biological information transfer. In cytochrome c (cyt c), ligation of Met80 to the heme iron is critical for the protein's electron-transfer (ET) function in oxidative phosphorylation and for suppressing its peroxidase activity in apoptosis. The hard base Lys is a better match for the hard ferric iron than the soft base Met is, suggesting the key role of the protein scaffold in favoring Met ligation. To probe the role of the protein structure in the maintenance of Met ligation, mutations T49V and Y67R/M80A were designed to disrupt hydrogen bonding and packing of the heme coordination loop, respectively. Electronic absorption, nuclear magnetic resonance, and electron paramagnetic resonance spectra reveal that ferric forms of both variants are Lys-ligated at neutral pH. A minor change in the tertiary contacts in T49V, away from the heme coordination loop, appears to be sufficient to execute a change in ligation, suggesting a cross-talk between the different regions of the protein structure and a possibility of built-in conformational switches in cyt c. Analyses of thermodynamic stability, kinetics of Lys binding and dissociation, and the pH-dependent changes in ligation provide a detailed characterization of the Lys coordination in these variants and relate these properties to the extent of structural perturbations. The findings emphasize the importance of the hydrogen-bonding network in controlling ligation of the native Met80 to the heme iron.

Increased dynamics in the 40-57 Ω-loop of the G41S variant of human cytochrome c promote its pro-apoptotic conformation

Thrombocytopenia 4 is an inherited autosomal dominant thrombocytopenia, which occurs due to mutations in the human gene for cytochrome c that results in enhanced mitochondrial apoptotic activity. The Gly41Ser mutation was the first to be reported. Here we report stopped-flow kinetic studies of azide binding to human ferricytochrome c and its Gly41Ser variant, together with backbone amide H/D exchange and ¹⁵N-relaxation dynamics using NMR spectroscopy, to show that alternative conformations are kinetically and thermodynamically more readily accessible for the Gly41Ser variant than for the wild-type protein. Our work reveals a direct conformational link between the 40–57 Ω-loop in which residue 41 resides and the dynamical properties of the axial ligand to the heme iron, Met80, such that the replacement of glycine by serine promotes the dissociation of the Met80 ligand, thereby increasing the population of a peroxidase active state, which is a key non-native conformational state in apoptosis.

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.

Structure-function relationships in human cytochrome c: The role of tyrosine 67

Spectroscopic and functional properties of human cytochrome c and its Tyr67 residue mutants (i.e., Tyr67His and Tyr67Arg) have been investigated. In the case of the Tyr67His mutant, we have observed only a very limited structural alteration of the heme pocket and of the Ω-loop involving, among others, the residue Met80 and its bond with the heme iron. Conversely, in the Tyr67Arg mutant the Fe-Met80 bond is cleaved; consequently, a much more extensive structural alteration of the Ω-loop can be envisaged. The structural, and thus the functional modifications, of the Tyr67Arg mutant are present in both the ferric [Fe(III)] and the ferrous [Fe(II)] forms, indicating that the structural changes are independent of the heme iron oxidation state, depending instead on the type of substituting residue. Furthermore, a significant peroxidase activity is evident for the Tyr67Arg mutant, highlighting the role of Arg as a basic, positively charged residue at pH7.0, located in the heme distal pocket, which may act as an acid to cleave the O-O bond in H2O2. As a whole, our results indicate that a delicate equilibrium is associated with the spatial arrangement of the Ω-loop. Clearly, Arg, but not His, is able to stabilize and polarize the negative charge on the Fe(III)-OOH complex during the formation of Compound I, with important consequences on cytochrome peroxidation activity and its role in the apoptotic process, which is somewhat different in yeast and mammals.

Interspecies Variation in the Functional Consequences of Mutation of Cytochrome c

The naturally occurring human cytochrome c variant (G41S) is associated with a mild autosomal dominant thrombocytopenia (Thrombocytopenia Cargeeg) caused by dysregulation of platelet production. The molecular basis of the platelet production defect is unknown. Despite high conservation of cytochrome c between human and mouse (91.4% identity), introducing the G41S mutation into mouse cytochrome c in a knockin mouse (Cycs^G41S/G41S) did not recapitulate the low platelet phenotype of Thrombocytopenia Cargeeg. While investigating the cause of this disparity we found a lack of conservation of the functional impact of cytochrome c mutations on caspase activation across species. Mutation of cytochrome c at residue 41 has distinct effects on the ability of cytochrome c to activate caspases depending on the species of both the cytochrome c and its binding partner Apaf-1. In contrast to our previous results showing the G41S mutation increases the ability of human cytochrome c to activate caspases, here we find this activity is decreased in mouse G41S cytochrome c. Additionally unlike wildtype human cytochrome c, G41S cytochrome c is unable to activate caspases in Xenopus embryo extracts. Taken together these results demonstrate a previously unreported species-specific component to the interaction of cytochrome c with Apaf-1. This suggests that the electrostatic interaction between cytochrome c and Apaf-1 is not the sole determinant of binding, with additional factors controlling binding specificity and affinity. These results have important implications for studies of the effects of cytochrome c mutations on the intrinsic apoptosis pathway.

Disruption of cytochrome c heme coordination is responsible for mitochondrial injury during ischemia

BACKGROUND

It was recently suggested that electron flow into cyt c, coupled with ROS generation, oxidizes cyt c Met(80) to Met(80) sulfoxide (Met-O) in isolated hearts after ischemia-reperfusion, and converts cyt c to a peroxidase. We hypothesize that ischemia disrupts Met(80)-Fe ligation of cyt c, forming pentacoordinated heme Fe(2+), which inhibits electron transport (ET) and promotes oxygenase activity.

METHODS

SS-20 (Phe-D-Arg-Phe-Lys-NH2) was used to demonstrate the role of Met(80)-Fe ligation in ischemia. Mitochondria were isolated from ischemic rat kidneys to determine sites of respiratory inhibition. Mitochondrial cyt c and cyt c Met-O were quantified by western blot, and cristae architecture was examined by electron microscopy.

RESULTS

Biochemical and structural studies showed that SS-20 selectively targets cardiolipin (CL) and protects Met(80)-Fe ligation in cyt c. Ischemic mitochondria showed 17-fold increase in Met-O cyt c, and dramatic cristaeolysis. Loss of cyt c was associated with proteolytic degradation of OPA1. Ischemia significantly inhibited ET initiated by direct reduction of cyt c and coupled respiration. All changes were prevented by SS-20.

CONCLUSION

Our results show that ischemia disrupts the Met(80)-Fe ligation of cyt c resulting in the formation of a globin-like pentacoordinated heme Fe(2+) that inhibits ET, and converts cyt c into an oxygenase to cause CL peroxidation and proteolytic degradation of OPA1, resulting in cyt c release.

GENERAL SIGNIFICANCE

Cyt c heme structure represents a novel target for minimizing ischemic injury. SS-20, which we show to selectively target CL and protect the Met(80)-Fe ligation, minimizes ischemic injury and promotes ATP recovery.

The hydrogen-peroxide-induced radical behaviour in human cytochrome c-phospholipid complexes: implications for the enhanced pro-apoptotic activity of the G41S mutant

We have investigated whether the pro-apoptotic properties of the G41S mutant of human cytochrome c can be explained by a higher than wild-type peroxidase activity triggered by phospholipid binding. A key complex in mitochondrial apoptosis involves cytochrome c and the phospholipid cardiolipin. In this complex cytochrome c has its native axial Met(80) ligand dissociated from the haem-iron, considerably augmenting the peroxidase capability of the haem group upon H2O2 binding. By EPR spectroscopy we reveal that the magnitude of changes in the paramagnetic haem states, as well as the yield of protein-bound free radical, is dependent on the phospholipid used and is considerably greater in the G41S mutant. A high-resolution X-ray crystal structure of human cytochrome c was determined and, in combination with the radical EPR signal analysis, two tyrosine residues, Tyr(46) and Tyr(48), have been rationalized to be putative radical sites. Subsequent single and double tyrosine-to-phenylalanine mutations revealed that the EPR signal of the radical, found to be similar in all variants, including G41S and wild-type, originates not from a single tyrosine residue, but is instead a superimposition of multiple EPR signals from different radical sites. We propose a mechanism of multiple radical formations in the cytochrome c-phospholipid complexes under H2O2 treatment, consistent with the stabilization of the radical in the G41S mutant, which elicits a greater peroxidase activity from cytochrome c and thus has implications in mitochondrial apoptosis.

Methionine ligand lability of homologous monoheme cytochromes c

Direct electrochemical analysis of adsorbed bacterial monoheme cytochromes c has revealed a phenomenological loss of the axial methionine when examined using pyrolytic "edge-plane" graphite (EPG) electrodes. While prior findings have reported that the Met-loss state may be quantitatively understood using the cytochrome c from Hydrogenobacter thermophilus as a model system, here we demonstrate that the formation of the Met-loss state upon EPG electrodes can be observed for a range of cytochrome orthologs. Through an electrochemical comparison of the wild-type proteins from organisms of varying growth temperature optima, we establish that Met-ligand losses at graphite surfaces have similar energetics to the "foldons" for known protein folding pathways. Furthermore, a downward shift in reduction potential to approximately -100 mV vs standard hydrogen electrode was observed, similar to that of the alkaline transition found in mitochondrial cytochromes c. Pourbaix diagrams for the Met-loss forms of each cytochrome, considered here in comparison to mutants where the Met-ligand has been substituted to His or Ala, suggest that the nature of the Met-loss state is distinct from either a His-/aquo- or a bis-His-ligated heme center, yet more closely matches the pKa values found for bis-His-ligated hemes., We find the propensity for adoption of the Met-loss state in bacterial monoheme cytochromes c scales with their overall thermal stability, though not with the specific stability of the Fe-Met bond.

Discovering co-occurring patterns and their biological significance in protein families

Background

The large influx of biological sequences poses the importance of identifying and correlating conserved regions in homologous sequences to acquire valuable biological knowledge. These conserved regions contain statistically significant residue associations as sequence patterns. Thus, patterns from two conserved regions co-occurring frequently on the same sequences are inferred to have joint functionality. A method for finding conserved regions in protein families with frequent co-occurrence patterns is proposed. The biological significance of the discovered clusters of conserved regions with co-occurrences patterns can be validated by their three-dimensional closeness of amino acids and the biological functionality found in those regions as supported by published work.

Methods

Using existing algorithms, we discovered statistically significant amino acid associations as sequence patterns. We then aligned and clustered them into Aligned Pattern Clusters (APCs) corresponding to conserved regions with amino acid conservation and variation. When one APC frequently co-occured with another APC, the two APCs have high co-occurrence. We then clustered APCs with high co-occurrence into what we refer to as Co-occurrence APC Clusters (Co-occurrence Clusters).

Results

Our results show that for Co-occurrence Clusters, the three-dimensional distance between their amino acids is closer than average amino acid distances. For the Co-occurrence Clusters of the ubiquitin and the cytochrome c families, we observed biological significance among the residing amino acids of the APCs within the same cluster. In ubiquitin, the residues are responsible for ubiquitination as well as conventional and unconventional ubiquitin-bindings. In cytochrome c, amino acids in the first co-occurrence cluster contribute to binding of other proteins in the electron transport chain, and amino acids in the second co-occurrence cluster contribute to the stability of the axial heme ligand.

Conclusions

Thus, our co-occurrence clustering algorithm can efficiently find and rank conserved regions that contain patterns that frequently co-occurring on the same proteins. Co-occurring patterns are biologically significant due to their three-dimensional closeness and other evidences reported in literature. These results play an important role in drug discovery as biologists can quickly identify the target for drugs to conduct detailed preclinical studies.

Enhancing the peroxidase activity of cytochrome c by mutation of residue 41: implications for the peroxidase mechanism and cytochrome c release

Enzymatic and mitochondrial release analyses of residue 41 variants of native and cardiolipin-bound mammalian cytochromes c implicate mobility of the 40–57 Ω loop in induction of peroxidase activity, but show that loss of axial Fe co-ordination is not essential.