TY - JOUR

T1 - Calculation of electric field gradients in Cd(II) model complexes of the CueR protein metal site

AU - O'Shea, Catriona A.

AU - Fromsejer, Rasmus

AU - Sauer, Stephan P. A.

AU - Mikkelsen, Kurt Valentin

AU - Hemmingsen, Lars Bo Stegeager

PY - 2023

Y1 - 2023

N2 - With this work we first test various DFT functionals against CCSD(T) for calculation of EFGs at the position of Cd(II) in a very small model system, Cd(SCH3)2. Moreover, the available basis sets in ADF are tested in terms of basis set convergence, and the effect of including relativistic effects using the scalar relativistic and spin orbit ZORA hamiltonians is explored. The results indicate that an error of up to around 10 % on the calculated EFG may be expected using spin-orbit ZORA and the BHandH functional with a locally dense basis set. Next, this method was applied to model systems of the CueR protein, aiming to interpret 111Ag-PAC spectroscopic data. Note that 111Ag decays to 111Cd on which the PAC data are recorded. Surprisingly, model systems truncated - as is often done - at the first C-C bond from the central Cd(II) are inadequate in size, and larger model systems must be employed to achieve reliable EFG calculations. The calculated EFGs agree well with experimental PAC data, and indicate that shortly after the nuclear decay the structure relaxes from linear two-coordinate AgS2 in the native protein, to a structure (or structures) where Cd(II) recruits additional ligands such as backbone carbonyl oxygens to achieve higher coordination number(s).

AB - With this work we first test various DFT functionals against CCSD(T) for calculation of EFGs at the position of Cd(II) in a very small model system, Cd(SCH3)2. Moreover, the available basis sets in ADF are tested in terms of basis set convergence, and the effect of including relativistic effects using the scalar relativistic and spin orbit ZORA hamiltonians is explored. The results indicate that an error of up to around 10 % on the calculated EFG may be expected using spin-orbit ZORA and the BHandH functional with a locally dense basis set. Next, this method was applied to model systems of the CueR protein, aiming to interpret 111Ag-PAC spectroscopic data. Note that 111Ag decays to 111Cd on which the PAC data are recorded. Surprisingly, model systems truncated - as is often done - at the first C-C bond from the central Cd(II) are inadequate in size, and larger model systems must be employed to achieve reliable EFG calculations. The calculated EFGs agree well with experimental PAC data, and indicate that shortly after the nuclear decay the structure relaxes from linear two-coordinate AgS2 in the native protein, to a structure (or structures) where Cd(II) recruits additional ligands such as backbone carbonyl oxygens to achieve higher coordination number(s).

U2 - 10.1039/D2CP05574K

DO - 10.1039/D2CP05574K

M3 - Journal article

VL - 25

SP - 12277

EP - 12283

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

ER -