Comparison of reduced and oxidized yeast iso-1-cytochrome c using proton paramagnetic shifts

Redox-dependent conformational changes in eukaryotic cytochromes revealed by paramagnetic NMR spectroscopy

Cytochrome c (Cc) is a soluble electron carrier protein, transferring reducing equivalents between Cc reductase and Cc oxidase in eukaryotes. In this work, we assessed the structural differences between reduced and oxidized Cc in solution by paramagnetic NMR spectroscopy. First, we have obtained nearly-complete backbone NMR resonance assignments for iso-1-yeast Cc and horse Cc in both oxidation states. These were further used to derive pseudocontact shifts (PCSs) arising from the paramagnetic haem group. Then, an extensive dataset comprising over 450 measured PCSs and high-resolution X-ray and solution NMR structures of both proteins were used to define the anisotropic magnetic susceptibility tensor, Δχ. For most nuclei, the PCSs back-calculated from the Δχ tensor are in excellent agreement with the experimental PCS values. However, several contiguous stretches-clustered around G41, N52, and A81-exhibit large deviations both in yeast and horse Cc. This behaviour is indicative of redox-dependent structural changes, the extent of which is likely conserved in the protein family. We propose that the observed discrepancies arise from the changes in protein dynamics and discuss possible functional implications.

Vibrational dynamics of iron in cytochrome C

Nuclear resonance vibrational spectroscopy (NRVS) and Raman spectroscopy on (54)Fe- and (57)Fe-enriched cytochrome c (cyt c) identify multiple bands involving vibrations of the heme Fe. Comparison with predictions from Fe isotope shifts reveals that 70% of the NRVS signal in the 300-450 cm(-1) frequency range corresponds to vibrations resolved in Soret-enhanced Raman spectra. This frequency range dominates the "stiffness", an effective force constant determined by the Fe vibrational density of states (VDOS), which measures the strength of nearest-neighbor interactions with Fe. The stiffness of the low-spin Fe environment in both oxidation states of cyt c significantly exceeds that for the high-spin Fe in deoxymyoglobin, where the 200-300 cm(-1) frequency range dominates the VDOS. This situation is reflected in the shorter Fe-ligand bond lengths in the former with respect to the latter. The longer Fe-S(Met80) in oxidized cyt c with respect to reduced cyt c leads to a decrease in the stiffness of the iron environment upon oxidation. Comparison with NRVS measurements allows us to assess assignments for vibrational modes resolved in this region of the heme Raman spectrum. We consider the possibility that the 372 cm(-1) band in reduced cyt c involves the Fe-S(Met80) bond.

Resilience of the iron environment in heme proteins

Conformational flexibility is essential to the functional behavior of proteins. We use an effective force constant introduced by Zaccai, the resilience, to quantify this flexibility. Site-selective experimental and computational methods allow us to determine the resilience of heme protein active sites. The vibrational density of states of the heme Fe determined using nuclear resonance vibrational spectroscopy provides a direct experimental measure of the resilience of the Fe environment, which we compare quantitatively with values derived from the temperature dependence of atomic mean-squared displacements in molecular dynamics simulations. Vibrational normal modes in the THz frequency range dominate the resilience. Both experimental and computational methods find a higher resilience for cytochrome c than for myoglobin, which we attribute to the increased number of covalent links to the peptide in the former protein. For myoglobin, the resilience of the iron environment is larger than the average resilience previously determined for hydrogen sites using neutron scattering. Experimental results suggest a slightly reduced resilience for cytochrome c upon oxidation, although the change is smaller than reported in previous Mössbauer investigations on a bacterial cytochrome c, and is not reproduced by the simulations. Oxidation state also has no significant influence on the compressibility calculated for cyt c, although a slightly larger compressibility is predicted for myoglobin.

Effects of axial methionine coordination on the in-plane asymmetry of the heme electronic structure of cytochrome c

The paramagnetic susceptibility ( chi) tensors of the oxidized forms of thermophile Hydrogenobacter thermophilus cytochrome c(552) (Ht cyt c(552)) and a quintuple mutant (F7A/V13 M/F34Y/E43Y/V78I; qm) of mesophile Pseudomonas aeruginosa cytochrome c(551) (Pa cyt c(551)) have been determined on the basis of the redox-dependent (1)H NMR shift changes of the main-chain NH and C(alpha)H proton resonances of non-coordinated amino acid residues and the NMR structures of the reduced forms of the corresponding proteins (J. Hasegawa, T. Yoshida, T. Yamazaki, Y. Sambongi, Y. Yu, Y. Igarashi, T. Kodama, K. Yamazaki, Y. Kyogoku, Y. Kobayashi (1998) Biochemistry 37:9641-9649; J. Hasegawa, S. Uchiyama, Y. Tanimoto, M. Mizutani, Y. Kobayashi, Y. Sambongi,Y. Igarashi (2000) J Biol Chem 275:37824-37828). From the chi tensors determined, we obtained the contact shifts for heme methyl proton resonances, which provided the heme electronic structures of the oxidized forms of Ht cyt c(552) and qm. We also characterized the heme electronic structure of the cyanide adducts of the proteins, where the axial Met was replaced by an exogenous cyanide ion, through the analysis of (1)H NMR spectra. The results indicated that the heme electronic structures of both the proteins in their oxidized forms with axial His and Met coordination are largely different to each other, while those in their cyanide adducts are similar to each other. These results demonstrated that the orientation of the axial Met sulfur lone pair, with respect to heme, predominantly contributes to the spin delocalization into the porphyrin-pi system of heme in the oxidized proteins with axial His and Met coordination.

The problem of proton transfer in membranes

Proton (H(+)) transfer has been examined in many molecular systems for more than 50 years. General mechanistic possibilities including tunnelling have been recognized for proton movement from local base-to-base centres. An especially fast mechanism over considerable distances, the Grotthus mechanism, has been described in water. Proton transfer over long distances in membranes which is now known to occur in many protein bio-energetic devices is not understood since ground state structures do not provide a continuous H-bond network. Here we consider the possible mechanisms and propose that the most likely pathway for protons in membranes uses an excited conformational state, equivalent to a partial denaturation.

Redox-related conformational changes in Rhodobacter capsulatus cytochrome c2

WEFT-NOESY and transfer WEFT-NOESY NMR spectra were used to determine the heme proton assignments for Rhodobacter capsulatus ferricytochrome c2. The Fermi contact and pseudo-contact contributions to the paramagnetic effect of the unpaired electron in the oxidized state were evaluated for the heme and ligand protons. The chemical shift assignments for the 1H and 15N NMR spectra were obtained by a combination of 1H-1H and 1H-15N two-dimensional NMR spectroscopy. The short-range nuclear Overhauser effect (NOE) data are consistent with the view that the secondary structure for the oxidized state of this protein closely approximates that of the reduced form, but with redox-related conformational changes between the two redox states. To understand the decrease in stability of the oxidized state of this cytochrome c2 compared to the reduced form, the structural difference between the two redox states were analyzed by the differences in the NOE intensities, pseudo-contact shifts and the hydrogen-deuterium exchange rates of the amide protons. We find that the major difference between redox states, although subtle, involve heme protein interactions, orientation of the heme ligands, differences in hydrogen bond networks and, possible alterations in the position of some internal water molecules. Thus, it appears that the general destabilization of cytochrome c2, which occurs on oxidation, is consistent with the alteration of hydrogen bonds that result in changes in the internal dynamics of the protein.

Circular dichroism and resonance raman comparative studies of wild type cytochrome c and F82H mutant

The UV-visible, circular dichroism (CD), and resonance Raman (RR) spectra of the wild type yeast iso-1-cytochrome c (WT) and its mutant F82H in which phenylalanine-82 (Phe-82) is substituted with His are measured and compared for oxidized and reduced forms. The CD spectra in the intrinsic and Soret spectral region, as well as RR spectra in high, middle, and low frequency regions, are discussed. From the analysis of the spectra, it is determined that in the oxidized F82H the two axial ligands to the heme iron are His-18 and His-82 whereas in the reduced form the sixth ligand switches from His-82 to Met-80 providing the coordination geometry similar to that of WT. Based on the spectroscopic data, the conclusion is that the porphyrin macrocycle is less distorted in the oxidized F82H compared to the oxidized WT. Similar distortions are present in the reduced form of the proteins. Frequency shifts of Raman bands, as well as the decrease of the alpha-helix content in the CD spectra, indicate more open conformation of the protein around the heme.

Solution structure of oxidized microsomal rabbit cytochrome b5. Factors determining the heterogeneous binding of the heme

Cytochrome b5 is heterogeneous in solution because of the presence of two isomers (A and B), differing in the rotation of the heme plane around the axis defined by the alpha and gamma meso protons. For rabbit cytochrome b5, the A/B ratio is 5 : 1. The solution structure of the major form of the oxidized soluble fragment of rabbit microsomal cytochrome b5 (94 amino acids) is here solved through NMR spectroscopy. From 1908 NOEs, of which 1469 were meaningful, there were 246 pseudocontact shifts and 18 3J couplings, a family of 40 energy-minimized conformers were obtained with average backbone rmsd (for residues 4-84) of 0.060 +/- 0.016 nm and average target function of 0.0078 nm2, no distance violations being larger than 0.03 nm. The structure was compared with the solution structures of the A (major) and B (minor) isomers of the rat cytochrome in the oxidized form. The A/B ratio for the rat cytochrome is 1.5 : 1, despite the very high sequence similarity (93%) to the rabbit protein. This comparison has provided insights into the factors determining the distribution in solution of the two isomers differing with respect to heme orientation. It appears that residues 23 and 74 are both important in determining this distribution, through interaction of their side chains with the prosthetic group. Hydrophobic and steric interactions are the key factors in determining the relative stability of one isomer with respect to the other.

Solution conformation of ferricytochrome c-551 from Pseudomonas stutzeri substrain ZoBell

The main chain protons and the majority of side chain protons have been assigned for the ferric form of Pseudomonas stutzeri substrain ZoBell (American Type Culture Collection 14405) cytochrome c-551. The chemical shifts were compared to those for the ferrous protein to determine the pseudocontact shift contribution. These observed values were compared to contributions calculated from the atomic coordinates of the ferrous cytochrome and an optimized effective room temperature g-tensor centered on the paramagnetic ferric iron. The agreement between observed and calculated values indicates that the conformations of the two forms are highly similar.

The structure of the complex of plastocyanin and cytochrome f, determined by paramagnetic NMR and restrained rigid-body molecular dynamics

BACKGROUND

The reduction of plastocyanin by cytochrome f is part of the chain of photosynthetic electron transfer reactions that links photosystems II and I. The reaction is rapid and is influenced by charged residues on both proteins. Previously determined structures show that the plastocyanin copper and cytochrome f haem redox centres are some distance apart from the relevant charged sidechains, and until now it was unclear how a transient electrostatic complex can be formed that brings the redox centres sufficiently close for a rapid reaction.

RESULTS

A new approach was used to determine the structure of the transient complex between cytochrome f and plastocyanin. Diamagnetic chemical shift changes and intermolecular pseudocontact shifts in the NMR spectrum of plastocyanin were used as input in restrained rigid-body molecular dynamics calculations. An ensemble of ten structures was obtained, in which the root mean square deviation of the plastocyanin position relative to cytochrome f is 1.0 A. Electrostatic interaction is maintained at the same time as the hydrophobic side of plastocyanin makes close contact with the haem area, thus providing a short electron transfer pathway (Fe-Cu distance 10.9 A) via residues Tyr1 or Phe4 (cytochrome f) and the copper ligand His87 (plastocyanin).

CONCLUSIONS

The combined use of diamagnetic and paramagnetic chemical shift changes makes it possible to obtain detailed information about the structure of a transient complex of redox proteins. The structure suggests that the electrostatic interactions 'guide' the partners into a position that is optimal for electron transfer, and which may be stabilised by short-range interactions.

Native tertiary structure in an A-state

The A-state is an equilibrium species that is thought to represent the molten globule, an on-pathway protein folding intermediate with native secondary structure and non-native, fluctuating tertiary structure. We used yeast iso-1-ferricytochrome c to test for an evolutionary-invariant tertiary interaction in its A-state. Thermal denaturation monitored by circular dichroism (CD)spectropolarimetry was used to determine A-state and native-state stabilities, delta GA reversible D and delta GN reversible D. We examined the wild-type protein, seven variants with substitutions at the interface between the N and C-terminal helices, and four control variants. The controls have the same amino acid changes as the interface variants, but the changes are close to, not at, the interface. We also examined the pH and sulfate concentration dependencies and found that while these factors affect the far-UV CD spectra of the least stable variants, they do not alter the difference in stability between the wild-type protein and the variants. A delta GA reversible D versus-delta GN reversible D plot for the interface variants has a slope near unity and the control variants have near-wild-type stability. These results show that the helix-helix interaction stabilizes the A-state and the native state to the same degree, confirming our preliminary report. We determined that the heat capacity change for A-state denaturation is approximately 60% of the value for native-state denaturation, indicating that the A-state interior is native-like. We discuss our results in relation to ferricytochrome c folding kinetics.

Solution structure of reduced microsomal rat cytochrome b5

The solution structure of the major form of the reduced soluble fragment of rat microsomal cytochrome b5 has been solved through 1H-NMR spectroscopy. The protein contains 98 amino acids. Proton assignment was available for residues 1-94, except 90 [Guiles, R. D., Basus, V. J., Kuntz, I. D. & Waskell, L. (1992) Biochemistry 31, 11,365-11,375] and has been confirmed. From 1722 NOEs, of which 1203 were found to be meaningful, a family of 40 energy-minimized structures has been obtained with average backbone rmsd (for residues 5-89) of 0.078 +/- 0.018 nm and average target function of 0.0045 nm2, no distance violations being larger than 0.029 nm. The structure has been compared with the X-ray structure of the oxidized rat mitochondrial isoenzyme and with that of the highly similar bovine microsomal isoenzyme in the oxidized form. The analysis of the elements of secondary structure is instructive in terms of their stability and of their occurrence in related structures, and of the capability of NMR and X-ray spectroscopy to observe them. Some detailed structural variations are noticed among the solved structures of the various isoenzymes and between solid and solution. The structural features in solution of the residues proposed to be involved in protein-protein recognition are found to be largely conserved with respect to the solid state.

Analysis of the structure and stability of omega loop A replacements in yeast iso-1-cytochrome c

Omega (omega)-loop A, residues 18-32 in wild-type yeast iso-1-cytochrome c, has been deleted and replaced with loop sequences from three other cytochromes c and one from esterase. Yeast expressing a partial loop deletion do not contain perceptible amounts of holoprotein as measured by low-temperature spectroscopy and cannot grow on nonfermentable media. Strains expressing loop replacement mutations accumulate holoprotein in vivo, but the protein function varies depending on the sequence and length of the replacement loop; in vivo expression levels do not correlate with their thermal denaturation temperatures. In vitro spectroscopic studies of the loop replacement proteins indicate that all fold into a native-like cytochrome c conformation, but are less stable than the wild-type protein. Decreases in thermal stability are caused by perturbation of loop C backbone in one case and a slight reorganization of the protein hydrophobic core in another case, rather than rearrangement of the loop A backbone. A single-site mutation in one of the replacement mutants designed to relieve inefficient hydrophobic core packing caused by the new loop recovers some, but not all, of the lost stability.

Analysis of the 1H-NMR chemical shifts of Cu(I)-, Cu(II)- and Cd-substituted pea plastocyanin. Metal-dependent differences in the hydrogen-bond network around the copper site

To compare cadmium-substituted plastocyanin with copper plastocyanin, the 1H-NMR spectra of CuI-, CuII- and Cd-plastocyanin from pea have been analyzed. Full assignments of the spectra of CuI- and Cd-plastocyanin indicate chemical shift differences up to 1 ppm. The affected protons are located in the four loops that surround the Cu site. The largest differences were found for protons in the hydrogen bond network which stabilizes this part of the protein. This suggests that the chemical shift differences are caused by very small but extensive structural changes in the network upon replacement of CuI by Cd. For CuII-plastocyanin the resonances of 72% of the protons observed in the CuI form have been identified. Protons within approximately 0.9 nm of the CuII were not observed due to fast paramagnetic relaxation. The protons between 0.9-1.7 nm from the CuII showed chemical shift differences up to 0.4 ppm compared to both CuI- and Cd-plastocyanin. These differences can be predicted assuming that they represent pseudocontact shifts. When corrected for the pseudocontact shift contribution, the CuII-plastocyanin chemical shifts were nearly all identical within error to those of the Cd form, but not of the CuI-plastocyanin, indicating that the CuII-plastocyanin structure, in as far as it can be observed, resembles Cd-rather than CuI-plastocyanin. In a single stretch of residues (64-69) chemical shift differences remained between all three forms after correction. The fact that pseudocontact shifts were observed for protons which were not broadened may be attributable to the weaker distance dependence of the pseudocontact shift effect compared to paramagnetic relaxation. This results in two shells around the Cu atom, an inner paramagnetic shell (0-0.9 nm), in which protons are not observed due to broadening, and an outer paramagnetic shell (0.9-1.7 nm), in which protons can be observed and show pseudocontact shifts. It is concluded that Cd-plastocyanin is a suitable redox-inactive substitute for Cu-plastocyanin.

The role of a conserved water molecule in the redox-dependent thermal stability of iso-1-cytochrome c

Eukaryotic cytochromes c contain a buried water molecule (Wat166) next to the heme that is associated through a network of hydrogen bonds to three invariant residues: tyrosine 67, asparagine 52, and threonine 78. Single-site mutations to two of these residues (Y67F, N52I, N52A) and the double-site mutation (Y67F/N52I) were introduced into Saccharomyces cerevisiae iso-1-cytochrome c to disrupt the hydrogen bonding network associated with Wat166. The N52I and Y67F/N52I mutations lead to a loss of Wat166 while N52A and Y67F modifications lead to the addition of a new water molecule (Wat166) at an adjacent site (Berghuis, A. M., Guillemette, J. G., McLendon, G., Sherman, F., Smith, M., and Brayer, G. D. (1994) J. Mol. Biol. 236, 786-799; Berghuis, A. M., Guillemette, J. G., Smith, M., and Brayer, G. D. (1994) J. Mol. Biol. 235, 1326-1341; Rafferty, S. P., Guillemette, J. G., Berghuis, A. M., Smith, M., Brayer, G. D., and Mauk, A. G. (1996) Biochemistry, 35, 10784-10792). We used differential scanning calorimetry (DSC) to determine the change in heat capacity (DeltaCp) and the temperature dependent enthalpy (DeltaHvH) for the thermal denaturation of both the oxidized and reduced forms of the iso-1 cytochrome c variants. The relative stabilities were expressed as the difference in the free energy of denaturation (DeltaGD) between the wild type and mutant proteins in both redox states. The disruption of the hydrogen bonding network results in increased stability for all of the mutant proteins in both redox states with the exception of the reduced Y67F variant which has approximately the same stability as the reduced wild type protein. For the oxidized proteins, DeltaGD values of 1.3, 4.1, 1.5, and 5.8 kcal/mol were determined for N52A, N52I, Y67F, and Y67F/N52I, respectively. The oxidized proteins were 8.2-11.5 kcal/mol less stable than the reduced proteins due to a redox-dependent increase in the entropy of unfolding.

Requirements for perpendicular helix pairing

Cassette mutagenesis was used to produce a library of mutations at the interface of the N- and C-terminal helices of Saccharomyces cerevisiae iso-1-cytochrome c. The library is random and comprises > 98% mutations. Over 11,000 candidates were assayed for function by selecting for the ability of yeast, with the mutated gene as their sole cytochrome c source, to grow on nonfermentable carbon sources. We estimate that approximately 0.5% of the 160,000 total amino acid combinations at these four residues result in a functional cytochrome c. Significant correlations are found between the phenotype of yeast harboring the alleles and both the Dayhoff mutation matrix and transfer free energies (cyclohexane-to-water and n-octanol-to-water). Similar correlations are observed with respect to growth rate. Finally, sequences that are consistent with function follow a binary amino acid pattern.

An optimized g-tensor for Rhodobacter capsulatus cytochrome c2 in solution: a structural comparison of the reduced and oxidized states

The optimized g-tensor parameters for the oxidized form of Rhodobacter capsulatus cytochrome c2 in solution were obtained using a set (50) of backbone amide protons. Dipolar shifts for more than 500 individual protons of R. capsulatus cytochrome c2 have been calculated by using the optimized g-tensor and the X-ray crystallographic coordinates of the reduced form of R. capsulatus cytochrome c2. The calculated results for dipolar shifts are compared with the observed paramagnetic shifts. The calculated and the observed data are in good agreement throughout the entire protein, but there are significant differences between calculated and experimental results localized to the regions in the immediate vicinity of the heme ligand and the region of the front crevice of the protein (residues 44-50, 53-57, and 61-68). The results not only indicate that the overall solution structures are very similar in both the reduced and oxidized states, but that these structures in solution are similar to the crystal structure. However, there are small structural changes near the heme and the rearrangement of certain residues that result in changes in their hydrogen bonding concomitant with the change in the oxidation states; this was also evident in the data for the NH exchange rate measurements for R. capsulatus cytochrome c2.

Possible Role of the Iron Coordination Sphere in Hemoprotein Electron Transfer Self-Exchange: (1)H NMR Study of the Cytochrome c-PMe(3) Complex

The rates of self-exchange electron transfer in the trimethylphosphine complex of cytochrome c have been measured by an NMR technique over a large range of ionic strengths. The rate constant is 1.56 x 10(4) M(-)(1) s(-)(1) at 23 degrees C (&mgr; = 0.34 M) at pH 6.9. Dependence on ionic strength of the rate constant is treated by van Leeuwen theory. Extrapolation of the rate constant to infinite ionic strength gives a rate constant of 3.9 x 10(5) M(-)(1) s(-)(1). This rate constant is compared with others reported for myoglobin and cytochrome b(5)(). The values for these systems range over 2 orders of magnitude with myoglobin-PMe(3) << cytochrome b(5)() < cytochrome c-PMe(3) < cytochrome c. Analysis of the data in terms of Marcus theory gives a reorganization energy, lambda, for self-exchange of 0.75 eV mol(-)(1) for cytochrome c-PMe(3). Substitution of Met-80 by PMe(3) appears to influence only weakly the rearrangement barrier to electron transfer.

Energised (entatic) states of groups and of secondary structures in proteins and metalloproteins

In this review I have examined the functional value of selected states of isolated groups in proteins energised away from their expected ground states whether they are observed with or without energy perturbation of larger parts of a protein structure. These energisations, found in the absence of substrates, are called 'entatic states of groups' [Vallee, B. L. & Williams, R. J. P. (1968) Proc. Natl Acad. Sci. USA 59, 498-505]. A group can be part of an amino acid or any bound metal ion or cofactor. In some particular cases the apoprotein, where the energised metal ion or cofactor has been removed, or the protein in which the energised amino acid has been replaced, has the same back-bone structure as the holoprotein and even side-chains are only slightly adjusted. This case is quite different from a condition of a group simultaneously energised with a protein fold, due to their combination, and which therefore involves conformational change in the protein and the group and which may adjust the group while the protein tightens. The final condition is again stable but removal of the group now must result in a reversed protein conformation change. This simultaneous energisation of both the adjustable protein and the group can be local, as in an induced fit or more extensive when it can be likened in some cases to a stretching by rack action, or may involve a change from an almost random to a structured protein when the group is energised in a very limited way. The various energisations must not be confused since they differ functionally. The first can give rise to optimal heightened catalytic (or other functional) potential of the local group but cannot be connected either to excitation of other parts of a protein as in induced fitting, or to a relay of energy (larger conformational change) in the protein. Clearly it restricts the rate of exchange of a group. Induced fit can also give rise to group activation, though to a somewhat reduced degree, while increasing exchange rate. A device such as a rack may rather give rise to a mechanical activity (message transmission), which is relayed a large distance into the protein, and can only give considerably lower activation of individual groups but exchange may now be fast. The final case involves very modest energisation of the group with gross rearrangement and energisation of the protein and may be associated with storage or carrier functions. The groups upon which I concentrate are metal ions since detailed electronic and structural knowledge of their ground states are well known, allowing energised states to be easily detected, but the ideas apply equally to organic side-chains of proteins as will be shown. A further energisation can arise from the addition of a substrate to each kind of protein. In fact all the ideas of energisation applicable to groups having cyclic activity in permanent features of protein structure are equally well applied to substrate binding or conversion of substrates through excited states to products.

A native tertiary interaction stabilizes the A state of cytochrome c

Certain kinetic intermediates in protein folding are similar to the molten globule, or A state, an equilibrium state of many proteins that is populated under high salt and low pH conditions. Many A states are nearly as compact as native proteins and have native-like secondary structure, but the extent to which nonlocal interactions stabilize the A state is unclear. In this study, thermal denaturation, monitored by circular dichroism, was used to determine the free energy of denaturation of the A state (delta GA<-->D) for Saccharomyces cerevisiae iso-1-ferricytochrome c. Specifically, we examined the wild-type protein, seven variants with amino acid substitutions at the interface between the N- and C-terminal helices, and two variants with mutations at a position close to, but not involved in, the interface. A plot of delta GA<-->D versus delta GN<-->D (the free energy of denaturation of the native state) has a slope near unity, showing that the evolutionarily conserved helix-helix interaction stabilizes the A state to the same degree that it stabilizes the native state.

Protein thermal denaturation, side-chain models, and evolution: amino acid substitutions at a conserved helix-helix interface

Random mutant libraries with substitutions at the interface between the N- and C-terminal helices of Saccharomyces cerevisiae iso-1-cytochrome c were screened. All residue combinations that have been identified in naturally occurring cytochrome c sequences are found in the libraries. Mutants with these combinations are biologically functional. Enthalpies, heat capacities, and midpoint temperatures of denaturation are used to determine the entropy and Gibbs free energy of denaturation (delta GD) for the ferri form of the wild-type protein and 13 interface variants. Changes in delta GD cannot be allocated solely to enthalpic or entropic effects, but there is no evidence of enthalpy-entropy compensation. The lack of additivity of delta GD values for single versus multiple amino acid substitutions indicates that the helices interact thermodynamically. Changes in delta GD are not in accord with helix propensities, indicating that interactions between the helices and the rest of the protein outweigh helix propensity. Comparison of delta GD values for the interface variants and nearly 90 non-cytochrome c variants to side-chain model data leads to several conclusions. First, hydrocarbon side chains react to burial-like transfer from water to cyclohexane, but even weakly polar side chains respond differently. Second, despite octanol being a poor model for protein interiors, octanol-to-water transfer free energies are useful stability predictors for changing large hydrocarbon side chains to smaller ones. Third, unlike cyclohexane and octanol, the Dayhoff mutation matrix predicts stability changes for a variety of substitutions, even at interacting sites.(ABSTRACT TRUNCATED AT 250 WORDS)

Protein structure refinement based on paramagnetic NMR shifts: applications to wild-type and mutant forms of cytochrome c

A new approach to NMR solution structure refinement is introduced that uses paramagnetic effects on nuclear chemical shifts as constraints in energy minimization or molecular dynamics calculations. Chemical shift differences between oxidized and reduced forms of horse cytochrome c for more than 300 protons were used as constraints to refine the structure of the wild-type protein in solution and to define the structural changes induced by a Leu 94 to Val mutation. A single round of constrained minimization, using the crystal structure as the starting point, converged to a low-energy structure with an RMS deviation between calculated and observed pseudo-contact shifts of 0.045 ppm, 7.5-fold lower than the starting structure. At the same time, the procedure provided stereospecific assignments for more than 45 pairs of methylene protons and methyl groups. Structural changes caused by the mutation were determined to a precision of better than 0.3 A. Structure determination based on dipolar paramagnetic (pseudocontact) shifts is applicable to molecules containing anisotropic paramagnetic centers with short electronic relaxation times, including numerous naturally occurring metalloproteins, as well as proteins or nucleic acids to which a paramagnetic metal ion or ligand may be attached. The long range of paramagnetic shift effects (up to 20 A from the iron in the case of cytochrome c) provides global structural constraints, which, in conjunction with conventional NMR distance and dihedral angle constraints, will enhance the precision of NMR solution structure determination.

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.

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.

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.

Novel heteronuclear methods of assignment transfer from a diamagnetic to a paramagnetic protein: application to rat cytochrome b5

15N and 1H resonance assignments for backbone and side-chain resonances of both equilibrium forms of rat ferricytochrome b5 have been obtained, using a combination of novel heteronuclear assignment transfer methods from the known assignments of the diamagnetic protein [Guiles, R. D., Basus, V. J., Kuntz, I. D., & Waskell, L. A. (1992) Biochemistry 31, 11365-11375] and computational methods which depend on an accurate determination of the orientation of the components of the susceptibility tensor. The transfer of amide proton resonance assignments takes advantage of the apparent insensitivity of amide 15N resonances to pseudocontact effects, evident in overlays of 15N-1H heteronuclear correlation spectra. Amide-proton resonance assignments tentatively transferred from the known diamagnetic assignments to the paramagnetic form of the protein were confirmed using conventional assignment strategies employing 600-MHz COSY, HOHAHA, and NOESY spectra of the oxidized protein. As was observed in rat ferrocytochrome b5, more than 40% of all residues exhibited NMR detectable heterogeneity due to the two different orientations of the heme. Complete assignment of both forms enabled accurate determination of the orientation of the susceptibility tensor for both conformations of the heme. The orientation of the z-component of the susceptibility tensors for the two forms are indistinguishable, while the in-plane components appear to differ by about 6 degrees. Differences in the orientation of the in-plane susceptibility components are undoubtedly due dominantly to the relative axial rotation of the heme of between 5 degrees and 10 degrees indicated by the NOESY contacts to the protein observed in the spectra of the ferrocytochrome [Guiles, R. D., Basus, V. J., Kuntz, I. D., & Waskell, L. A. (1992) Biochemistry 31, 11365-11375; Pochapsky, T. C., Sligar, S. G., McLachlan, S. J., & LaMar, G. N. (1990) J. Am. Chem. Soc. 112, 5258-5263].

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.

The history of proton-driven ATP formation

This article sets down the beginnings of some thoughts in bio-energetics. It illustrates how difficult it is in science as elsewhere to know how a new idea is generated. The literature needs very careful examination and separation from personalities.

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.

Redox-dependent changes in beta-extended chain and turn structures of cytochrome c in water solution determined by second derivative amide I infrared spectra

The redox-dependent changes in secondary structure of cytochromes c from horse, cow, and dog hearts in water at 20 degrees C have been determined by amide I infrared spectroscopy. Second derivative amide I spectra were obtained by use of a procedure that includes a convenient method for the effective subtraction of the spectrum of water vapor in the system. The band at 1657 cm-1 representing the helix structure was unaffected by a change in redox state whereas changes in bands due to turns at 1680, 1672, and 1666 cm-1, unordered structure at 1650 cm-1, and beta-structures at 1632 and 1627 cm-1 occurred. About one-fourth of the beta-extended chain spectral region and one-fifth of the beta-turn region (involving a total of approximately 9-13 residues) were sensitive to the oxidation state of heme iron. No significant changes in the secondary structure of either the reduced or oxidized protein due to changes in ionic strength were detected. The localized structural rearrangements triggered by the changes in oxidation state of heme iron are consistent with differences in the binding of heme iron to a histidine imidazole nitrogen and a methionine sulfur atom from the beta-extended chain. The demonstrated ability to obtain highly reproducible second derivative amide I infrared spectra confirms the unique utility of such spectral measurements for localization of subtle changes in secondary structure within a protein, especially for changes among the multiple turns and beta-structures.

Temperature-sensitive variants of Saccharomyces cerevisiae iso-1-cytochrome c produced by random mutagenesis of codons 43 to 54

In vitro random mutagenesis within the CYC1 gene from the yeast Saccharomyces cerevisiae was used to produce a library of mutants encompassing codons 43 to 54 of iso-1-cytochrome c. This region consists of an evolutionarily conserved structure within an evolutionarily diverse sequence. The library, on a low-copy-number yeast shuttle phagemid, was introduced into a yeast strain lacking cytochrome c. The ability of transformants harboring a functional cytochrome c to grow on the non-fermentable carbon source glycerol at 30 degrees C and 37 degrees C was used to determine the phenotype of nearly 1000 transformants. Approximately 90% of the missense mutants present in the library give rise to the wild-type phenotype, 7% result in the temperature-sensitive (Cycts) phenotype, and 3% give rise to the non-functional (Cyc-) phenotype. Phagemids from 20 Cycts and 30 Cyc- transformants were subjected to DNA sequence analysis. All the mutations occur within the targeted region. One-third of the mutants from Cyc- transformants and all the mutants from Cycts transformants are missense mutants. The remaining mutants from Cyc- transformants are nonsense or frame-shift mutants. Missense mutations within the codons for Gly45, Tyr46, Thr49, Asn52 or Ile53 alone are sufficient to produce temperature-sensitive behavior both in vivo and in the variant proteins. The deduced amino acid substitutions correlate remarkably well with side-chain dynamics, secondary structure and tertiary structure of the wild-type protein.

Constraints on amino acid substitutions in the N-terminal helix of cytochrome c explored by random mutagenesis

The interaction of the N- and C-terminal helices is a hallmark of the cytochrome c family. Oligodeoxyribonucleotide-directed random mutagenesis within the gene encoding the C102T protein variant of Saccharomyces cerevisiae iso-1-cytochrome c was used to generate a library of mutations at the evolutionary invariant residues Gly-6 and Phe-10 in the N-terminal helix. Transformation of this library (contained on a low-copy-number yeast shuttle phagemid) into a yeast strain lacking a functional cytochrome c, followed by selection for cytochrome c function, reveals that 4-10% of the 400 possible amino acid substitutions are compatible with function. DNA sequence analysis of phagemids isolated from transformants exhibiting the functional phenotype elucidates the requirements for a stable helical interface. Basic residues are not tolerated at position 6 or 10. There is a broad volume constraint for amino acids at position 6. The amino acid substitutions observed to be compatible with function at Phe-10 show that the hydrophobic effect alone is sufficient to promote helical association. There are severe constraints that limit the combinations consistent with function, but the number of functionally consistent combinations observed exemplifies the plasticity of proteins.

The role of the internal hydrogen bond network in first-order protein electron transfer between Saccharomyces cerevisiae iso-1-cytochrome c and bovine microsomal cytochrome b5

An internal water molecule (designated WAT166) is found in iso-1-cytochrome c which is part of a redox-state-dependent hydrogen bond network. The position of this water molecule with respect to the polypeptide fold can be altered or even displaced by site-directed mutagenesis leading to structural perturbations and associated changes in redox potential. Using saturation transfer 1H-NMR methods, this study measures changes in the electron transfer reactivity for three variants of yeast iso-1-cytochromes c in which the position of this water molecule is altered. In particular, the reverse electron transfer rate is measured within a complex formed between either wild-type or variant yeast iso-1-cytochromes c and the tryptic fragment of bovine liver microsomal cytochrome b5. For three variants of yeast iso-1-cytochrome c the rate constants measured by saturation transfer are wild-type (Asn52, E0 = 270 mV, kex = 0.3 s-1), Asn52----Ala (E0 = 240 mV, kex = 0.6 s-1), Asn52----Ile (E0 = 220 mV, kex = 1.0 s-1). The first-order rates are compared with that of a fourth variant Phe82----Gly which has been measured previously (E0 = 220 mV, kex = 0.7 s-1). An analysis of the variation in the observed cross exchange rate using Marcus theory shows that these changes can be predicted quantitatively by the shift in redox potential that accompanies mutagenesis. So, although the perturbation of the internal water molecule by mutagenesis alters both the structure and redox potential of cytochrome c, surprisingly it does not significantly influence the intrinsic electron transfer reactivity of the protein. Studies of the activation parameters suggests that a variation of temperature changes both delta G* and also the prefactor. These data are discussed in terms of models involving dynamic molecular recognition between proteins.

Proton nuclear magnetic resonance as a probe of differences in structure between the C102T and F82S,C102T variants of iso-1-cytochrome c from the yeast Saccharomyces cerevisiae

Differences in chemical shifts and in nuclear Overhauser effects between the C102T and F82S,C102T variants of Saccharomyces cerevisiae iso-1-cytochrome c in both the reduced and oxidized forms are reported and analyzed. There is evidence for small conformational differences in both oxidation states of the double variant near position 82. Differences in structure are more evident in the oxidized forms of the variants. These differences extend to distant parts of the protein. It is concluded that the oxidized double variant has undergone a small rearrangement of several regions of the protein that are linked by a hydrogen-bond network. It is shown that the rearrangement involves hydrogen bonds associated with the two heme propionates and associated water molecules. The deductions from nuclear magnetic resonance data are compared with the differences in the crystal structures of the reduced forms of wild-type protein and the F82S variant [Louie, G. V., Pielak, G. J., Smith, M., & Brayer, G. D. (1988) Biochemistry 27, 7870-7876].