Premium
Influence of the protein and DFT method on the broken‐symmetry and spin states in nitrogenase
Author(s) -
Cao Lili,
Ryde Ulf
Publication year - 2018
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/qua.25627
Subject(s) - basis set , state (computer science) , chemistry , protonation , hybrid functional , ground state , van der waals force , density functional theory , spin states , atomic physics , spin (aerodynamics) , electronic correlation , computational chemistry , crystallography , physics , thermodynamics , molecule , ion , organic chemistry , algorithm , computer science
The enzyme nitrogenase contains a complicated MoFe 7 CS 9 cofactor with 35 possible broken‐symmetry (BS) states. We have studied how the energies of these states depend on the geometry, the surrounding protein, the DFT functional and the basis set, studying the resting state, a one‐electron reduced state and a protonated state. We find that the effect of the basis set is small, up to 11 kJ/mol. Likewise, the effect of the surrounding protein is restricted, up to 10 and 7 kJ/mol for the electrostatic and van der Waals energy terms. Single‐point energies calculated on a single geometry give a good correlation ( R 2 = 0.92‐0.98) to energies calculated after geometry optimization, but some BS states may be disfavored by up to 37 kJ/mol. A change from the pure TPSS functional to the hybrid B3LYP functional may change the relative energies by up to 58 kJ/mol and the correlation between the two results is only 0.57‐0.72. Both functionals agree that BS7 is the most stable BS state and that the ground spin state is the quartet for the resting state and the quintet for the reduced state. With the TPSS functional, the BS6 state is the second most stable state, always at least 21 kJ/mol less stable than the BS7 state. However, with the B3LYP functional, BS10 is the second most stable state and for the protonated state it comes close in energy. Based on these results, we suggest a procedure how to consider the 35 BS states in future investigations of the nitrogenase reaction mechanism.