We describe electron transfer through the NrfHA nitrite reductase heterodimer using a thermodynamic integration scheme based upon molecular dynamics simulations. From the simulation data, we estimate two of the characteristic energies of electron transfer, the thermodynamic driving forces, ΔG, and the reorganization energies, λ. Using a thermodynamic network analysis, the statistical accuracy of the ΔG values can be enhanced significantly. Although the reaction free energies and activation barriers are hardly affected by protein aggregation, the complete reaction mechanism only emerges from the simulations of the dimer rather than focussing on the individual protein chains: it involves an equienergetic transprotein element of electron storage and conductivity.

No abstract

We simulate electron transfer within a fragment of the NADH:ubiquinone oxidoreductase (respiratory complex I) of the hyperthermophilic bacterium Aquifex aeolicus. We apply molecular dynamics simulations, thermodynamic integration, and a thermodynamic network least squares analysis to compute two key parameters of Marcus' theory of charge transfer, the thermodynamic driving force and the reorganization energy. Intramolecular contributions to the Gibbs free energy differences of electron and hydrogen transfer processes, ΔG, are accessed by calibrating against experimental redox titration data. This approach permits the computation of the interactions between the species NAD+, FMNH2, N1a-, and N3-, and the construction of a free energy surface for the flow of electrons within the fragment. We find NAD+ to be a strong candidate for the regulation of charge transfer.