dePolo Gwen E, Kaliakin Danil S, Varganov Sergey A
Department of Chemistry, University of Nevada , Reno, 1664 North Virginia Street, Reno, Nevada 89557-0216, United States.
J Phys Chem A. 2016 Nov 3;120(43):8691-8698. doi: 10.1021/acs.jpca.6b07717. Epub 2016 Oct 25.
Rubredoxin is a small iron-sulfur protein involved in biological electron transfer, which is accomplished by changing the oxidation state of the iron atom in the active site. We investigate the possibility of spin-forbidden transitions between the lowest energy electronic states with different spin multiplicities in the rubredoxin active site models [Fe(SCH)] (n = 2-, 1-, 0) using nonadiabatic transition state theory (NA-TST). The equilibrium structures, minimum energy crossing point structures and Hessians were obtained with density functional theory. The spin-orbit coupling (SOC) was calculated with the complete active space configuration interaction method using the two-electron spin-orbit Breit-Pauli Hamiltonian. We found several crossings between the lowest energy spin states associated with the changes in Fe coordination. However, only triplet/quintet crossings in [Fe(SCH)] and [Fe(SCH)], as well as a quartet/sextet crossing in [Fe(SCH)] are characterized by nonzero first-order SOC responsible for transitions between these spin states. The rates of spin-forbidden transitions in the [Fe(SCH)] complex are 1 and 2 orders of magnitude higher than the rates in the [Fe(SCH)] and [Fe(SCH)] complexes, respectively. These rate differences are related to a large variation of the SOC between the complexes with different charges, which in turn comes from different molecular orbitals involved in the spin-flip transitions. Finally, we demonstrate that the differences between the NA-TST rates and the rates calculated under the assumption of completely spin-allowed transitions could be as large as 4 orders of magnitude. This means that even in qualitative discussions of the reaction mechanisms involving changes in spin states the partially spin-forbidden nature of the transitions between these states must be taken into account.
红素氧还蛋白是一种参与生物电子转移的小铁硫蛋白,它通过改变活性位点中铁原子的氧化态来实现电子转移。我们使用非绝热过渡态理论(NA-TST)研究了红素氧还蛋白活性位点模型[Fe(SCH)](n = 2-、1-、0)中具有不同自旋多重性的最低能量电子态之间自旋禁阻跃迁的可能性。利用密度泛函理论获得了平衡结构、最低能量交叉点结构和海森矩阵。使用双电子自旋轨道布赖特-泡利哈密顿量,通过完全活性空间组态相互作用方法计算了自旋轨道耦合(SOC)。我们发现与铁配位变化相关的最低能量自旋态之间存在几个交叉点。然而,只有[Fe(SCH)]和[Fe(SCH)]中的三重态/五重态交叉点,以及[Fe(SCH)]中的四重态/六重态交叉点,其特征是负责这些自旋态之间跃迁的非零一阶SOC。[Fe(SCH)]配合物中自旋禁阻跃迁的速率分别比[Fe(SCH)]和[Fe(SCH)]配合物中的速率高1个和2个数量级。这些速率差异与不同电荷配合物之间SOC的大幅变化有关,而这又源于自旋翻转跃迁中涉及的不同分子轨道。最后,我们证明NA-TST速率与在完全自旋允许跃迁假设下计算的速率之间的差异可能高达4个数量级。这意味着即使在定性讨论涉及自旋态变化的反应机理时,也必须考虑这些态之间跃迁的部分自旋禁阻性质。