Ab initio molecular dynamics of heme in cytochrome c

Out of plane distortions of the heme b of Escherichia coli succinate dehydrogenase

The role of the heme b in Escherichia coli succinate dehydrogenase is highly ambiguous and its role in catalysis is questionable. To examine whether heme reduction is an essential step of the catalytic mechanism, we generated a series of site-directed mutations around the heme binding pocket, creating a library of variants with a stepwise decrease in the midpoint potential of the heme from the wild-type value of +20 mV down to −80 mV. This difference in midpoint potential is enough to alter the reactivity of the heme towards succinate and thus its redox state under turnover conditions. Our results show both the steady state succinate oxidase and fumarate reductase catalytic activity of the enzyme are not a function of the redox potential of the heme. As well, lower heme potential did not cause an increase in the rate of superoxide production both in vitro and in vivo. The electron paramagnetic resonance (EPR) spectrum of the heme in the wild-type enzyme is a combination of two distinct signals. We link EPR spectra to structure, showing that one of the signals likely arises from an out-of-plane distortion of the heme, a saddled conformation, while the second signal originates from a more planar orientation of the porphyrin ring.

NMR and DFT investigation of heme ruffling: functional implications for cytochrome c

Out-of-plane (OOP) deformations of the heme cofactor are found in numerous heme-containing proteins and the type of deformation tends to be conserved within functionally related classes of heme proteins. We demonstrate correlations between the heme ruffling OOP deformation and the (13)C and (1)H nuclear magnetic resonance (NMR) hyperfine shifts of heme aided by density functional theory (DFT) calculations. The degree of ruffling in the heme cofactor of Hydrogenobacter thermophilus cytochrome c(552) has been modified by a single amino acid mutation in the second coordination sphere of the cofactor. The (13)C and (1)H resonances of the cofactor have been assigned using one- and two-dimensional NMR spectroscopy aided by selective (13)C-enrichment of the heme. DFT has been used to predict the NMR hyperfine shifts and electron paramagnetic resonance (EPR) g-tensor at several points along the ruffling deformation coordinate. The DFT-predicted NMR and EPR parameters agree with the experimental observations, confirming that an accurate theoretical model of the electronic structure and its response to ruffling has been established. As the degree of ruffling increases, the heme methyl (1)H resonances move upfield while the heme methyl and meso (13)C resonances move downfield. These changes are a consequence of altered overlap of the Fe 3d and porphyrin pi orbitals, which destabilizes all three occupied Fe 3d-based molecular orbitals and decreases the positive and negative spin density on the beta-pyrrole and meso carbons, respectively. Consequently, the heme ruffling deformation decreases the electronic coupling of the cofactor with external redox partners and lowers the reduction potential of heme.

Free energies for biological electron transfer from QM/MM calculation: method, application and critical assessment

Computer simulations of biological electron transfer reactions are reviewed with a focus on the calculation of reaction free energy (driving force) and reorganization free energy. Then a mixed quantum mechanical/molecular mechanical (QM/MM) approach is described which is designed for computation of these quantities for pure electron transfer reactions with large donor-acceptor separation distances. The method is applied to intra-protein electron transfer in Ru(bpy)(2)(im)His33 cytochrome c and the results compared to experimental data. Several modeling aspects which are important for successful calculation of free energies with QM/MM are discussed in detail.