Heightened Dynamics of the Oxidized Y48H Variant of Human Cytochrome c Increases Its Peroxidatic Activity

A computational study of the fold and stability of cytochrome c with implications for disease

Cytochrome c (Cyt-c), encoded by the CYCS gene, is crucial for electron transport, peroxidase activity, and apoptosis. Mutations in CYCS cause thrombocytopenia 4 (THC4), a disorder with low platelet counts. We have, for instance, recently described six Italian families with five different heterozygous missense CYCS variants. These mutations likely enhance peroxidase and apoptotic activities, yet the mechanisms causing reduced platelet production and increased apoptosis are unclear. This study investigates clinically-related Cyt-c variants using an integrated bioinformatics approach. Our findings reveal that all variants are at evolutionarily conserved sites, potentially disrupting Cyt-c function and contributing to disease phenotypes. Specific variants are predicted to affect phosphorylation (T20I, V21G, Y49H), and ubiquitination (G42S, A52T, A52V, T103I). Molecular dynamics simulations (500 ns) revealed significant structural deviations from the wild-type protein, with mutants showing reduced stability and increased unfolding and flexibility, particularly in the Ω-loops. These changes result in the displacement of the Ω-loops away from the heme iron, weakening critical hydrogen bonds and consequently opening the heme active site. This open conformation may enhance accessibility to small molecules such as H₂O₂, thereby promoting peroxidase activity, which may enhance apoptosis and likely impact megakaryopoiesis and platelet homeostasis in THC4.

Phosphorylation of cytochrome c at tyrosine 48 finely regulates its binding to the histone chaperone SET/TAF-Iβ in the nucleus

Post-translational modifications (PTMs) of proteins are ubiquitous processes present in all life kingdoms, involved in the regulation of protein stability, subcellular location and activity. In this context, cytochrome c (Cc) is an excellent case study to analyze the structural and functional changes induced by PTMS as Cc is a small, moonlighting protein playing different roles in different cell compartments at different cell-cycle stages. Cc is actually a key component of the mitochondrial electron transport chain (ETC) under homeostatic conditions but is translocated to the cytoplasm and even the nucleus under apoptotic conditions and/or DNA damage. Phosphorylation does specifically alter the Cc redox activity in the mitochondria and the Cc non-redox interaction with apoptosis-related targets in the cytoplasm. However, little is known on how phosphorylation alters the interaction of Cc with histone chaperones in the nucleus. Here, we report the effect of Cc Tyr48 phosphorylation by examining the protein interaction with SET/TAF-Iβ in the nuclear compartment using a combination of molecular dynamics simulations, biophysical and structural approaches such as isothermal titration calorimetry (ITC) and nuclear magnetic resonance (NMR) and in cell proximity ligation assays. From these experiments, we infer that Tyr48 phosphorylation allows a fine-tuning of the Cc-mediated inhibition of SET/TAF-Iβ histone chaperone activity in vitro. Our findings likewise reveal that phosphorylation impacts the nuclear, stress-responsive functions of Cc, and provide an experimental framework to explore novel aspects of Cc post-translational regulation in the nucleus.

Human cytochrome C natural variants: Studying the membrane binding properties of G41S and Y48H by fluorescence energy transfer and molecular dynamics

Cytochrome C (cyt C), the protein involved in oxidative phosphorylation, plays several other crucial roles necessary for both cell life and death. Studying natural variants of cyt C offers the possibility to better characterize the structure-to-function relationship that modulates the different activities of this protein. Naturally mutations in human cyt C (G41S and Y48H) occur in the protein central Ω-loop and cause thrombocytopenia 4. In this study, we have investigated the binding of such variants and of wild type (wt) cyt C to synthetic cardiolipin-containing vesicles. The mutants have a lower propensity in membrane binding, displaying higher dissociation constants with respect to the wt protein. Compressibility measurements reveal that both variants are more flexible than the wt, suggesting that the native central Ω-loop is important for the interaction with membranes. Such hypothesis is supported by molecular dynamics simulations. A minimal distance analysis indicates that in the presence of cardiolipin the central Ω-loop of the mutants is no more in contact with the membrane, as it happens instead in the case of wt cyt C. Such finding might provide a hint for the reduced membrane binding capacity of the variants and their enhanced peroxidase activity in vivo.

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.

The Y49H cytochrome c variant enhances megakaryocytic maturation of K-562 cells

Five pathogenic variants in the gene encoding cytochrome c (CYCS) associated with mild autosomal dominant thrombocytopenia have been reported. Previous studies of peripheral blood CD34+ or CD45+ cells from subjects with the G42S CYCS variant showed an acceleration in megakaryopoiesis compared to wild-type (WT) cells. To determine whether this result reflects a common feature of the CYCS variants, the c.145T>C mutation (Y49H variant) was introduced into the endogenous CYCS locus in K-562 cells, which undergo megakaryocytic maturation in response to treatment with a phorbol ester. The c.145T>C (Y49H) variant enhanced the megakaryocyte maturation of the K-562 cells, and this effect was seen when the cells were cultured at both 18 % and 5 % oxygen. Thus, alteration of megakaryopoiesis is common to both the G42S and Y49H CYCS variants and may contribute to the low platelet phenotype. The Y49H CYCS variant has previously been reported to impair mitochondrial respiratory chain function in vitro, however using extracellular flux analysis the c.145T>C (Y49H) variant does not alter mitochondrial bioenergetics of the K-562 cells, consistent with the lack of a phenotype characteristic of mitochondrial diseases in CYCS variant families. The Y49H variant has also been reported to enhance the ability of cytochrome c to trigger caspase activation in the intrinsic apoptosis pathway. However, as seen in peripheral blood cells from G42S CYCS variant carriers, the presence of Y49H cytochrome c in K-562 cells did not significantly change their response to an apoptotic stimulus.

Phosphorylations and Acetylations of Cytochrome c Control Mitochondrial Respiration, Mitochondrial Membrane Potential, Energy, ROS, and Apoptosis

Cytochrome c (Cytc) has both life-sustaining and cellular death-related functions, depending on subcellular localization. Within mitochondria, Cytc acts as a single electron carrier as part of the electron transport chain (ETC). When released into the cytosol after cellular insult, Cytc triggers the assembly of the apoptosome, committing the cell to intrinsic apoptosis. Due to these dual natures, Cytc requires strong regulation by the cell, including post-translational modifications, such as phosphorylation and acetylation. Six phosphorylation sites and three acetylation sites have been detected on Cytc in vivo. Phosphorylations at T28, S47, Y48, T49, T58, and Y97 tend to be present under basal conditions in a tissue-specific manner. In contrast, the acetylations at K8, K39, and K53 tend to be present in specific pathophysiological conditions. All of the phosphorylation sites and two of the three acetylation sites partially inhibit respiration, which we propose serves to maintain an optimal, intermediate mitochondrial membrane potential (ΔΨm) to minimize reactive oxygen species (ROS) production. Cytc phosphorylations are lost during ischemia, which drives ETC hyperactivity and ΔΨm hyperpolarization, resulting in exponential ROS production thus causing reperfusion injury following ischemia. One of the acetylation sites, K39, shows a unique behavior in that it is gained during ischemia, stimulating respiration while blocking apoptosis, demonstrating that skeletal muscle, which is particularly resilient to ischemia-reperfusion injury compared to other organs, possesses a different metabolic strategy to handle ischemic stress. The regulation of Cytc by these post-translational modifications underscores the importance of Cytc for the ETC, ΔΨm, ROS production, apoptosis, and the cell as a whole.

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.

Structural Dynamics of the Precatalytic State of Human Cytochrome c upon T28C, G34C, and A50C Mutations: A Molecular Dynamics Simulation Perspective

The native structure of cytochrome c (cytc) contains hexacoordinate heme iron with His18 and Met80 residues ligated at the axial sites. Mutations of cytc at Ω-loops have been investigated in modulating the peroxidase activity and, hence, related to the initiation of the apoptotic pathway. Our previous experimental data reported on the peroxidase activity of the cysteine-directed mutants at different parts of the Ω-loop of human cytc (hCytc), that is, T28C, G34C, and A50C. In this work, we performed 1 μs molecular dynamics (MD) simulations to elucidate the detailed structural and dynamic changes upon these mutations, particularly at the proximal Ω-loop. The structures of hCytc were modeled in the hexacoordinated form, which was referred to as the "precatalytic state". The results showed that the structural features of the G34C mutant were more distinctive than those of other mutants. G34C mutation caused local destabilization and flexibility at the proximal Ω-loop (residues 12-28) and an extended distance between this Ω-loop region and heme iron. Besides, analysis of the orientation of the Arg38 side chain of the G34C mutant revealed the Arg38 conformer facing away from the heme iron. The obtained MD results also suggested structural diversity of the precatalytic states for the three hCytc mutants, specifically the effect of G34C mutation on the flexibility of the proximal Ω-loops. Therefore, our MD simulations combined with previous experimental data provide detailed insights into the structural basis of hCytc that could contribute to its pro-apoptotic function.

Alkaline State of the Domain-Swapped Dimer of Human Cytochrome c: A Conformational Switch for Apoptotic Peroxidase Activity

A 2.08 Å structure of an alkaline conformer of the domain-swapped dimer of K72A human cytochrome c (Cytc) crystallized at pH 9.9 is presented. In the structure, Lys79 is ligated to the heme. All other domain-swapped dimer structures of Cytc have water bound to this coordination site. Part of Ω-loop D (residues 70-85) forms a flexible linker between the subunits in other Cytc domain-swapped dimer structures but instead converts to a helix in the alkaline conformer of the dimer combining with the C-terminal helix to form two 26-residue helices that bracket both sides of the dimer. The alkaline transition of the K72A human dimer monitored at both 625 nm (high spin heme) and 695 nm (Met80 ligation) yields midpoint pH values of 6.6 and 7.6, respectively, showing that the Met80 → Lys79 and high spin to low spin transitions are distinct. The dimer peroxidase activity increases rapidly below pH 7, suggesting that population of the high spin form of the heme is what promotes peroxidase activity. Comparison of the structures of the alkaline dimer and the neutral pH dimer shows that the neutral pH conformer has a better electrostatic surface for binding to a cardiolipin-containing membrane and provides better access for small molecules to the heme iron. Given that the pH of mitochondrial cristae ranges from 6.9 to 7.2, the alkaline transition of the Cytc dimer could provide a conformational switch to tune the peroxidase activity of Cytc that oxygenates cardiolipin in the early stages of apoptosis.

Electric field-induced functional changes in electrode-immobilized mutant species of human cytochrome c

Post-translational modifications and naturally occurring mutations of cytochrome c have been recognized as a regulatory mechanism to control its biology. In this work, we investigate the effect of such in vivo chemical modifications of human cytochrome c on its redox properties in the adsorbed state onto an electrode. In particular, tyrosines 48 and 97 have been replaced by the non-canonical amino acid p-carboxymethyl-L-phenylalanine (pCMF), thus mimicking tyrosine phosphorylation. Additionally, tyrosine 48 has been replaced by a histidine producing the natural Y48H pathogenic mutant. Thermodynamics and kinetics of the interfacial electron transfer of wild-type cytochrome c and herein produced variants, adsorbed electrostatically under different local interfacial electric fields, were determined by means of variable temperature cyclic film voltammetry. It is shown that non-native cytochrome c variants immobilized under a low interfacial electric field display redox thermodynamics and kinetics similar to those of wild-type cytochrome c. However, upon increasing the strength of the electric field, the redox thermodynamics and kinetics of the modified proteins markedly differ from those of the wild-type species. The mutations promote stabilization of the oxidized form and a significant increase in the activation enthalpy values that can be ascribed to a subtle distortion of the heme cofactor and/or difference of the amino acid rearrangements rather than to a coarse protein structural change. Overall, these results point to a combined effect of the single point mutations at positions 48 and 97 and the strength of electrostatic binding on the regulatory mechanism of mitochondrial membrane activity, when acting as a redox shuttle protein.

Converting cytochrome c into a DyP-like metalloenzyme

Dye-decolorizing peroxidase (DyP), which can degrade anthraquinone dyes using H₂O₂, is an attractive prospect for potential biotechnological applications for environmental purification. We previously designed an artificial DyP with an optimal pH for reactive blue 19 (RB19) degradation shifting from pH 4.5 to 6.5. We then attempted to degrade RB19 using Escherichia coli expressing this mutant, but RB19 was degraded equally compared with bacteria expressing wild-type (WT) DyP because most DyP was expressed in a heme-free form. In this study, we attempted to design an artificial peroxidase based on cytochrome c (cyt c), whose heme is covalently bound to the protein. We found that cyt c can degrade RB19, but its ability at pH 7.0 was ∼60% of that of DyP from Vibrio cholerae at pH 4.5. To enhance this activity we constructed several mutants using three approaches. Initially, to improve reactivity with H₂O₂, Met80 was replaced with a noncoordinating residue, Ala or Val, but catalytic efficiency (k_cat/K_m) was increased by only ∼1.5-fold. To enhance the substrate binding affinity we introduced an additional Trp by replacing Pro76 (P76W). The catalytic efficiency of this mutant was ∼3-fold greater than that of WT cyt c. Finally, to form a hydrogen bond to axial histidine Gly29 was replaced with Asp (G29D). This mutant exhibited an ∼80-fold greater dye-decolorizing activity. Escherichia coli expressing the G29D mutant was unable to degrade RB19 in solution due to degradation of heme itself, but this study provides new insights into the design of artificial DyPs.

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

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

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.

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.

How to Turn an Electron Transfer Protein into a Redox Enzyme for Biosensing

Cytochrome c is a small globular protein whose main physiological role is to shuttle electrons within the mitochondrial electron transport chain. This protein has been widely investigated, especially as a paradigmatic system for understanding the fundamental aspects of biological electron transfer and protein folding. Nevertheless, cytochrome c can also be endowed with a non-native catalytic activity and be immobilized on an electrode surface for the development of third generation biosensors. Here, an overview is offered of the most significant examples of such a functional transformation, carried out by either point mutation(s) or controlled unfolding. The latter can be induced chemically or upon protein immobilization on hydrophobic self-assembled monolayers. We critically discuss the potential held by these systems as core constituents of amperometric biosensors, along with the issues that need to be addressed to optimize their applicability and response.

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.

In Situ Spectroelectrochemical Investigations of Electrode-Confined Electron-Transferring Proteins and Redox Enzymes

This perspective analyzes recent advances in the spectroelectrochemical investigation of redox proteins and enzymes immobilized on biocompatible or biomimetic electrode surfaces. Specifically, the article highlights new insights obtained by surface-enhanced resonance Raman (SERR), surface-enhanced infrared absorption (SEIRA), protein film infrared electrochemistry (PFIRE), polarization modulation infrared reflection-absorption spectroscopy (PMIRRAS), Förster resonance energy transfer (FRET), X-ray absorption spectroscopy (XAS), electron paramagnetic resonance (EPR), and differential electrochemical mass spectrometry (DMES)-based spectroelectrochemical methods on the structure, orientation, dynamics, and reaction mechanisms for a variety of immobilized species. This includes small heme and copper electron shuttling proteins, large respiratory complexes, hydrogenases, multicopper oxidases, alcohol dehydrogenases, endonucleases, NO-reductases, and dye decolorizing peroxidases, among other enzymes. Finally, I discuss the challenges and foreseeable future developments toward a better understanding of the functioning of these complex macromolecules and their exploitation in technological devices.

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.

Exploring protein phosphorylation by combining computational approaches and biochemical methods

Post-translational modifications of proteins expand their functional diversity, regulating the response of cells to a variety of stimuli. Among these modifications, phosphorylation is the most ubiquitous and plays a prominent role in cell signaling. The addition of a phosphate often affects the function of a protein by altering its structure and dynamics. However, these alterations are often difficult to study and the functional and structural implications remain unresolved. New approaches are emerging to overcome common obstacles related to the production and manipulation of these samples. Here, we summarize the available methods for phosphoprotein purification and phosphomimetic engineering, highlighting the advantages and disadvantages of each. We propose a general workflow for protein phosphorylation analysis combining computational and biochemical approaches, building on recent advances that enable user-friendly and easy-to-access Molecular Dynamics simulations. We hope this innovative workflow will inform the best experimental approach to explore such post-translational modifications. We have applied this workflow to two different human protein models: the hemeprotein cytochrome c and the RNA binding protein HuR. Our results illustrate the usefulness of Molecular Dynamics as a decision-making tool to design the most appropriate phosphomimetic variant.

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.

Wheel and Deal in the Mitochondrial Inner Membranes: The Tale of Cytochrome c and Cardiolipin

Cardiolipin oxidation and degradation by different factors under severe cell stress serve as a trigger for genetically encoded cell death programs. In this context, the interplay between cardiolipin and another mitochondrial factor—cytochrome c—is a key process in the early stages of apoptosis, and it is a matter of intense research. Cytochrome c interacts with lipid membranes by electrostatic interactions, hydrogen bonds, and hydrophobic effects. Experimental conditions (including pH, lipid composition, and post-translational modifications) determine which specific amino acid residues are involved in the interaction and influence the heme iron coordination state. In fact, up to four binding sites (A, C, N, and L), driven by different interactions, have been reported. Nevertheless, key aspects of the mechanism for cardiolipin oxidation by the hemeprotein are well established. First, cytochrome c acts as a pseudoperoxidase, a process orchestrated by tyrosine residues which are crucial for peroxygenase activity and sensitivity towards oxidation caused by protein self-degradation. Second, flexibility of two weakest folding units of the hemeprotein correlates with its peroxidase activity and the stability of the iron coordination sphere. Third, the diversity of the mode of interaction parallels a broad diversity in the specific reaction pathway. Thus, current knowledge has already enabled the design of novel drugs designed to successfully inhibit cardiolipin oxidation.

Cytochrome c phosphorylation: Control of mitochondrial electron transport chain flux and apoptosis

Cytochrome c (Cytc)¹is a cellular life and death decision molecule that regulates cellular energy supply and apoptosis through tissue specific post-translational modifications. Cytc is an electron carrier in the mitochondrial electron transport chain (ETC) and thus central for aerobic energy production. Under conditions of cellular stress, Cytc release from the mitochondria is a committing step for apoptosis, leading to apoptosome formation, caspase activation, and cell death. Recently, Cytc was shown to be a target of cellular signaling pathways that regulate the functions of Cytc by tissue-specific phosphorylations. So far five phosphorylation sites of Cytc have been mapped and functionally characterized, Tyr97, Tyr48, Thr28, Ser47, and Thr58. All five phosphorylations partially inhibit respiration, which we propose results in optimal intermediate mitochondrial membrane potentials and low ROS production under normal conditions. Four of the phosphorylations result in inhibition of the apoptotic functions of Cytc, suggesting a cytoprotective role for phosphorylated Cytc. Interestingly, these phosphorylations are lost during stress conditions such as ischemia. This results in maximal ETC flux during reperfusion, mitochondrial membrane potential hyperpolarization, excessive ROS generation, and apoptosis. We here present a new model proposing that the electron transfer from Cytc to cytochrome c oxidase is the rate-limiting step of the ETC, which is regulated via post-translational modifications of Cytc. This regulation may be dysfunctional in disease conditions such as ischemia-reperfusion injury and neurodegenerative disorders through increased ROS, or cancer, where post-translational modifications on Cytc may provide a mechanism to evade apoptosis.

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.

Tissue-specific regulation of cytochrome c by post-translational modifications: respiration, the mitochondrial membrane potential, ROS, and apoptosis

Cytochrome c (Cyt c) plays a vital role in the mitochondrial electron transport chain (ETC). In addition, it is a key regulator of apoptosis. Cyt c has multiple other functions including ROS production and scavenging, cardiolipin peroxidation, and mitochondrial protein import. Cyt c is tightly regulated by allosteric mechanisms, tissue-specific isoforms, and post-translational modifications (PTMs). Distinct residues of Cyt c are modified by PTMs, primarily phosphorylations, in a highly tissue-specific manner. These modifications downregulate mitochondrial ETC flux and adjust the mitochondrial membrane potential (ΔΨm), to minimize reactive oxygen species (ROS) production under normal conditions. In pathologic and acute stress conditions, such as ischemia-reperfusion, phosphorylations are lost, leading to maximum ETC flux, ΔΨm hyperpolarization, excessive ROS generation, and the release of Cyt c. It is also the dephosphorylated form of the protein that leads to maximum caspase activation. We discuss the complex regulation of Cyt c and propose that it is a central regulatory step of the mammalian ETC that can be rate limiting in normal conditions. This regulation is important because it maintains optimal intermediate ΔΨm, limiting ROS generation. We examine the role of Cyt c PTMs, including phosphorylation, acetylation, methylation, nitration, nitrosylation, and sulfoxidation and consider their potential biological significance by evaluating their stoichiometry.-Kalpage, H. A., Bazylianska, V., Recanati, M. A., Fite, A., Liu, J., Wan, J., Mantena, N., Malek, M. H., Podgorski, I., Heath, E. I., Vaishnav, A., Edwards, B. F., Grossman, L. I., Sanderson, T. H., Lee, I., Hüttemann, M. Tissue-specific regulation of cytochrome c by post-translational modifications: respiration, the mitochondrial membrane potential, ROS, and apoptosis.

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.

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.