Free energy of transition for the individual alkaline conformers of yeast iso-1-cytochrome c

Direct protein electrochemistry was used to obtain the thermodynamic parameters of transition from the native (state III) to the alkaline (state IV) conformer for untrimethylated Saccharomyces cerevisiae iso-1-cytochrome c expressed in E. coli and its single and multiple lysine-depleted variants. In these variants, one or more of the lysine residues involved in axial Met substitution (Lys72, Lys73, and Lys79) was mutated to alanine. The aim of this work is to determine the thermodynamic affinity of each of the substituting lysines for the heme iron and evaluate the interplay of enthalpic and entropic factors. The equilibrium constants for the deprotonation reaction of Lys72, 73, and 79 were computed for the minimized MD average structures of the wild-type and mutated proteins, applying a modified Tanford-Kirkwood calculation. Solvent accessibility calculations for the substituting lysines in all variants were also performed. The transition enthalpy and entropy values within the protein series show a compensatory behavior, typical of a process involving extensive solvent reorganization effects. The experimental and theoretical data indicate that Lys72 most readily deprotonates and replaces M80 as the axial heme iron ligand, whereas Lys73 and Lys79 show comparably higher pK(a) values and larger transition free energies. A good correlation is found within the series between the lowest calculated Lys pK(a) value and the corresponding experimental pK(a) value, which can be interpreted as indicative of the deprotonating lysine itself acting as the triggering group for the conformational transition. The triple Lys to Ala mutant, in which no lysine residues are available for heme iron binding, features transition thermodynamics consistent with a hydroxide ion replacing the axial methionine ligand.