Spectroscopic and Computational Studies of the Azide-Adduct of

Gayan B. Wijeratne , Briana Corzine , Victor W. Day , and Timothy A. Jackson ..... Understanding the influence of the protein environment on the Mn(II...
0 downloads 0 Views 549KB Size
Published on Web 09/14/2004

Spectroscopic and Computational Studies of the Azide-Adduct of Manganese Superoxide Dismutase: Definitive Assignment of the Ligand Responsible for the Low-Temperature Thermochromism Timothy A. Jackson,† Anush Karapetian,‡ Anne-Frances Miller,‡ and Thomas C. Brunold*,† Contribution from the Department of Chemistry, UniVersity of Wisconsin-Madison, Madison Wisconsin 53706 and Department of Chemistry, UniVersity of Kentucky, Lexington, Kentucky 40506 Received March 26, 2004; E-mail: [email protected]

Abstract: A variety of spectroscopic and computational techniques have been used to examine the thermochromic transition previously reported for the oxidized state of Mn-dependent superoxide dismutase from E. coli in the presence of substrate analog azide (N3-Mn3+SOD).[Whittaker, M. M.; Whittaker, J. W. Biochemistry 1996, 35, 6762-6770.] Although previous spectroscopic studies had shown that this thermochromic event corresponds to a change in coordination number of the active-site Mn3+ ion from 6 to 5 as temperature is increased, the ligand that dissociates in this conversion had yet to be identified. Through the use of electronic absorption, circular dichroism (CD), and magnetic CD (MCD) spectroscopies, both dfd and ligand-to-metal charge-transfer (LMCT) transition energies have been determined for native Mn3+SOD (possessing a five-coordinate Mn3+ center) and Y34F N3-Mn3+SOD (forming a six-coordinate N3-Mn3+ adduct at all temperatures). These two systems provide well-defined reference points from which to analyze the absorption and CD data obtained for N3-Mn3+SOD at room temperature (RT). Comparison of excited-state spectroscopic data reveals that Mn3+SOD and RT N3-Mn3+SOD exhibit virtually identical dfd transition energies, suggesting that these two species possess similar geometric and electronic structures and, thus, that azide does not actually coordinate to the active-site Mn3+ ion at RT. However, resonance Raman spectra of both N3-Mn3+SOD and Y34F N3-Mn3+SOD at 0 °C exhibit azide-related vibrations, indicating that azide does interact with the active site of the native enzyme at this temperature. To gain further insight into the nature of the azide/Mn3+ interaction in RT N3-Mn3+SOD, several viable active-site models designed to promote either dissociation of coordinated solvent, Asp167, or azide were generated using DFT computations. By utilizing the time-dependent DFT method to predict absorption spectra for these models of RT N3-Mn3+SOD, we demonstrate that only azide dissociation is consistent with experimental data. Collectively, our spectroscopic and computational data provide evidence that the active site of N3-Mn3+SOD at RT exists in a dynamic equilibrium, with the azide molecule either hydrogenbonded to the second-sphere Tyr34 residue or coordinated to the Mn3+ ion. These results further highlight the role that second-sphere residues, especially Tyr34, play in tuning substrate (analog)/metal ion interactions.

1. Introduction

Superoxide dismutases (SODs)1 are a family of enzymes that utilize either an Fe, Mn, Cu/Zn, or Ni cofactor to disproportion† ‡

Department of Chemistry, University of Wisconsin-Madison. Department of Chemistry, University of Kentucky.

(1) Abbreviations: ADF, Amsterdam density functional; CD, circular dichroism; CT, charge transfer; DFT, density functional theory; DOMO, doubly occupied molecular orbital; EDDM, electron density difference map; EPR, electron paramagnetic resonance; HOMO, highest occupied molecular orbital; INDO/S-CI, intermediate neglect of differential overlap/spectroscopic parametrization - configuration interaction; LMCT, ligand-to-metal charge transfer; LT, low temperature; LUMO, lowest-unoccupied molecular orbital; MCD, magnetic circular dichroism; MO, molecular orbital; NIR, near-infrared; PDB, protein data bank; RR, resonance Raman; RT, room temperature; SDS-PAGE, sodium dodecyl sulfate - polyacrylamide gel electrophoresis; SOD, superoxide dismutase; SOMO, singly occupied molecular orbital; TD, time-dependent; UV, ultraviolet; VTVH, variabletemperature, variable-field; WT, wild type. 10.1021/ja0482583 CCC: $27.50 © 2004 American Chemical Society

ate the toxic superoxide radical anion, O2.-, that is produced as a byproduct during photosynthesis and aerobic respiration and also plays a major role in reperfusion injury.2,3 The Mndependent SODs (MnSODs) are present in a variety of prokaryotes and eukaryotes. They are highly homologous to FeSODs but structurally unrelated to Cu/Zn- and NiSODs. X-ray crystallography has been used to determine the structures of MnSODs from a variety of diverse organisms, including humans,4,5 bacteria (i.e., Escherichia coli and the thermophilic (2) (a) Packer, L., Ed.; Methods in Enzymology: Superoxide Dismutase; Academic Press: San Diego, 2002; Vol. 349. (b) Miller, A.-F.; Sorkin, D. L. Comments Mol. Cell. Biophys. 1997, 9, 1-48. (c) Miller, A.-F. In ComprehensiVe Coordination Chemistry II; McCleverty, J. A., Meyer, T. J., Eds.; Elsevier Ltd.: Oxford, UK, 2004; Vol. 8, pp 479-506. (d) Miller, A.-F. Curr. Opin. Chem. Biol. 2004, 8, 162-168. (3) Valentine, J. S.; Wertz, D. L.; Lyons, T. J.; Liou, L.; Goto, J. J.; Gralla, E. B. Curr. Opin. Chem. Biol. 1998, 2, 253-262. J. AM. CHEM. SOC. 2004, 126, 12477-12491

9

12477

Jackson et al.

ARTICLES

Figure 1. Active-site structures of E. coli Mn3+SOD (left) and T. thermophilus N3-Mn3+SOD (right) based on coordinates from the protein data bank files 1VEW and 1MNG-L, respectively. Dashed lines indicate important H-bonding interactions. Note that the E. coli amino acid numbering scheme is used for both species.

bacterium Thermus thermophilus),6-14 and a fungus that has been implicated in certain lung diseases.15 These studies have invariably revealed a distorted trigonal bipyramidal ligand environment around the manganese ion, consisting of two histidines (His81 and His171, using the E. coli MnSOD numbering scheme)16 and an aspartate residue (Asp167) in the equatorial plane and a third histidine residue (His26) and a solvent molecule in the axial positions (Figure 1, left). In addition to these active-site ligands, several second-sphere residues such as Gln146 and Tyr34, which are part of a hydrogen-bonding network involving coordinated solvent and Asp167 (Figure 1, dashed lines), are highly conserved among MnSODs. Experimental and computational studies have suggested that the protonation state of the solvent ligand is coupled to the oxidation state of the metal ion, with the oxidized state coordinating a hydroxide ion and the reduced state binding a water molecule (see eq 1).17-20 Thus, coordinated solvent may provide one of the two protons required in the two-step catalytic mechanism originally elucidated by Fridovich and co-workers based on pulse radiolysis studies (eq 2).21

Mn3+SODOH + e- + H+ f Mn2+SODH2O Mn3+SODOH + O2•- + H+ f Mn2+SODH2O + O2

(1) (2a)

Mn2+SODH2O + O2•- + H+ f Mn3+SODOH + H2O2 (2b) (4) Børgstahl, G. E. O.; Parge, H. E.; Hickey, M. J.; Beyer, W. F.; Hallewell, R. A.; Tainer, J. A. Cell 1992, 71, 107-118. (5) Wagner, U. G.; Pattridge, K. A.; Ludwig, M. L.; Stallings, W. C.; Werber, M. M.; Oefner, C.; Frolow, F.; Sussman, J. L. Protein Sci. 1993, 2, 814825. (6) Edwards, R. A.; Baker, H. M.; Whittaker, M. M.; Whittaker, J. W.; Jameson, G. B.; Baker, E. N. J. Biol. Inorg. Chem. 1998, 3, 161-171. (7) Børgstahl, G. E. O.; Pokross, M.; Chehab, R.; Sekher, A.; Snell, E. H. J. Mol. Biol. 2000, 296, 951-959. (8) Schmidt, M.; Meier, B.; Parak, F. J. Biol. Inorg. Chem. 1996, 1, 532541. (9) Liao, J.; Li, M.; Liu, M.-Y.; Chang, T.; Le Gall, J.; Gui, L.-L.; Zhang, J.-P.; Liang, D.-C.; Chang, W.-R. Biochem. Biophys. Res. Comm. 2002, 294, 60-62. (10) Liao, J.; Liu, M.-Y.; Chang, T.; Li, M.; Le Gall, J.; Gui, L.-L.; Zhang, J.-P.; Jiang, T.; Liang, D.-C.; Chang, W.-R. J. Struct. Biol. 2002, 139, 171180. (11) Ludwig, M. L.; Metzger, A. L.; Pattridge, K. A.; Stallings, W. C. J. Mol. Biol. 1991, 219, 335-358. (12) Lah, M. S.; Dixon, M. M.; Pattridge, K. A.; Stallings, W. C.; Fee, J. A.; Ludwig, M. L. Biochemistry 1995, 34, 1646-1660. (13) Parker, M. W.; Blake, C. C. J. Mol. Biol. 1988, 199, 649-171. (14) Atzenhofer, W.; Regelsberger, G.; Jacob, U.; Peschek, G. A.; Furtmu¨ller, P. G.; Huber, R.; Obinger, C. J. Mol. Biol. 2002, 321, 479-489. (15) Flu¨ckiger, S.; Mittl, P. R. E.; Scapozza, L.; Fijten, H.; Folkers, G.; Gru¨tter, M. G.; Blaser, K.; Crameri, R. J. Immunol. 2002, 168, 1267-1272. 12478 J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004

Small anions such as azide (N3-) and fluoride (F-) have been frequentlyusedtomodelsubstrate/metalioninteractions.11,12,17,22-25 Studies of Mn3+SODs in the presence of azide (yielding N3Mn3+SOD) have shown that these systems exhibit spectroscopic properties that are temperature dependent.23,24 While a room temperature (RT) X-ray structure of N3-Mn3+SOD from T. thermophilus revealed that N3- binding to Mn3+SOD increases the coordination number of the Mn3+ center from 5 to 6, generating a distorted octahedral geometry (Figure 1, right),12 previous studies utilizing RT absorption and circular dichroism (CD) spectroscopies had suggested that azide binding to E. coli Mn3+SOD did not expand the coordination sphere of the Mn3+ ion.22 Subsequent variable-temperature spectroscopic studies of N3-Mn3+SODs from E. coli and T. thermophilus23,24 indicated that these species are thermochromic, exhibiting distinct absorption and CD spectra above and below an organism-specific transition temperature (225 and 305 K for the E. coli and T. thermophilus enzymes, respectively).23,24 Due to the fact that (i) the low temperature (LT) species exhibits a “red-band” in the absorption spectrum at ∼11 000 cm-1, characteristic of sixcoordinate Mn3+ systems,26 whereas the lowest-energy band for the high temperature form occurs at ∼16 500 cm-1, and (ii) the ∆H for the thermochromic conversion (∆H ≈ 5 kcal/mol) is in the range of bond dissociation enthalpies for Mn3+ complexes, it was proposed that the thermochromic event corresponds to a conversion from a six- to a five-coordinate Mn3+ ion as temperature is increased.23,24 Intriguingly, mutagenesis studies of E. coli Mn3+SOD have shown that replacement of the second-sphere Tyr34 residue (Figure 1) with Phe followed by treatment of the corresponding variant with azide results in a species (Y34F N3-Mn3+SOD) that exhibits a temperature-independent absorption spectrum (Figure 2) resembling that of LT N3-Mn3+SOD,27 thus suggesting that azide binding to Y34F Mn3+SOD results in the formation of a sixcoordinate Mn3+ ion at all temperatures. However, the mechanism by which this second-sphere substitution prevents ligand dissociation at RT is currently unclear, and the identity of the ligand that dissociates upon the LT to RT conversion of the native enzyme also remains unknown. It has generally been assumed that azide dissociation is not responsible for the LT to RT conversion of N3-Mn3+SOD because the absorption spectrum of the RT form (Figure 2) appears to exhibit N3- f Mn3+ charge-transfer (CT) transitions in the near-UV,23 though detailed assignments of near-UV spectral features for Mn3+SOD with or without azide have not yet been reported. Consequently, coordinated solvent or Asp167 (16) The E. coli MnSOD amino acid numbering scheme is used in this study unless otherwise indicated. (17) Stallings, W. C.; Metzger, A. L.; Pattridge, K. A.; Fee, J. A.; Ludwig, M. L. Free Radic. Res. Commun. 1991, 12-13, 259-268. (18) Tierney, D. L.; Fee, J. A.; Ludwig, M. L.; Penner-Hahn, J. E. Biochemistry 1995, 34, 1661-1668. (19) Han, W.-G.; Lovell, T.; Noodleman, L. Inorg. Chem. 2002, 41, 205-218. (20) Miller, A. F.; Padmakumar, F.; Sorkin, D.; Karapetian, A.; Vance, C. K. J. Inorg. Biochem. 2003, 93, 71-83. (21) Pick, M.; Rabani, J.; Yost, F.; Fridovich, I. J. Am. Chem. Soc. 1974, 96, 7329-7333. (22) Whittaker, J.; Whittaker, M. J. Am. Chem. Soc. 1991, 113, 5528-5540. (23) Whittaker, M. M.; Whittaker, J. W. Biochemistry 1996, 35, 6762-6770. (24) Whittaker, M. M.; Whittaker, J. W. J. Biol. Inorg. Chem. 1997, 2, 667671. (25) Un, S.; Dorlet, P.; Voyard, G.; Tabares, L. C.; Cortez, N. J. Am. Chem. Soc. 2001, 123, 10 123-10 124. (26) Lever, A. B. P. Inorganic Electronic Spectroscopy, 2nd ed.; Elsevier: Amsterdam; New York, 1984. (27) Whittaker, M. M.; Whittaker, J. W. Biochemistry 1997, 36, 8923-8931.

Azide-Adduct of Manganese Superoxide Dismutase

Figure 2. Room-temperature absorption spectra of Mn3+SOD (solid line), N3-Mn3+SOD (dotted line), and Y34F N3-Mn3+SOD (dashed line).

have been viewed as the most likely candidates for the dissociating ligand upon conversion from LT to RT N3-Mn3+SOD. Support for a model involving Asp167 dissociation is provided by the partial bleaching of the prominent absorption band at ∼21 000 cm-1 when azide is added to Mn3+SOD at RT (Figure 2). As the d f d transitions of the Mn3+ center, which produce the dominant contribution to this band, appear to gain their intensities through mixing with Asp- f Mn3+ CT transitions,28 dissociation of Asp167 should drastically reduce the intensities of these transitions. However, we have recently shown that deprotonation of Tyr34 (Figure 1) also results in slight optical bleaching,29,30 suggesting that changes in the dielectric environment of the Mn3+ ion, which could be caused by the presence of a nearby, but noncoordinating, N3- or OHmolecule, are sufficient to cause a loss of absorption intensity. Furthermore, while the dissociation of Asp167 would essentially eliminate mixing between Mn3+ d f d and Asp- f Mn3+ CT transitions, the ligation of N3- should yield N3- f Mn3+ CT transitions in the near-UV that would provide an alternative intensity source for d f d transitions. Thus, while several models of the five-coordinate form of N3-Mn3+SOD have been proposed, the rigorous evaluation of these models on the basis of experimental or computational evidence has been lacking. In this study, a variety of spectroscopic tools were used to explore the electronic properties of resting Mn3+SOD, RT N3Mn3+SOD, and Y34F N3-Mn3+SOD, the last of which provides a six-coordinate N3-Mn3+ adduct that can be studied at room temperature. Additionally, DFT and time-dependent DFT (TDDFT) computations were employed to assist in the analysis of these spectroscopic data and to allow for the evaluation of possible models of RT N3-Mn3+SOD on the basis of experimental data. Our results indicate that in RT N3-Mn3+SOD the azide molecule alternates between two binding sites: (i) H-bonding strongly to Tyr34 and remote from the Mn3+ ion, and (ii) directly coordinated to the Mn3+ ion, as in LT N3Mn3+SOD. Possible implications of these results with respect to the catalytic function of MnSOD are discussed. 2. Materials and Methods 2.1 Sample Preparation. The variant protein Y34F MnSOD was generated by introducing the Y34F mutation in the SodA gene through site-directed mutagenesis using the primer 5′-TTGTTTACAAAGGTCTGATG-3′. Wild-type (WT) MnSOD and Y34F MnSOD from E. coli were purified using an improved version31 of the normal (28) Whittaker, M. M.; Ekberg, C. A.; Edwards, R. A.; Baker, E. N.; Jameson, G. B.; Whittaker, J. W. J. Phys. Chem. B 1998, 102, 4668-4677. (29) Jackson, T. A.; Xie, J.; Yikilmaz, E.; Miller, A.-F.; Brunold, T. C. J. Am. Chem. Soc. 2002, 124, 10 833-10 845. (30) Maliekal, J.; Karapetian, A.; Vance, C.; Yikilmaz, E.; Wu, Q.; Jackson, T.; Brunold, T. C.; Spiro, T. G.; Miller, A.-F. J. Am. Chem. Soc. 2002, 124, 15 064-15 075. (31) Karapetian, A.; Miller, A.-F., unpublished results.

ARTICLES procedure.30,32,33 In brief, the ammonium sulfate concentration used in previous procedures was restricted to only 45% saturation and the step that involves heating the crude extract to 60 °C was eliminated. Fractions from all columns were assayed for SOD activity,34 and only active fractions were carried forward, in contrast to the previous practice of assessing the presence of SOD on the basis of SOD protein content as revealed by SDS-PAGE. Moreover, an additional (third) column was used to achieve higher purity (sephadex A-50, using 5 mM potassium phosphate at pH 7.8 both to preequilibrate and develop the column). This modified procedure resulted in purified WT and Y34F MnSODs with specific activities of ∼12 000 units/mg‚min and ∼9000 units/mg‚min, respectively, a drastic improvement over those reported previously (∼7000 units/mg‚min and ∼4800 units/mg‚min for WT and Y34F Mn3+SOD, respectively).30,32,35 Both MnSOD and Y34F MnSOD were isolated predominantly in the oxidized state and were thus not treated with chemical oxidants. The Mn3+SOD species used in this study were ∼0.5-5 mM in 50 mM phosphate buffer (pH 7.0). N3-Mn3+SOD and Y34F N3-Mn3+SOD were generated by incubating the resting oxidized enzyme with excess azide (for electronic spectroscopic experiments) or a concentrated solution of azide to ensure complexation of ∼95-100% Mn3+SOD active sites with azide (for resonance Raman experiments, as determined by the Kd values for these azide adducts: Kd ) 7.2 and 1.6 mM for WT Mn3+SOD and Y34F Mn3+SOD, respectively).22,27 Azide was not added in excess for resonance Raman experiments in order to avoid possible interference between the signals associated with free and coordinated azide. The exact concentrations of Mn3+SOD and azide used for these experiments are given in the appropriate figure captions. 2.2 Absorption, CD, and MCD Experiments. RT absorption spectra were recorded on a Cary 5E spectrophotometer equipped with near-infrared, visible, and ultraviolet detection capabilities. LT absorption as well as variable-temperature CD and MCD spectra were collected on a spectropolarimeter (Jasco J-715) in conjunction with a magnetocryostat (Oxford Instruments SM-4000) capable of fields up to 8 T. Glycerol, used as a glassing agent in LT spectroscopic experiments, had no effect on the RT CD spectra of Mn3+SOD and Y34F N3-Mn3+SOD, indicating that this glassing agent does not interact with the Mn3+ active site. Typically, LT samples were prepared in 55% (v/v) glycerol and 50 mM phosphate buffer (pH 7.0). 2.3 Resonance Raman Experiments. Resonance Raman (RR) spectra were collected upon excitation with either an Ar+ (Coherent I-305) or Kr+ (Coherent I-90) ion laser with incident power in the 5-20 mW range for high-temperature samples and 20-60 mW for LT samples, using a ∼135° backscattering arrangement. Samples were kept in NMR tubes that were placed in an EPR dewar. To obtain data at 0 °C (273 K), the EPR dewar was packed with an ice-water mixture, whereas for 77 K experiments the EPR dewar was filled with liquid N2. The scattered light was dispersed by a triple monochromator (Acton Research, equipped with 300, 1200, and 2400 g/mm gratings) and detected with a back-illuminated CCD camera (Princeton Instruments, 1340 × 100 pixels). Note that while previous RR experiments were hindered by a strong fluorescence background that prevented detection of vibrational modes associated with N3-Mn3+SOD, the modified purification procedure used for this study31 resulted not only in higher specific activities but also drastically lowered fluorescence. A resonance Raman excitation profile was obtained by monitoring the intensity of the peak associated with the symmetric N3- stretching mode, νs(N3), at ∼1340 cm-1 relative to the intensity of the 984 cm-1 mode of sulfate anion that was added to ∼25-50 mM in some samples as an internal standard. For 0 °C experiments, sample degradation was monitored via absorption spectroscopy. In the case of visible laser excitation (676-407 nm), absorption data collected before and after (32) Vance, C. K.; Miller, A.-F. Biochemistry 1998, 37, 5518-5527. (33) Slykhouse, T. O.; Fee, J. A. J. Biol. Chem. 1976, 251, 5472-5477. (34) Nishikimi, M.; Rao, N.; Jagy, K. Biochem. Biophys. Res. Comm. 1972, 46, 846-856. J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004 12479

Jackson et al.

ARTICLES resonance Raman experiments were similar, though the latter were often partially bleached to a varying extent, probably due to photoreduction of the Mn3+ center. Experiments performed with UV excitation, on the other hand, did result in sample degradation, as evidenced by the generation of a dark precipitate in the protein solution. Thus, only data collected using visible excitation are presented. 2.3 Computations. (A) Active-Site Models. The active-site model of resting Mn3+SOD and hypothetical active-site models of RT N3Mn3+SOD were generated using the PDB file of Mn3+SOD from E. coli (1VEW)6 and N3-Mn3+SOD from T. thermophilus (1MNG-L),12 respectively. For all models, the Asp and His ligands (Figure 1) were replaced by acetate and 4-methylimidazole, respectively, and the secondsphere Tyr and Gln residues were modeled by 4-methylphenol and acetamide, respectively. As all first sphere amino acid residues were truncated at the Cβ positions, the CR-Cβ bond was replaced by an H-Cβ bond of length 1.113 Å; thus, the terminal methyl groups of all first sphere ligands used in our active-site models are referred to as the H3Cβ groups. Hypothetical active-site models for the RT form of N3-Mn3+SOD (referred to as N3-Mn3+SOD(-OH), N3-Mn3+SOD(-H2O), N3-Mn3+SODOH(-Asp), N3-Mn3+SODH2O(-Asp), and N3-Mn3+SODOH(-N3)) also included N-methylacetamide to model the backbone amide group between Val168 and Trp169 that appears to donate a H-bond to Asp167.12 For N3-Mn3+SOD(-OH) and N3-Mn3+SOD(-H2O), the axial solvent ligand was positioned such that the oxygen atom of this ligand, Osol, was initially 2.4 Å from the Mn3+ ion. For N3-Mn3+SODOH(-Asp) and N3-Mn3+SODH2O(-Asp) the noncoordinating oxygen atom of the acetate group used to model Asp167, ONC(Asp167), was protonated, and the C-O(Asp167)36 and C-ONC(Asp167) bond lengths were altered to model CdO double and C-O single bonds, respectively. Additionally, the Mn3+ ion was moved along the Mn-O(Asp167) axis in order to increase the Mn-O(Asp167) distance to 2.4 Å. For N3-Mn3+SOD(-OH) and N3-Mn3+SODOH(Asp), the H-bond network was altered by moving a proton initially bound to Gln146 halfway between N(Gln146) and Osol and by positioning the phenolic H+ of Tyr34 in a manner allowing for H-bond donation to N(Gln146). N3-Mn3+SODOH(-N3) was generated by combining the first coordination sphere of resting Mn3+SOD (1VEW) with the second-sphere Tyr34 and Gln146 residues from N3-Mn3+SOD (1MNG-L). The positions of the second-sphere residues were determined by superimposing the Mn-Osol and Mn-O(Asp167) vectors for 1VEW and 1MNG-L. The azide ligand was positioned within H-bonding distance of the Tyr34 residue (∼1.5 Å) and separated ∼3.4 Å from the Mn3+ ion. The initial coordinates of all hypothetical activesite models for RT N3-Mn3+SOD are included as Supporting Information (Tables S2-S6). (B) DFT Geometry Optimizations. Geometry optimizations were performed by fixing the coordinates of the H3Cβ groups of the firstsphere amino acid residues and all atoms of second-sphere residues not involved in the H-bond network (see Figure 1, dotted line). A separate study of LT N3-Mn3+SOD has shown that this approach allows for a realistic degree of conformational flexibility, while suitably accounting for geometric constraints imposed on the active site by the protein matrix.37 Density functional theory (DFT) energy minimizations were performed using the Amsterdam Density Functional (ADF) 2003.01 software package.38-40 These computations were carried out on a cluster (35) Vance, C. K.; Miller, A.-F. J. Am. Chem. Soc. 1998, 120, 461-467. (36) O(Asp167) and ONC(Asp167) refer to the coordinating and noncoordinating oxygen atoms, respectively, of the acetate group used to model Asp167. (37) Jackson, T. A.; Karapetian, A.; Miller, A.-F.; Brunold, T. C., manuscript in preparation. (38) te Velde, G.; Bickelhaupt, F. M.; van Gisbergen, S. J. A.; Guerra, C. F.; Baerends, E. J.; Snijders, J. G.; Ziegler, T. J. Comput. Chem. 2001, 22, 931-967. (39) Guerra, C. F.; Snijders, J. G.; te Velde, G.; Baerends, E. J. Theor. Chem. Acc. 1998, 99, 391-403. (40) ADF2003.01, SCM, Theoretical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands, http://www.scm.com. 12480 J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004

consisting of Intel Xeon processors (Ace Computers) using uncontracted triple-ζ Slater-type orbitals including a single set of polarization functions as the basis set (ADF basis set IV), an integration constant of 3.0 (for geometry optimizations) or 4.0 (for single-point calculations), and the local density approximation of Vosko, Wilk, and Nusair41 with the nonlocal gradient corrections of Becke42 and Perdew.43 Core orbitals were frozen through 1s (O, N, C) and 2p (Mn). All calculations were spin-unrestricted and converged to the S ) 2 spin state. The coordinates of all DFT energy minimized models presented in this study are included as Supporting Information (Tables S1 and S7-S11). (C) Calculation of Ground and Excited-State Properties. The ORCA 2.2 software package, developed by Dr. Frank Neese,44 was used to perform semiempirical INDO/S-CI, single point DFT, and timedependent (TD) DFT calculations on our geometry optimized activesite models. For INDO/S-CI calculations, the ORCA program employs the model of Zerner and co-workers,45,46 the valence shell ionization potentials and Slater-Condon parameters listed by Bacon and Zerner,47 the standard interaction factors fpσpσ ) 1.266 and fpπpπ ) 0.585, and the following spin-orbit coupling constants: ζ3d(Mn) ) 300 cm-1, ζ4p(Mn) ) 334 cm-1, ζ2p(N) ) 76 cm-1, ζ2p(O) ) 150 cm-1. Restricted open-shell Hartree-Fock self-consistent field calculations for our Mn3+SOD model were tightly converged on the quintet ground state which was used as the reference state for configuration interaction (CI) calculations. Transition energies and intensities were calculated by including single excitations within the 79 highest energy doubly occupied MOs (DOMOs), the 4 singly occupied MOs (SOMOs), and the 80 lowest energy virtual MOs for the quintet states. Ground-state spin Hamiltonian parameters were calculated by including single excitations within the 40 highest energy DOMOs, the 4 SOMOs, and the 26 lowest energy virtual MOs, together with double excitations within the highest 35 DOMOs, the SOMOs, and the 7 lowest energy virtual MOs for the quintet states as well as single electron excitations within the 4 SOMOs for the triplet states. In each case, the completeness of the active space was verified by the insensitivity of the calculated zero-field splitting parameters or excited-state energies to the inclusion of additional one- and two-electron excitations. For single-point DFT calculations, the SV(P) (Ahlrichs polarized split valence) basis48,49 with the SV/C auxiliary basis50 were used for all atoms except Mn, where the larger TZVP (Ahlrichs polarized triple-ζ valence) basis49,51 in conjunction with the TZV/J auxiliary basis were employed. These computations employed Becke’s three-parameter hybrid functional for exchange52,53 combined with the Lee-YangParr correlation functional54 (B3LYP/G), using the default 20% Hartree-Fock exchange. The size of the integration grid used for all calculations was 3 (Lebedev 194 points). Isosurface plots of natural orbitals were generated by the gOpenMol program developed by Laaksonen55,56 using an isodensity value of 0.03 b-3. Electronic excitation energies and absorption intensities for all Mn3+ active-site models used in this study were also determined by the TD(41) (42) (43) (44) (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56)

Vosko, S. H.; Wilk, L.; Nusair, M. Can. J. Phys. 1980, 58, 1200. Becke, A. D. J. Chem. Phys. 1986, 84, 4524-4529. Perdew, J. P. Phys. ReV. B 1986, 33, 8822-8824. Neese, F. ORCA - an ab initio, Density Functional and Semiempirical Program Package, Version 2.2, Max Planck Institut fu¨r Strahlenchemie, Mu¨lheim, 2002. Ridley, J.; Zerner, M. C. Theor. Chim. Acta 1973, 32, 111. Zerner, M. C.; Loew, G. H.; Kirchner, R. F.; Mueller-Westerhof, U. T. J. Am. Chem. Soc. 1980, 102, 589. Bacon, A. D.; Zerner, M. C. Theor. Chim. Acta 1979, 53, 21. Scha¨fer, A.; Horn, H.; Ahlrichs, R. J. Chem. Phys. 1992, 97, 2571-2577. Ahlrichs, R., unpublished results. The Ahlrichs auxiliary basis sets were obtained from the TurboMole basis set library under ftp.chemie.uni-karlsruhe.de/pub/cbasen. Weigend, F.; Ha¨ser, M. Theor. Chem. Acc. 1997, 97, 331-340. Scha¨fer, G.; Horn, H.; Ahlrichs, R. J. Chem. Phys. 1994, 100, 5829-5835. Becke, A. D. J. Chem. Phys. 1993, 98, 1372-1377. Becke, A. D. J. Chem. Phys. 1993, 98, 5648-5652. Lee, C.; Yang, W.; Parr, R. G. Phys. ReV. B 1988, 37, 785-789. Laaksonen, L. J. Mol. Graphics 1992, 10, 33-34. Bergman, D. L.; Laaksonen, L.; Laaksonen, A. J. Mol. Graphics & Modelling 1997, 15, 301-306.

Azide-Adduct of Manganese Superoxide Dismutase

Figure 3. 10 K electronic absorption (top), 4.5 K CD (center), and 4.5 K, 7 T MCD (bottom) spectra of Mn3+SOD. Individual Gaussian curves (dotted lines) and their sums (dashed lines) obtained from iterative fits of these data sets are displayed on their respective spectra. Conditions: 3 mM Mn3+SOD in 50 mM phosphate buffer (pH 7.0) and 50% (v/v) glycerol.

DFT57-59 method within the Tamm-Dancoff approximation,60,61 employing the B3LYP/G hybrid functional and the basis sets described above. The resolution of the identity approximation in calculating the Coulomb term62 was used to ensure convergence of all TD-DFT calculations. In each case, 40 excited states were calculated by including all one-electron excitations within an energy window of ( 3 hartrees with respect to the HOMO/LUMO energies. Isosurface plots of electron density difference maps (EDDMs) were generated with the gOpenMol program55,56 using an isodensity value of 0.003 b-3.

3. Results and Analysis

3.1 Electronic Spectroscopy. (A). WT Mn3+SOD. Our LT absorption, CD, and MCD data collected for Mn3+SOD shown in Figure 3 are virtually identical to those published previously.22,23,28,29,63 Using the d f d transition energies determined by polarized single-crystal absorption studies of RT Mn3+SOD as a starting point,28 iterative Gaussian fitting of our absorption, CD, and MCD spectra reveals that at least 13 transitions contribute to the spectral region between 14 000 and 33 000 cm-1 (Figure 3, Table S12). On the basis of their energetic positions and relatively moderate widths, and supported by previous assignments,22,28 bands 2, 4, 7, and 8, which give rise to the broad absorption envelope in the visible region, are assigned as the four spin-allowed Mn3+ d f d transitions expected for this 3d4 ion. As these d f d transitions are formally symmetry forbidden, their moderate intensities likely reflect significant mixing with intense ligand-to-metal charge-transfer (LMCT) transitions that occur at higher energies. Because the d f d transitions are polarized roughly along the MnO(Asp167) axis, it has been proposed that they gain their absorption intensities primarily through mixing with Asp- f (57) Bauernschmitt, R.; Ahlrichs, R. Chem. Phys. Lett. 1996, 256, 454-464. (58) Casida, E. M.; Jamorski, C.; Casida, K. C.; Salahub, D. R. J. Chem. Phys. 1998, 108, 4439-4449. (59) Stratman, R. E.; Scuseria, G. E.; Frisch, M. J. J. Chem. Phys. 1998, 109, 8218-8224. (60) Hirata, S.; Head-Gordon, M. Chem. Phys. Lett. 1999, 302, 375-382. (61) Hirata, S.; Head-Gordon, M. Chem. Phys. Lett. 1999, 314, 291-299. (62) Neese, F.; Olbrich, G. Chem. Phys. Lett. 2002, 362, 170-178. (63) The increase in absorption intensity upon cooling, (20 920 cm-1) ) 850 M-1cm-1 at RT versus (20 920 cm-1) ) 1080 M-1cm-1 at 10 K, is consistent with previous VT absorption studies performed for this system (see ref 22).

ARTICLES

Figure 4. Room-temperature electronic absorption (top), CD (center), and 4.5 K, 7 T MCD (bottom) spectra of Y34F N3-Mn3+SOD. The near-IR region of the MCD spectrum is shown on an expanded scale in the inset. Individual Gaussian curves (dotted lines) and their sums (dashed lines) obtained from iterative fits of these data sets are displayed on their respective spectra. Conditions: 2.0 mM Y34F Mn3+SOD and 180 mM NaN3 in 50 mM phosphate buffer (pH 7.0) and, for MCD data, 55% (v/v) glycerol.

Mn3+ CT transitions.28 Bands 1, 3, 5, 6, and 11, which contribute negligibly to the absorption spectrum and are significantly narrower than the bands associated with the spin-allowed d f d transitions, are assigned as spin-forbidden quintet to triplet d f d (i.e., spin-flip) transitions of the Mn3+ ion. Bands 9, 10, 12, and 13 are generally broader and less intense in the CD and MCD spectra than the d f d bands and are thus assigned as LMCT transitions. The presence of LMCT transitions in the visible to near-UV region is not unusual for a Mn3+ species26 and reflects the moderately high reduction potential of the Mn3+ ion. While in Mn3+SOD the reduction potential of the Mn3+ ion is greatly depressed from its free ion value by the protein matrix,35,64 this redox-tuning is largely achieved by disfavoring the protonation of the solvent ligand upon reduction.65 As the formal reduction of the Mn3+ ion accompanying a LMCT excitation is too rapid to allow for protonation of the solvent ligand, the LMCT energies of Mn3+SOD should indeed reflect a typical Mn3+ species. Overall, our analysis of absorption, CD, and MCD spectra obtained for Mn3+SOD is consistent with, but more complete than, previous analyses that focused mainly on the d f d region. Importantly, our analysis provides both d f d and LMCT transition energies and intensities, providing an excellent framework within which to interpret spectroscopic data of RT N3-Mn3+SOD. (B). Y34F N3-Mn3+SOD. While the absorption, CD, and MCD spectra of Y34F Mn3+SOD in the absence of azide are virtually identical to those of WT Mn3+SOD (cf., Figures 3 and S1), the addition of azide greatly perturbs the spectroscopic properties of this variant protein at both RT and LT (Figure 4). Most significantly, the LT MCD spectra of Y34F N3-Mn3+SOD exhibits a weak, negative band in the NIR (∼11 000 cm-1) that has no counterpart in the MCD spectrum of Mn3+SOD, where the lowest energy feature is observed at ∼16 000 cm-1. (64) Vance, C. K.; Miller, A. F. Biochemistry 2001, 40, 13 079-13 087. (65) (a) Yikilmaz, E.; Xie, J.; Brunold, T. C.; Miller, A.-F. J. Am. Chem. Soc. 2002, 124, 3482-3483. (b) Schwartz, A. L.; Yikilmaz, E.; Vance, C. K.; Vathyam, S.; Miller, A.-F. J. Inorg. Biochem. 2000, 80, 247-256. J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004 12481

Jackson et al.

ARTICLES

These data suggest that the electronic and, thus, the geometric structuresof Y34F N3-Mn3+SOD are greatly different from those of resting Mn3+SOD. A Gaussian deconvolution of the absorption, CD, and MCD spectra of Y34F N3-Mn3+SOD reveals the presence of at least 11 transitions that span the NIR to the near-UV spectral region (Figure 4). On the basis of their absorption and MCD intensities as well as their moderate widths, bands 1, 2, 6, and 7 are assigned as the four spin-allowed d f d transitions. The drastic red-shifting of the lowest-energy, spin-allowed d f d transition of Y34F N3-Mn3+SOD (band 1) relative to Mn3+SOD (band 2, see Figure 3) is indicative of the formation of a six-coordinate Mn3+ center upon the addition of azide to Y34F Mn3+SOD. Bands 3-5, based on their sharpness and lack of appreciable absorption intensity, are assigned as spin-forbidden, spin-flip transitions of the Mn3+ center, whereas bands 8 through 11, which carry substantial absorption intensity but relatively minor MCD intensity, are assigned as LMCT transitions. By analogy to N3-Fe3+SOD66 and the LT form of N3-Mn3+SOD,37 the ligand most likely involved in these LMCT transitions is azide. As both the spin-allowed and spin-forbidden d f d transition energies observed for Y34F N3-Mn3+SOD are very similar to those determined for LT N3-Mn3+SOD37 (Table S13), the ligand environment of the Mn3+ center must be nearly identical in these two species. Moreover, the LMCT transition energies are also rather similar, although the two transitions responsible for bands 9 and 10 in the LT N3-Mn3+SOD spectra (Table S13) are unresolved in the spectra of Y34F N3-Mn3+SOD. This small difference in the LMCT region between these two species is not surprising considering that the substitution of Tyr34 to Phe eliminates the putative H-bond between the second-sphere Tyr34 residue and the azide ligand (Figure 1, right). The elimination of this H-bond is expected to alter the energies of the azide frontier orbitals and, therefore, the N3- f Mn3+ CT transition energies. Nonetheless, our analysis clearly indicates that Y34F N3-Mn3+SOD and LT N3-Mn3+SOD have very similar spectroscopic properties and thus likely possess virtually identical active-site structures, as proposed previously on the basis of qualitative considerations.27 Consequently, Y34F N3Mn3+SOD provides a suitable model of the six-coordinate form of N3-Mn3+SOD that, however, may be studied at RT. (C) Five-Coordinate N3-Mn3+SOD. RT absorption and CD data collected for N3-Mn3+SOD are shown in Figure 5. The absorption spectrum is strikingly similar to that of resting Mn3+SOD (Figure 3), consisting of a broad feature in the visible region with a shoulder on its lower energy side and more intense features in the near-UV region (∼30 000 cm-1). While previous absorption and CD studies suggested that the broad absorption envelope can be fit with four Gaussian bands corresponding to the four spin-allowed d f d transitions of the Mn3+ ion,23 our analysis of the visible region of the CD spectrum (∼15 00027 000 cm-1) reveals that at least five transitions actually contribute to this region (see Figure 5, bottom and Table S14).67 Collectively, our Gaussian deconvolutions of the RT absorption and CD spectra of N3-Mn3+SOD (Table S14) reveal the (66) Xie, J.; Yikilmaz, E.; Miller, A. F.; Brunold, T. J. Am. Chem. Soc. 2002, 124, 3769-3774. (67) While in a previous study the absorption spectrum of RT N3-Mn3+SOD was fit without a band at ∼21 500 cm-1 (see ref 23), our analysis of the CD spectrum of RT N3-Mn3+SOD clearly requires the presence of a positive Gaussian band at ∼21 500 cm-1 (band 3) to achieve satisfactory agreement between the simulated and experimental spectra. 12482 J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004

Figure 5. Room-temperature electronic absorption (top) and CD (bottom) spectra of N3-Mn3+SOD. Individual Gaussian curves (dotted lines) and their sums (dashed lines) obtained from iterative fits of these data sets are displayed on their respective spectra. Conditions: 2.0 mM Mn3+SOD and 180 mM NaN3 in 50 mM phosphate buffer (pH 7.0).

Figure 6. Experimentally determined energies of d f d excited states for Mn3+SOD, RT N3-Mn3+SOD, and six-coordinate Y34F N3-Mn3+SOD (states are designated using the hole formalism). The coordinate system relates to the D-tensor orientation obtained from INDO/S-CI computations for Mn3+SOD and six-coordinate N3-Mn3+SOD.

presence of four bands of comparable widths in the visible region (bands 1-4) and at least three broader features in the near-UV (bands 5-7). As previously noted by Whittaker and co-workers,23,24 the d f d transition energies of N3-Mn3+SOD at RT are suggestive of a five-coordinate active site, since six-coordinate Mn3+ complexes invariably exhibit one d f d transition in the nearIR region.26 Indeed, inspection of our experimentally derived d f d excited-state correlation diagram for resting Mn3+SOD, RT N3-Mn3+SOD, and the temperature-independent Y34F N3Mn3+SOD azide-adduct (Figure 6) reveals remarkable similarities between the former two but drastic differences between the two azide-treated proteins. The relative energies of the d f d excited states shift by at most ∼300 cm-1 upon the addition of azide to Mn3+SOD at RT, suggesting that the coordination geometry of RT N3-MnSOD is, in fact, remarkably similar to the trigonal bipyramidal geometry of resting Mn3+SOD. Others have reported that azide must coordinate to the Mn3+ ion at RT, as the near-UV region changes upon the addition of azide due to the appearance of N3- f Mn3+ CT transitions in this spectral region.23 However, our Gaussian deconvolution of the RT absorption spectrum of resting Mn3+SOD (Figure S2, Table

Azide-Adduct of Manganese Superoxide Dismutase

ARTICLES Table 1. Peak Positions (cm-1) of the Azide-Related RR Features of Wild-Type (WT) and Y34F N3-Mn3+SOD at 273 K (0 °C) and 77 K (-196° C) N3−MnSOD species

WT (273 K) WT (77 K) Y34F (273 K) Y34F (77 K)

Figure 7. Resonance Raman spectra at 0 °C of N3-Mn3+SOD and Y34F N3-Mn3+SOD using 466 nm (21 459 cm-1) excitation. Inset: Expanded view of the νs(N3) feature for samples prepared with 14N3 (solid line) and 15N-14N (dotted line). Conditions for N -Mn3+SOD: 2.6 mM Mn3+SOD 2 3 and 56 mM NaN3 in 50 mM phosphate buffer (pH 7.0). Conditions for Y34F N3-Mn3+SOD: 4.3 mM Y34F Mn3+SOD and 25 mM NaN3 in 50 mM phosphate buffer (pH 7.0).

S15)68 reveals that the intensities of the near-UV CT transitions of Mn3+SOD (fexp ≈ 11 × 10-3 for band 12) and N3-Mn3+SOD (fexp ≈ 9.1 × 10-3 for band 7) are actually very similar. Collectively, the minor spectral changes in the NIR/visible/nearUV regions that accompany conversion of Mn3+SOD to N3Mn3+SOD at RT suggest that azide does not coordinate directly to the active-site metal ion at RT. 3.2 Resonance Raman Spectroscopy. To further characterize the active site of RT N3-Mn3+SOD, resonance Raman (RR) spectroscopy was utilized as a selective probe of metal-ligand vibrations. For these experiments the N3-Mn3+SOD sample was cooled with ice water (∼0 °C, 273 K) during data collection so as to minimize photoinduced sample degradation while ensuring that the majority of Mn3+ active sites were in the five-coordinate form (at 273 K, Keq ) 6.3 in favor of the RT species23). To interpret RR data obtained for this structurally ill-defined species, parallel experiments were performed on Y34F N3Mn3+SOD, a model of the LT, six-coordinate form of N3-Mn3+SOD that can be studied at RT.69 RR spectra of N3-Mn3+SOD and Y34F N3-Mn3+SOD obtained at 0 °C with 466 nm (21 459 cm-1) excitation reveal the presence of five prominent features (Figure 7 and Table 1). Based on RR studies of N3-Mn3+ porphyrins,70 as many as four features associated with the N3-Mn3+ unit are anticipated in RR spectra of N3-Mn3+ systems; namely, the Mn-N(N2) stretch, ν(Mn-N), at ∼350 cm-1, the N3- bend, δ(N3), at ∼630 cm-1, the symmetric azide stretch, νs(N3), at ∼1320 cm-1, and the antisymmetric azide stretch, νas(N3), at ∼2040 cm-1. By analogy, the RR features of N3-Mn3+SOD and Y34F N3Mn3+SOD that occur at ∼350 and 1340 cm-1 are therefore assigned as ν(Mn-N) and νs(N3), respectively. Consistent with (68) Here, we have used the RT data for Mn3+SOD (Figure S2) rather than the LT data presented in Figure 3 to permit a meaningful comparison of oscillator strengths, which are weakly temperature dependent, between resting Mn3+SOD and RT N3-Mn3+SOD. (69) All RR data were collected using excitation energies below 25 000 cm-1, where LT N3-Mn3+SOD and Y34F N3-Mn3+SOD are spectroscopically indistinguishable (see Table S13). (70) Czernuszewicz, R. S.; Wagner, W.-D.; Ray, G. B.; Nakamoto, K. J. Mol. Struct. 1991, 242, 99-117.

azide isotope

ν(Mn−N3)

νs(N3)

14N 3 15N-14N 2 14N 3 14N 3 15N-14N 2 14N 3

315; 354 ND ND 318; 358 ND ND

1342 1323; 1343 1342 1344 1320; 1342 1343

the assignment of the feature at ∼1340 cm-1 as νs(N3), this peak broadens and shifts to lower energy when terminally labeled 15N-14N2 is used instead of 14N3 (Figure 7, insets and Table 1). While the ν(Mn-N) vibration is also expected to show slight isotopic sensitivity (a shift of ∼3 cm-1),70 the low intensity of the RR feature at ∼350 cm-1 prevented a conclusive analysis of its isotopic behavior. However, as both the energy and width of this low frequency band are sensitive to substitution of Tyr34 (cf. Figure 7, top and bottom), a residue proposed to interact with the azide ligand (Figure 1, right), this feature can safely be attributed to the N3-Mn3+ unit. Similarly, as the width and energy of the band at ∼320 cm-1 are altered by substitution of Tyr34, this band is assigned as a further component of ν(MnN). The presence of more than one feature associated with the Mn-N(N2) stretching mode suggests that this vibration is coupled with other vibrational modes. It is interesting to note that such coupling was also observed for N3-Fe3+SOD, where a total of three closely spaced bands were observed between 355 and 380 cm-1.66 Despite the fact that features at ∼630 or 2040 cm-1, corresponding to δ(N3) or νas(N3), respectively, could potentially also appear in the RR spectra of N3-Mn3+SOD and Y34F N3-Mn3+SOD, no such features were actually observed, suggesting that excited-state distortions along these coordinates are negligible. The remaining RR features at 1002 and ∼1445 cm-1 are characteristic of amino acid vibrations71,72 and are thus not believed to be directly associated with the protein active site. In support of this conclusion, RR excitation profile data show that the band at ∼1002 cm-1 is not resonance enhanced in the visible spectral region (Figure S3). It is interesting to note that all vibrations associated with the N3-Mn3+ moiety are considerably sharper for N3-Mn3+SOD than for Y34F N3-Mn3+SOD (Figure 7). This difference possibly reflects the loss of a H-bond between Tyr34 and azide in the Y34F variant protein. In the WT protein, this H-bond may stabilize the N3-Mn3+ unit in a single configuration, thus producing relatively sharp RR features. Alternatively, in the Y34F variant, the N3-Mn3+ unit may sample more conformational space, giving rise to a sizable frequency distribution for vibrations associated with the N3-Mn3+ moiety and thus broader RR features. RR experiments for N3-Mn3+SOD and Y34F N3-Mn3+SOD were also performed at 77 K. For these experiments, however, laser excitation typically generated a strong fluorescence background, such that no vibrations associated with the protein samples were detectable at excitation wavelengths >413 nm (25 000 cm-1 resulted in sample degradation, thereby preventing RRassisted band assignments for near-UV transitions. Thus, while our RR data clearly demonstrate that excitation into the visible absorption band of N3-Mn3+SOD at 0 °C gives rise to sizable enhancement of vibrational modes associated with the N3-Mn3+ unit, they do not provide evidence for or against the presence of N3- f Mn3+ CT transitions in the near-UV region of the RT N3-Mn3+SOD or Y34F N3-Mn3+SOD absorption spectra. In summary, while absorption and CD data of RT N3-Mn3+SOD suggest that azide does not directly coordinate to the active-site Mn3+ ion, our RR data provide compelling evidence that azide does in fact associate with the Mn3+ center at 273 K. Taken together, these seemingly disparate results do not permit direct identification of the ligand that dissociates upon conversion from LT to RT N3-Mn3+SOD. Nonetheless, these data provide an extensive experimental framework within which this issue can be explored computationally. 12484 J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004

3.3 Computations. (A) Resting Mn3+SOD. An active-site model of resting Mn3+SOD was generated using the X-ray structure coordinates (PDB file 1VEW)6 of the first shell ligands, Tyr34, and Gln146 (Figure 1, left). This model was subsequently refined by performing a constrained DFT energy minimization in which the positions of the H3Cβ groups of the first-sphere amino acid residues as well as all coordinates of second-sphere residues except for the H atoms involved in the H-bond network were held fixed. In this optimization process, all metal-ligand bond distances shortened, particularly the Mn-Osol bond length, which decreased from 2.2 to 1.8 Å. This finding is consistent with previous computational studies by Whittaker et al.28,73 as well as Noodleman and co-workers,19,74 further validating the assumption that differences in experimental and calculated bond distances may be attributed to redox heterogeneity of the MnSOD crystals used for X-ray structure determination.28,74 The shortening of the Mn-Osol and Mn-N(His26) bond lengths leads to a compressed Jahn-Teller axis for this system, greatly destabilizing the Mn 3dz2-derived MO relative to the other four Mn 3d-based MOs and, thus, giving rise to a 3dz2 orbital ground state (using the generally adopted hole formalism). To validate our Mn3+SOD model, INDO/S-CI calculations, shown to be very reliable and reasonably accurate at predicting zero-field splitting parameters for inorganic systems,75 were performed to predict ground state properties that could be compared directly with experimental results. These computations predict D ) 1.46 cm-1 and E/D ) 0.019, in reasonable agreement with EPR data (D ) 2.11 cm-1, E/D ) 0.11).76 Consistent with the distorted trigonal bipyramidal active-site geometry, INDO/S-CI calculations predict that the z-axis of the D-tensor coincides with the Osol-Mn-N(His26) axis, while the y-axis is oriented roughly along the Mn-N(His81) bond vector (Figure 9, bottom inset). Knowledge of the D-tensor orientation for this system is important as it allows transition polarizations obtained from an analysis of VTVH MCD data37 to be related to the molecular structure (vide infra). As a means of further validating our model of Mn3+SOD, INDO/S-CI calculations were performed to predict the corresponding electronic absorption spectrum. To facilitate comparison with spectroscopic data, the computed transition energies and oscillator strengths were used to simulate absorption spectra.77 Comparison of experimental and INDO/S-CI calculated absorption spectra for Mn3+SOD (Figure 9) reveals that the computed d f d transition energies and intensities are fairly accurate, consistent with our previous work on Mn3+SOD.29 However, the INDO/S-CI calculation fails to predict LMCT transitions below ∼39 000 cm-1, in clear disagreement with experimental data (Figure 9). This discrepancy between experimental and predicted transition energies using the INDO/S-CI method most likely reflects the propensity of this computational method to overestimate the ionic contributions to metal-ligand bonding, placing the metal d orbitals too high in energy relative (73) Whittaker, J. Int. J. Quantum Chem. 2002, 90, 15 629-15 635. (74) Li, J.; Fisher, C. L.; Konecny, R.; Lovell, T.; Bashford, D.; Noodleman, L. Inorg. Chem. 1999, 38, 929-939. (75) Neese, F. Curr. Opin. Chem. Biol. 2003, 7, 125-135. (76) Campbell, K. A.; Yikilmaz, E.; Grant, C. V.; Gregor, W.; Miller, A.-F.; Britt, R. D. J. Am. Chem. Soc. 1999, 121, 4714-4715. (77) For these simulations it was assumed that each electronic transition may be described by a Gaussian band with a full width at half-maximum (ν1/2) of 3 000 cm-1, in which case the oscillator strength (fcalc) of the electronic transition is related to its molar extinction coefficient max (M-1cm-1) by fcalc ) 4.61 × 10-9 max ν1/2.

Azide-Adduct of Manganese Superoxide Dismutase

ARTICLES Table 2. Metal-Ligand Bond Lengths for Six-Coordinate N3-Mn3+SOD Active Sites as Determined by X-ray Crystallography (PDB file 1MNG-L) and Partial DFT Energy Minimizations

bond lengths (Å) Mn-N(N2) Mn-Osol Mn-O(Asp167) Mn-N(His171) Mn-N(His26) Mn-N(His81)

1MNG-L

N3−MnSODOH

N3−MnSODH2O

2.22 1.95 2.25 2.12 2.12 2.07

2.38 1.80 2.14 2.10 2.12 2.11

2.00 2.28 2.00 2.18 2.26 2.13

to the ligand orbitals. Thus, while d f d transition energies are accurately predicted by this method, LMCT transition energies, which are particularly sensitive to the energy separation between ligand-centered and metal d-based MOs, are in poor agreement with experimental data. To find a better method for predicting electronic transition energies, time-dependent (TD) DFT computations using the B3LYP hybrid functional were performed on our active-site model of Mn3+SOD.78 While TD-DFT computations have previously been performed on hypothetical active-site models of the superoxide adducts of Mn2+SOD and Fe2+SOD to explore the mechanism of antiferromagnetic coupling between the superoxide ligand and the active-site metal ions,79 to our knowledge, this method has not previously been applied to any Mn3+SOD species for predicting d f d and LMCT transition energies, intensities, and polarizations. TD-DFT calculations were initially performed using our entire Mn3+SOD model; however, nearly identical results were obtained for a model lacking the second-sphere Tyr34 and Gln146 residues. Thus, these residues were omitted in subsequent computations unless otherwise noted.80 TD-DFT computations performed on our active-site model of Mn3+SOD predict four d f d transitions between 17 000 and 24 000 cm-1 to give rise to a broad absorption envelope ( ≈ 1 mM-1cm-1), in excellent agreement with experimental data (Figure 9). In contrast to the results of

INDO/S-CI computations, TD-DFT computations also predict several LMCT transitions near 30 000 cm-1, in much better agreement with experimental data, where the onset of CT transitions is observed at ∼25 000 cm-1. Moreover, TD-DFT computations performed on our model of Mn3+SOD aligned according to the INDO/S-CI derived D-tensor orientation shown in Figure 9 predict the lowest energy d f d transition, which, on the basis of its energetic location and intensity, corresponds to band 2 in the experimental Mn3+SOD absorption spectrum (Figure 3), to be 14% x-, 85% y-, and 1% z-polarized. This prediction is in excellent agreement with the experimental polarization of this band obtained from an analysis of corresponding VTVH MCD data; i.e., 21% x-, 73% y-, and 7% z-polarization.37 As the TD-DFT method is successful in describing the excited-state properties of Mn3+SOD, it is reasonable to assign bands 9, 10, 12, and 13 in the experimental absorption spectrum as His f Mn3+ and Asp- f Mn3+ CT transitions on the basis of these computations. Specifically, TDDFT computations predict two LMCT transitions from the two equatorial His residues to the Mn 3dz2-based LUMO at ∼32 000 cm-1 and Asp- f Mn3+ and His26 f Mn3+ CT transitions at slightly higher energies. Overall, TD-DFT computations prove to be considerably better suited than the INDO/S-CI method for predicting excited-state properties and have been used throughout the remainder of this study. (B) Y34F N3-Mn3+SOD. A detailed spectroscopic and computational analysis of the geometric and electronic structures of the LT form of N3-Mn3+SOD is beyond the scope of this study and will therefore be presented in a separate paper.37 Nonetheless, key results from this analysis will be briefly summarized here to facilitate comparison of active-site models of LT (six-coordinate) and RT (five-coordinate) N3-Mn3+SOD81 as well as to aid in the interpretation of spectroscopic data obtained for Y34F N3-Mn3+SOD. A model of LT N3-Mn3+SOD was generated by using the X-ray structure coordinates (PDB file 1MNG-L) as a starting point for DFT geometry optimizations in which the H3Cβ groups of all first-sphere amino acids as well as all second-sphere residues not involving in the H-bond network were held fixed. Two models of LT N3-Mn3+SOD were considered that possessed either an axial OH- (N3-Mn3+SODOH) or H2O ligand (N3-Mn3+SODH2O); relevant structural parameters for these two

(78) TD-DFT computations performed using the Perdew-Wang pure functional yielded d f d and LMCT transition energies in poor agreement with experimental data. (79) Carrasco, R.; Morgenstern-Badarau, I.; Cano, J. Chem. Commun. 2003, 436-437. (80) In TD-DFT computations performed for putative active-site models of N3Mn3+SOD, the H-bond donation from Tyr34 to the azide ligand, which may be important in determining N3-f Mn3+ CT transition energies, was accounted for by replacing the Tyr34 residue with a H2O molecule oriented in a manner identical to the phenolic -COH group of Tyr34.

(81) Based on published thermodynamic parameters (see ref 23), there is still a fraction of six-coordinate N3-Mn3+SOD present at RT; thus, it may initially seem inappropriate to equate the RT form of N3-Mn3+SOD with fivecoordinate N3-Mn3+SOD. However, the fraction of six-coordinate N3Mn3+SOD present at RT is predicted to be very small (∼8%), as exemplified by the fact that the RT spectroscopic properties of N3-Mn3+SOD are suggestive of a five-coordinate Mn3+ center, thus reflecting the vast majority of five-coordinate active sites present at this temperature. Therefore, our computational models of N3-Mn3+SOD possessing a five-coordinate Mn3+ center are referred to as models of RT N3-Mn3+SOD.

Figure 9. Top: Experimental absorption spectrum of Mn3+SOD at 10 K (solid line) and corresponding Gaussian deconvolution (dotted lines). Center and bottom: Simulated absorption spectra based on INDO/S-CI and TDDFT computations, respectively, performed on our active-site model of Mn3+SOD. The four spin-allowed d f d transitions are indicated by vertical sticks. Inset: Truncated form of the Mn3+SOD active-site model used for the computations. The axis system relates to the D-tensor orientation obtained from INDO/S-CI computations.

J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004 12485

Jackson et al.

ARTICLES

Figure 10. Top: Experimental absorption spectrum of Y34F N3-Mn3+SOD at RT (solid line) and corresponding Gaussian deconvolution (dotted lines). Bottom: TD-DFT computed absorption spectrum for the N3-Mn3+SODOH model. The four spin-allowed Mn3+ d f d transitions are indicated by vertical sticks. Inset: EDDMs of the two highest energy d f d transitions predicted for N3-Mn3+SODOH. Red and gray denote gain and loss of electron density, respectively.

models after geometry optimization are given in Table 2. To estimate the acidity of the axial solvent ligand of six-coordinate N3-Mn3+SOD, a previously employed method was used to compute a ∆pK value between N3-Mn3+SOD and Mn2+SOD, the latter of which was chosen as a reference model because it possesses the same overall charge as N3-Mn3+SOD.37 These computations revealed that the axial solvent ligand of N3-Mn3+SOD is more acidic than that of Mn2+SOD by ∼21 pH units. As previous DFT/electrostatics calculations by Noodleman and co-workers led to the proposal that the absolute pK value of the axial solvent ligand of Mn2+SOD is ∼10-18, our results strongly suggest that N3-Mn3+SOD binds OH- rather than H2O in the axial position at neutral pH. In further support of this hypothesis, the TD-DFT computed absorption spectrum of N3Mn3+SODOH bears a striking resemblance to the experimental spectrum of six-coordinate Y34F N3-Mn3+SOD (Figure 10). In contrast, the TD-DFT computed N3-Mn3+SODH2O absorption spectrum is dramatically different, exhibiting an intense ( ≈ 4 mM-1 cm-1) feature at ∼17 500 cm-1 that has no counterpart in the experimental spectrum. In light of the striking spectral similarities between LT N3Mn3+SOD and Y34F N3-Mn3+SOD (vide supra), it is reasonable to use the TD-DFT results obtained for N3-Mn3+SODOH as a basis for interpreting the spectroscopic properties of Y34F N3-Mn3+SOD. Consequently, band 2 in the experimental spectrum of Y34F N3-Mn3+SOD (Figure 10), corresponding to feature b in the TD-DFT computed spectrum, may be assigned as the Mn 3dxz f 3dz2 transition that is considerably mixed with an N3- f Mn3+ CT transition, thus rationalizing the enhancement of νs(N3) in this region (Figure 8). Bands 6 and 7 are assigned as Mn 3dxy f 3dz2 and 3dyz f 3dz2 transitions, respectively, where both donor MOs contain significant contributions from the N3- π nonbonding (πnb) orbitals. Examination of electron density difference maps (EDDMs) for these transitions (see Figure 10, bottom inset) reveals that in addition 12486 J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004

to a generalredistribution of electron density at the Mn3+ center (reflecting the d f d character of these transitions), the corresponding excited states are stabilized by transfer of electron density from the N3- πnb-based MOs to the Mn3+ 3d-derived donor MOs. Importantly, these EDDMs reveal symmetric loss of electron density from the terminal nitrogens of the N3- ligand, suggesting that excitation in resonance with bands 6 and 7 should give rise to symmetric distortions of this ligand, thus providing an intuitively appealing explanation for the selective RR enhancement of νs(N3) over νas(N3) in this region (Figure 8). TD-DFT computations performed on N3-Mn3+SODOH also reveal that multiple N3- f Mn3+ CT transitions, possessing varying degrees of His f Mn3+ and Asp- f Mn3+ CT character, give rise to bands d and e. Specifically, bands d and e (Figure 10, bottom) contain contributions from N3- f Mn3+ CT transitions where the acceptor MOs are the Mn3+ t2g- and eg-derived MOs (using symmetry labels of the parent Oh point group), respectively, and LMCT transitions that give rise to band e are predicted to be roughly twice as intense as those giving rise to band d. On the basis of these results, bands 8 and 9 of Y34F N3-Mn3+SOD, which exhibit moderate absorption intensities, are assigned as N3- f Mn3+ t2g-type CT transitions, whereas the more intense bands 10 and 11 are attributed to N3f Mn3+ eg-type CT transitions. These assignments are fully consistent with our absorption, CD, and MCD studies of N3Mn3+SOD and Y34F N3-Mn3+SOD (vide supra), which also suggest that all near-UV features above ∼25 000 cm-1 are due to N3- f Mn3+ CT transitions. (C) RT N3-Mn3+SOD. Whittaker and Whittaker have shown that the RT state of N3-Mn3+SOD represents a local enthalpy minimum that is ∼5 kcal/mol above the global minimum defined by the LT form of N3-Mn3+SOD (specifically, ∆H ) 5 and 6.3 kcal/mol for N3-Mn3+SOD from E. coli and T. thermophilus, respectively).23,24 This enthalpy difference provided one criterion that we have used to evaluate hypothetical active-site models of RT N3-Mn3+SOD; i.e., for a model to be deemed reasonable, its total energy must be ∼5 kcal/mol higher than that of our experimentally validated N3-Mn3+SODOH model of LT N3-Mn3+SOD37 described above. The second criterion used to evaluate these five-coordinate N3-Mn3+SOD models was that their TD-DFT computed absorption spectra (in particular the splitting pattern of the d f d excited states) must resemble the experimental spectrum of RT N3-Mn3+SOD. While our spectroscopic and computational data indicate that the axial solvent ligand of LT N3-Mn3+SOD is OH- rather than H2O,37 it may be argued that the RT form of N3-Mn3+SOD, which has a lower coordination number of the Mn3+ ion, possesses an elevated pK. Therefore, in generating hypothetical active-site models for RT N3-Mn3+SOD, we have considered both protonation states of the axial solvent ligand. Models of RT N3-Mn3+SOD were generated by distorting the geometry of the six-coordinate N3-Mn3+SOD active site provided by X-ray structure data to promote either solvent (N3-Mn3+SOD(-OH) and N3-Mn3+SOD(-H2O)), Asp167 (N3-Mn3+SODOH(-Asp) and N3-Mn3+SODH2O(-Asp)), or N3- (N3-Mn3+SODOH(-N3)) dissociation from the Mn3+ center. These models were subsequently refined using constrained DFT energy minimizations (Figure 11, right and Table 3), where the H3Cβ groups of the first-sphere residues as well as all second-sphere

Azide-Adduct of Manganese Superoxide Dismutase

ARTICLES

Figure 11. Experimental absorption spectrum of RT N3-Mn3+SOD (top left) and predicted absorption spectra (left) based on TD-DFT computations performed on the geometry-optimized hypothetical active-site models of RT N3-Mn3+SOD shown on the far right. The structures of these models, shown before and after geometry optimization, have been truncated so as to show only the coordinating N atoms of the His residues and the H atoms involved in the conserved H-bond network. Relevant bond lengths before and after geometry optimization are given. Table 3. Metal-Ligand Bond Lengths for DFT Geometry-Optimized Active-Site Models of RT N3-Mn3+SOD

bond lengths (Å) Mn-N(N2) Mn-Osol Mn-O(Asp167) Mn-N(His171) Mn-N(His26) Mn-N(His81)

N3−MnSOD(−OH)

N3−MnSODOH(−Asp)

N3−MnSOD(−H2O)

N3−MnSODH2O(−Asp)

N3−MnSODOH(−N3)

2.30 1.81 2.11 2.13 2.13 2.12

2.37 1.81 2.14 2.09 2.12 2.10

2.12 2.26 2.14 2.13 2.20 2.10

2.08 2.33 2.07 2.17 2.20 2.13

3.58 1.81 1.98 2.07 2.06 2.04

residues not involved in the H-bond network (Figure 1) were held fixed. Note that N3-Mn3+SOD(-OH), N3-Mn3+SODOH(-Asp), and N3-Mn3+SODOH(-N3) contain the same overall protonation state as N3-Mn3+SODOH, whereas N3-Mn3+SOD(-H2O) and N3-Mn3+SODH2O(-Asp) possess an additional proton. As computed energies can be compared only for systems

that possess the same number and type of atoms, the ∆H values for N3-Mn3+SOD(-H2O) and N3-Mn3+SODH2O(-Asp) have been determined using N3-Mn3+SODH2O (Table 2) as a reference point. N3-Mn3+SOD(-OH) and N3-Mn3+SOD(-H2O) converged to geometries that are 2.5 and 13.3 kcal/mol higher in energy J. AM. CHEM. SOC.

9

VOL. 126, NO. 39, 2004 12487

Jackson et al.

ARTICLES

than N3-Mn3+SODOH and N3-Mn3+SODH2O, respectively. As the relative energy of N3-Mn3+SOD(-H2O) is much greater than that of N3-Mn3+SOD(-OH), the pK of N3-Mn3+SOD(-H2O) must be significantly lower than that of N3-Mn3+SODH2O, and, thus, the axial solvent ligand of N3-Mn3+SOD(H2O) is predicted to be even more acidic (pK < 0). Table 3 and Figure 11 reveal that after geometry optimization, the axial solvent ligand of both the N3-Mn3+SOD(-OH) and the N3Mn3+SOD(-H2O) models reassociates with the Mn3+ ion to generate a six-coordinate active site, suggesting that dissociation of the solvent ligand is not likely to be the source of the thermochromism observed for N3-Mn3+SOD. In support of this conclusion, the corresponding TD-DFT computed absorption spectra exhibit a “red-band” at ∼10 000 cm-1, in clear disagreement with the experimental absorption spectrum of RT N3Mn3+SOD (Figure 11) that lacks any discernible features at energies