Mechanism of water oxidation catalyzed by vitamin B12:Redox non-innocent nature of corrin ligand and crucial role of phosphate

Ying-Ying Li,Rong-Zhen Liao

Key Laboratory of Material Chemistry for Energy Conversion and Storage,Ministry of Education,Hubei Key Laboratory of Bioinorganic Chemistry and Materia Medica,Hubei Key Laboratory of Materials Chemistry and Service Failure,School of Chemistry and Chemical Engineering,Huazhong University of Science and Technology,Wuhan 430074,China

ABSTRACT Vitamin B12(macrocyclic cobalamin)has been recently reported to be capable of electrochemically catalyzing water oxidation in a neutral phosphate buffer solution.In this work,density functional calculations were employed to elucidate the water oxidation mechanism catalyzed by vitamin B12.The calculations showed that the catalytic cycle starts from the L?-CoII-OH2 complex 1.A proton-coupled electron transfer process then leads to the formation of a L?-CoIII-OH complex 2,followed by another protoncoupled electron transfer event to afford a corrin ligand radical cation intermediate 3(L?-CoIII-O?).The redox non-innocent nature of the corrin ligand plays an essential role in the oxidation process.3 is capable of triggering the O-O bond formation via a water nucleophilic attack mechanism,in which a hydrophosphate dianion functions as a base to accept a proton from the water nucleophile.A dioxygen molecule is released after the oxidation of the CoIII-OOH intermediate.The rate-determining step was calculated to be the O-O bond formation with a total barrier of 16.5 kcal/mol.While the use of water molecules as the proton acceptor was found to be less feasible for the O-O bond formation,with a barrier of 31.2 kcal/mol,further highlighting the crucial of phosphate in water oxidation.

Keywords:Water oxidation DFT calculations Mechanism Vitamin B12 Cobalt

Nowadays,the energy crisis is one of the urgent issues faced by us.Sunlight-driven water splitting in artificial photosystem provides a promising strategy for sustainable and clean hydrogen fuel production[1–5].Water oxidation is a half and crucial step in water splitting,consisting of releasing four protons and four electrons and dissociating the molecular dioxygen.During the last few decades,extraordinary efforts have been devoted to this challenging reaction[6–12].Ru-based water oxidation catalysts(WOCs)are most studied among homogeneous water oxidation communities[13–25].However,considering the cost,low toxicity and availability,first-row transition metal-based complexes are believed to be better candidates for WOCs.Even though significant progress has been made in this field[24–32],the ultimate goal of application in artificial photosystem still needs more efficient,stable and robust WOCs.To achieve this goal,a thorough mechanistic understanding by combing experiments and computation is essential.Quantum chemical calculations are playing more and more critical roles in this area[33–35].

The commercially available vitamin B12has been reported to be an environmentally friendly catalyst for chemical reactions,such as dehalogenation,CO2electrochemical reduction,hydrogenation of double bonds[36–39].Because of the nontoxic and rich redox chemical property,vitamin B12has attracted more and more attention in the catalytic community.Recently,Verpoort and coworkers reported that vitamin B12could function as a homogeneous WOC electrochemically in pH 7.0 phosphate buffer(NaPi)solution with an overpotential of 0.58 V[40].UV-vis,ESI-MS and FTIR experiments suggested that the catalyst remains its molecular structure up to 11 h[40].The onset potential was observed at 1.2 Vvs.Ag/AgCl,together with a more significant current density of water oxidation at around 1.5 Vvs.Ag/AgCl.The linear relationship between the catalytic current and the catalyst concentration was also established,indicating a unimolecular mechanism.Furthermore,the increase of the concentration of NaPi leads to the linear growth of the current density in the working conditions,which suggested that phosphate ions may play an essential role in water oxidation[40].

Fig.1.Optimized structures of 1.Distances are shown in ?A in red.Mulliken spin densities on selected atoms are shown in black.The multiplicity is shown in superscript.Hydrogen atoms of methyl are not shown for clarity.

In the present work,DFT calculations were used to investigate the mechanism of water oxidation catalyzed by vitamin B12and address the role of phosphate ion in this chemical process.The methods used here are to a large extent similar to those in our previous water oxidation studies[33,41,42],and the details are shown in Supporting information.To construct the whole energy diagram,a reference potential of 1.72 Vvs.SHE(corresponding to 1.5 Vvs.Ag/AgCl)was used.

We first investigated the penta-coordinated complex 1(CoIIIOH2,Fig.1).The corrin ring has four nitrogen atoms coordinating with the cobalt center in the equatorial positions.An imidazole nitrogen atom is interacting with the metal center in the axial direction.The axial water molecule does not coordinate with the metal,instead interacts with the complexviatwo hydrogen bonds with the amide group of the ligand(Fig.1).These two hydrogen bond distances are 2.03 and 1.90 ?A,respectively.1 prefers to be a triplet,while the singlet lies at +5.5 kcal/mol relative to the triplet.As shown in Fig.1 and Fig.S2(Supporting information),the spin density in 1 is partially delocalized on the ligand(total spin density of 0.91).Even though the formal oxidation state of the metal is CoIII,the electronic structure of 1(L?-CoII)would be better described as a low-spin CoIIcenter(SCo=1/2)interacting with a corrin ligand radical cation(SL=1/2).The coordination of water to cobalt to form a hexa-coordinated complex 1’(Fig.S3 in Supporting information)is 3.6 kcal/mol higher in energy.

The oxidation of 1 to form 2(formally CoIV-OH,Fig.S4 in Supporting information)is coupled with the release of one proton from the water ligand,since the pKaof 2pt(the protonated form of 2,Fig.S2)was calculated to be -5.9.2 is a doublet with the quartet state lying at +16.3 kcal/mol.The redox potential was calculated to be 0.99 V.Taking 1.72 V as a reference,this step is exergonic by 16.8 kcal/mol(Fig.2).The spin delocalization on the ligand(spin density of 1.25)in 2 indicates that its electronic structure(L?-CoIIIOH)can be described as a low-spin CoIII(SCo=0)with a ligand radical cation(SL=1/2)(Fig.S2).The Co-O1 bond is 1.88 ?A in 2.

Next,a proton-coupled electron transfer(PCET)process of 2 leads to the generation of 3(formally CoV=O,Fig.S3)with a redox potential of 1.69 V.The calculated value is fully consistent with the experimental observation of 1.72 V[40].The pKaof 3pt(the protonated form of 3)was calculated to be 1.8(Fig.S2,3pt).The Mulliken spin population analysis(1.0 on the corrin and 1.0 on O1)suggested that the electronic structure of 3(L?-CoIII-O?)can be interpreted as a low-spin CoIII(SCo=0)interacting with an oxyl radical(SO=1/2)as well as a ligand radical(SL=1/2)ferromagnetically,which is reminiscent of previously reported,cobalt hangmancorroles WOC[43,44].3 is a triplet,while the broken-symmetry singlet is slightly higher in energy,locating at +0.8 kcal/mol above.The Co-O1 distance in 3 is 1.81 ?A which is somewhat shorter than that in 2.

The formally CoV=O species has been widely proposed to trigger the O-O bond formation for cobalt-based WOCs[34,44–47].Here,two possible mechanistic scenarios for the O-O bond formation were explored.First,a water molecule performs the nucleophilic attack on the oxyl moiety of 3 to form the O-O bond,during which one proton is transferred from the attacking water molecule to phosphate ion,either HPO42?or H2PO4?in working conditions(vide infra).The use of phosphate as a base has been suggested by the experiment and other metal-based water oxidation processes[34,46–50].Second,water molecules in the aqueous solution could also function as a proton acceptor.

The transition state(TS1)for the water nucleophilic attack(WNA)pathway with the HPO42?behaving as a proton acceptor is shown in Fig.2.TS1 prefers to be a broken-symmetry singlet,and a spin crossing from the triplet to the broken-symmetry singlet is required to prompt the O-O bond formation.In TS1,the atomic spin density on the oxyl radical is 0.49α,while the ligand has the opposite spin density of 0.53β,which provides the most efficient electronic state for the O-O bond formation.(Fig.2)During the OO bond formation,the attacking water molecule provides a pair of electrons to form the O-O bond[33,51].Theα-electron from the nucleophilic water molecule transfers to the ligand,while the proton is transferred to the phosphate.This can be considered an uncommon PCET process in which proton and electron are delivered to different acceptors.Theβ-electron left could easily combine with the oxyl radial to afford the O-O bond,which is similar to our previous work on the iron WOC[41].TS1 was confirmed to have only one imaginary frequency of 1237.2icm?1,which corresponds to the O-O bond formation and proton transfer from water to phosphate.At TS1,the nascent O1-O2 bond distance is 1.96 ?A.TS1 has a barrier of +16.5 kcal/mol above the 3.The triplet state lies at +9.5 kcal/mol relative to the broken-symmetry singlet.Downhill from TS1,the singlet intermediate Int1 with a hydroperoxide coordinated to CoIII(Fig.S4)is generated.In Int1,the O1-O2 bond length is 1.45 ?A,and the Co-O1 distance increases from 1.81 ?A in 3 to 1.83 ?A in Int1.Int1 lies at +4.6 kcal/mol relative to 3.When H2PO4?was used as the proton acceptor,the barrier was calculated to be 17.9 kcal/mol(TS1’,Fig.2),which is 1.4 kcal/mol higher than the HPO42?assisted pathway.The NBO calculations were carried out for the transition states,and the NBO charge for HPO42?moiety in TS1 is ?1.75,while in TS1’the charge is calculated to be ?0.84 for H2PO4?moiety(Table S2 in Supporting information).This explains the favorability for the former pathway for WNA.In TS1’,the O1-O2 distance is 0.03 ?A shorter than that in TS1.In the working condition(pH 7.0,0.1 mol/L NaPi),the concertation of H2PO4?and HPO42?are 0.062 mol/L and 0.038 mol/L,respectively.In any case,the HPO42?assisted pathway is preferred.

Alternatively,the transition state(TS1’’)with water molecules featuring as the proton acceptor is shown in Fig.2.In this configuration,three more water molecules were added to solvate the generated hydronium explicitly[41,52,53].TS1’’was confirmed as a real transition state with only one imaginary frequency of 610.7icm?1.TS1’’prefers to be a broken-symmetry singlet,in which the oxyl radical and the corrin ligand have opposite spin populations(0.44 on the oxyl,while -0.44 on the ligand).The crucial O1-O2 distance is 1.89 ?A at TS1’’.The Co- O1 length increases to 1.86 ?A from 1.81 ?A in 3.The O2-H1 and O4-H1 distances are 1.02 ?A and 1.58 ?A,respectively.The barrier for this pathway was calculated to be 31.2 kcal/mol,which is much higher than the phosphate pathway(Fig.2).The triplet TS1’’was calculated to be 8.5 kcal/mol higher than the broken singlet.This WNA process also generated the hydro-peroxide intermediate Int1,which joins the previous pathway.

Fig.2.Gibbs energy profile for water oxidation at the B3LYP-D3 level.Optimized TS1(top right),TS1’and TS1’’(bottom left)in broken singlet are also shown.The imaginary frequencies of TS1,TS1’and TS1’’were 1237.2 i cm?1,1164.0 i cm?1 and 610.7 i cm?1,respectively.Some of the unimportant atoms are not shown for clarity.Distances are shown in ?A in red.Mulliken spin densities on selected atoms are shown in black.

From Int1,a PECT process leads to the formation of Int2(Fig.S6 in Supporting information)with a redox potential of only 0.13 V.This demonstrated that this oxidation process is quite facile,with an exergonicity of 36.7 kcal/mol(Fig.2).Int2 is a quartet state,while the doublet state is 3.0 kcal/mol higher.In Int2,the triplet dioxygen is already generated since the spin densities on O1 and O2 are 0.95 and 0.99,respectively.Meanwhile,the O2dissociates from the cobalt center with the calculated Co-O1 distance of 3.26 ?A.Furthermore,during the oxidation process,one electron is transferred from the hydro-peroxide moiety to the Co center with the formation of a low spin CoII,whose spin density is 0.99(Fig.S4).The release of O2from Int2 is followed by the binding of a water molecule to form the doublet intermediate Int3(Fig.S4),with the quartet 12.1 kcal/mol higher in energy.This ligand exchange process is exergonic by 4.0 kcal/mol(Fig.2).Finally,the last oneelectron oxidation process from Int3 regenerates 1,with a calculated redox potential of 0.32 V.It should be mentioned that the ligand is involved in this oxidation process,as the spin is partially delocalized on the ligand in 1(L?-CoII).

From Fig.2,it can be seen that the O-O bond formation is the rate-determining step with a calculated total barrier of 16.5 kcal/mol at the B3LYP-D3 level.It should be mentioned that the use of different applied potentials may change the total barrier.When using a reference potential lower than 1.69 V(oxidation potential of 2 to 3),the total barrier will increase somewhat as the formation of 3 becomes endergonic.For example,when a reference potential of 1.4 V was used to set up the energy diagram,the total barrier increases to 23.2 kcal/mol(Table S1 in Supporting information).If the reference potential is higher than 1.69 V,the total barrier does not change as the formation of 3 is exergonic.3 is the resting state for the O-O bond formation and its electronic structure is L?-CoIII-O?.Therefore,the catalytic activity of water oxidation is assisted by the redox non-innocent corrin ligand,which avoids the access of the high-valent metal center.The use of non-innocent ligands has been disclosed in other metal-based WOCs,a beneficial strategy in the catalytic water oxidation process[11,41,48,54].Furthermore,single-point calculations were performed using the MN15 and MN15L functionals[55,56].The calculated redox potentials and total barriers are summarized in Table 1.For the first oxidation process,the MN15 gives a relatively higher redox potential value than B3LYP-D3 and MN15L.Importantly,all functionals give similar results for the O-O bond formation.When HPO42?was used to accept the proton,the barrier is 14.7 kcal/mol lower than water receiving the proton at the B3LYP-D3 level.Thebarrier differences at the MN15 and MN15L levels are 16.9 kcal/mol and 18.9 kcal/mol,respectively.This demonstrated that HPO42?is a better base for the O-O bond formation in the water oxidation process.

Table 1 Comparison of redox potentials and total barriers at different levels.a

To conclude,we provide insights into the mechanism of water oxidation catalyzed by commercially available cobalamin(vitamin B12)through density functional calculations.The suggested catalytic cycle is summarized in Scheme 1,and the corresponding energy diagram is shown in Fig.2.The starting point of the catalytic cycle is the formal CoIII-OH2complex 1.Two sequential PCET processes result in the formation of a formal CoV=O intermediate 3,which can be described as L?-CoIII-O?.The redox non-innocent character of the corrin ligand in the water oxidation was illuminated.3 initiates the O-O bond formationviaWNA with the production of a hydro-peroxide intermediate.Importantly,HPO42?in the phosphate buffer was suggested to be the proton acceptor during the O-O bond formation.This step was suggested to be ratedetermining with a total barrier of 16.5 kcal/mol.The alternative pathway with water accepting the proton during the O-O bond formation was calculated to be kinetically unfavorable since the barrier is 14.7 kcal/mol higher.The triplet dioxygen molecule dissociated from the cobalt center spontaneously after one PCET process from the CoIIIhydro-peroxide intermediate.Finally,the second water molecule enters after the release of O2,followed by another one-electron oxidation process to regenerate 1.

Scheme 1.Suggested water oxidation mechanism.Inset:The electron flow during the O-O bond formation is described.

Declaration of competing interest

All authors declare that no conflict of interest exists.

Acknowledgments

This work was supported by the National Natural Science Foundation of China(No.21873031)and the National Key R&D Program of China(No.2018YFA0903500).

Supplementary materials

Supplementary material associated with this article can be found,in the online version,at doi:10.1016/j.cclet.2021.06.028.