Sep 28, 1983 - *Department of Chemistry, University of Manchester Institute of Science and Technology, P.O. Box 88,. Manchester M60 IQD, U.K., tDaresbury ...
985
Biochem. J. (1984) 219, 985-990 Printed in Great Britain
An extended-X-ray-absorption-fine-structure study of bovine erythrocyte superoxide dismutase in aqueous solution Direct evidence for three-co-ordinate Cu(I) in reduced enzyme
Ninian J. BLACKBURN,* S. Samar HASNAIN,f Norman BINSTED,4 Gregory P. DIAKUN,f C. David GARNER: and Peter F. KNOWLES§ *Department of Chemistry, University of Manchester Institute of Science and Technology, P.O. Box 88, Manchester M60 IQD, U.K., tDaresbury Laboratory, Science and Engineering Research Council, Warrington WA4 4AD, Cheshire, U.K., IDepartment of Chemistry, University of Manchester, Manchester M13 9PL, U.K., and §Department of Biophysics, University of Leeds, Leeds LS2 9JT, U.K.
(Received 28 September 1983/Accepted 26 January 1984) Copper and zinc K-edge e.x.a.f.s. (extended X-ray-absorption fine structures) were measured for the metal sites of oxidized and reduced bovine superoxide dismutase in aqueous solution. Detailed analysis of the spectra indicates that the copper site of the enzyme changes on reduction and is most probably co-ordinated to three imidazole groups at a shorter distance Cu-N(a)=0.194nm (1.94A) in the reduced form compared with a co-ordination of four imidazole groups at 0. 199nm (1.99 A) and an oxygen atom from solvent water at 0.224nm (2.24A) in the oxidized form. Examination of the edge, near-edge structure and e.x.a.f.s. of the zinc sites indicates that the stereochemical changes at copper that accompany reduction introduce minimal perturbation on the stereochemistry at zinc.
Cu,Zn superoxide dismutases (superoxide: superoxide oxidoreductase, EC 1.15.1.1) are widely distributed in both plant and animal kingdoms and are found in eukaryotic cytosols. Although some controversy surrounds the true catalytic function of these proteins (Fee, 1982), they exhibit catalytic activity towards the dismutation of superoxide anion according to the equation: 202- + 2H+ -O2+ H2O2 The proximity of the metal centres and the early availability of a crystal structure to 0.3nm resolution (Richardson et al., 1975) stimulated interest in the co-ordination chemistry of the copper and zinc sites, which have since been investigated by a wide variety of spectroscopic techniques (Fee & Valentine, 1977; Fee, 1981; Valentine & Pantoliano, 1981). Many of the conclusions deduced from the spectroscopy have been confirmed by the recent crystallographic refinement of the structure to 0.2nm (Tainer et al., 1982), ii which the copper ligands (nitrogen atoms of imidazole groups of His44, -46, -61 and -118) show an uneven tetrahedral Abbreviations used: e.x.a.f.s., extended X-ray-absorption fine structure; x.a.n.e.s., X-ray-absorption nearedge structure.
Vol. 219
distortion from a square plane, with the fifth axial co-ordination position occupied by solvent (Fee & Gaber, 1972; Boden et al., 1979). His-61 forms a bridge to the zinc, the other ligands of which are N atoms of His-69 and -78 and an oxygen atom of Asp-81. The geometry around zinc is tetrahedral, with a strong distortion towards a trigonal pyramid having the buried Asp-81 at the apex. Tainer et al. (1982) make the point that the presence of non-crystallographic symmetry relating the four copper and zinc sites of the two superoxide dismutase dimers per asymmetric unit introduces ambiguity as to whether the observed inequivalence in metal-ligand bond lengths reflects structural differences between the four copper and four zinc sites or results from lack of resolution in the refinement. They argue in favour of the latter and quote Cu-N and Zn-N bond lengths of about 0.21 nm and a shorter Zn-0 bond length from Asp-81 of nearer 0.20nm for bond lengths averaged over all the metal sites. E.x.a.f.s. can in principle determine bond lengths to a precision of + 0.002nm in metalloproteins (Cramer & Hodgson, 1979; Teo, 1981), and therefore offers a method for improving the structural characterization of the local environ-
N. J. Blackburn and others
986 ment of the metal sites by providing precise bond lengths. We have recently reported an e.x.a.f.s. study of the copper(II) and zinc(II) sites in freezedried bovine superoxide dismutase (Blackburn et al., 1983), but ambiguities arising from structural changes that accompany freeze-drying (Strothkamp & Lippard, 1982) preclude relating that study to the crystallographic data. However, e.p.r. studies (Lieberman et al., 1981) have indicated that the crystal structure is most probably maintained in aqueous buffered solution. In the present paper we report a detailed e.x.a.f.s. study of bovine superoxide dismutase in aqueous 10mM-phosphate buffer that, in addition to precise determination of metal-protein contacts, has allowed location of a solvent molecule co-ordinated to copper at 0.224nm. In marked contrast with the native Cu(II),Zn(II) protein, structural information on the Cu(I),Zn(II) form is limited. A vacant co-ordination position is known to be available on Cu(I) (Fee & Ward, 1976), and the pH-dependence of the redox potential (Fee & DiCorletto, 1973; Lawrence & Sawyer, 1979) has been taken to imply loss of a coordinated imidazole group on reduction of the copper. We have undertaken an e.x.a.f.s. study of the reduced form of bovine superoxide dismutase in aqueous solution, and report here the first direct structural evidence that copper is three-co-ordinate in this form of the enzyme. Materials and methods Superoxide dismutase was freshly prepared from bovine erythrocytes by the method of McCord & Fridovich (1969) with the modifications and additions described previously (Blackburn et al., 1983). The enzyme samples used for X-rayabsorption measurements [5 mm in Cu(II) and Zn(II)] were dissolved in 10mM-sodium phosphate buffer, pH 7.4. Reduction of the enzyme was performed by addition of a 2-fold excess of solid Na2S204, followed by rapid transfer to a sealed e.x.a.f.s. cell that was itself bathed in a stream of N2 during the run. For the copper site, two spectra were recorded; individual spectra were carefully examined for changes in edge position, edge structure or e.x.a.f.s. that might result from radiation damage during data collection. In particular, the intensity of the sharp peak on the absorption edge characteristic of the Cu(I) form of the enzyme (Fig. 1, and Blumberg et al., 1978) was used to monitor the occurrence of any photo-redox reactions during the measurements. No such changes were found for either oxidized or reduced forms of the enzyme. For zinc, only a single scan was made. X-ray-absorption spectra were recorded in the
fluorescence mode as fluorescence excitation spectra with the recently commissioned fluorescencedetection system (Hasnain et al., 1983) at the Daresbury Synchrotron Radiation Source operating at an energy of 2 GeV with an average current of 150mA. Harmonic contamination was minimized by using the double-crystal Si 220 monochromator (Greaves et al., 1983), so that e.x.a.f.s. amplitudes can be accurately measured. However, to minimize the monochromator instability, a delay of 2.0s per energy step was found to be necessary. Data analysis utilized the single scattering spherical-wave method for calculating e.x.a.f.s. with phase shifts derived from 'ab-initio' calculations as described previously (Lee & Pendry, 1975; Perutz et al., 1982). The quality of fit and refinement was carried out as described previously (Blackburn et al., 1983). Results Absorption edges and near-edge structure (x.a.n.e.s.) Fig. 1 shows the Cu K-absorption edge and x.a.n.e.s. for oxidized and reduced superoxide dismutase. The edge features marked with arrows 0.3 -
.A. *1
t:
.
0-
\*-J
(a)
-0.3 -
1.
.
t
t
-0.6-
I (b)
*I: -1.0-
..
.
.
.
.
..0
:-
':
-0.6-
-1.0 -
8900
8950
9000
9050
9100
Energy (eV) Fig. 1. Cu K-absorption edges and near-edge structure for oxidized and reduced bovine superoxide dismutase (a) Oxidized enzyme; (b) dithionite-reduced enzyme. Arrows indicate edge features that occur at the same energy in both spectra. The Cu K-edge of metallic copper foil was used to calibrate the monochromator.
1984
E.x.a.f.s. of superoxide dismutase
987 0.2-
0.3-
(a)
0.1 s; .
*
~~~~~~~(c) 0-
'A0''
[441
.
.AI=r rl r?p
T~~~~~~~~~f
4.
-0.3 *+ ' ', ,'
*
~~~~~~(b)
-0.1
~~~~~~~~(a)
-0.2
-0.6 _. .
-1.0
.z
-o.3
0-v. .11
O..2
-1.0 -
9600
(b)
O..1
9650
9700
9750
9800
Energy (eV)
Fig. 2. Zn K-absorption edge and near-edge structure for oxidized and reduced bovine superoxide dismutase (a) Oxidized enzyme; (b) dithionite-reduced enzyme; (c) freeze-dried oxidized enzyme. The Zn Kedge of metallic zinc foil was used to calibrate the monochromator.
-0..1
-0.
/
.2-i
-0.3
4
6
5
8
7
9
10
k(A-') and also the half-height of the absorption edges appear at the same energy in both cases, but show clear differences in relative intensity, suggestive of stereochemical change accompanying reduction. Fig. 2 shows a similar comparison of zinc K-edge and x.a.n.e.s. region. In marked contrast with the copper sites, the close correspondence between all prominent x.a.n.e.s. features requires that reduction of the copper site causes no substantial change in co-ordination number, bond lengths or bond angles at zinc.
E.x.a.f.s. region Figs. 3(a) and 3(b) show the e.x.a.f.s. spectra of the oxidized and the reduced bovine superoxide dismutase in solution. The spectra are typical of imidazole co-ordination (Blackburn et al., 1983), but the differences in local structure apparent in the edge region are also reflected in the e.x.a.f.s. of the copper sites. The main differences between the two e.x.a.f.s. spectra occur in the first two (low-k) beat regions, k = 4-5 A - and k = 5.5-7 A- 1. (k is quoted in A-l, which is equivalent to I0nm-1.) The overall amplitude of the reduced-enzyme Vol. 219
Fig. 3. Cu K-edge e.x.a.s. ofoxidized and reduced bovine superoxide dismutase E.x.a.f.s. [x(k)] is presented as kX(k) plotted versus k, where k is in A- and is equivalent to 10nm- 1. Continuous lines represent the experimental spectra, broken lines the simulated spectra. (a) Oxidized enzyme, fit-index parameter 0.31; (b) reduced enzyme, fit-index parameter 0.30. Samples were dissolved in 10mM-phosphate buffer, pH7.4. I
spectrum (Fig. 3b) is seen to have decreased to approx. 80% of the oxidized-enzyme spectrum, suggesting a decrease in co-ordination number.
Analysis of e.x.a.f.s. of the oxidized copper site Simulation of the e.x.a.f.s. spectrum of the site of oxidized superoxide dismutase in solution initially utilized 'best-fit' parameters for the freeze-dried enzyme as being the closest model for Cu(II),Zn(II) superoxide dismutase in solution. These parameters gave a fit index defined as [k{x,,p.(k) Xthr.(k)}2]
copper
-
k
of 0.62 and
a very poor
fit in the k=7-9A-
988
N. J. Blackburn and others
region. Least-squares refinement resulted in a lengthening of one set of Cu-C(p) distances, and the fit index was lowered to 0.46. However, the fit was visually worse in the first e.x.a.f.s. beat (k = 45 A - 1) and did not reproduce the low-k shoulder of the second beat. Inclusion of an oxygen atom refining at 0.224nm corrected these defects and resulted in a further improvement of fit index to 0.31. The final fit including this oxygen atom is shown in Fig. 3(a). We note that, in our simulations of the enzyme solution spectra, the Debye-Waller factor a2 was left unaltered from that used for the freeze-dried enzyme. This is reasonable, since we do not expect the static disorder to change significantly within the same enzyme system. The structural parameters as obtained from the simulations are given in Table 1 and are compared with those reported for the freeze-dried enzyme. Analysis of e.x.a.fs. of reduced copper site The simulation of the e.x.a.f.s. spectrum of the reduced enzyme was also carried out by using starting parameters derived from the freeze-dried oxidized enzyme. This immediately showed that the phase of the first-shell oscillation was incorrect, suggesting a shortening of Cu-N(a) distance, and that the amplitude of the experimental e.x.a.f.s. was significantly less. An overall fit index of 0.85 was obtained with this simulation. An attempt to refine resulted in a shortening of the first-shell CuN(a) distance to 0.194nm and an immediate decrease of fit index to 0.46 due to matching of the phase of experimental and simulated spectra. The quality of fit further improved when the coordination was decreased to three imidazole groups. The fit index improved to 0.30 and the simulated and experimental e.x.a.f.s. matched well in amplitude. The addition of an oxygen atom that refined at 0.218nm gave a slightly poorer fit, the fit index deteriorating to 0.33. Therefore the final
simulation shown in Fig. 3(b) is based on the coordination of the copper by three imidazole groups, Table 1 gives the details of the parameters used. Zinc near-edge structure and e.x.a.f.s. The near-edge region (including the initial portion of the e.x.a.f.s.) shown in Fig. 2, which indicates close similarities in stereochemistry at zinc for oxidized and reduced superoxide dismutase, are virtually identical with the x.a.n.e.s. of the zinc site of freeze-dried enzyme, also shown in Fig. 2 for comparison. This observation allows us to predict that zinc-imidazole and zinc-aspartate bond lengths are insensitive to changes in oxidation state or co-ordination number at copper, and are identical with those derived for freeze-dried superoxide dismutase. The fluorescence e.x.a.f.s. data for the zinc sites were of insufficient quality to allow meaningful simulations to be carried out, but a theoretical fit that used bond lengths and DebyeWaller factors obtained for the freeze-dried zinc site was found to be consistent with the experimental e.x.a.f.s. of both oxidized and reduced superoxide dismutase solutions. Discussion The interpretation of e.x.a.f.s. data for the copper site of bovine erythrocyte superoxide dismutase in aqueous 10mM-phosphate buffer, pH 7.4, indicates a co-ordinate structure in close agreement with the most recent crystallographic results, reported by Tainer et al. (1982), namely four imidazole groups with Cu-N(s) bond length 0.199 nm and one solvent water molecule with Cu0 bond length 0.224nm as ligands to copper. The Cu-N(a) bond lengths are shorter than the crystallographic value of 0.21 nm, determined as the average over the four copper sites of the asymmet-
Table 1. Parameters used to simulate the e.x.af s. associated with the copper K-edge of oxidized and reduced superoxide dismutase in aqueous solution Results for freeze-dried native enzyme are also included for comparison. Note that Debye-Waller factors (a2) are quoted in nm2. Solution Freeze-dried
Atom N(a) 0
Number 4
a2 (nm)2
Oxidized R (nm) 0.200
Number 4 1 4 0.0001 4 0.292 C(p) 4 4 0.0001 0.305 C(p) 4 4 0.0001 0.383 N(y)* 4 4 0.0001 0.384 C(y)* * These distances are expected to be underestimated by
0.0001
a2 (nm)2 0.0001 0.0001 0.0001 0.0001 0.0001
Reduced
a2 (nm)2 Number R (nm) 0.199 0.0001 3 0.224 0.296 3 0.0001 0.0001 0.316 3 0.390 3 0.0001 0.0001 0.0001 3 0.380 approx. 0.02nm (Perutz et al., 1982).
R (nm) 0.194
0.298 0.308 0.390
0.380
1984
E.x.a.f.s. of superoxide dismutase ric unit. However, the e.x.a.f.s. value agrees closely with Cu-N(,) distances reported for similar Cu(II)-imidazole complexes, e.g. [Cu(imid)4](NO3)2 (McFadden et al., 1976) and [Cu(imid)3H20]S04 (Fransson & Lundberg, 1974), where Cu-N(a) bond lengths are 0.201 and 0.200nm respectively. The e.x.a.f.s. value is also within 0.001 nm of that determined by e.x.a.f.s. for the freeze-dried derivative of superoxide dismutase (Table 1). We have previously shown (Blackburn et al., 1983) that the data-analysis routines used here can reproduce the crystallographic Cu-N(,) distance for [Cu(imid)4](NO3)2 to within 0.003nm, and this may be taken as an indication of the likely error in the bond lengths we report in the present paper for the solution superoxide dismutase structure. The e.x.a.f.s. interpretation also locates an oxygen atom at 0.224nm. Both crystallography (Tainer et al., 1982) and water proton relaxation studies (Fee & Gaber, 1972; Boden et al., 1979) provide strong evidence for the presence of solvent in an axial-type position. The distance of 0.224nm found for this oxygen atom as the fifth ligand is consistent with an axially co-ordinated water molecule. Examples of crystallographically characterized tetragonal complexes containing axial water include the Tutton salts [(NH4)2Cu(H20)6](SO4)2 (Duggan et al., 1979) and [Rb2Cu(H20)6j(SO4)2 (Smith et al., 1975). The crystal structures of these compounds have been determined at 77K, where the tetragonal distortion becomes static, thereby allowing the true axial bond lengths to be determined under non-fluxional conditions (Hathaway, 1981). Axial bond lengths of 0.2278 and 0.2317nm are found in these complexes respectively. Likewise, e.x.a.f.s. studies on [Cu(H2O)6]2+ ions in aqueous solution, where fluxional motion is slow on the e.x.a.f.s. time scale, have identified two Cu-O bond lengths at 0.21 and 0.23nm for equatorial and axial co-ordination respectively (I. Ross, personal communication). For the reduced form of the enzyme, three-coordinate Cu(I) involving only imidazole co-ordination has produced the lowest fit index. The DebyeWaller term (a2) was fixed, throughout the entire fitting procedure for the reduced enzyme, at the value obtained for the oxidized native enzyme in solution, where the co-ordination number (N) is known (Tainer et al., 1982). This procedure helps to minimize the uncertainty in co-ordination numbers of equidistant ligands that results from floating two highly correlated parameters a2 and N in fitting routines, and amounts to using native superoxide dismutase as a model compound for the unknown reduced structure. The 33% decrease in fit index that resulted when N was decreased from four to three while keeping a2 fixed is therefore
Vol. 219
989 strong evidence in favour of a three-imidazole structure in the reduced enzyme. However, inclusion of an oxygen atom refining at 0.218nm resulted in only slight deterioration in quality of fit, so that the presence of a water molecule as the fourth ligand cannot be excluded on the basis of e.x.a.f.s. amplitudes alone. The value of the Cu-N(a) bond lengths of 0. 194nm in the reduced form is itself indicative of a co-ordination number lower than four, since reduction of the copper is expected to lengthen the Cu-N(