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 -