The role of histidines 26 and 33 in the structural stabilization of cytochrome c

Screen-printed graphite macroelectrodes for the direct electron transfer of cytochrome c: a deeper study of the effect of pH on the conformational states, immobilization and peroxidase activity

The direct electron transfer of cytochrome c has been studied at screen-printed graphite macroelectrodes without recourse to mediators or the need for any electrode pre-treatment as is commonly employed within the literature. A wide range of pH values from 2.0 to 11.0 have been explored upon the electrochemical response of cytochrome c and different voltammetric signatures have been observed. The direct electron transfer of the alkaline transition of cytochrome c was found impeded within alkaline media leading to either an irreversible redox process or even no voltammetric responses. In acidic aqueous media the electrochemical process is observed to undergo a mixed diffusion and adsorption controlled process rather than a purely diffusional process of the native conformation as observed at pH 7.0. Interestingly, at pH 3.5 a new conformational state is revealed in cooperation with the native conformation. The immobilization of the protein was satisfactorily obtained using a simple method by cycling the protein at specific solution pH values allowing amperometric responses to be obtained and gives rise to useful pseudo-peroxidase activity for sensing H2O2. Apparent Michaelis-Menten constant values (Km) were calculated via the Lineweaver-Burk method with deduced values of 25 ± 4, 98 ± 12 and 230 ± 30 mM, respectively for pH values of 2.0, 3.0 and 7.0. Such work is important for those utilising cytochrome c in bio-electrochemical and related applications.

Compressing the free energy range of substructure stabilities in iso-1-cytochrome c

Evolutionary conservation of substructure architecture between yeast iso-1-cytochrome c and the well-characterized horse cytochrome c is studied with limited proteolysis, the alkaline conformational transition and global unfolding with guanidine-HCl. Mass spectral analysis of limited proteolysis cleavage products for iso-1-cytochrome c show that its least stable substructure is the same as horse cytochrome c. The limited proteolysis data yield a free energy of 3.8 +/- 0.4 kcal mol(-1) to unfold the least stable substructure compared with 5.05 +/- 0.30 kcal mol(-1) for global unfolding of iso-1-cytochrome c. Thus, substructure stabilities of iso-1-cytochrome c span only approximately 1.2 kcal mol(-1) compared with approximately 8 kcal mol(-1) for horse cytochrome c. Consistent with the less cooperative folding thus expected for the horse protein, the guanidine-HCl m-values are approximately 3 kcal mol(-1)M(-1) versus approximately 4.5 kcal mol(-1)M(-1) for horse versus yeast cytochrome c. The tight free energy spacing of the yeast cytochrome c substructures suggests that its folding has more branch points than for horse cytochrome c. Studies on a variant of iso-1-cytochrome c with an H26N mutation indicate that the least and most stable substructures unfold sequentially and the two least stable substructures unfold independently as for horse cytochrome c. Thus, important aspects of the substructure architecture of horse cytochrome c, albeit compressed energetically, are preserved evolutionally in yeast iso-1-cytochrome c.

Characterization of the polyanion-induced molten globule-like state of cytochromec

Cytochrome c (cyt c) undergoes a poly(vinylsulphate) (PVS)-induced transition at slightly acidic pH into a molten globule-like state that resembles the effect that negatively charged membrane surfaces have on this protein. In this work, the thermodynamic properties of the molten globule-like state of cyt c in complex with PVS are studied using differential scanning calorimetry, circular dichroism, fluorescence, and absorbance spectroscopy. The temperature-induced transition of the molten globule-like state of cyt c in the complex with PVS is characterized by a significantly lower calorimetric enthalpy than in the "typical" molten globule state of cyt c, i.e. free protein at pH 2.0 in high ionic strength. Moreover, the thermally-denatured state of cyt c in the complex at pH < 6 contains nearly 50% of the native secondary structure. The dependence of the transition temperature on the pH indicates a role for histidine residues in the destabilization of the cyt c structure in the PVS complex and in stabilization of the denatured state with the residual secondary structure. A comparison of the effects of small anions and polyanions demonstrates the importance of cooperativity among the anions in the destabilization of cyt c. Predictably, other hydrophilic flexible polyanions such as heparin, polyglutamate, and polyadenylate also have a destabilizing effect on the structure of cyt c. However, a correlation between the properties of the polyanions and their effect on the protein stability is still unclear.

Cooperative omega loops in cytochrome c: role in folding and function

Hydrogen exchange experiments under slow exchange conditions show that an omega loop in cytochrome c (residues 40-57) acts as a cooperative unfolding/refolding unit under native conditions. This unit behavior accounts for an initial step on the unfolding pathway, a final step in refolding, and a number of other structural, functional and evolutionary properties.

Denatured state thermodynamics: residual structure, chain stiffness and scaling factors

A set of nine variants of yeast iso-1-cytochrome c with zero or one surface histidine have been engineered such that the N-terminal amino group is acetylated in vivo. N-terminal acetylation has been confirmed by mass spectral analysis of intact and proteolytically digested protein. The histidine-heme loop-forming equilibrium, under denaturing conditions (3 M guanidine hydrochloride), has been measured by pH titration providing an observed pK(a), pK(a)(obs), for each variant. N-terminal acetylation prevents the N-terminal amino group-heme binding equilibrium from interfering with measurements of histidine-heme affinity. Significant deviation is observed from the linear dependence of pK(a)(obs) on the log of the number of monomers in the loop formed, expected for a random coil denatured state. The maximum histidine-heme affinity occurs for a loop size of 37 monomers. For loop sizes of 37-83 monomers, histidine-heme pK(a)(obs) values are consistent with a scaling factor of -4.2+/-0.3. This value is much larger than the scaling factor of -1.5 for a freely jointed random coil, which is commonly used to represent the conformational properties of protein denatured states. For loop sizes of nine to 22 monomers, chain stiffness is likely responsible for the decreases in histidine-heme affinity relative to a loop size of 37. The results are discussed in terms of residual structure and sequence composition effects on the conformational properties of the denatured states of proteins.

The heme-containing N-fragment (residues 1-56) of cytochrome c is a bis-histidine functional system

The structural and redox properties of a heme-containing fragment (1-56 residues) of cytochrome c have been investigated by spectroscopic (circular dichroism, electronic absorption, and EPR) and voltammetric techniques. The results indicate that the N-fragment lacks ordered secondary structure and has two histidines axially bound to the heme-iron (the native His18 and a misligated His26 or His33). Despite the absence of ordered secondary structure, the peptide chain shields the heme group from solvent, as shown by (i) the pK(a) of protonation of the nonnative histidine ligand (5.18 +/- 0.05), lower than that of the bis-histidine guanidine-unfolded cytochrome c (5.58 +/- 0.05), and (ii) the redox potential, E(o) = 0 +/- 5 mV versus NHE, close to that of bis-histidine cytochrome c mutants but less negative than that of bis-histidine complexes of microperoxidase with short peptides. The electroactive N-fragment may be taken as a "minichrome c" model, with interesting potential for application to biosensor technology; further, the system provides useful information for a deeper understanding of cytochrome c folding and structural/functional organization.

Measuring denatured state energetics: deviations from random coil behavior and implications for the folding of iso-1-cytochrome c

The changes in the free energy of the denatured state of a set of yeast iso-1-cytochrome c variants with single surface histidine residues have been measured in 3 M guanidine hydrochloride. The thermodynamics of unfolding by guanidine hydrochloride is also reported. All variants have decreased stability relative to the wild-type protein. The free energy of the denatured state was determined in 3 M guanidine hydrochloride by evaluating the strength of heme-histidine ligation through determination of the pK(a) for loss of histidine binding to the heme. The data are corrected for the presence of the N-terminal amino group which also ligates to the heme under similar solution conditions. Significant deviations from random coil behavior are observed. Relative to a variant with a single histidine at position 26, residual structure of the order of -1.0 to -2.5 kcal/mol is seen for the other variants studied. The data explain the slower folding of yeast iso-1-cytochrome c relative to the horse protein. The greater number of histidines and the greater strength of ligation are expected to slow conversion of the histidine-misligated forms to the obligatory aquo-heme intermediate during the ligand exchange phase of folding. The particularly strong association of histidine residues at positions 54 and 89 may indicate regions of the protein with strong energetic propensities to collapse against the heme during early folding events, consistent with available data in the literature on early folding events for cytochrome c.

Formation of a molten-globule-like state of cytochrome c induced by high concentrations of glycerol

The effect of glycerol on the structure of cytochrome c was investigated by circular dichroism, absorbance and EPR spectroscopy. The results obtained show that an increasing concentration of the organic solvent (70-99.2%, v/v) in aqueous-polyalcohol mixtures converts native cytochrome c into a new, low spin form through a fully reversible, two-state transition. The glycerol-stabilized form (that we call here the G state) retains native-like amounts of alpha-helix structure while rigid tertiary structure and native Fe(III)-Met(80) axial bond are lost. Analysis of data suggests a molten globule character of the G state; support to this view is afforded by the striking similarities between the spectroscopic (and, thus, structural) properties of the G state with those of the acidic molten globule of the protein (A state).

The source of heterogeneity in the heme vicinity of ferricytochrome c

Heterogeneity in the heme vicinity of ferricytochrome c was reported to be detectable by a split of the NMR signal of the heme methyl 3 group [P.D. Burns and G.N. La Mar, J. Am. Chem. Soc. 101 (1979) 5844]. Using cytochrome c mutants and computer simulations of the native and mutated cytochromes, the source of this heterogeneity is found to originate from the His-33 residue motions. The H33F mutation abolished the NMR split and computer simulations of the H33F mutant revealed a narrower distribution of fluctuations of the radius of gyration, suggesting a more rigid structure due to the mutation. The stabilization of the mutant was further demonstrated by a reduction in the H33F mutant of 4 Kcal/mol in the calculated interaction energy between residue 33 and the rest of the cytochrome, in keeping with known experimental results.

The interaction of borate ions with cytochrome c surface sites: a molecular dynamics study

Ionic interactions of cytochrome c play an important role in the electron transfer process. Molecular dynamics simulations of the binding of borate ion, which serves as a model ion, at three different cytochrome c surface sites are performed. This work is motivated by previous NMR studies of cytochrome c in borate solution, which indicate the existence of two types of binding sites, a slow exchange site and a fast exchange site. These two types of binding behavior were observed in the dynamic simulations, offering a molecular interpretation of "loose" and "tight" binding. At the "loose" binding sites (near Lys25/Lys27 and Lys55/Lys73) the ion forms two to three hydrogen bonds to the nearest lysine residue. This binding is transient on the time scale of the simulation, demonstrating the feasibility of fast exchange. At the "tight" binding site (near Lys13/Lys86), on the other hand, the ion becomes integrated into the protein hydrogen bond network and remains there for the duration of the simulation (exemplifying slow exchange). Binding simulations of the ion at the "tight" site of H26Q mutant cytochrome c also showed integration of the ion into the protein's hydrogen bond network. However, this integration differs in details from the binding of the ion to the native protein, in agreement with previous NMR observations.

Mutagenesis of histidine 26 demonstrates the importance of loop-loop and loop-protein interactions for the function of iso-1-cytochrome c

In yeast iso-1-cytochrome c, the side chain of histidine 26 (His26) attaches omega loop A to the main body of the protein by forming a hydrogen bond to the backbone atom carbonyl of glutamic acid 44. The His26 side chain also forms a stabilizing intra-loop interaction through a hydrogen bond to the backbone amide of asparagine 31. To investigate the importance of loop-protein attachment and intra-loop interactions to the structure and function of this protein, a series of site-directed and random-directed mutations were produced at His26. Yeast strains expressing these variant proteins were analyzed for their ability to grow on non-fermentable carbon sources and for their intracellular production of cytochrome c. While the data show that mutations at His26 lead to slightly decreased intracellular amounts of cytochrome c, the level of cytochrome c function is decreased more. The data suggest that cytochrome c reductase binding is affected more than cytochrome c oxidase or lactate dehydrogenase binding. We propose that mutations at this residue increase loop mobility, which, in turn, decreases the protein's ability to bind redox partners.

A 1H NMR study of structurally relevant inter-segmental hydrogen bond in cytochrome c

NMR signal arising from His 26 N(epsilon)H proton in horse and tuna ferrocytochromes c has been assigned. This His residue is highly conserved in most mitochondrial cytochromes c and X-ray crystallographic studies strongly suggested that its side-chain imidazole participates in an internal hydrogen bond network which is relevant to the stability of the non-helical protein folding near the heme active site. The shift and line width of the assigned signal indicated that this NH hydrogen is indeed involved in an internal hydrogen bond. On the basis of the X-ray crystal structures, the carbonyl oxygen of the residue at 44 is thought to act as a proton-acceptor for this hydrogen. The observation of nuclear Overhauser effect correlation between His 26 C(epsilon)H and Asn 31 main-chain amide NH proton signals in the present proteins also demonstrated the formation of the hydrogen bond between these residues. Consequently, the presence of a unique triad hydrogen bond network in these cytochromes c in solution has been confirmed. Taking advantage of the sensitivity of His 26 N(epsilon)H proton signal to the structural properties of this hydrogen bond network, influences of the presence of high concentration of salt or various concentrations of denaturant on the protein folding were inferred from the analysis of the NMR spectral parameters of the signal.

A 1H NMR comparative study of the structure of the critical packing interfaces between helix and non-helical region in various ligation states of sperm whale myoglobin

NMR signals arising from HisB5 N(delta)H and HisEF5 N(epsilon)H protons in sperm whale met-aquo, met-azido, met-imidazole and deoxy myoglobins have been assigned. HisB5 N(delta)H proton and N(epsilon) atom are hydrogen-bonded to peptide carbonyl oxygen of AspB1 and HisGH1 N(epsilon)H proton, respectively. This unique triad hydrogen bond stabilizes the interface between the B helix and the GH comer, which is critical to the stability of its protein folding. HisEF5 N(epsilon)H proton is hydrogen-bonded to one of the side-chain carboxyl oxygen atoms of AspH18 and this internal hydrogen bond stabilizes the orientation of the EF corner, the hinge region between the heme binding E and F helices, relative to the H helix. The assigned His NH proton signals were used as probes to investigate possible structural differences of these critical interfaces among various ligation states of myoglobin. Shift differences of the signals among the proteins strongly suggested that the AspB1-HisB5-HisGH1 and HisEF5-AspH18 hydrogen bonds are affected by heme ligation. The present study demonstrates the presence of a structural correlation between the active site and the critical interfaces in myoglobin.

Introduction of a proline residue into position 31 of the loop of the dimeric 4-alpha-helical protein ROP causes a drastic destabilization

The exchange of an alanine with a proline residue in position 31 of the loop region of the dimeric 4-alpha-helical-bundle protein ROP causes a reduction in the alpha-helix content of 7% and a reduction in stability of about 40% compared to the wild type parameters. The Gibbs energy of unfolding by denaturants extrapolated linearly to zero denaturant concentration, delta G0D (buffer, 25 degrees C), has been determined to be 43 kJ (mol dimer)-1. The corresponding ROPwt value is 72 kJ (mol dimer)-1 (Steif et al., 1993). The extrapolated delta G0D values obtained from urea and GdmHCI un- and refolding studies are identical within error limits. Deconvolution of the stability values into enthalpy and entropy terms resulted in the following parameters. At T1/2 = 43 degrees C (Cprotein = 0.05 mg.ml-1) the ROP A31P mutant is characterized by delta Hv.H.0 = 272 kJ (mol dimer)-1, delta Cp = 7.2 kJ (mol dimer)-1 K-1, delta S0 = 762 J (mol dimer)-1 K-1. These parameters are only approximately 50% as large as the corresponding values of ROPwt. We assume that the significant reduction in stability reflects the absence of at least one hydrogen bond as well as deformation of the protein structure. This interpretation is supported by the reduction in the change in heat capacity observed for the A31P mutant relative to ROPwt, by the increased aggregation tendency of the mutant and by the reduced specific CD absorption at 222 nm. All results support the view that in the case of ROP protein the loop region plays a significant role in the maintenance of native structure and conformational stability.

Fast folding of cytochrome c

Native iso-2 cytochrome c contains two residues (His 18, Met 80) coordinated to the covalently attached heme. On unfolding of iso-2, the His 18 ligand remains coordinated to the heme iron, whereas Met 80 is displaced by a non-native heme ligand, His 33 or His 39. To test whether non-native His-heme ligation slows folding, we have constructed a double mutant protein in which the non-native ligands are replaced by asparagine and lysine, respectively (H33N,H39K iso-2). The double mutant protein, which cannot form non-native histidine-heme coordinate bonds, folds significantly faster than normal iso-2 cytochrome c: gamma = 14-26 ms for H33N,H39K iso-2 versus gamma = 200-1,100 ms for iso-2. These results with iso-2 cytochrome c strongly support the hypothesis that non-native His-heme ligation results in a kinetic barrier to fast folding of cytochrome c. Assuming that the maximum rate of a conformational search is about 10(11) s-1, the results imply that the direct folding pathway of iso-2 involves passage through on the order of 10(9) or fewer partially folded conformers.