Analysis of the electrochemistry of hemes with Ems spanning 800 mV

Authors:
Publisher website
Google Scholar
Abstract
The free energy of heme reduction in different proteins is found to vary over more than 18 kcal/mol. It is a challenge to determine how proteins manage to achieve this enormous range of Ems with a single type of redox cofactor. Proteins containing 141 unique hemes of a‐, b‐, and c‐type, with bis‐His, His‐Met, and aquo‐His ligation were calculated using Multi‐Conformation Continuum Electrostatics (MCCE). The experimental Ems range over 800 mV from −350 mV in cytochrome c3 to 450 mV in cytochrome c peroxidase (vs. SHE). The quantitative analysis of the factors that modulate heme electrochemistry includes the interactions of the heme with its ligands, the solvent, the protein backbone, and sidechains. MCCE calculated Ems are in good agreement with measured values. Using no free parameters the slope of the line comparing calculated and experimental Ems is 0.73 (R2 = 0.90), showing the method accounts for 73% of the observed Em range. Adding a +160 mV correction to the His‐Met c‐type hemes yields a slope of 0.97 (R2 = 0.93). With the correction 65% of the hemes have an absolute error smaller than 60 mV and 92% are within 120 mV. The overview of heme proteins with known structures and Ems shows both the lowest and highest potential hemes are c‐type, whereas the b‐type hemes are found in the middle Em range. In solution, bis‐His ligation lowers the Em by ≈205 mV relative to hemes with His‐Met ligands. The bis‐His, aquo‐His, and His‐Met ligated b‐type hemes all cluster about Ems which are ≈200 mV more positive in protein than in water. In contrast, the low potential bis‐His c‐type hemes are shifted little from in solution, whereas the high potential His‐Met c‐type hemes are raised by ≈300 mV from solution. The analysis shows that no single type of interaction can be identified as the most important in setting heme electrochemistry in proteins. For example, the loss of solvation (reaction field) energy, which raises the Em, has been suggested to be a major factor in tuning in situ Ems. However, the calculated solvation energy vs. experimental Em shows a slope of 0.2 and R2 of 0.5 thus correlates weakly with Ems. All other individual interactions show even less correlation with Em. However the sum of these terms does reproduce the range of observed Ems. Therefore, different proteins use different aspects of their structures to modulate the in situ heme electrochemistry. This study also shows that the calculated Ems are relatively insensitive to different heme partial charges and to the protein dielectric constant used in the simulation. Proteins 2009.

This publication has 131 references indexed in Scilit: