Calculating standard reduction potentials of [4Fe-4S] proteins.

The oxidation-reduction potentials of electron transfer proteins determine the driving forces for their electron transfer reactions. Although the type of redox site determines the intrinsic energy required to add or remove an electron, the electrostatic interaction energy between the redox site and its surrounding environment can greatly shift the redox potentials. Here, a method for calculating the reduction potential versus the standard hydrogen electrode, E°, of a metalloprotein using a combination of density functional theory and continuum electrostatics is presented. This work focuses on the methodology for the continuum electrostatics calculations, including various factors that may affect the accuracy. The calculations are demonstrated using crystal structures of six homologous HiPIPs, which give E° that are in excellent agreement with experimental results.

Characterizing the effects of the protein environment on the reduction potentials of metalloproteins.

The reduction potentials of electron transfer proteins are critically determined by the degree of burial of the redox site within the protein and the degree of permanent polarization of the polypeptide around the redox site. Although continuum electrostatics calculations of protein structures can predict the net effect of these factors, quantifying each individual contribution is a difficult task. Here, the burial of the redox site is characterized by a dielectric radius R(p) (a Born-type radius for the protein), the polarization of the polypeptide is characterized by an electret potential ϕ(p) (the average electrostatic potential at the metal atoms), and an electret-dielectric spheres (EDS) model of the entire protein is then defined in terms of R(p) and ϕ(p). The EDS model shows that for a protein with a redox site of charge Q, the dielectric response free energy is a function of Q(2), while the electret energy is a function of Q. In addition, R(p) and ϕ(p) are shown to be characteristics of the fold of a protein and are predictive of the most likely redox couple for redox sites that undergo different redox couples.

Identifying sequence determinants of reduction potentials of metalloproteins.

The reduction potential of an electron transfer protein is one of its most important functional characteristics. Although the type of redox site and the protein fold are the major determinants of the reduction potential of a redox-active protein, its amino acid sequence may tune the reduction potential as well. Thus, homologous proteins can often be divided into different classes, with each class characterized by a biological function and a reduction potential. Site-specific mutagenesis of the sequence determinants of the differences in the reduction potential between classes should change the reduction potential of a protein in one class to that of the other class. Here, a procedure is presented that combines energetic and bioinformatic analysis of homologous proteins to identify sequence determinants that are also good candidates for site-specific mutations, using the [4Fe-4S] ferredoxins and the [4Fe-4S] high-potential iron-sulfur proteins as examples. This procedure is designed to guide site-specific mutations or more computationally expensive studies, such as molecular dynamics simulations. To make the procedure more accessible to the general scientific community, it is being implemented into CHARMMing, a Web-based portal, with a library of density functional theory results for the redox site that are used in the setting up of Poisson-Boltzmann continuum electrostatics calculations for the protein energetics.

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.