Oxidation state-dependent conformational changes in cytochrome c

Thermodynamics of the equilibrium unfolding of oxidized and reduced Saccharomyces cerevisiae iso-1-cytochromes c

We report thermodynamic data for the chemical denaturation of iso-1-cytochromes c from Saccharomyces cerevisiae having amino acid substitutions R38A, N52I, and F82S in all possible combinations. The guanidine hydrochloride denaturation of isolated proteins was monitored by fluorescence measurements. The redox potentials, Eo', for both the folded and unfolded conformations have been measured. Free energy changes of chemical unfolding together with direct electrochemical measurement of the free energy changes of reduction for both the native and unfolded proteins yield a complete thermodynamic cycle, which includes four states of cytochrome c: oxidized folded, oxidized unfolded, reduced folded, and reduced unfolded. Completed cycles illustrate that the stability of cytochrome c to denaturing conditions is different for each amino acid substitution by an amount that depends on the heme oxidation state. Thus, the differential protein stability cannot be interpreted simply in terms of a hydrophobic effect, without also considering coupled Coulombic effects.

Effects of mutating Asn-52 to isoleucine on the haem-linked properties of cytochrome c

Asn-52 of rat cytochrome c and baker's yeast iso-1-cytochrome c was changed to isoleucine by site-directed mutagenesis and the mutated proteins expressed in and purified from cultures of transformed yeast. This mutation affected the affinity of the haem iron for the Met-80 sulphur in the ferric state and the reduction potential of the molecule. The yeast protein, in which the sulphur-iron bond is distinctly weaker than in vertebrate cytochromes c, became very similar to the latter: the pKa of the alkaline ionization rose from 8.3 to 9.4 and that of the acidic ionization decreased from 3.4 to 2.8. The rates of binding and dissociation of cyanide became markedly lower, and the affinity was lowered by half an order of magnitude. In the ferrous state the dissociation of cyanide from the variant yeast cytochrome c was three times slower than in the wild-type. The same mutation had analogous but less pronounced effects on rat cytochrome c: it did not alter the alkaline ionization pKa nor its affinity for cyanide, but it lowered its acidic ionization pKa from 2.8 to 2.2. These results indicate that the mutation of Asn-52 to isoleucine increases the stability of the cytochrome c closed-haem crevice as observed earlier for the mutation of Tyr-67 to phenylalanine [Luntz, Schejter, Garber and Margoliash (1989) Proc. Natl. Acad. Sci. U.S.A. 86, 3524-3528], because of either its effects on the hydrogen-bonding of an interior water molecule or a general increase in the hydrophobicity of the protein in the domain occupied by the mutated residues. The reduction potentials were affected in different ways; the Eo of rat cytochrome c rose by 14 mV whereas that of the yeast iso-1 cychrome c was 30 mV lower as a result of the change of Asn-52 to isoleucine.

Stability of yeast iso-1-ferricytochrome c as a function of pH and temperature

Absorbance-detected thermal denaturation studies of the C102T variant of Saccharomyces cerevisiae iso-1-ferricytochrome c were performed between pH 3 and 5. Thermal denaturation in this pH range is reversible, shows no concentration dependence, and is consistent with a 2-state model. Values for free energy (delta GD), enthalpy (delta HD), and entropy (delta SD) of denaturation were determined as functions of pH and temperature. The value of delta GD at 300 K, pH 4.6, is 5.1 +/- 0.3 kcal mol-1. The change in molar heat capacity upon denaturation (delta Cp), determined by the temperature dependence of delta HD as a function of pH (1.37 +/- 0.06 kcal mol-1 K-1), agrees with the value determined by differential scanning calorimetry. pH-dependent changes in the Soret region indicate that a group or groups in the heme environment of the denatured protein, probably 1 or both heme propionates, ionize with a pK near 4. The C102T variant exhibits both enthalpy and entropy convergence with a delta HD of 1.30 kcal mol-1 residue-1 at 373.6 K and a delta SD of 4.24 cal mol-1 K-1 residue-1 at 385.2 K. These values agree with those for other single-domain, globular proteins.

Mutations of iso-1-cytochrome c at positions 13 and 90. Separate effects on physical and functional properties

Residues at positions 13 (lysine or arginine) and 90 (glutamate or aspartate) of eukaryotic cytochromes c have been conserved during evolution; Cys102, however, is found only in yeast cytochrome c. The positively charged residue at position 13 and the negatively charged residue at position 90 are close together in those cytochromes c for which three-dimensional structures are available. We have replaced the amino acids at these two positions by cysteine in Saccharomyces cerevisiae iso-1-cytochrome c; in an earlier study, Cys102 was replaced by threonine without negatively influencing the physical or enzymic properties of the protein. The mutated proteins [R13C, C102T]cytochrome c (iso-1-cytochrome c containing Arg13-->Cys and Cys102-->Thr mutations), [D90C, C102T]cytochrome c (iso-1-cytochrome c containing Asp90-->Cys and Cys102-->Thr mutations) and [R13C, D90C, C102T]cytochrome c (iso-1-cytochrome c containing Arg13-->Cys, Asp90-->Cys, and Cys102-->Thr mutations) are functional in vivo. Free sulfhydryl titration shows that the doubly mutated forms each contain one sulfhydryl group while the triple mutant contains two sulfhydryl groups. The stability of mutant [R13C, C102T]cytochrome c resembles that of [C102T] cytochrome c, whereas the stability of [D90C, C102T]cytochrome c resembles the stability of [R13C, D90C, C102T]cytochrome c. The activity of cytochrome-c oxidase using cytochrome c was monitored polarographically. Compared to the wild-type or [C102T]cytochrome c, which shows two kinetic phases with cytochrome-c oxidase, [D90C, C102T]cytochrome c has much the same profile; [R13C, C102T]cytochrome c and [R13C, D90C, C102T]cytochrome c exhibit one kinetic phase with decreased activity. Electron-transfer activity of the mutant cytochromes c is inhibited by Hg2+. The inhibition is highest for the triple mutant, less for [R13C, C102T]cytochrome c, even less for [D90C, C102T]cytochrome c and insignificant for the wild type. It would appear as though the stability of the triple mutant follows the changes that result from the Asp90-->Cys mutation while the activity changes follow those of the Arg13-->Cys mutation.

Structural water in oxidized and reduced horse heart cytochrome c

The existence of structural water in the interior of both oxidized and reduced horse-heart cytochrome c in solution is demonstrated using nuclear magnetic resonance spectroscopy. Six water molecules have been located in ferrocytochrome c and five in ferricytochrome c, with residence times greater than a few hundred picoseconds. Two water molecules are located in the haem crevice, one of which is found to undergo a large change in position with a change of oxidation state. Both of these observations indicate that buried structural waters in the haem crevice have, by microscopic dielectric effects, significant roles in the setting of the solvent reorganization energy associated with electron transfer.

The significance of denaturant titrations of protein stability: a comparison of rat and baker's yeast cytochrome c and their site-directed asparagine-52-to-isoleucine mutants

The residue asparagine-52 of rat cytochrome c and baker's yeast iso-1-cytochrome c was mutated to isoleucine by site-directed mutagenesis, and the unfolding of the wild-type and mutant proteins in urea or guanidinium chloride solutions was studied. Whereas the yeast mutant cytochrome unfolded in 4-7 M urea with a rate constant (k) of 1.7 x 10(-2) s-1, the rat mutant protein unfolded with k = 5.0 x 10(-2) s-1, followed by a slow partial refolding with k = 5.0 x 10(-4) s-1. Denaturant titrations indicated that the mutation increased the stability of the yeast cytochrome by 6.3 kJ (1.5 kcal)/mol, while it decreased that of the rat protein by 11.7 kJ (2.8 kcal)/mol. These results probably reflect structural differences between yeast iso-1 and vertebrate cytochromes c in the vicinity of the Asn-52 side chain.

Laue diffraction study on the structure of cytochrome c peroxidase compound I

BACKGROUND

Cytochrome c peroxidase from yeast is a soluble haem-containing protein found in the mitochondrial electron transport chain where it probably protects against toxic peroxides. The aim of this study was to obtain a reliable structure for the doubly oxidized transient intermediate (termed compound I) in the reaction of cytochrome c peroxidase with hydrogen peroxide. This intermediate contains a semistable free radical on Trp191, and an oxyferryl haem group.

RESULTS

Compound I was produced in crystals of yeast cytochrome c peroxidase by reacting the crystalline enzyme with hydrogen peroxide in a flow cell. The reaction was monitored by microspectrophotometry and Laue crystallography in separate experiments. A nearly complete conversion to compound I was achieved within two minutes of the addition of hydrogen peroxide, and the concentration of the intermediate remained at similar levels for an additional half an hour. The structure of the intermediate was determined by Laue diffraction. The refined Laue structure for compound I shows clear structural changes at the peroxide-binding site but no significant changes at the radical site. The photographs were processed with a new software package (LEAP), overcoming many of the former problems encountered in extracting structural information from Laue exposures.

CONCLUSIONS

The geometry of the haem environment in this protein allows structural changes to be extremely small, similar in magnitude to those observed for the Fe2+/Fe3+ transition in cytochrome c. The results suggest that these molecules have evolved to transfer electrons with a minimal need for structural adjustment.

Analysis of the bimolecular reduction of ferricytochrome c by ferrocytochrome b5 through mutagenesis and molecular modelling

Site-directed mutagenesis has been used to produce variants of cytochrome c in which selected structural or functional properties of this protein are altered that have been implicated previously in contributing to the rate at which ferricytochrome c is reduced by ferrocytochrome b5. In total, 18 variants have been studied by kinetics and electrochemical methods to assess the contributions of thermodynamic driving force, surface charge and hydrophobic interactions, and redox-linked structural reorganization of the protein to the rate of electron transfer between these two proteins under conditions where the reaction is bimolecular. While some variants (those at position-38) appear to affect primarily the driving force of the reaction, others appear to influence the rearrangement barrier to electron transfer (those at positions-67 and -52) while the interface between electron donor and acceptor centers is the principal effect of substitutions for a conserved aromatic heme contact residue at the surface of the protein (position-82). Interpretation of these results has been facilitated through the use of energy minimization calculations to refine the hypothetical models previously suggested for the cytochrome c- cytochrome b5 precursor complex on the basis of Brownian dynamics simulations of the bimolecular encounter event.

Interaction between cytochrome c and the photosynthetic reaction center of purple bacteria: behaviour at low temperature

In purple photosynthetic bacteria the electron donor to the special pair, after its oxidation by a light-induced reaction, is a c-type cytochrome: either a soluble monoheme cytochrome which forms a transitory complex with the reaction center, or a tetraheme cytochrome which remains permanently bound to the reaction center. The effects of low temperatures on electron transfer in the complex are presented and discussed. They provide estimates for the reorganization energy. The most prominent effect of low temperature is that a dominant fast phase of electron transfer becomes impossible at a temperature of around 250 K (monoheme cytochrome) or located between 250 K and 80 K according to the redox state (tetraheme cytochrome). This inhibition is attributed to a freezing-like transition of pools of water molecules which blocks structural changes of the protein which are normally associated with the cytochrome oxidation.

Involvement of electrostatic interactions in cytochrome c complex formations

Structural studies on various electron transfer complexes involving the tetrahemic cytochrome c3 provided evidence that one of the hemes (heme 4) is the interacting site on the molecule. The reactivity of this particular heme is allocated to the positive charges found around the heme group which are strongly involved in the electrostatic interaction processes. Electrostatic and hydrophobic effects in complex formation are considered on the basis of two electron transfer complex examples: the soluble cytochrome c-cytochrome c peroxidase and the membrane bound photosynthetic reaction center.

Probing weakly polar interactions in cytochrome c

Theoretical, statistical, and model studies suggest that proteins are stabilized by weakly polar attractions between sulfur atoms and properly oriented aromatic rings. The two sulfur-containing amino acids, methionine and cysteine, occur frequently among functional alleles in random mutant libraries of Saccharomyces cerevisiae iso-1-cytochrome c genes at positions that form a weakly polar aromatic-aromatic interaction, the wild-type protein. To determine if a weakly polar sulfur-aromatic interaction replaced the aromatic-aromatic interaction, the structure and stability of two variants were examined. Phenylalanine 10, which interacts with tyrosine 97, was replaced by methionine and cysteine. The cysteine was modified to form the methionine and cysteine analog, S-methyl cysteine (CysSMe). Proton NMR studies indicate that changing Phe 10 to Met or CysSMe affects only local structure and that the structures of sulfur-containing variants are nearly identical. Analysis of chemical shifts and nuclear Overhauser effect data indicates that both sulfur-containing side chains are in position to form a weakly polar interaction with Tyr 97. The F10M and F10CSMe variants are 2-3 kcal mol-1 less stable than iso-1-cytochrome c at 300 K. Comparison of the stabilities of the F10M and F10CSMe variants allows evaluation of the potential weakly polar interaction between the additional sulfur atom of F10CSMe and the aromatic moiety of Tyr 97. The F10CSMe;C102T variant is 0.7 +/- 0.3 kcal mol-1 more stable than the F10M;C102T protein. The increased stability is explained by the difference in hydrophobicity of the sulfur-containing side chains. We conclude that any weakly polar interaction between the additional sulfur and the aromatic ring is too weak to detect or is masked by destabilizing contributions to the free energy of denaturation.

Amide proton exchange rates of oxidized and reduced Saccharomyces cerevisiae iso-1-cytochrome c

Proton NMR spectroscopy was used to determine the rate constant, kobs, for exchange of labile protons in both oxidized (Fe(III)) and reduced (Fe(II)) iso-1-cytochrome c. We find that slowly exchanging backbone amide protons tend to lack solvent-accessible surface area, possess backbone hydrogen bonds, and are present in regions of regular secondary structure as well as in omega-loops. Furthermore, there is no correlation between kobs and the distance from a backbone amide nitrogen to the nearest solvent-accessible atom. These observations are consistent with the local unfolding model. Comparisons of the free energy change for denaturation, delta Gd, at 298 K to the free energy change for local unfolding, delta Gop, at 298 K for the oxidized protein suggest that certain conformations possessing higher free energy than the denatured state are detected at equilibrium. Reduction of the protein results in a general increase in delta Gop. Comparisons of delta Gd to delta Gop for the reduced protein show that the most open states of the reduced protein possess more structure than its chemically denatured form. This persistent structure in high-energy conformations of the reduced form appears to involve the axially coordinated heme.

Changes in global stability and local structure of cytochrome c upon substituting phenylalanine-82 with tyrosine

We have examined the F82Y;C102T variant of Saccharomyces cerevisiae iso-1-cytochrome c using high-resolution proton nuclear magnetic resonance spectroscopy, chemical denaturation, and differential scanning calorimetry. Comparison of proton chemical shifts, paramagnetic shifts, and nuclear Overhauser effects indicates structural changes are localized to the vicinity of position 82. One alteration involves the rearrangement of the side chain of leucine-85. Using many more proton assignments than were available in the initial report [G. J. Pielak, R. A. Atkinson, J. Boyd, and R. J. P. Williams, Eur. J. Biochem. 177, 179-185 (1988)], a second alteration involving an interaction between arginine-13 and tyrosine-82 is observed. The interaction appears to involve a hydrogen bond with the eta-protons of arginine's guanido group acting as donor and tyrosine's phenolic eta-oxygen as acceptor. In spite of this potentially-stabilizing interaction, the free energy of denaturation decreases by approximately 2.4 kcal mol-1. Results are discussed with respect to alterations in the native and denatured states.

Importance of protein rearrangement in the electron-transfer reaction between the physiological partners cytochrome f and plastocyanin

Cytochrome f from turnip and plastocyanin from French bean were noninvasively cross-linked in the presence of the carbodiimide EDC so that the exposed heme edge in the former protein abuts the acidic patch remote from the copper site in the latter [Morand, L.Z., Frame, M.K., Colvert, K.K., Johnson, D.A., Krogmann, D.W., & Davis, D.J. (1989) Biochemistry 28, 8039]. The molecular mass, reduction potentials, and UV-visible and ESR spectra of the covalent complex were consistent with the composition cyt/pc and with a lack of noticeable structural perturbations of the protein molecules. Isoelectric focusing showed the presence of N-acylurea groups, byproducts of the cross-linking reaction [Zhou, J.S., Brothers, H.M. II, Neddersen, J.P., Peerey, L.M., Cotton, T.M., & Kostić, N.M. (1992) Bioconjugate Chem. 3, 382]. Laser flash spectroscopy, with riboflavin semiquinone as the reductant, showed that the electrontransfer reaction within the covalent complex cyt(II)/pc(II) is either undetectably slow or reversible. The question was resolved by monitoring, during redox titrations, the 1H NMR line widths of the heme methyl groups in free ferricytochrome f and in this protein cross-linked to plastocyanin.(ABSTRACT TRUNCATED AT 250 WORDS)

Evaluation of 13C and 1H Fermi contact shifts in horse cytochrome c. The origin of the anti-Curie effect

Many ferricytochromes c exhibit a peculiar effect in which the 1H chemical shifts of the haem methyl groups appear in pairs and, although the paramagnetic shifts of the two groups with the larger shifts decrease with temperature, those of the pair with the smaller shifts actually increase. Recent NMR studies [Santos, H. and Turner, D. L. (1992) Eur. J. Biochem. 206, 721-728] gave 1H and 13C assignments for most of the haem substituents and the axial ligands in horse cytochrome c at 30 degrees C and 50 degrees C in both oxidation states. These data are used together with an empirically determined magnetic susceptibility tensor to evaluate the Fermi contact contribution to the paramagnetic shift and hence map the delocalization of the unpaired electron. The anti-Curie effect is explained by a Boltzmann distribution between partially filled porphyrin 3e(pi) molecular orbitals with an energy difference of 3 kJ/mol. The fact that the energy gap is small with respect to the energy of binding to the electron transfer partners calls into question the significance of the asymmetry of the electron distribution in the electron transfer process.

1H- and 13C-NMR investigation of redox-state-dependent and temperature-dependent conformation changes in horse cytochrome c

The redox-state dependent changes in chemical shift, which have been measured for almost 100 CHn groups in the 13C-NMR spectra of horse cytochrome c [Santos, H., and Turner, D. L. (1992) Eur. J. Biochem. 206, 721-728], have been used to investigate the nature of the redox-related change in conformation. Apart from the haem and its axial ligands, the shifts are found to be dominated by the electron-nuclear dipolar coupling in the oxidised form, as was the case in 1H-NMR studies. These pseudocontact shifts are well described by using an empirically determined magnetic susceptibility tensor in conjunction with atomic coordinates for the horse cytochrome c. The groups which fit least well are located in the vicinity of Trp59. Comparison between 1H and 13C shifts and their temperature dependence shows that the differences from expectation based on a single structure for both oxidation states are caused largely by changes in the diamagnetic contribution to the chemical shifts. Since these are different for 1H and 13C resonances they indicate, independently from crystal structure data, some redox-related movement of the protein under the haem. The significance of these results for understanding electron transfer pathways is discussed. Finally, the temperature dependence of the pseudocontact shifts in the range 30-50 degrees C is shown to be anomalous. Approximately half of the anomalous effect may be attributed to Zeeman mixing of the electronic wavefunctions with a spin-orbit coupling constant lambda = 241 cm-1, while the other half is attributed to thermal expansion of the protein.

Exploring the interface between the N- and C-terminal helices of cytochrome c by random mutagenesis within the C-terminal helix

Buried within cytochrome c lies a highly-conserved helix-helix interface formed by the perpendicular packing of the C-terminal helix against the N-terminal helix. This interface involves a peg-in-hole interaction between Gly-6 and Leu-94 and an aromatic-aromatic interaction between Phe-10 and Tyr-97. To gain insight into protein design, we investigated the relationship between the sequence of the interface and the physiological function of yeast iso-1-cytochrome c. A library of mutants at positions 94 and 97 of the C-terminal helix was created to examine the effect of novel amino acid combinations. We isolated 45 of the 400 possible amino acid combinations, 32 of which result in a functional cytochrome c. Contrary to evolutionary conservation of the peg-in-hole and aromatic-aromatic interactions, we find that side-chain volume and conservation of aromatic residues do not play an essential role in determining function. Additionally, we find negatively-charged residues within the interface that result in a functional cytochrome c. Examination of the 45 missense mutants indicates that approximately 120 unique combinations are compatible with function. These results show that the interface is flexible. However, truncation of the C-terminal helix at position 94 abolishes function, suggesting that the interface is essential. The correlation observed between our library of mutants and the mutation matrix compiled by Gonnet et al. [Gonnet, G. H., Cohen, M. A., & Benner, S. A. (1992) Science 256, 1443-1445] demonstrates the potential use of the matrix to predict the effect of sequence changes on natural proteins and to optimize the design of novel proteins.

Expression of recombinant cytochromes c from various species in Saccharomyces cerevisiae: post-translational modifications

A complete protocol for the expression of recombinant cytochrome c genes from yeast, Drosophila melanogaster, and rat in a yeast strain, GM-3C-2, which does not express its own cytochromes c is described. The construction of the expression vectors, transformation and large-scale growth of the yeast, and preparation and purification of the recombinant cytochromes c are described. It was found that, contrary to the way yeast modifies its own cytochromes c, the recombinant proteins were partially acetylated at their N-terminus, except for the drosophila protein, which remained entirely unblocked. Furthermore, the yeast and rat proteins were close to fully trimethylated at lysine 72, while the drosophila protein could be separated chromatographically into forms containing tri-, di-, mono-, and unmethylated lysine 72 showing corresponding resonances in the NMR spectrum. These observations emphasize that, in employing expression procedures to obtain native or mutant forms of cytochrome c, it is essential to identify the variety and extent of post-translational modifications and to separate the preparation into pure monomolecular species. Otherwise, it may become impossible to distinguish between the influence of a site-directed mutation and unexamined post-translational modifications.

Structure determination and analysis of yeast iso-2-cytochrome c and a composite mutant protein

As part of a study of protein folding and stability, the three-dimensional structures of yeast iso-2-cytochrome c and a composite protein (B-2036) composed of primary sequences of both iso-1 and iso-2-cytochromes c have been solved to 1.9 A and 1.95 A resolutions, respectively, using X-ray diffraction techniques. The sequences of iso-1 and iso-2-cytochrome c share approximately 84% identity and the B-2036 composite protein has residues 15 to 63 from iso-2-cytochrome c with the rest being derived form the iso-1 protein. Comparison of these structures reveals that amino acid substitutions result in alterations in the details of intramolecular interactions. Specifically, the substitution Leu98Met results in the filling of an internal cavity present in iso-1-cytochrome c. Further substitutions of Val20Ile and Cys102Ala alter the packing of secondary structure elements in the iso-2 protein. Blending the isozymic amino acid sequences in this latter area results in the expansion of the volume of an internal cavity in the B-2036 structure to relieve a steric clash between Ile20 and Cys102. Modification of hydrogen bonding and protein packing without disrupting the protein fold is illustrated by the His26Asn and Asn63Ser substitutions between iso-1 and iso-2-cytochromes c. Alternatively, a change in main-chain fold is observed at Gly37 apparently due to a remote amino acid substitution. Further structural changes occur at Phe82 and the amino terminus where a four residue extension is present in yeast iso-2-cytochrome c. An additional comparison with all other eukaryotic cytochrome c structures determined to date is presented, along with an analysis of conserved water molecules. Also determined are the midpoint reduction potentials of iso-2 and B-2036 cytochromes c using direct electrochemistry. The values obtained are 286 and 288 mV, respectively, indicating that the amino acid substitutions present have had only a small impact on the heme reduction potential in comparison to iso-1-cytochrome c, which has a reduction potential of 290 mV.