Theoretical Study on the Mechanism of Photosynthetic Oxygen Evolution by ABEEM/MM/MD and BS-DFT
Guo Yu,Yao Yuan,Li Hui,He Lanlan,Zhu Zunwei,Yang Zhongzhi,Gong Lidong,Liu Cui,Zhao Dongxia
DOI: https://doi.org/10.6023/a17040151
2017-01-01
Acta Chimica Sinica
Abstract:Charge parameters of atom-bond electronegativity equalization method (ABEEMsp) for the oxygen-evolving complex (OEC) in the S-2 state were established, which were applied to molecular dynamic (MD) simulation based on ABEEM/MM polarizable force field in order to study the mechanism of photosynthetic oxygen evolution, in combination with broken-symmetry density functional theory (BS-DFT). Charge fitting results at HF/STO-3G level (the reason why the basis set is adopted is shown in the cited literatures) show good linear correlation, proving the rationality and efficiency of the ABEEMsp model in calculating charge distributions. It can be seen from MD simulations that bidirectional isomerizations of the Mn4CaO5 cluster accompany by the transfer of Ca-bound water molecule W3 to pentacoordinate Mn1(III)/Mn4(III). For the g = 2 to g = 4.1 form, W3 leaves for Mn4(III), while W3 moves to Mn1(III) for the inverse course. Both processes involve motions to the original positions of W4 to W3, W588 to W4, and local rearrangements of the water environment, which may indicate the importance of hydrogen bond network to biocatalysis. The observation of W3 coordination to the vacant site of Mn(III) as the sixth ligand proximal to O-5 may imply W3 could be the fast-exchanging substrate water (W-f) in the S-2 state, which makes O-2 with O-5 in the S-4 state. Based on the inference, we investigate O-O bond formation in all the possible spin states for the two isomeric structures under the framework of oxo-oxyl radical coupling mechanism. It is demonstrated from BS-DFT calculations that O-2 formation activity is significantly advantageous for the open-cubane structure than the closed-cubane form, i.e. the differences of barriers and driving forces are beyond 20 kcal/mol and 25 kcal/mol, respectively. For the open-cubane structure, the antiferromagnetic coupling of Mn1(IV)-O-center dot stabilizes the reactants, and the spin-parallel feature between O-center dot radical and Mn4(IV) lowers the barriers, and ferromagnetic coupling of Mn1-Mn3 ensures the release of triplet O-2. For the closed-cubane structure, the antiferromagnetic coupled Mn4(IV)-O-center dot changes to more stabilized Mn4(III)-oxo, losing the radical character of the ligand oxygen, which greatly increases the barriers. The ferromagnetic coupled Mn4(IV)-O-center dot does not belong to the relatively stable electronic configuration of reactants, and cannot be formed by spin frustration from the stable spin states, for their large energy differences above 20 kcal/mol. Thus it cannot be accepted as the accesses of effective reaction channels. Our work expounds the exclusive role of the open-cubane OEC in oxo-oxyl coupling mechanism of O-2 creation. It also reveals that the oxygen evolution reactivity is extremely dependent on the spin coupling ways of manganeses and oxygen radical, meanwhile, the structural flexibility of the OEC is essential to the S-state cycle and photosynthetic water splitting.