Fold versus sequence effects on the driving force for protein-mediated electron transfer

Electron transport chains composed of electron transfer reactions mainly between proteins provide fast efficient flow of energy in a variety of metabolic pathways. Reduction potentials are essential characteristics of the proteins because they determine the driving forces for the electron transfers. As both polar and charged groups from the backbone and side chains define the electrostatic environment, both the fold and the sequence will contribute. However, although the role of a specific sequence may be determined by experimental mutagenesis studies of reduction potentials, understanding the role of the fold by experiment is much more difficult. Here, continuum electrostatics and density functional theory calculations are used to analyze reduction potentials in [4Fe-4S] proteins. A key feature is that multiple homologous proteins in three different folds are compared: six high potential iron-sulfur proteins, four bacterial ferredoxins, and four nitrogenase iron proteins. Calculated absolute reduction potentials are shown to be in quantitative agreement with electrochemical reduction potentials. Calculations further demonstrate that the contribution of the backbone is larger than that of the side chains and is consistent for homologous proteins but differs for nonhomologous proteins, indicating that the fold is the major protein factor determining the reduction potential, whereas the specific amino acid sequence tunes the reduction potential for a given fold. Moreover, the fold contribution is determined mainly by the proximity of the redox site to the protein surface and the orientation of the dipoles of backbone near the redox site.