INADEQUACIES OF THE POINT-DIPOLE APPROXIMATION FOR DESCRIBING ELECTRON-NUCLEAR INTERACTIONS IN PARAMAGNETIC PROTEINS - HYBRID DENSITY-FUNCTIONAL CALCULATIONS AND THE ANALYSIS OF NMR RELAXATION OF HIGH-SPIN IRON(III) RUBREDOXIN
Sj. Wilkens et al., INADEQUACIES OF THE POINT-DIPOLE APPROXIMATION FOR DESCRIBING ELECTRON-NUCLEAR INTERACTIONS IN PARAMAGNETIC PROTEINS - HYBRID DENSITY-FUNCTIONAL CALCULATIONS AND THE ANALYSIS OF NMR RELAXATION OF HIGH-SPIN IRON(III) RUBREDOXIN, JOURNAL OF PHYSICAL CHEMISTRY B, 102(42), 1998, pp. 8300-8305
High-level, all-electron, density functional calculations have been us
ed, in conjunction with high-resolution X-ray structural data, to pred
ict, and to compare with experiment, the contribution of unpaired elec
trons to the relaxation times for N-15 nuclei in oxidized Clostridium
pasteurianum rubredosin. Published X-ray structures for the iron(III)
rubredoxin from C. pasteurianum were employed to construct a 104-atom
model for the iron center that included all atoms shown to have strong
electronic interactions with the unpaired iron electrons. The remaind
er of the amide nitrogen resonances in the protein, which show no appa
rent Fermi contact contribution to the chemical shift, are represented
in the model by ghost atoms (atoms with no charge or basis functions)
. This model served as a starting point for quantum mechanical calcula
tions at the B3LYP/6-311G* level, which, in turn, yielded calculated
values for eigenvalues of the spin-differential field gradient tensor,
which finally yielded expectation values for effective distances betw
een nuclei and the delocalized spin-density. We report here that using
effective distances, which are calculated from the spin-differential
field gradient tensor, in the Solomon-Bloembergen equation in place of
distances measured from the crystal structures greatly improves the c
orrelation for a plot of experimental relaxation rates versus r(-6) fo
r N-15 resonances in C. pasteurianum iron(III) rubredoxin. With increa
ses in the speed of computers and algorithms, iterative quantum chemic
al optimization of paramagnetic center geometries based on NMR-derived
distance and angular constraints from paramagnetic interactions shoul
d lead to significant improvements in the determination of the structu
res of paramagnetic centers in proteins by NMR spectroscopy.