Jafari S, Tavares Santos YA, Bergmann J, Irani M, Ryde U. Benchmark Study of Redox Potential Calculations for Iron-Sulfur Clusters in Proteins.
Inorg Chem 2022;
61:5991-6007. [PMID:
35403427 PMCID:
PMC9044450 DOI:
10.1021/acs.inorgchem.1c03422]
[Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Redox potentials
have been calculated for 12 different iron–sulfur
sites of 6 different types with 1–4 iron ions. Structures were
optimized with combined quantum mechanical and molecular mechanical
(QM/MM) methods, and the redox potentials were calculated using the
QM/MM energies, single-point QM methods in a continuum solvent or
by QM/MM thermodynamic cycle perturbations. We show that the best
results are obtained with a large QM system (∼300 atoms, but
a smaller QM system, ∼150 atoms, can be used for the QM/MM
geometry optimization) and a large value of the dielectric constant
(80). For absolute redox potentials, the B3LYP density functional
method gives better results than TPSS, and the results are improved
with a larger basis set. However, for relative redox potentials, the
opposite is true. The results are insensitive to the force field (charges
of the surroundings) used for the QM/MM calculations or whether the
protein and solvent outside the QM system are relaxed or kept fixed
at the crystal structure. With the best approach for relative potentials,
mean absolute and maximum deviations of 0.17 and 0.44 V, respectively,
are obtained after removing a systematic error of −0.55 V.
Such an approach can be used to identify the correct oxidation states
involved in a certain redox reaction.
We
have studied redox potentials of 12 iron−sulfur
sites of 6 types with 1−4 iron ions. Structures were optimized
with combined quantum mechanical and molecular mechanical (QM/MM)
methods, and the redox potentials were calculated with QM/MM, QM calculations
in a continuum solvent or by QM/MM thermodynamic cycle perturbations.
The best results are obtained with the second approach using ∼300
atoms in the QM model and a large dielectric constant.
Collapse