Coordination of Cu(I1) to Lipophilic Bis-hydroxamate Binders As

Orientation-selective electron spin echo envelope modulation (ESEEM) experiments were carried out at 8.45 and. 9.15 GHz. The spectra obtained were ...
0 downloads 0 Views 1MB Size
J. Am. Chem. SOC. 1995,117, 12771-12778

12771

Coordination of Cu(I1) to Lipophilic Bis-hydroxamate Binders As Studied by One- and Two-Dimensional Electron Spin Echo Spectroscopy V. Kofman,t J. J. ShaneJ S. A. Dikanov,”s M. K. Bowman,#J. Libman: A. Shanzer,l and D. Goldfarb*gt Contribution from the Departments of Chemical Physics and Organic Chemistry, Weizmann Institute of Science, 76100 Rehovot, Israel, Institute of Chemical Kinetics and Combustion, Novosibirsk, 630090, Russian Federation, and Macromolecular Structure and Dynamics, Environmental Molecular Science Laboratory, Pacific Northwest Laboratories, Richland, Washington 99352 Received June 19, 1995@

Abstract: The binding of Cu(II) to a lipophilic bis-hydroxamate binder, RL252, and its parent RL239, was investigated by pulsed EPR. The binders have two arms, each terminated with a hydroxamate group which serves as a donor. The major difference between RL252 and RL239 is the absence of the leucine amino acid bridge in RL239. Orientation-selective electron spin echo envelope modulation (ESEEM) experiments were carried out at 8.45 and 9.15 GHz. The spectra obtained were exceptionally well resolved and indicated that the cancellation condition, which requires that the hyperfine coupling is approximately twice the nuclear Larmor frequency, is met at 8-9 GHz. The spectra of both Cu-RL252 and Cu-RL239 showed the nuclear quadrupole resonance (NQR) frequencies, yo, v-, and v+, of the I4N of the hydroxamate groups at 1.6-1.7,2.15-2.25, and 3.8-3.9 MHz. A peak corresponding to the double quantum transition in the “noncanceled manifold”, VDQ, was observed at 5-5.2 MHz, and at some magnetic fields a peak corresponding to a single quantum transition, VSQ, and its combination harmonic appeared as well. The assignment of all the ESEEM frequencies was achieved by the application of the two-dimensional hyperfine sublevel correlation (HYSCORE) experiment. Following the assignment, simulations of the orientation-selective ESEEM spectra were performed, yielding the magnitude and orientation of the quadrupole and hyperfine tensors of the coupled nitrogens. While only one best fit set of quadrupolar parameters was found, two such sets were obtained for the hyperfine interaction. Analysis of the orientation of the quadrupole tensor with respect to the g-tensor showed that in both Cu-RL252 and Cu-RL239 the binding site is close to coplanar. Very subtle differences in the spin Hamiltonian parameters were observed between Cu-RL252 and Cu-RL239. The I4N quadrupolar parameters and the anisotropic hyperfine component were slightly larger in Cu-RL239. The relatively small I4N hyperfine coupling is attributed to a node in the molecular orbital, occupied by the unpaired electron, at the nitrogen.

Introduction Lipophilic carriers of metal ions are of great interest for biological applications in view of their capability to penetrate biological membranes. Thanks to these properties, such compounds may function either as scavengers of intracellular metal ions, or as carriers for the delivery of metal ions needed in the cell. One example is the use of lipophilic ion carriers for the removal of iron from red blood Such carriers have been shown to provide potential antimalarial agents by inhibiting intracellular parasite growth. Another example is the use of lipophilic bis-hydroxamate carriers for the delivery of vanadyl ions to fat c e l k 4 Vanadyl ions have been shown to simulate most of the properties of insulin, and their effective intracellular delivery has been demonstrated to significantly augment their insulinomimetic properties. Recently, a lipophilic bis-hydrox-

’ Department of Chemical Physics.

Institute of Chemical Kinetics and Combustion. Macromolecular Structure and Dynamics. Department of Organic Chemistry. Abstract published in Aduance ACS Abstracts, December 1, 1995. (1) Shanzer, A.; Libman, J.; Lytton, S . ; Glickstein, 2.;Cabantchik, 2.I. Proc. Natl. Acad. Sci. 1991, 88, 6585. (2)Lytton, S . D.; Mester, B.; Dayan, I.; Glickstein, H.; Libman, J.; Shanzer, A,; Cabantchik, 2. I. Blood 1993, 81, 214. (3)Lytton, S . D.; Mester, B.; Landau, I.; Libman, J.; Shanzer, A,; Cabantchik, Z. I. Am. J . Hemntol. 1993, 43, 217. (4) Shechter, Y . ; Shisheva, A,; Lazar, R.; Libman, J.; Shanzer, A. Biochemistry 1992, 31, 2063.

amate for Cu(I1) (termed RL252, Figure 1) was designed and synthesized. The Cu(II) replaces the two protons, thus forming a neutral complex. It is a highly selective Cu(I1)-extractant, and it selectively transports Cu(1I) across artificial membranes in the presence of five competing metal ions, namely Co(II), Ni(II), Zn(II), Cd(II), and Pb(II).5 RL252 proved to bind Cu(11)with an exceptionally high binding constant, log 8 2 = 13.99’ and thereby to significantly exceed the binding efficiency of a single acetohydroxamate m ~ i e t y .The ~ unique Cu(I1)-binding characteristics of RL252, and the lack of any crystallographic data on Cu(I1)-bis-hydroxamates (except for some powder data8), prompted us to examine the Cu(I1)-bis-hydroxamate configuration of Cu-RL252 as a prototype of Cu(1I)-bishydroxamate. Encouraged by the recent application of electron spin echo envelope modulation (ESEEM) spectroscopy for the characterization of the vanadyl-bis-hydroxamate~,~orientationselective three-pulse ESEEM along with two-dimensional (2D)

5

@

0002-7863/95/1517-12771$09.00/0

(5) Bromberg, L.; Levin, G.; Libman, L.; Shanzer, A. J . Membrane Sci. 1992, 69, 143.

(6) Farkas, E. Unpublished results. (7) Schwartenbach, G.; Schwartenbach, K. Helu. Chim. Acta 1963, 46, 1390. (8) Hutchinson, B.; Bostick, Y . ;Crowder, Y . ;Eversdyk, D.; Hurley, D.; Jett, D.; Olbricht, S . ; Sample, S . ; Thompson, L. Inorg. Chim. Acta 1983, 74, 29. (9) Kofman, V.; Dikanov, S . A,; Haran, A.; Libman, J.; Shanzer, A,; Goldfarb, D. J . Am. Chem. SOC. 1995,117, 383.

0 1995 American Chemical Society

Kofman et al.

12772 J. Am. Chem. SOC., Vol. 117, No. 51, 1995

of this experiment over the standard ESEEM approaches. The structure of the Cu(II) binding site in Cu-RL252 and CuRL239 was found to be close to square coplanar, and only small differences were observed in the spectroscopic properties of Cu-RL252 and Cu-RL239.

Experimental Section

RL239

RL252 ( R=iBu, m=2 ) RL261 ( R=Me, m=l )

Figure 1. Schematic structures of the bis-hydroxamate binders.

hyperfine sublevel correlation (HYSCORE) expenmentsl0 were applied to study the binding geometry of Cu-RL252 and CuRL239 (Figure 1). The ESEEM method was particularly effective for investigating the structure of the binding site of the VO(I1) complexes due to the so-called “exact cancellation condition”‘ which for the hydroxamate nitrogens is met at X-band EPR frequency (-9 G H Z ) . ~The cancellation condition requires that the hyperfine coupling is approximately twice the Larmor frequency of the nitrogen ( 2 v ~GZ aiso),leading to an effective field close to zero in one of the electron spin manifolds due to the cancellation of the nuclear Zeeman field by the hyperfine field. The modulation frequencies corresponding to nuclear transitions within this manifold are, therefore, the nuclear quadrupole resonance (NQR) frequencies, yo, v - , Y + . ~ ’ , ’ ~The Fourier-transformed (FT) ESEEM spectrum of an orientationally disordered system typically consists of the three peaks at the NQR frequencies and a fourth line corresponding to the Am1 = 2 transition (VDQ) in the other electron spin manifold. The other two signals, corresponding to single quantum transitions (YSQ,Am1 = 1) within this manifold, are usually broadened beyond detection.’ When the EPR spectrum is governed by a large anisotropic interaction it is possible to perform orientation-selectiveESEEM experiment^'^ which provide the anisotropic component of the hyperfine interactionI4 and the orientation of the hyperfine and quadrupole tensors. 14,’ As in the VO(II)-bis-hydroxamate c o m p l e ~ e sthe , ~ hyperfine coupling of the hydroxamate nitrogens in the Cu(II) complexes fulfills the cancellation condition at -9 GHz, hence providing the means to determine their hyperfine and quadrupole interactions and the geometry of the binding site. The orientationselective FT-ESEEM spectra of the Cu(I1) complexes are exceptionally well resolved, and in some spectra VSQ can be observed. Complete assignment of the spectra was achieved by the application of the 2D HYSCORE experiment which generates spectra with cross peaks between ESEEM frequencies belonging to different M S manifolds.’0,’6 So far, only a few examples of the effectiveness of HYSCORE in peak assignment have been reportedI6 and this work demonstrates the advantage (10) Hofer, P.; Grupp, A,; Nebenfiihr, H.; Mehring, M. Chem. Phys. Lett. 1986, 132, 279. (11)Mims, W. B.; Peisach, J. J . Chem. Phys. 1978, 69, 4921. (121 Flanaean. H. L.: Sineel. D. J. J . Chem. Phvs. 1987. 87. 5606. (13) HurstrG. C.; Henderson, T. pi.; Kreilick, R: W. J . Am. Chem. SOC. 1985, 107, 7294. (14) Goldfarb, D.; Fauth, J.-M.; Tor, Y.; Shanzer, A. J . Am. Chem. SOC. 1991, 113, 1941. (15) Flanagan, H. L.; Gerfen, G. L.; Lai, A.; Singel, D. J. J . Chem. Phys. 1988, 88, 2162.

(16) Shane, J. J.; van der Heijden, P. A. A. W.; Reijerse, E. J.; de Boer, E. Appl. Magn. Res. 1994, 6, 427. Reijerse, E. J.; Shane, J. J.; de Boer, E.; Hbfer, P.; Collison, D. In Electron Magnetic Resonance of Disordered Systems N; Yordanov, N. D., Ed.; World Scientific: Singapore, 1991; p 253.

Ligand Synthesis. The hydroxamate ligands were prepared as described p r e v i o u ~ l y . The ~ ~ preparation involved essentially three steps: (i) preparation of the bis-pentachlorophenolates (CH&(CHzOCHzCOOC&l&, (ii) preparation of the aminohydroxamate H2NCHiBuCONOHCHj, and (iii) coupling of the phenolates with N-hydroxylamine to provide RL239 and with aminohydroxamate to provide RL252. Both compounds were purified by column chromatography and were fully characterized by their spectroscopic features, which were in agreement with the assigned structures. Sample Preparations. 50 p L of a methanol or acetonitrile solution of the ligand (5 or 10 mM) were mixed with 10 p L of an aqueous solution of CuSO4 (20 mM). The resulting solution was then transferred into an EPR quartz tube (3 mm 0.d.) and sealed. Spectroscopic Measurements. Continuous wave (CW) EPR measurements were performed at 140-150 K on a Varian E-12 spectrometer operating at -9.3 GHz. The X-band ESEEM and 2D HYSCORE experiments were carried out at 4- 10 K at frequencies around 9.3 and 8.5 GHz, using a home built spectrometer, described e1~ewhere.l~ The magnetic field and spectrometer frequency were measured with a Bruker NMR gaussmeter (ER-035M) and a HP-5350B frequency counter. The ESEEM waveforms were recorded using the three-pulse sequence, n12t-nl2-T-nl2-z-echo, where the echo intensity is measured as a function of the time interval T. Measurements were carried out at a repetition rate of 100 Hz and t-values of 180, 220, and 300 ns. Typically, the n12 pulse length was 20 ns and T was incremented in steps of 20 ns. 2D HYSCORE spectra were recorded using the sequence nI2-t-nl2-tl -n-tz-n/2-s-echo, where the echo is measured as a function of tl and t*.lfl The duration of the n12 and n pulses was 15 ns and the amplitude of the n pulse was twice that of the n12 pulses. From 100 to 140 points were collected in each dimension, and tl and tz were incremented in steps of 40 ns from their initial value of 90 ns. The appropriate phase cycles, eliminating unwanted echoes, were employed in all experiment^.^^,^^ Data Manipulation. (a) 1D ESEEM. To avoid distortions of the FT-ESEEM spectrum due to spectrometer dead time, which for the three-pulse experiment was typically t 40 ns, the missing data points were reconstructed using the linear prediction singular value decomposition (LPSVD) method?’ Prior to Fourier transformation the background decay was subtracted either by a polynomial or exponential fit or by subtracting the appropriate zero frequency component as obtained from the LPSVD procedure, the data were convoluted with an exponential decay, and zero-filling was performed. All spectra shown are the cosine FT-ESEEM. (b) 2D ESEEM. The background decay in both tl and 12 dimensions was subtracted using a third-order polynomial fit and after convolution with a cosine roll-off function and zero filling to 256 points in each dimension Fourier transformation was carried out in the two dimensions. The spectra shown are contour plots in the magnitude mode with logarithmic scaling of the contour intervals. Basic Equations and Simulations Procedure. The spin Hamiltonian for Cu(I1) (S = 112, I = 312) interacting with two 14N, I = 1, nuclei is:

+

The last term includes the I4N nuclear interactions, which are small relative to the other interactions and in our case do not generate resolved (17) Shanzer. A.: Libman. J.: Lifson. S. US Patent 5.101.066. 1992. and US‘Patent 5.149.845, 1992. (18) Fauth, J.-M.; Schweiger, A.; Braunschweiler, L.; Forrer, J.; Emst, R. R. J . Magn. Res. 1986, 66, 74. (19) Gemperle, C.; Aebli, G.; Schweiger, A.; Emst, R. R. J . Magn. Res. 1990, 88, 241. (20) Barkhuijsen, H.; de Beer, R.; Bovee, W. M. M. J.; van Ormondt, D. J . Magn. Res. 1985, 61,465.

Coordination of Cu(II) to Lipophilic Bishydroxamate Binders

J. Am. Chem. SOC., Vol. 117, No. 51, I995 12113

a

b f'

Figure 2. Schematic structures of the Cu(I1) binding site and an arbitrary presentation of the coordinate systems. (a) Transoid configuration. (b) Cisoid configuration. The X and Y axes were placed arbitrarily and the same coordination systems applies for both configurations. splittings in the EPR spectrum. Therefore, simulations of the EPR spectra were done taking into account only the first two terms, treating the hyperfine term by second-order perturbation theory?I neglecting the 63@Cuquadrupole interaction, and assuming that the principal axes of g and A ( 6 3 W u )coincide. The ESEEM waveform can be conveniently calculated using Mims general expressions22once the spin Hamiltonian in eq 1 is brought into a block diagonal form, thus defining a nuclear Hamiltonian for each M S manifold. Diagonalization of the electron Zeeman interaction, neglecting the Cu hyperfine interaction yields the two nuclear Hamiltonians

~

2600

2800

3000

3200

3400

MAGNETIC FIELD, H (G)

Figure 3. EPR spectra of frozen solutions (150 K) of (top) Cu-RL252 in acetonitrile (9.32 GHz) and (bottom) Cu-RL239 in MeOH (9.32 GHz). The power used was 9 mW, and the modulation amplitude was 6.3 G. The dashed traces represent computer simulations obtained with the parameters given in Table 1. The asterisk denotes spectral features of a minor Cu(I1) species. The solid arrows specifies the g-values at which orientation-selective ESEEM experiments were performed and the dotted arrows correspond to positions where HYSCORE was carried out.

where V N is the nuclear Larmor frequency. The expression for YDQ was derived under the assumption that the anisotropic part of the hyperfine interaction can be neglected. Equations 4 and 5 were used to estimate e2qQlh, q, and Aiso from the experimental positions of the and the values were then refined by simulations. A detailed description of this Hamiltonian is given e l ~ e w h e r e . ' ~ . ~ ~ . peaks, ~~ The three-pulse echo intensity E(T,t) was calculated according t0'4.27 The index i denotes the two different nitrogens and AI = (Ai,l Ai)g.I4 The principal components of Ai, which is assumed to be axially symmetric and the same for both nitrogens, are -Al, -Al, and 2Al. For a small g anisotropy, aiso= Aidis,, and a1 = A s i s , . The principal values of the I4N quadrupole tensor, Q , are

+

(3)

IQLl

where q is the asymmetry parameter andjQLl 2 2 IQLI. The orientation of the external magnetic field, H,,with respect to principal is given by 8, and #,,. The principal axis system axis of g (X,Y , Z), of A, x', y', z', is transformed into that of g using the polar angles 8 and # and that of Q, x", y", z", by the Euler angles a,fi, y. Figure 2 shows a schematic representation of the binding site and the coordinate systems used. Simulated spectra were generated by numerical diagonalization of the nuclear Hamiltonians described in eq 2, and the eigenvalues and eigenvectors so obtained were used to calculate the modulation patterns using the Mims expressions.22 At exact cancellation, Y,, Y-, and Y+, and YDQ are given

(4) (21) Ovchinnikov, I. V.; Konstantinov, V. N. J . Magn. Res. 1978, 32, 179. (22) Mims, W. B. Phys. Rev. B 1972, 5, 2409; 1973, 6, 3543. (23) Gerfen, G. J.; Singel, D. J. J. Chem. Phys. 1990, 93, 4571. (24) Rowan, L. G.; Hahn, E. L.; Mims, W. B. Phys. Rev. A 1965, 137, 61. (25) Goldfarb, D.; Kevan, L. J . Chem. Phys. 1987, 87, 6323.

where A@,,&) is a weighting function accounting for the inhomogeneous EPR line width and the pulse band width. The ranges of Bo and 4" that contribute to the echo at each magnetic field within the EPR line shape were determined graphically from plots of the resonant field us e,, taking into account the different Cu(M1) hyperfine component^.'^ These were generated by simulations of the EPR spectra in which both the Cu hyperfine and the g-anisotropies were taken into account. For each Bo a range of 0-180" for 4" was taken. After the zero frequency component was subtracted, the simulated time-domain pattems were convoluted with an exponential decay function with a time-constant of typically 4 p, to account for relaxation effects, and then Fourier transformation was performed.

Results The EPR spectra of frozen solutions of Cu-RL252 and CuRL239, shown in Figure 3, are similar, and the same spectra were observed for methanol or acetonitrile solutions. While the spectrum of Cu-RL252 is indicative of the existence of only one Cu(I1) species the spectrum of Cu-FU239 shows the presence of an additional minor species. Coarse simulations of the EPR spectra were carried out in order to obtain the g and Cu hyperfine matrices from which the orientations con(26) Dikanov, S. A.; Tsvetkov, Yu. D.; Bowman, M. K.; Astashkin, A. V. Chem. Phys. Lett. 1982, 90, 149. (27) Dikanov, S. A.; Yudanov, V. F.; Tsvetkov, Yu. D. J . Magn. Res. 1979, 34, 631.

12774 J. Am. Chem. SOC., Vol. 117, No. 51, 1995

Kofman et al.

Table 1. The EPR Parameters of Cu-RL252 and Cu-RL239 Cu-RL252 Cu-RL239

All

RI I

RI

2.274 2.272

2.061 2.061

cm-I) 173 179

1

2648

Frequency, MHz

00-200

' 0 14'

0 2 4 6 8

~

1

~

,

~

1

~

1

0 2 4 6 8 10121416

k--

0 2 4 6 8 10121416

C

21 15

4-

2666

b

a

cm-I)

A1

0 2 4 6 8

Frequency,MHz

Figure 4. (a) Experimental three-pulse FT-ESEEM spectra of CuRL252 in CH3CN recorded at different magnetic fields within the EPR powder pattern (t = 0.22 ps, v = 9.15 GHz, and T = 5 K, for the traces recorded at 2648-2750 G, v = 9.32 GHz, and for 2900, 3149 G , v = 9.14 GHz). (b, c) Simulated spectra obtained with parameter sets A and B, respectively, as listed in Table 3 and the same t values and spectrometer frequency as in (a). The angles listed on the right hand side correspond to the range of eoused in the simulations.

tributing to the echo at each magnetic field within the EPR powder pattern were determined. The principal components of these matrices are given in Table 1. These values are typical for tetragonal Cu(I1) with four oxygen ligandsSz8The simulations were too coarse to detect any biaxility in g or A. Three-pulse FT-ESEEM spectra of Cu-RL252 recorded at 9.13 GHz with t = 0.22 ps and t = 0.30 ps at several magnetic fields within the EPR powder pattem (see Figure 3) are presented in Figures 4a and 5a. A signal at the proton Larmor frequency and three peaks at 1.66-1.78, 2.15, and 3.85 MHz are visible in all spectra. These peaks are relatively narrow, and their frequencies are field independent within the range 2648-3150 G. This indicates that at these fields the cancellation condition is met and we thus assign these peaks to the NQR frequencies yo, v-, and v+, respectively. Simulations have shown that the cancellation condition is met when IYN - (ai,d 2)IhN < 0.25.29 Over the range 2648-3150 G, YN varies by about 16% which is still within the cancellation range. Low intensity peaks appear at 3.27 MHz in the high field spectra and at -7.3 MHz in the low field spectra. These are tentatively assigned to the sum combinations of the 1.66 and 3.7 MHz peaks, respectively. Spectra recorded with t = 0.22 ps are very similar to those for which t = 0.3 ps (Figure 5a), the main difference being the reduced intensity of the v+ peak for t = 0.22 ps, due to the suppression effect of t.22For both values of t,a weak, narrow signal at -5.1 MHz appears in spectra recorded at magnetic fields within the region of gl. Orientation-selectivethree-pulse ESEEM spectra recorded at 8.46 GHz (with t = 0.22 ps and 0.3 ps) exhibit the same peaks as observed at 9.15 GHz. A slight shift of the 5.1 MHz peak to 5.0 MHz supports its assignment to VDQ. Similar spectra were obtained for CuRL252 in either methanol or acetonitrile solutions. Usually, the YDQ peak is broad compared to the NQR peaks (28) Peisach, J.; Blumberg, W. E. Arch. Eiochem. Eiophys. 1974, 185, 691. (29) Reijerse, E. J.; Keijzers, C. P. J . Mugn. Reson. 1987, 71, 83.

Frequency,MHr

~

-1

~

0 2

1

~

1

4 6 8

~

1

0 2 4 6 8

Frequency,MHz

Figure 5. (a) Experimental three-pulse FT-ESEEM spectra of CuRL252 in CH3CN recorded at different magnetic fields within the EPR powder pattem (t = 0.30 ps, v = 9.15 MHz, and T = 5 K, for the traces recorded at 2648, 2666 G, v = 9.32 GHz, for 2900-3150 G, v = 9.14 GHz). (b, c) Simulated spectra obtained with parameter sets A and B, respectively, as listed in Table 3 and the same t values and spectrometer frequency as in (a). The angles listed on the right hand side correspond to the range of Bo used in the simulations.

and shows a significant field d e p e n d e n ~ e . ' ~Since ! ~ ~ this behavior is not observed in our case, the VDQ peak cannot be assigned on the basis of the 1D ESEEM spectra alone. To unambiguously assign this peak, 2D HYSCORE experiments were carried out. The 2D HYSCORE spectra exhibit cross peaks only between frequencies belonging to different Ms manifolds.I0J6 Therefore, peaks at (YDQ, yo), (YDQ, v-1, and (YDQ, v+) are expected. The 2D HYSCORE spectrum of CuRL252 in CH3CN recorded at 3236 G, which corresponds to the resonant magnetic field at which the echo amplitude reaches maximum, is presented in Figure 6b. The spectrum exhibits an intense cross peak at (5.1, 3.8(v+)) MHz and a weaker one at (5.1, 2.15(v-)) MHz. The appearance of these cross peaks provides an unambiguous assignment of the 5.1 MHz peak to VDQ. The spectrum also shows a cross peak in the form of a ridge at (2-4, 1.66) MHz attributed to cross peaks between YSQ and vo. A similar spectrum was obtained at €€, = 3256 G and at other t values. A 2D HYSCORE spectrum recorded close to the gll position (2786 G, 9.33 GHz) is shown in Figure 6a. In this spectrum, the cross peak ( Y ~ Q yo) , appears as a narrow peak rather than as a ridge, clearly determining the frequency VSQ, at 3.6 MHz. The frequency of the second VSQ can now be determined according to VDQ - VSQ, = Y ~ Qyielding ~ , VSQ= ~ 1.4- 1.5 MHz. A close look at the 1D spectrum recorded at 2648 G (Figures 4a and 5a) shows a low frequency shoulder at m 1.4 MHz, which can now be attributed to v ~ Q The ~ . assignment of VSQ, indicates that the peak at 3.7-3.8 MHz, observed at low fields in the 1D spectra (Figures 4-5) is a superposition of YSQ, and v+. The low field 2D HYSCORE spectrum also shows cross peaks at (v+, YDQ) and (YO, 2 % ~ ) . The assignment of the 3.7-3.8 MHz peak observed at low fields to a superposition of V S Q ~and v+ was confirmed by a careful tracking of changes in its width at low fields as shown in Figure 4a. At the lowest field, corresponding to the gll feature of M1(Cu) = -3/2, the peak is narrow but as the field increases toward the gll position of the M1(Cu) = -1/2 state, where more orientations contribute to the echo, a broadening due to shifts of Y S Q , is observed. At higher fields, the range of orientations ~

~~~

(3O)iang, F.; Peisach, J.; Ming, L.-J.; Que, L., Jr.; Chen, V. J. Biochemistry 1991, 30, 11437.

Coordination of Cu(II) to Lipophilic Bishydroxamate Binders -5

-10

1

'

g6k

h

"

'

5

0

~

"

'

'

1

'

'

"

~

J. Am. Chem. SOC., Vol. 117, No. 51, 1995 12775

'

'

b

a

10

"

l

C

tu-RL252 H=2786G

6

a

I

I.

800-900

4

2

t

0

0 0

0

2

4

6

8 10121416

Frequency, MHz

0

2

4

6

8

0 2 4 6 8

Frequency, MHz

Figure 7. (a) Experimental three-pulse FT-ESEEM spectra of CuRL239 in MeOH recorded at 12 K and t = 0.3 ps. The spectrum recorded at 3262 G was obtained with t = 0.22 ps. The spectrum recorded at 2787 G was obtained at 9.26 GHz whereas all the other at 9.47 GHz. (b, c) Simulated spectra obtained with parameter sets A and B for Cu-RL239, respectively, as listed in Table 3. The angles listed on the right side correspond to the range of 0, used in the simulations.

/

8-

'

,

'

I

I

"

c

-10

"

-8

-5

5

0

10

Frequency (MHz) Figure 6. 2D HYSCORE spectra (contour plots) of (a) Cu-RL252 in CH3CN recorded at 5 K, t = 0.35 ps, v = 9.32 GHz, Ho= 2786 G; (b) same but H, = 3236 G (max echo position). (c) Cu-RL239 in CH3CN recorded at 12 K, t = 0.3 ps, v = 9.47 GHz, Ho= 3288 G (max echo position). Table 2. The ESEEM Frequencies Observed for Cu-RL252 and Cu-RL239 and Their Assignments Cu-RL252 (MHz)

Cu-RL239 (MHz)

1.66-1.71 2.15 3.3 3.5-3.8 1.4- 1.5 3.86-3.91 5.0 5.5 7.4

1.71-1.8 2.2 3.32 3.46 1.27 3.96 5.1-5.2

assignment

vsQ~

7.7

2vSQ

comment

vo VvSQ

vSQz

only at -gll only at -gll

Vt VOQ yo

+ v+

contributing to the echo increases significantly and only the v+ component of the peak is detected. The field dependence of V S Q ~is well illustrated by following its sum combination signal, 2VSQI,at 7-8 MHz (see Figures 4a and sa). To summarize, Table 2 lists all the observed peaks and their assignments. From the frequencies and eqs 4 and 5, we estimate e2qQlh

4.04 MHz, q 0.85, and ais0 % 1.2-1.5 MHz for the hydroxamate nitrogens in Cu-RL252. The 1D ESEEM spectra of Cu-RL239 are very similar to those of Cu-RL252. This is not surprising because the EPR spectra are not indicative of a radically different structure. Orientation-selective three-pulse FT-ESEEM spectra of CuRL239 in methanol are shown in Figure 7a. The v,, v-, and v+ peaks appear at 1.71-1.8, 2.2, and 3.9 MHz, respectively, and VDQ is found at 5.1-5.2 MHz. From these values, e2qQ/h = 4.1 MHz, q 0.83, and aisoG 1.5 MHz are calculated. These values are only slightly different from those of Cu-RL252. As in Cu-RL252, the spectrum recorded close to the gll edge of MI(CU)= -312 shows a strong peak at 3.9 MHz, assigned to VSQ,. In all spectra, the intensity of the v+ peak is somewhat reduced in comparison to the spectra of Cu-E252 (see Figure 4a). The 2D HYSCORE spectrum of Cu-RL239 recorded at 3288 G (g = 2.058) and t = 0.3 ,us (Figure 6c) confirms the assignment of the 5.1 MHz peak to VDQ and of the peak at 3.9 MHz to a superposition of v+ and VSQ, based on the cross peak at (5.1 MHz, v+) and the ridge at (2-3.9 MHz, v,). Simulations. In order to determine the anisotropic part of the hyperfine interaction and the orientation of the principal axis systems of the quadrupole and hyperfine interactions with respect to the principal axis system of g, computer simulations of the orientation-selective 1D ESEEM spectra were performed. The simulations were carried out for two chemically equivalent nitrogens related either by Cz (trans) or axz(cis) symmetry (see Figure 2).9 Three-pulse FT-ESEEM spectra calculated for the two structures showed no significant differences. The appearance of a rather narrow V S Q , peak when the field is set to gll (M*(Cu) = -3/2) indicates that at this orientation the inhomogeneous broadening due to the quadrupole interaction is minimal. Although a minimum range of the complex orientations (8,)is excited when the magnetic field is set to this position, a disorder with respect to the I4N hyperfine and quadrupole interactions is still possible, depending on the relative orientation of the principal axis system of these tensors and that of the g matrix. While the range of 8, is small, a full integration over 4, is still required, and the latter is the source of possible disorder. When 7 % 1, minimum disorder is achieved when either the y" or z" principal axes of the quadrupole tensor are approximately parallel to the gli direction, thus /3 0" or 90".

12776 J. Am. Chem. Soc., Vol. 117,No. 51, 1995

Kofman et al.

b

a

I

A i = 0.1

U

I

.

I

.

1

*

I

.

l

.

W

b

a H=2940G

I

h

H=2400G

I ~ ' ~ ~ ~ ' ~ ~ ~ ' ~ ~ ~ 4 ~ ' d . d

0 1 2 3 4 5 6 7 6 9 1 0

0 1 2 3 4 5 6 7 6 9 1 0

Frequency,MHz

Frequency, MHz

Figure 8. Simulated three-pulse FT-ESEEM spectra calculated with H,= 2400 G, v = 8.46 GHz, t = 0.22 ,us, 0, = 0-24" and B = y = 9 = 0". (a) The dependence on 8 for A1 = 0.1 MHz. The values of e2qQ/h,7, and A,,, were, from top to bottom, 3.99 MHz, 0.815, 0.82 MHz; 3.98 MHz, 0.835, 0.74 MHz; 3.98 MHz, 0.85, 0.70-MHz; 3.98 MHz, 0.85,0.68 MHz; 3.99 MHz, 0.85,0.68 MHz. (b) The dependence on A1 for 0 = 30". The values for e2qQ/h,7, and A,,, were, from top to bottom, 4 MHz, 0.85, 0.66 MHz; 3.98 MHz, 0.85, 0.68 MHz; 4 MHz, 0.845, 0.72 MHz; 4 MHz, 0.845, 0.74 MHz.

These orientations and the estimated values of e2qQ/h, 7, and a,,, obtained from the peaks positions were used as a starting point in the simulation procedure where we searched for values of a, 8, and A1 and refined the above parameters to reproduce the experimental positions and relative intensities of all peaks. This was done separately for the two cases of ,L? = 0" and 90". While many sets of parameters could reproduce the spectrum at the g l region, simulating YSQ] at the gll(MI(Cu) = -312) region was more difficult. Among the well-resolved lines in these spectra VSQ, was found to be most sensitive to 8, A l , and A,,,, and it was therefore used as the main feature for the detennination of 8 and A1 and the refinement of A,,,. For ,L? = O", the position of VSQ] could be reproduced with small values of 8, A,,, > 0, and A1 > 0 (or, alternatively, with A,,, C 0 and A1 C 0). This is shown in Figure 8a, where simulated spectra, calculated for the low field edge and A1 = 0.1 MHz, are depicted as a function of 8. In Figure 8b, 8 is kept constant at 30" and A1 is varied. Rather than keeping the parameters e2qQ/h, 7, and A,,, constant they were slightly refined for each spectrum in order to reproduce the NQR and YDQ frequencies and show how 8 and A1 were determined from the simulations. The dependence of the -7.3 MHz signal on A,,, indicates that it is due to sum combination of V S Q ] rather than of v+. This also explains why this peak is observed in the experimental spectra only at low fields (see Figures 4a and 5a) where the intensity of VSQ]is large. Comparison of the experimental spectra (Figure 4a) and the calculated spectra of Figure 8 shows that the best agreement is achieved for 8 15'. Nevertheless spectra calculated with these parameters for fields within the center of the EPR powder pattern (H,, = 2900 G) exhibit a too intense v- peak for 't = 0.3 ps. The frequency and intensity of VSQ] can be reproduced also with A,,, > 0, A1 < 0, and 8 = 90". Figure 9 shows examples of simulated spectra for A1 = -(0.1 - 0.4) MHz, calculated for the low field gll edge and for the gl position. The best agreement for both fields is achieved when A l = -0.2 MHz.

0 1 2 3 4 5 6 7 8 910

0 1 2 3 4 5 6 7 8 910

Frequency, MHz

Frequency, MHz

Figure 9. Simulated three-pulse FT-ESEEM spectra calculated for /3 = a = 0 = 90"; y = q5 = 90" and t = 0.22 ,us as a function of A1 for (a) H,= 2940 G; (b) H, = 2400 G. The values of e2qQ/h,17, and Aiso were, from top to bottom, 4.0 MHz, 0.86, 0.67 MHz; 4.0 MHz, 0.86, 0.61 MHz; 4.0 MHz, 0.88, 0.649 MHz; 4.0 MHz, 0.88, 0.40 MHz.

The same spectra were obtained for Ai,, C 0 and A l > 0. In general, the spectra were not very sensitive to Q, and y, and they were therefore arbitrarily set to 0. The simulated spectra shown in Figures 8 and 9 were calculated for v = 8.46 GHz which is lower than the frequency at which the experimental spectra shown in Figures 4 and 5 were measured (v = 9.15 GHz). The trends observed, however, are similar at both frequencies. Spectra calculated for the low field edge with a = 8 = ,L? = 90" also exhibit a V S Q ] peak at 3.71 MHz for A1 = -(0.1-0.2) MHz. The relative intensities of the v, and v- peaks, however, are not satisfactory. Moreover, spectra calculated with the same parameters for the g l position did not reproduce the experimental spectra well, the major problem being the intensity and line shape of v, and v- peaks. The case of B = 90" and a = 0" was excluded since the sum combination harmonic of V S Q ) at 7.3 MHz could not be reproduced, indicating that the signal appearing at 3.8 MHz is v+ and not V S Q ] . We failed to find any other sets of parameters and values of /3 between 0-90" that reproduced the experimental spectra satisfactorily. Finally, two best fit sets of parameters, listed in Table 3, were found. The quadrupole tensor parameters in both sets are similar, but the sets differ in A l and 8. Note that both sets yield very similar values for the hyperfine couplings along gll. Using these parameters, all orientation-selective three-pulse FTESEEM spectra in Figures 4a and 5a were simulated (Figures 4c,d and 5c,d). Both sets exhibit a good agreement with the experimental spectra, and we cannot prefer one set over the other only on the basis of their agreement with the experimental data. This demonstrates the limitations of the orientation selective experiment, where only at one field position a single crystallike spectrum is observed. Comparison of simulated spectra with the best fit parameters for one and two I4N nuclei show that while the peak at 7.3 MHz (observed at & = 2648-2700 G) is a sum combination of V S Q ] , the signal at 3.27 MHz (H = 2900-3300 G) is not a sum combination of v, but rather corresponds to a singularity of one of the single quantum powder patterns. The orientation-selective ESEEM spectra of Cu-RL239 could also be simulated satisfactorily with two sets of Hamil-

J. Am. Chem. SOC.,Vol. 11 7, No. 51, 1995 12777

Coordination of Cu(II) to Lipophilic Bishydroxamate Binders

Table 3. Best-Fit Spin Hamiltonian Parameters (Used in Figures 4, 5, and 7) As Determined from the Simulations of the Orientation-Selective Three-Pulse ESEEM Spectra

CU-RL252 Cu-RL252 CU-FU239 Cu-RL239 a

a, Y

set

le2qQ/hl (MHz 5 0.02)

q f 0.01

(deg f 10)

(MHz f 0.1")

(MHz f 0.05O)

A B A B

3.9 3.9 4.0 4.01

0.85 0.85 0.86 0.86

0 0 0 0

f 1.28 i1.28

F0.2

f1.25 51.36

F0.26 10.12

P

ais0

3

The relatively small uncertainty is due to the high sensitivity of

Y ~ to Q

tonian parameters, differing from each other mainly in the amplitude and orientation of the hyperfine interaction as listed in Table 3. The simulated spectra are presented in Figure 7b,c. As for Cu-RL252, the two options could not be distinguished. The best fit sets of parameters for Cu-RL239 and Cu-RL252 are very similar; the major differences are the slightly larger quadrupole coupling constant and anisotropic hyperfine component in Cu-RL239. This is manifested mainly in the subtle shifts of vo and YDQ and the apparent broadening of v-, which is probably a consequence of overlap with a YSQ signal. In order to differentiate between cis and trans configurations of the complexes we compared the modulation depth of the ESEEM waveforms calculated with the best fit parameters relating the parameters of the two nitrogens according to the symmetry of the binding site.9 Unfortunately, the differences observed were not significant enough to draw unambiguous conclusions regarding the configuration. We note, however, that in the case of Cu-RL239 a trans configuration is very unlikely due to steric constraints.

Discussion Ab-initio SCF-MO calculations, performed by Ruiqin et aL3' on several substituted hydroxamic acids, and previous ESEEM results on VO-RL26l9 (see Figure 1) show that the principal axis of the I4N quadrupole tensor, z", is along the nitrogen lone pair, perpendicular to the C-N-0 plane. Therefore, /.? = 0" indicates that the hydroxamate planes are perpendicular to the gll direction, which is to a good approximation along the symmetry axis of the complex, as found for Cu(1I) acetyl acetonate ( a ~ a c ) .Hence, ~~ the complex has either a square coplanar structure or an octahedral structure with an elongated tetragonal distortion. Comparison of the best fit parameters obtained for the Cu(I1) and the VO(I1) bis-hydroxamate (Q = 3.8, 4.0 MHz; 7 = 0.84, 0.91, respectively) shows that the quadrupole coupling constant and the asymmetry parameters are very similar. This implies that the quadrupole coupling constant is relatively invariant to the type of ion bound. Nevertheless, the quadrupole coupling constant in the complexes is significantly reduced compared to that of substituted hydroxamic acids (4.6-5.2 MHz).~' The reduction in the quadrupole coupling constant is attributed to some delocalization of the nitrogen lone pair over the 0-C-N-0 ~tructure.~ The comparison between similar Cu(II) and VO(II) complexes is of interest due to their complementary configurations, d9 vs d'. Superhyperfine ligand couplings are usually larger in Cu(11) complexes than in VO(I1) complexes. For example, while the isotropic hyperfine constant of directly bound nitrogens is usually in the range of 6-7 MHz for V0(II),33-37 it is 30-40 (31) Ruiqin, W.; Xiaolan, Y.; Zhenye, F.; Haq, M. M. I.; Khurshid, M. M. P.; Payner, T. J.; Smith, A. S.; Palmer, M. H. J . Am. Chem. Soc. 1989, 1 1 1 . 114. (32)Maki, A. S.; McGarvey, B. R. J . Chem. Phys. 1958, 29, 31. (33) Dikanov, S. A.; Tyryshkin, A. M.; Huttermann, J.; Bogumil, R.; Witzel, H. J . Am. Chem. SOC.1995,117,4976,Table 2 and reference therein. (34) Reijerse, E. J.; Shane, J. J.; de Boer, E.; Collison, D. In Electron Magnetic Resonance of Disordered Systems; Yordanov, N. D., Ed.; World Scientific: Singapore, 1989; p 189.

Al

FO.1

e

4

a1 (MHz f 0.1")

(deg f 10)

(deg & 10)

'i0.43 F0.21 F0.55 f0.26

90 15 90 15

0 0 0 0

the hyperfine parameters.

MHz for The large difference stems from their different electronic structures. In Cu(II), in a tetragonal symmetry, which is ubiquitous for Cu(II), the unpaired electron resides in the bl, molecular orbital, the lobes of which are directed toward the ligands as it contains contributions from both the copper dX2-).2orbital and the u orbitals of the ligands.32 The mixing of the ligand's u orbitals leads to a significant superhyperfine interaction. For vanadyl in a C4" symmetry the unpaired electron resides in a bzg orbital which is a linear combination of the vanadium dq orbital and in-plane p orbitals of the ligands?' The lobes of this molecular orbital are pointing between the ligands, and therefore, the nitrogen hyperfine is mostly a consequence of mixing in the x plane which is usually small. Hence, in VO(I1) the unpaired electron is almost completely localized on the vanadium dq orbital, and the ligand hyperfine couplings are The comparable values of I4N hyperfine interactions in the copper and vanadyl bis-hydroxamate complexes are thus unexpected. In the VO(I1) complexes also = 1.8 MHz and a1 = 0.6-0.7 M H z , ~and in the Cu(I1) complexes the values are somewhat lower, aiso=Z 1.3 MHz and a1 = 0.2-0.6 MHz. Comparison of the tetra-imidazole complexes of Cu(I1) and VO(I1) shows that the hyperfine coupling of both the directly bound nitrogen and the remote nitrogens are significantly larger in the Cu(I1) complex (39.442 and 1.7543MHz, respectively) than in the VO(I1) complexes (6-733 and 0.344 MHz, respectively). The isotropic hyperfine constant of the remote nitrogen, which is three bonds away from the metal ion, is comparable to that of the nitrogen in Cu(I1) bis(N,N-di-n-butyldithiocarbamate), aiso= 1.25 MHz, which is also three bonds away from C U ( I I ) . ~In~ this complex, aisoof the directly bound 33Sis 34.8 MHZ,"~which corresponds to a value of 32.7 MHz for I4N, considering only the different gyromagnetic ratios. This value is in the order of those found in imidazole complexes. The bis-hydroxamate binders, which are bidentate chelates, consist of two hydroxamate groups each, and the donors are the two oxygens (see Figure 1). Therefore, the ground state of the Cu(II) complexes should also be mostly of a dq character (35) Astashkin. A. V.: Dikanov. S. A.: Tsvetkov. Yu. D. J . Struct. Chem. 1985,26, 363. (36) Goslar. J.; Piekara-Sadv. L.: Kisoert. L. In Handbook of Electron Spin Resonance; Poole, C. P.,-Jr., Faralh, H. F., Eds.; AIP Press: New York, 1994; p 360. (37) Kivelson, D.; Lee, S.-K. J . Chem. Phys. 1964, 41, 1896. (38) Werst, M. M.; Davoust, C. E.; Hoffman, B. M. J . Am. Chem. SOC. 1991, 113, 1533. (39) Harrison, S. E.; Assour, J. M. J . Chem. Phys. 1965, 43, 1746. (40) Kivelson, D.; Neiman, R. J . Chem. Phys. 1961, 35, 149. (41) Ballhausen, C. J.; Gray, H. B. Inorg. Chem. 1962, 1, 111. (42) Sholl, H.-J.; Huttermann, J. J . Phys. Chem. 1992, 96, 9684. (43) Mondovi, B.; Graziani, M. T.; Mims, W. B.; Oltzik, R.; Peisach, J. Biochemistry 1977, 16, 4198. McCracken, J.; Pember, S.; Benkovic, S. J.; Villafranca, J. J.; Miller, R. J.; Peisach, J. J . Am. Chem. SOC.1988, 110, 1069. McCracken, J.; Peisach, J.; Dooley, D. M. J . Am. Chem. SOC.1987, 109, 4064. (44) Dikanov, S.A,; Burgard, C.; Huttermann, J. Chem. Phys. Lett. 1993, 212, 493. (45) Reijerse, E. J.; Paulissen, M. L. H.; Keijzers, C. P. J . Magn. Reson. 1984, 60, 66.

12778 J. Am. Chem. SOC., Vol. 117, No. 51, 1995

Kofman et al.

as in C ~ ( a c a c ) 2 ~ and ~ - ~in * Cu(I1) bis-dithi~carbamate.~~ EPR single crystal studies of Cu(I1) ketoenolato complexes showed that the ,g and g, directions, where the x and y axes were defined along the directions of the out of plane n orbitals d,, and dyz,are in between the oxygen ligands, in the nodes of the ground state orbital.46 Hence, the unpaired electron resides in a molecular orbital with a dq character, the lobes of which are directed mostly toward the similar to many cases where the ground state has a d+2 character. It is therefore expected that also in the bis-hydroxamate copper complexes, as in the imidazoles and thiocarbamates, there will be a significant spin density on the directly bound nuclei (oxygens) due to mixing with 0 orbitals. The relatively small couplings found for the nitrogen, which is similar to that found on nitrogens three bonds away and to that of the VO(I1) complex, may thus indicate that the ground state molecular orbital has a node at the nitrogen. Whether this behavior is unique to the hydroxamate binders or represents a general trend in Cu(I1) coordination to chelates is an open question due to the lack of more experimental examples and detailed MO calculations. We note, however, that a similar situation has been observed in Cu(11) bis(cis- 1,2-dicyanoethylenedithiolate),where the hyperfine couplings of the I3Cl,2situated two bonds away from the Cu(11) (qs0= -4.0, -3.81 MHz) are smaller than those of the I3C3,4which are three bonds away (also= 5.83, 5.56 MHZ).~O In this case, extended Hiickel molecular orbital calculations were performed, showing that the coefficients of the atomic orbitals of C1.2 in the wave function of the unpaired electron are smaller than those of C3.4. Using the values obtained for I3Cl,2,a rough estimate of = l MHz is obtained for a I4N two bonds away. This is the order of the values found in Cu-RL252 and CuRL239. One, however, has to remember that the comparison with the dithiolene complex is rather coarse since our system involves negatively charged oxygen donors rather than neutral sulfurs. The I4N isotropic hyperfine coupling of the VO(I1) complexes, 1.8 MHz, is between the values obtained for directly bound nitrogen and nitrogen situated three bonds away. This value is comparable to other hyperfine couplings with nuclei situated 2 bonds away from the vanadium in VO(II) complexes such as values of 5.4 and 5.0 MHz for I3C in VO(glicinate)2 and V0([2-13C]pyruvate)2.51These values scale down to 1.54 and 1.42 MHz for I4N. For both Cu-RL252 and Cu-RL239, the simulations generated two possible sets of parameters for the hyperfine interaction (see Table 3). In the first set (A), the principal axis z’ is approximately perpendicular to the 0-C-N-0 plane (8 15”), whereas in the second set (B), it is in the plane (0 90”). Set A implies that the major contribution to the anisotropic hyperfine interaction comes from spin density in the pn orbital of the nitrogen. The mechanism leading to a finite spin density in the nitrogen pn orbital could be exchange polarization from the dV orbital to the dxzand/or dyzorbitals involved in the n bonding. In this case, however, the anisotropic hyperfine interaction should not have axial symmetry as was initially assumed. The second option suggests that the anisotropic hyperfine interaction is dominated mainly by the dipole-dipole interaction between the unpaired electron on the metal ion and the nitrogen. The value of lull = 0.43 MHz corresponds to a

- -

(46) Hitchman, M. A.; Olson, C. D.; Belford, R. L. J . Chem. Phys. 1%9, 50, 1195.

(47) Funck, L. L.; Ortolano, T. R. Inorg. Chem. 1968, 7, 507. (48) Hathaway, B. J.; Billing, D. E. Coord. Chem. Rev. 1970, 5, 143. (49) Keijzers, C. P.; Smaathrorst. Chem. Phys. Lett. 1980, 69, 2. (50) Kirmse, R.; Stach, J.; Abram, U.; Dietzch, W.; Bottcher, R.; Gribnau, M. C. M.; Keijzers, C. P. Inorg. Chem. 1984, 23, 3333. (51) Tipton, P. A.; McCracken, J.; Cornelius, J. B.; Peisach, J. Biochemisfry 1988, 28, 5720.

Cu-N distance of 2.35 A. Rough calculations of the Cu-N distance taking the Cu-0 distance as 2 A, the N-0 bond as 1.4 A, and a CU-0-N angle of 120” yields a Cu-N distance around 2.8 A. This is significantly larger than the experimental value and indicates that the point-dipole approximation is not valid and that there are additional contributions to the anisotropic hyperfine coupling. This is not surprising as in Cu(II) complexes with substituted imidazole a~ of the remote nitrogen was found to depend on the ~ u b s t i t u e n t . ~ ~ Finally, the differences found between the Cu-RL239 and Cu-RL252 systems are small, manifested mainly in the magnitude of a l , which for both sets A and B was larger in Cu-RL239. Interestingly, a l is also larger in VO-RL239 than in VO-RL261 .9 While both RL252 and RL261 contain amino acid bridges, alanine and leucine, respectively, this bridge is absent in RL239. The different properties of the amide and the 0-CH2 groups apparently affect the electron distribution in the binding sites, leading to small changes in al. Inspection of molecular models suggests that the square planar Cu(I1) complex of RL239 is only feasible in a cisoid, rather than a transoid arrangement. The remarkably sharp signals of CuRL252 are suggestive of a highly defined complex that adopts a single conformation. Whether the preferred configuration for the more flexible Cu-RL252 complex is cisoid or transoid is not possible to determine from the experimental results.

Conclusions The hyperfine coupling of the nitrogens in the hydroxamate groups falls within the “cancellation condition” range at 8-9 GHz leading to an exceptionally well-resolved ESEEM spectrum consisting of the NQR frequencies, VDQ and VSQ. The latter is rarely detected in orientationally disordered systems. Complete assignment of all signals in the orientation-selective ESEEM spectra was achieved using 2D HYSCORE experiments. Computer simulations of the orientation-selective 1D ESEEM spectra gave the 14N hyperfine and quadrupolar parameters which indicated that the structure of the binding site of both Cu-RL252 and the parent Cu-RL239 is close square coplanar. The appearance of Y ~ Qprovided excellent constraints on the simulations, leading to a more accurate determination of the Hamiltonian parameters. The relatively small I4N hyperfine coupling, found for both complexes, is attributed to a node in the molecular orbital, occupied by the unpaired elec‘tron, at the nitrogen. Only minor differences in the Hamiltonian parameters of Cu-RL252 and Cu-RL239 were detected, the quadrupole parameters and the anisotropic hyperfine component being slightly larger in Cu-RL239. These subtle differences indicate that the amido bridge practically does not affect the binding site geometry. Acknowledgments. This research was supported by the Israel Science Foundation administered by the Israel Academy of Sciences and Humanities and partially by the US-Israel Binational Science Foundation (Grant 90-00128). We thank Dr. Ed. Reijerse from the University of Nijmegen for performing C-band (4-5 GHz) ESEEM measurements and Dr. D. van Ormondt for giving us the LPSVD computer program. S.A.D. thanks the Department of Chemical Physics at the Weizmann Institute of Science for supporting his stay at this Institute and for the hospitality shown. Part of this work was also sponsored by the U S . Department of Energy under contract DE-AC0676 RLO and by Associated Westem Universities, Inc., Northwest Division, under Grant DE-FG06-89 ER-75522 or DEFG06-92RI-12451 with the U.S. Department of Energy. JA95 19795 (52) Jiang, F.; McCracken, J.; Peisach, J. J . Am. Chem. Soc. 1990, 112, 9035.