Amide Ligation in Nickel Superoxide

Article Cite This: Inorg. Chem. XXXX, XXX, XXX−XXX

pubs.acs.org/IC

The Role of Mixed Amine/Amide Ligation in Nickel Superoxide Dismutase Hsin-Ting Huang,† Stephanie Dillon,‡ Kelly C. Ryan,† Julius O. Campecino,† Olivia E. Watkins,‡ Diane E. Cabelli,§ Thomas C. Brunold,‡ and Michael J. Maroney*,†

Inorg. Chem. Downloaded from pubs.acs.org by UNIV OF SUNDERLAND on 10/03/18. For personal use only.

†

Department of Chemistry, University of Massachusetts at Amherst, 104 Lederle Graduate Research Tower A, 710 North Pleasant Street, Amherst, Massachusetts 01003, United States ‡ Department of Chemistry, University of Wisconsin−Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States § Department of Chemistry, Building 555A, Brookhaven National Laboratory, P.O. Box 5000, Upton, New York 11973, United States S Supporting Information *

ABSTRACT: Superoxide dismutases (SODs) utilize a pingpong mechanism in which a redox-active metal cycles between oxidized and reduced forms that differ by one electron to catalyze the disproportionation of superoxide to dioxygen and hydrogen peroxide. Nickel-dependent SOD (NiSOD) is a unique biological solution for controlling superoxide levels. This enzyme relies on the use of cysteinate ligands to bring the Ni(III/II) redox couple into the range required for catalysis (∼300 mV vs. NHE). The use of cysteine thiolates, which are not found in any other SOD, is a curious choice because of their well-known oxidation by peroxide and dioxygen. The NiSOD active site cysteinate ligands are resistant to oxidation, and prior studies of synthetic and computational models point to the backbone N-donors in the active site (the N-terminal amine and the amide N atom of Cys2) as being involved in stabilizing the cysteines to oxidation. To test the role of the backbone Ndonors, we have constructed a variant of NiSOD wherein an alanine residue was added to the N-terminus (Ala0-NiSOD), effectively altering the amine ligand to an amide. X-ray absorption, electronic absorption, and magnetic circular dichroism (MCD) spectroscopic analyses of as-isolated Ala0-NiSOD coupled with density functional theory (DFT) geometry optimized models that were evaluated on the basis of the spectroscopic data within the framework of DFT and time-dependent DFT computations are consistent with a diamagnetic Ni(II) site with two cysteinate, one His1 amide, and one Cys2 amidate ligands. The variant protein is catalytically inactive, has an altered electronic absorption spectrum associated with the nickel site, and is sensitive to oxidation. Mass spectrometric analysis of the protein exposed to air shows the presence of a mixture of oxidation products, the principal ones being a disulfide, a bis-sulfenate, and a bis-sulfinate derived from the active site cysteine ligands. Details of the electronic structure of the Ni(III) site available from the DFT calculations point to subtle changes in the unpaired spin density on the S-donors as being responsible for the altered sensitivity of Ala0-NiSOD to O2.

■

2O2•− + 2H+ → O2 + H 2O2

INTRODUCTION Superoxide dismutases (SODs) are a family of metalloenzymes that protect aerobic organisms from oxidative stress caused by the superoxide radical anion (O2•−) and reactive oxygen species (ROS) that are derived from O2•−. Specifically, SODs catalyze the disproportionation of superoxide to molecular oxygen (O2) and hydrogen peroxide (H2O2) at rates that are near the diffusion limit. Therefore, SODs are regarded as the first line of defense against this harmful reactive oxygen species.1 During catalysis, the redox-active metal cofactor cycles between oxidized and reduced states that differ by one electron (eqs 1−3). M (n + 1) + O2•− → M n + + O2

(3)

SODs can be categorized into three distinct classes based on amino acid sequence homology: copper- and zinc-containing SOD (CuZnSOD), iron- and manganese-containing SOD (FeSOD/MnSOD), and nickel-containing SOD (NiSOD).1 Much like the Fe- and Mn-SODs, the more recently discovered NiSOD enzyme catalyzes the disproportionation of O2•− by cycling between M(III) and M(II) oxidation states. However, NiSOD has a unique active site structure that is necessary to support redox catalysis by nickel. The first two residues of the NiSOD protein are responsible for donating four of the five ligands involved in coordination of the nickel cofactor. In the resting reduced state of NiSOD, the Ni(II) cofactor adopts a

(1)

M n + + O2•− + 2H+ → M (n + 1) + H 2O2

Received: May 30, 2018 (2)

© XXXX American Chemical Society

A

DOI: 10.1021/acs.inorgchem.8b01499 Inorg. Chem. XXXX, XXX, XXX−XXX

Article

Inorganic Chemistry

amide linkage. This small structural change has a large effect on the electronic structure and results in dramatic changes in the redox properties, catalytic activity, and oxygen sensitivity of the active site. These changes are interpreted in the framework of density functional theory (DFT) calculations that were validated on the basis of spectroscopic data.

four-coordinate, planar NiN2S2 geometry involving coordination by the N-terminal amine, a deprotonated amide from the Cys2 backbone, and the side-chains of Cys2 and Cys6. In the oxidized state of NiSOD, the Ni(III) ion additionally binds the imidazole side chain of His1, which results in a five-coordinate pyramidal geometry around the Ni(III) ion (Figure 1).2−4

■

MATERIALS AND METHODS

The single Ala0-NiSOD variant was prepared by polymerase chain reaction (PCR) of a WT-NiSOD pET-30 Xa/LIC plasmid with mutagenic 5′- and 3′-primers using the QuikChange XL Site-Directed Mutagenesis Kit (Table S1). NovaBlue (Novagen) competent cells were transformed with the PCR product, and transformants were selected by antibiotic resistance. The plasmid containing mutagenic sodN was isolated using a plasmid mini prep kit from Qiagen, and DNA sequencing was performed to confirm the expected base sequence for the mutant (W. M. Keck DNA facility at Yale University, New Haven, CT). BL21 (DE3) pLysS (Novagen) competent cells were transformed with the plasmid containing the desired insertion. Transformants were selected by antibiotic resistance; the mutant gene was expressed, and the variant protein was purified as previously described.17−19 Expression, Purification, Processing, and Reconstitution of Variant Protein. Single colonies were grown overnight at 37 °C with shaking in 10 mL of Luria−Bertani broth supplemented with chloramphenicol (cam) and kanamycin (kan) for selection. These cultures were added to 1 L of prewarmed fresh media and grown to an OD600 of 0.6 and then induced with 0.8 mM isopropyl β-D-1thiogalactopyranoside (IPTG) for 3−5 h. Cells were harvested by centrifugation, resuspended in 40 mL of Ni-NTA binding buffer (10 mM imidazole, 50 mM tris(hydroxymethyl)aminomethane (Tris), 300 mM sodium chloride pH 8.0), and then frozen at −80 °C to lyse the cells. The cell harvests were thawed and treated with 100 μL of DNase I solution (10 mg/mL DNase I, 10 mM magnesium chloride, 20 mM Tris pH 7.5, 40% glycerol) at 37 °C until the viscosity of the solution was significantly reduced. The cell lysate was centrifuged for 5 min at 8,000 rpm, and the supernatant was used for chromatographic protein purification. All chromatographic purifications employed an AKTA-FPLC (Amersham Biosciences). The cell lysate supernatant was loaded onto a column (Pharmacia HR10) containing Ni-NTA HisBind Superflow resin (Novagen) at 3 mL/min with Ni-NTA binding buffer. Once the absorbance at 280 nm returned to the baseline, the buffer was changed to 12% elute buffer (250 mM imidazole, 50 mM Tris, 300 mM sodium chloride pH 8.0) in one step, and the column was washed with 7 column volumes of 12% elute buffer. The fusion protein was then eluted from the column using 70% elute buffer. Electrospray ionization mass spectrometry (ESI-MS) was used to confirm the molecular weight (MW) of the expected fusion protein (vide infra). For imidazole to be removed from the protein sample, the fusion protein was loaded onto a column (Pharmicia HR10) containing QSepharose resin at 3 mL/min with bind buffer (50 mM Tris pH 8.50). Once the absorbance at 280 nm returned to the baseline, the buffer

Figure 1. Active site of NiSOD (PDB ID: 1T6U) in the (a) reduced and (b) oxidized states.

Ligation by an N-terminal amine is rare in metalloproteins and to date has only been observed in NiSOD, HypA, E. coli HypB, RcnR, and albumins, of which NiSOD is the only enzyme.5−9 On the basis of the results obtained from computational and synthetic models, mixed amine/amide ligation in NiSOD has been suggested to stabilize the Ni(III) state, making the oxidation more metal-centered and thus also protecting the S-donors from oxidation by substrate O2•− and products H2O2 and O2.10−14 It has been demonstrated for several NiN2S2 complexes featuring bisamine ligation that the thiolate donors are sensitive to oxidation, whereas bisamide complexes tend to be more stable to sulfur oxidation.10,15 However, recent studies employing synthetic peptide models reached the opposite conclusion, that “all NiSOD mimics are highly unstable, showing sulfur oxidation through the formation of disulfide-bridges between monomeric peptide units and sulfur oxygenation”.16 Further, it was proposed that “the ability to degrade superoxide by the model peptides is coupled to an increased sensitivity for oxidation, which results from the asymmetric amine/amide Ni(II) coordination”. These trends have not been studied directly in the NiSOD enzyme. In this study, we have produced an insertion mutation involving addition of an alanine residue at the N-terminus adjacent to the His1 residue (Ala0-NiSOD), which extends the N-terminus by one residue. Our results suggest that extension of the N-terminus causes changes in the structure of the reduced enzyme consistent with coordination of the His1 Table 1. Characterization of WT- and Ala0-NiSOD NiSOD sample

MW (ESI-MS) (calcd value)

quaternary structure (SEC)

Ni/ subunit

WT

18171.2a (18169.6)

hexamer

0.88

Ala0

anaerobic: 13271.6 (13272.0) air oxidized: 13269.9 (disulfide) 13303.1 (Ala0NiSOD + 2 O = 13304), 13335.4 (Ala0-NiSOD + 4 O = 13336)

hexamer

0.89

electron paramagnetic resonance as isolated g = 2.30, 2.23, 2.01; Azz = 24.9G silent

kinetics kcat @ pH 7.5 (×109 M−1 s−1)

Tm (°C)

0.7

84.8,b 65.2c

0.001

73.8b

ref 19 this work

a

Indicates fusion protein. bHeat capacity peak maximum. cHeat capacity shoulder. B


Article

Inorganic Chemistry was changed to 40% elute buffer (50 mM Tris, 500 mM sodium chloride pH 8.5) in one step, and the fusion protein was then eluted from the column and collected. The eluted sample buffer was then adjusted to the conditions used for N-terminal processing using factor Xa (5 mM CaCl2, 50 mM Tris, 100 mM NaCl, pH 8.0). The purified fusion protein was cleaved to yield NiSOD with Ala0 at the Nterminus using factor Xa. The concentration of fusion protein in factor Xa cleavage buffer was determined using a bicinchoninic acid (BCA) assay. The assay was performed using the Enhanced Test Tube Protocol outlined in Pierce’s BCA Protein Assay Kit instruction manual. Factor Xa was then added (30 μL of 2.1 U/μL factor Xa/50 mg of fusion protein), and the mixture was incubated at 4 °C. The extent of cleavage was monitored using 16% SDS-PAGE. Processed protein was reduced with a 5-fold excess of dithiothreitol (DTT) and reconstituted with a 3-fold excess of NiCl2 in an anaerobic glovebox (Coy Laboratory Products, Inc.). The hexameric quaternary structure of the holo-protein was confirmed by size exclusion chromatography (Figure S1). Mass Spectrometry. To determine the molecular weight (MW) of the NiSOD expression products, ESI-MS analysis was performed using a Bruker Esquire mass spectrometer equipped with an HPHPLC. Analysis was performed under protein denaturing conditions, 0.1% formic acid, and skimmer voltage set to 120 V, and samples were prepared to a final concentration of 20 μM monomer (Table 1). Mass spectrometry was also used to identify the species produced by air-oxidation of Ala0-NiSOD, which resulted in a mixture of Cys-S oxygenates and a disulfide. The presence of a disulfide bond in a portion of the air-oxidized Ala0-NiSOD sample mixture was identified by ESI-MS analysis before and after incubating the sample with 5 mM tris(2-carboxyethyl)phosphine (TCEP) for 3 h. ABI-SCIEX QSTARXL equipped with a turbo-spray ESI source in sequence with an Agilent 1100 HPLC system with a BioBasic-8 (Thermo Scientific) column was used, and the protein was eluted with a 0.1% formic acid/ 5% aqueous acetonitrile solution. As a control, ESI-MS was also used on a similarly treated sample of WT-NiSOD, which showed no evidence of S-oxygenation (see Figure S2). For the positions of oxygen modifications to be identified, peptide analysis employing a Waters Synapt G2-Si mass spectrometer equipped with a nanoACQUITY UPLC system was used. The sample was run through a 2.1 × 30 mm pepsin-immobilized Enymate column (Waters), and the peptides were separated using an analytical HSS T3 column (Waters) before injection into the mass spectrometer. The data were analyzed using Waters ProteinLynx Global Server 3.0.1 (PLGS). Metal Analysis. A PerkinElmer Optima 4300 DV inductively coupled plasma-optical emission spectrometer (ICP-OES) was used to quantify the nickel content of the reconstituted variant protein. This instrument is equipped with a 40 MHz free-running generator and a segmented-array charge-coupled device (SCD) detector. The sample introduction system consisted of a concentric nebulizer with a cyclonic spray chamber. The concentration of nickel in the sample was determined at λ = 231.604 nm. Quaternary Structure and Stability. For the oligomeric state of the proteins to be determined, size-exclusion chromatography was performed using a Superdex 75 15/300 GL (GE Life Sciences) column. The column was standardized with albumin (67 kDa), ovalbumin (43 kDa), chymotrypsinogen (25 kDa), and ribonuclease A (13.7 kDa). A standard curve was constructed by plotting Ve/V0 vs MW (where Ve = elution volume of the peak, and V0 = elution volume of blue dextran), and the data were fit with a second order polynomial. WT- and Ala0-NiSOD were injected onto the column at a concentration of 500 μM (monomer) in 20 mM Tris, 100 mM NaCl pH 8.0, and the retention volume of the peak(s) (Ve/V0) from the chromatogram was analyzed using the standard curve to determine the molecular weight of the eluted protein (see Figure S1). Melting temperature (Tm) was measured with a Microcal VP-DSC with a 0.5 mL sample and reference cells. The sample was concentrated to 100−140 μM monomer in 50 mM Tris buffer. The protein sample and blank (50 mM Tris buffer) were degassed under a vacuum for 10 min and syringed into the cells using a pulsing motion

to expel air bubbles. The protein sample was run at 30 psi over 25− 100 °C at a scan rate of 30 °C/hour. Baseline correction and normalization were performed with the Microcal interface to the Origin graphing program, and the Tm was taken to be the peak maximum of the thermogram. Kinetics. Enzyme kinetics were studied using pulse radiolytic generation of O2•− and were carried out using a 2 MeV van de Graaff accelerator at Brookhaven National Laboratory. Superoxide radicals were generated upon pulse radiolysis of an aqueous, air/O2-saturated solution containing 10 mM phosphate, 30 mM formate, and 5 μM ethylenediaminetetraacetic acid (EDTA) (eqs 6−10). Catalytic rate constants were obtained by monitoring the disappearance of O2•− OH• + HCO2 ‐ → H 2O + CO2•−

(6)

CO2•− + O2 CO2 + O2•−

(7)

eaq − + O2 → O2•−

(8)

H• + O2 → HO2

(9)

HO2 F O2•− + H+

(10)

at 260 nm in the presence of micromolar concentrations of WT- and Ala0-NiSOD. The path length of the quartz cell used was 2.0 cm and a 200−1900 ns pulse width was chosen, resulting in the generation of ∼6−45 μM O2•− per pulse. The disappearance of O2•− was followed in identical pH 8.5 solutions where one contained 10 μM Ala0NiSOD and one only buffer. The bimolecular rate for the disappearance of O2•− was 3 × 104 M−1 s−1. In the presence of enzyme, there was some evidence for a first order rate of 1 s−1 leading to a possible bimolecular rate for catalytic activity of the enzyme of ≤1 × 105 M−1 s−1. This bimolecular rate was calculated using the nickel concentration with the assumption that all nickel ions are specifically bound and equally contributing to SOD activity. Electronic Spectroscopy. Variable-temperature absorption and magnetic circular dichroism (MCD) spectra of Ala0-NiSOD were obtained on a Jasco J-715 spectropolarimeter in conjunction with an Oxford Instruments SM-4000 8T magnetocryostat. To eliminate contributions from the CD background and glass strain to the MCD signal, the difference between MCD data collected with the magnetic field aligned parallel and antiparallel to the light propagation axis was taken. The final protein concentration used for samples of Ala0NiSOD prepared in an anaerobic environment was 0.655 mM, and the final concentration of Ala0-NiSOD after long-term exposure to air was 0.720 mM. All samples contained at least 55% (v/v) of the glassing agent glycerol. X-ray Absorption Spectroscopy. X-ray absorption spectroscopy (XAS) data collection and analysis were performed as previously described.20 Nickel K-edge XAS data were collected on beamline X9B at the National Synchrotron Light Source (Brookhaven National Laboratory). Samples of frozen protein solutions (1 mM, based on nickel content, in 20 mM Tris·HCl, pH 8.0) were placed in polycarbonate holders inserted into aluminum blocks and held near 50 K using a He displex cryostat. The ring conditions for data collection were 2.8 GeV and 120−300 mA. A sagittally focusing Si(111) double crystal monochromator and a 13-channel Ge fluorescence detector (Canberra) were used for data collection. Xray absorption near-edge spectroscopy (XANES) data were collected from ±200 eV relative to the nickel K-edge. The edge energy reported was taken to be the maximum of the first-derivative of the XANES spectrum. Extended X-ray absorption fine structure (EXAFS) data were collected to 9307 eV (k = 16 Å−1). The X-ray energy for the Kedge of nickel was internally calibrated to 8331.6 eV using transmission data from a nickel foil. The data shown are the average of 6 scans. Sixpack software was used to remove bad detector signals, calibrate the edge energy of Ni foil, and average the data. Athena software was used for data reduction and normalization. The EXAFS123 software package was used for XANES analysis.21 EXAFS data were analyzed using the Artemis software as previously described.22 Scattering parameters were generated using the FEFF 8 C


Article

Inorganic Chemistry software package23 with the starting point for refinement derived from the atomic coordinates of the geometry-optimized computational model (vide infra). The k3-weighted data were fitted by setting the FT window to 2.0−12.5 Å and fitting from r = 1.0−4.0 Å in r-space with an S0 value of 0.9. Scattering from atoms in the first, second, and third (where noted) coordination shells was accounted for in the fits, including multiple-scattering pathways that contribute more than 15% to the EXAFS. Imidazole ligands were fit as rigid rings with one adjustable Ni−N distance and one value of σ2 for all the atoms in the ring. For different models to be compared to the same data set, ifeffit was used to obtain three goodness of fit parameters, χ2 (eq 11), reduced χ2, and the R-factor (eq 12), where Nidp is the number of independent data points, Nε2 is the number of uncertainties to minimize, Re(f i) is the real part of the EXAFS fitting function, and Im(f i) is the imaginary part of the EXAFS fitting function. Reduced χ2 = χ2/(Nind − Nvarys) (where Nvarys is the number of refining parameters) and represents the degrees of freedom in the fit. χ2 =

Nidp Nε

2

are coordinated to the Ni(II)-ion. For comparison, a DFT geometry optimization was also performed using previously published coordinates of an active site model of a Cys → Ser NiSOD variant (C2S/C6S6‑c)25 using the same functional and basis sets as listed below. All DFT geometry optimizations, as well as single-point and timedependent (TD-DFT) calculations, were performed using Ahlrich’s polarized split and auxiliary basis sets26 (SV(P) and SV/C, respectively) for all atoms except for nickel and the atoms directly coordinating to it, for which Ahlrich’s valence triple-ζ basis set (TZV(P)) was employed.26 Becke’s three-parameter hybrid functional27 in conjunction with the Lee−Yang−Parr functional28 (B3LYP) were employed in all computations. Spin-restricted and -unrestricted formalisms were used for calculations involving diamagnetic, S = 0 Ni(II) (i.e., WTred, Ala0red, and Ala0-depred) and paramagnetic, S = 1 Ni(II) (i.e., 6C−O and 6C−S) or S = 1/2 Ni(III) (i.e., WTox and Ala0ox) complexes, respectively. Isosurface plots of molecular orbitals (MOs), electron density difference maps, and total spin density plots were generated using Pymol29 with isodensity values of 0.05, 0.005, and 0.005 au, respectively. All MO assignments were made by visual inspection.

N

∑ {[Re(fi )]2 + [Im(fi )]2 } i=1

(11)

■

Ifeffit also calculates the R-factor for the fit, which is given by eq 12 and is scaled to the magnitude of the data, making it proportional to χ2. For different models (fits) to be compared, the R-factor and reduced χ2 parameters can be evaluated to determine which model provides the best fit, in which case both parameters should be minimized. Although the R-factor will always improve with an increasing number of adjustable parameters, reduced χ2 will go through a minimum and then increase, indicating that the model is overfitting the data.

RESULTS Protein Characterization. The purified Ala0-NiSOD fusion protein had the MW expected for the single alanine insertion (Table 1). Upon N-terminal processing with Factor Xa to remove the fusion peptide and nickelation, the Ala0NiSOD variant bound nickel with a stoichiometric ratio of ∼1:1 (Table 1). Unlike the native enzyme, which is goldenbrown in color, the Ala0-NiSOD variant was a deep reddishbrown. This color change corresponds to an apparent red-shift of the dominant absorption feature from the UV in the reduced WT enzyme to 457 nm in Ala0-NiSOD (Figure 2).

N

R=

∑i = 1 {[Re(fi )]2 + [Im(fi )]2 } N

∑i = 1 {[Re(χ ̃ data i)]2 + [Im(χ ̃ data i)]2 }

(12)

Computational Methods. All DFT geometry optimizations and single-point calculations were carried out using the ORCA 3.0 program developed by Dr. Frank Neese.24 X-ray crystallographic data of the active site from subunit A of WT-NiSOD (PDB ID: 1T6U)2,3 served as the starting point for the generation of the different computational models investigated. In each case, a constrained geometry optimization was carried out whereby all atoms were allowed to move except for those indicated in boldface type below. All initial structures included the primary coordination sphere residues (His1, Cys2, and Cys6) of NiSOD. The carbonyl group of the amide of Cys2 was transformed to a methyl (-CH3), and Cys6 was truncated at the α carbon and modeled as CH3CH2S−. These atoms made up the active site models of reduced and oxidized WT-NiSOD (WTred and WTox, respectively). For models of the Ala0-NiSOD nickel site, alanine was built as a new fragment using the PyMOL program and was incorporated via a peptide bond to the His1 terminal amine to account for possible amide coordination to the Ni(II) center (Ala0red and Ala0ox). The side chain was first truncated and modeled as a glycine to determine whether the terminal amine could also coordinate to the nickel and form a six-coordinate complex. Geometry optimizations of this model led to the dissociation of the terminal amine from the nickel. For this reason, alanine was truncated at the αcarbon and methylated (-COCH2CH3) in all subsequent geometry optimizations of Ala0-NiSOD active site models. For the protonation state of the His1 backbone amide to be explored, an additional model, denoted as Ala0-dep red, was constructed using the same initial coordinates as for the other models except that the proton was removed from the His1 backbone N atom. Six-coordinate, high-spin Ni(II) models were also generated to evaluate possible Cys thiolate oxidation products. In these models, Ala0 was modeled as a glycine by removing the -CH3 side chain and capping the α-carbon with a hydrogen atom. The terminal amine was also retained. The thiolates from Cys2 and Cys6 were modified to represent sulfinates by the addition of two oxygen atoms to each sulfur. Two models were considered, one in which one of the sulfinate oxygens (6C−O) and another in which the sulfinate sulfurs (6C−S)

Figure 2. Room temperature UV−vis of anaerobic and air-exposed Ala0-NiSOD.

The variant was unstable during storage in air; over a period of ∼2 weeks, the 457 nm peak intensity slowly decreased, ultimately resulting in a light tan-colored sample (Figure 2), suggesting oxidation of the Ni site. This change in the electronic absorption spectrum was accompanied by a change in the ESI mass spectrum that is consistent with the formation of a mixture of protein products that feature oxidation of the two cysteinate ligands to form a disulfide as well as the addition of oxygen atoms to form oxygenates (Figure 3a and Figure S3). The mass associated with the major peak of the airexposed Ala0-NiSOD mass spectrum is 13269.9 Da, which was 2 Da less than the expected unmodified mass. When 5 mM D


Article

Inorganic Chemistry

Figure 3. Mass spectral data and analysis for Ala0-NiSODox. (a) ESI-MS spectrum of Ala0-NiSODox before the addition of 5 mM TCEP. (b) Overlaid spectra of Ala0-NiSODox before and after the addition of 5 mM TCEP.

spectrum correspond to Ala0-NiSOD oxygenates with +2 O and +4 O atoms. Features that might correspond to trace amounts of other products (e.g., +1 O and +3 O) are comparable to the noise level of the data. There is no evidence

TCEP was added to the solution, the mass associated with the major peak increased to 13271.6 Da, indicating that the mass difference was caused by the presence of a disulfide bond (Figure 3b). Two other major peaks in the deconvoluted E


Article

Inorganic Chemistry for unmodified Ala0-NiSOD, and thus, the air-exposed sample is 100% oxidized but features a mixture of several protein species. On the basis of the chemistry of model Ni thiolate complexes, the sulfoxygenates likely represent sulfenate and/ or sulfinate complexes.30−32 Pepsin digestion shows that the oxygen incorporation is associated with peptides that contain Cys2 and Cys6 (Figures S2). MS/MS analyses of N-terminal peptides containing the Cys residues show modifications of both Cys residues by one O atom in the +2O product, consistent with a bis-sulfenate, and modification of both Cys residues by 2 O atoms in the +4 O product, consistent with a bis-sulfinate. Other oxidation products were not present in amounts sufficient for MS/MS analysis. Size-exclusion chromatography of as-isolated Ala0-NiSOD gives rise to an elution peak with a retention volume corresponding to the molecular weight appropriate for a homohexamer (Figure S1). The protein thermal stability was minimally altered by the Ala0 N-terminal extension, as revealed by the DSC data (Figure 4), melting temperatures (Tm, Table

Figure 5. Disappearance of O2− monitored at 260 nm in the presence and absence of 10 μM Ala0-NiSOD. Conditions: pulse-radiolytic generation of 30 μM (bottom curves) or 45 μM O2−(top curves) in 10 mM phosphate buffer pH 8.5, 30 mM formate, and 5 μM EDTA.

nickel K-edge XAS spectrum provide information about the nickel center coordination number and geometry.33 The XANES spectrum obtained for Ala0-NiSOD (Figure 6) has a

Figure 4. DSC thermogram of Ala0-NiSOD compared with WTNiSOD.

1), and stability of the hexamer. The DSC thermogram associated with Ala0-NiSOD is similar in general appearance to that of WT-NiSOD but reveals a lower melting temperature (Figure 4, Table 1). The DSC thermograms of Ala0- and WTNiSOD both exhibit an asymmetric curve with a tail to the low temperature side that suggests a complex thermal denaturation of these proteins, perhaps reflecting dissociation of the hexamer followed by unfolding of monomers.19 Kinetics. The rate of superoxide dismutation in the presence of Ala0-NiSOD was determined by monitoring the disappearance of the optical absorbance of pulse radiolytically generated O2•− at 260 nm in the presence and absence of enzyme. The catalytic rate observed for WT-NiSOD at pH 7.5 is 0.71 × 109 M−1 s−1.18 In contrast, in the presence of Ala0NiSOD, virtually no increase in the rate of disappearance of superoxide above the spontaneous rate of superoxide dismutation under these conditions (∼2 × 105 M−1 s−1) is observed (Table 1). At pH 8.5, there is only a trace of catalytic activity of ≤1 × 105 M−1 s−1 (Figure 5). Thus, the activity of NiSOD is essentially abolished by the Ala N-terminal extension. Geometric Structure of the Ala0-NiSOD Ni(II) Site. XAS. Ni K-edge XAS was used to examine the structure of the Ni site in Ala0-NiSOD prepared anaerobically. Features in the X-ray absorption near-edge structure (XANES) region of the

Figure 6. Ni K-edge XANES spectra of anaerobic Ala0-NiSOD (blue) and dithionite-reduced WT-NiSOD (black) from ref 4.

peak near 8331 eV that is associated with a 1s → 3d electronic transition, the area of which (Table 2) is most consistent with the presence of a five-coordinate nickel site. Although a highly distorted four-coordinate structure cannot be completely ruled out, such as a C2v idealized geometry with dihedral angles greater than 0° (the value of a planar complex),34 this is not supported by our computational model (vide infra). The presence of a shoulder near 8338 eV, rather than a resolved maximum, in the Ala0-NiSOD XANES spectrum is associated with the Ni 1s → 4pz electronic transition and is typical of F


Article

Inorganic Chemistry Table 2. Summary of XAS Analyses for Ala0-NiSODa edge energy (eV)

1s → 3d peak area (×102 eV)

geometry

8345.6

5.4(6)

Spy

r (Å)

σ2 (×103 Å−2)

ΔE0 (eV)

red. χ2

% residual

1N 1N 2S

1.83(5) 1.97(9) 2.26(3)

0(8) 0(8) 4(2)

−5(4)

1066

25.7

1N 1N 2S 2C 2C

1.88(3) 2.06(4) 2.26(2) 2.72(9) 2.83(9)

0(4) 0(4) 8(2) 0(12) 0(12)

−3(1)

325

9.3

1N 1N 2S 2C 2C 2O

1.86(2) 2.05(3) 2.24(2) 2.82(3) 2.70(3) 3.31(3)

0(3) 0(3) 7(2) 0(3) 0(3) 0(3)

−4(1)

236

3.9

1.88(3) 2.10(6) 2.22(2) 2.83(5) 2.72(4) 3.30(3) 2.61(8) [3.54(8)] [3.66(8)] [4.82(8)] [4.75(8)] [3.54(8)]

2(3) 2(3) 6(3) 0(3) 0(3) 0(3) 2(9) [2(9)] [2(9)] [2(9)] [2(9)] [2(9)]

−5(2)

208

3.5

N

1N 1N 2S 2C 2C 2O 1N(Im) 1C(Im) 1C(Im) 1C(Im) 1N(Im) 1C(Im) a

Numbers in [ ] correspond to distances to second and third shell imidazole atoms that were not refined independently.

Table 3. Comparison of Bond Distances (in Å) for Optimized Structures of the S = 0, Ni(II)-Bound Models of the WT- and Ala0-NiSOD Active Sites and Experimentally Derived Distances for Reduced WT-NiSOD and Ala0-NiSOD DFT bond Ni−NHis1,ax Ni−NHis1,eq Ni−NCys2 Ni−SCys2 Ni−SCys6

red

WT

5.25 1.99 1.91 2.21 2.25

Ala0

red

4.12 2.13 1.90 2.20 2.27

X-ray structure (1T6U) red

Ala0-dep

Ni(II)SOD

5.24 2.00 1.90 2.27 2.32

3.98 1.88 1.93 2.16 2.20

Ni(II) complexes with a five-coordinate pyramidal geometry.33 This is distinct from dithionite-reduced WT-NiSOD (Figure 6), where the XANES spectrum exhibits a much more resolved maximum near 8338 eV that provides unambiguous evidence for a Ni(II) site with four-coordinate planar geometry,33 which has been confirmed by crystallographic analysis.2,3 The comparison of the XANES spectra for Ala0-NiSOD and reduced WT-NiSOD also reveals a more intense maximum (white line) for Ala0-NiSOD, indicating an increase in the N:S donor ratio relative to that of reduced WT-NiSOD.33 The EXAFS region of the XAS spectrum provides information regarding the metal ligands, including the donoratom types, numbers, and bond distances. Analysis of the EXAFS (Table 3 and Figure 7) was performed by first modeling the primary coordination sphere provided by ligands other than the His1 imidazole. This gave a best fit for N2S2 coordination with two resolved Ni−N distances and a single average Ni−S distance (Table 2). Attempts to split the S-

2

EXAFS WT DT-red ND 1.91(1) 1.91(1) 2.16(4) 2.16(4)

4

Ala0 2.61(8) 2.10(6) 1.88(3) 2.22(2) 2.22(2)

scattering atoms into separate shells gave Ni−S distance differences that were below the resolution of the data (∼0.15 Å) (Table S2). Next, C atoms from chelate rings formed by coordination of backbone N atoms (His1 amide N and Cys2 amide N atoms and by coordination of the Cys2 amide N and thiolate S atoms) were added to the model, and the distances were refined, improving the fit (Table 2). Addition of twothird-coordination sphere O atoms to model the amide carbonyl O atoms that should be ordered by the fivemembered chelate rings also improved the fit. This fourcoordinate fit is compared with the computational model structure in Table 3. For testing whether the His1 imidazole is coordinated, an imidazole ring was added to the planar model. This gave a very modest improvement to the fit over the four-coordinate model and gave a Ni−N(Im) distance (Ni−N = 2.61(8) Å) that is in reasonable agreement with a crystallographically determined Ni−N(Im) distance (2.70 Å) in oxidized WT-NiSOD.3 This G


Article

Inorganic Chemistry

crystallographic coordinates of the oxidized active site of WT NiSOD2 and optimized using DFT. Because our MCD data indicate that Ala0-NiSOD possesses a diamagnetic ground state (vide infra), only S = 0, Ni(II) models were considered. The optimized structures of these models are compared to the analogous active site model of reduced WT-NiSOD (WTred) in Figure 8 and Table 3. From the optimized bond distances for Ala0red and Ala0depred (Table 3), these active site models are not predicted to maintain the five-coordinate geometry initially imposed (chosen on the basis of XAS data, vide supra). Instead, the optimized bond distances for the geometry-optimized WTred model closely match those determined by X-ray crystallography and EXAFS analysis for the reduced WT-NiSOD enzyme (Table 3) and Ala0-NiSOD (note, however, that a weak axial Ni−N interaction as suggested by our XAS data cannot be ruled out on the basis of these computations as the WTred model did not include outer-sphere interactions that may limit the conformational freedom of the His1 side chain). The Ni(II) ion in all reduced models resides in a fourcoordinate, distorted planar ligand environment with coordination from the His1 backbone NHis1,eq, the Cys2 backbone amidate, and the Cys2 and Cys6 thiolates. Below, the novel electronic absorption features exhibited by Ala0-NiSOD will be used to determine which of the two Ala0-NiSOD computational models more faithfully reproduces the Ni(II) coordination environment in this variant. Electronic Structure of Ala0-NiSOD. Spectroscopy. The Ala0 insertion dramatically affects the electronic properties of the nickel center. X-band electron paramagnetic resonance (EPR) spectroscopy was used to characterize the redox state of the nickel center in the as-isolated protein. Resting (asisolated) native and recombinant WT-NiSOD exhibit a rhombic EPR spectrum that arises from a five-coordinate, low-spin (S = 1/2) Ni(III) center with the unpaired electron residing in the Ni dz2-based molecular orbital.18 Anaerobic Ala0-NiSOD and air-oxidized Ala0-NiSODox samples lack this EPR signal; both were EPR-silent (data not shown), indicating that only Ni(II) was present.

Figure 7. Ni K-edge EXAFS spectra of Ala0-NiSOD (blue) and best fit from Table 2 (red).

result is also in agreement with the XANES analysis (vide supra) and raises the possibility that, under the conditions used for XAS data collection, the complex is a low-spin fivecoordinate Ni(II) complex with an apical imidazole ligand. However, a five-coordinate Ni(II) site in Ala0-NiSOD is at odds with the computational structures, all of which are consistent with a four-coordinate planar site and indicate that the imidazole side chain of His1 is not a ligand (vide infra). This finding suggests that, although the Ni(II) center favors a four-coordinate, planar coordination environment, outersphere interactions not accounted for in our computational models may prevent complete dissociation of the His1 side chain by limiting the Ni−N(Im) distance to ∼2.6 Å (note that this weak axial Ni−N interaction is not expected to affect the electronic structure of the Ni(II) site of Ala0-NiSOD). In any event, the XAS data are consistent with coordination of the two amide N-donors and two Cys S-donors predicted by the computational models. Computational Models. To elucidate the Ni coordination environment of Ala0-NiSOD under anaerobic conditions, active site models were generated starting from the published

Figure 8. Geometry optimized active site models of S = 0, Ni(II)-bound WT and Ala0 NiSOD: (a) WTred, (b) Ala0red, and (c) Ala0-depred. Asterisks (*) denote the atoms that were frozen during geometry optimizations. H


Article

Inorganic Chemistry The low temperature electronic absorption spectrum of Ala0-NiSOD under anaerobic conditions (Figure 9, top)

Figure 10. Low temperature absorption (top) and MCD (bottom) spectra of Ala0-NiSOD following long-term exposure to air.

mixture by mass spectrometry (vide supra), the spectra were not analyzed in further detail. DFT Computations. To assess which geometry-optimized Ala0-NiSOD model is most likely representative of the coordination environment of Ala0-NiSOD, TD-DFT was used to compute absorption spectra for Ala0red, Ala0-depred, and WTred. The appearance of an intense feature at 22,300 cm−1 in the experimental Ala0-NiSOD absorption spectrum that is absent in the absorption spectrum of WT-NiSOD served as a marker for evaluating the computational models. The TD-DFT computed absorption spectra for all models are shown in Figure 11. Interestingly, even though all models possess similar ligand environments, the TD-DFT computed absorption spectrum for Ala0red exhibits striking differences from those predicted for WTred and Ala0-depred. The presence of an intense feature for this model, absent in the calculated spectra for WTred and Ala0-depred, is consistent with the experimental absorption spectrum of Ala0-NiSOD prepared under anaerobic conditions (Figure 9). The energy of this feature is overestimated compared to that of the experimental data, which is common for TD-DFT calculated spectra.35 For a direct comparison with the experimental data to be facilitated, the TD-DFT-computed absorption spectra of all models were red-shifted by 10,000 cm−1. The predicted absorption spectrum for Ala0red is dominated by an intense feature at ∼23,400 cm−1 and two higher energy bands at ∼26,600 and 29,700 cm −1 , which have no counterparts in the computed WTred and Ala0-depred spectra. An electron density difference map (EDDM) reveals the LMCT character of the transition associated with the ∼23,400 cm−1 feature (Figure 11). Specifically, this transition involves electronic excitation from the deprotonated amide (N/O) of the Cys2 backbone to the Ni(II) ion. These results indicate

Figure 9. Low temperature absorption (top) and MCD (bottom) spectra of Ala0-NiSOD prepared under anaerobic conditions.

exhibits an intense feature centered at ∼22,300 cm−1 along with several higher energy bands at ∼29,000 and 31,300 cm−1 and a lower energy shoulder at ∼19,000 cm−1. This spectrum is vastly different from that obtained for reduced WTNiSOD,10 most notably in terms of the appearance of the intense band at 22,300 cm−1, indicating that Ala0-NiSOD possesses a unique Ni(II) site. The corresponding MCD spectrum is essentially featureless (Figure 9, bottom). The lack of temperature-dependent features coinciding with the dominant absorption bands provides compelling evidence that Ala0-NiSOD contains a low-spin, diamagnetic (S = 0) Ni(II) center like reduced WT-NiSOD.10 After long-term exposure of Ala0-NiSOD to air, the absorption feature at 22,300 cm−1 is no longer present (Figures 1, 10), indicating a change in the electronic structure of the active site possibly due to thiolate oxidation. Consistent with this hypothesis, the corresponding MCD spectrum (Figure 10, bottom) is distinctly different from that of the oxidized form of WT-NiSOD. 10 Most notable is the appearance of temperature-dependent, negatively signed features located at ∼16,800 and ∼28,000 cm−1, indicating the presence of a paramagnetic species in the mixture of protein oxidation products. The Ala0-NiSODox MCD spectrum closely resembles the MCD spectrum obtained for active site Cys → Ser NiSOD variants, which have no thiolate ligands and feature high-spin, S = 1, Ni(II) centers. This result is consistent with oxidation of Cys2 and Cys6 in Ala0NiSODox.19,25 Because Ala0-NiSODox was shown to be a I


Article

Inorganic Chemistry

occupied molecular orbital (HOMO) is ligand based with large contributions from the S 3p-orbitals of both Cys2 and Cys6 (64.6%) that engage in π-antibonding interactions with the Ni dxz-based orbital (26.2% character). The HOMO-2 is predominantly Ni dz2-based (79.4% 3d-character), whereas the remaining Ni 3d-based MOs have large ligand contributions and are difficult to unambiguously identify. Electronic Structure of Oxidized Ala0-NiSOD. To further examine the consequences of the Ala0 modification, the oneelectron oxidized derivative of the Ala0red model (Ala0ox) was generated computationally and compared with WTox (Figure 13 and Table 3). In agreement with X-ray crystal structures of the active site of oxidized Ni(III)SOD,2,3 both the WTox and Ala0ox models retain a five-coordinate, distorted square pyramidal geometry with additional ligation from the imidazole ring of His1. The most dramatic difference between the Ala0ox and WTox models involves the Ni−NHis1,eq bond distance consistent with the difference in amine versus amide coordination. The DFT-computed MO diagram for Ala0ox is shown in Figure 14. Because of spin-polarization and the consequent mixing between Ni and ligand-based orbitals, only the spindown (β)-MOs are shown. Plots of the β-spin MOs and the unpaired spin density for Ala0ox (Figure 14, inset) reveal that the formally singly occupied MO has significant orbital contributions from Ni (69.5% 3d-orbital character) and SCys2 (12.9% 3p-orbital character), which participate in a πantibonding interaction. This sizable percentage of sulfur orbital character (along with other factors, such as an altered ability of the His1 side chain to coordinate to the Ni(III) ion) could be partially responsible for the instability of Ala0-NiSOD in air, as evidenced by the changes in the experimental absorption spectrum and mass spectrometric characterization (Figures 2 and 3). To higher energy, β-LUMO+1 has significant contributions from the S ligand atoms (23.9% 3p orbital character) as well as Ni (51.3% 3d-orbital character), which engage in σ-antibonding interactions. The remaining Ni 3d-based MOs can be identified at lower energies. Ligand orbitals produce the dominant contributions to the MOs in the energy window between the β-LUMO and the filled Ni 3dbased β-MOs.

Figure 11. TD-DFT-computed absorption spectra for (a) WTred, (b) Ala0red, and (c) Ala0-depred. All spectra were red-shifted by 10,000 cm−1 to facilitate a direct comparison with the experimental absorption spectrum. The EDDM for the key transition of Ala0red is shown at the top of the corresponding spectrum, where blue and white represent areas of loss and gain of electron density.

that replacing a coordinating amine with an amide results in dramatic changes to the interaction between Ni(II) and the Cys2 amidate ligand. Note that for WT-NiSOD, a similar amidate → Ni(II) LMCT transition is predicted at ∼31,100 cm−1. The slightly decreased S 3p orbital contribution from Cys6 (which is trans to the Cys2 backbone amide) to the LUMO of the variant likely contributes to the red-shift of this transition from WT-NiSOD to Ala0-NiSOD. Alternatively, the computed absorption spectrum for Ala0depred only exhibits slight differences from the WTred spectrum (Figure 11). Two features, located at ∼21,500 and 23,100 cm−1, have intensities that are reasonably similar to that of the prominent band at 22,300 cm−1 in the experimental Ala0NiSOD absorption spectrum (Figure 9). However, nearly identical features are present in the computed absorption spectrum for WTred, which is inconsistent with the drastic changes to the experimental absorption spectrum in response to the Ala0 substitution. Thus, we conclude that the Ala0red model more faithfully reproduces the active site environment of the Ala0-NiSOD variant under anaerobic conditions. Bonding Description. Our results suggest that amide coordination to the Ni(II) center by the His1 backbone dramatically changes the electronic structure of the active site, the effects of which are most evident in the appearance of an intense feature in the absorption spectrum of Ala0red that has no counterpart in the WTred spectrum. The results from a DFT single point calculation for Ala0red reveal extensive mixing between the Ni and ligand frontier orbitals (Figure 12). The lowest unoccupied molecular orbital (LUMO) is Ni dx2-y2 based (45.9% 3d-character), consistent with the LUMO composition of WTred and the d-orbital splitting of a four coordinate, square planar Ni(II) complex.10 The highest

■

DISCUSSION

Nickel binding in NiSOD is achieved utilizing only three invariant amino acid residues (His1, Cys2, and Cys6),36 where, in addition to the side chains, the N-terminal amine of His1 and the backbone amidate of Cys2 provide ligands to the metal binding site (Figure 1).2,3 A remarkable characteristic of the NiSOD active site is the presence of metal−thiolate bonds that remain intact during the catalytic cycle despite the propensity of thiolates to oxidize in the presence of reactive oxygen species (e.g., peroxide37,38) or even O2.30−32 This remarkable stability, as well as the Ni-centered rather than S-centered oxidation, have in part been attributed to the presence of mixed amine/amidate N-donors in the active site.10−15 In fact, a structural survey of redox active Ni-containing enzymes, with the exception of methyl coenzyme M reductase, reveals a trend: Ni enzymes featuring amidate ligands, such as acetylcoenzyme synthase39 and Ni superoxide dismutase,2,3 are oxygen-tolerant whereas redox active Ni enzymes that do not have any amidate ligation, such as NiFe-hydrogenase40 and carbon monoxide dehydrogenase,39 are oxygen-sensitive.41 J


Article

Inorganic Chemistry

Figure 12. DFT-calculated MO diagram and isosurface plots of relevant filled (blue/white) and empty (red/white) MOs of Ala0red. Note that the HOMO/LUMO gap is not drawn to scale.

properties and are also stable in air. One example is (Ni(BEAAM))−, which is a planar N2S2 nickel complex with ligand distances similar to those observed for reduced WTNiSOD. In MeCN, this complex features a diamagnetic nickel center, and the S-donors are stable toward oxygen for weeks.12 DFT studies have suggested an explanation for this trend: The S-orbital character in the “redox-active” MO (RAMO) increases as the N-donors are systematically changed from amide to amine: A bis-amide has the lowest, amine-amide (e.g., Ni(BEAAM)−) has intermediate, and bis-amine (e.g., Ni(bmmp-dmed)) has the highest S-orbital character in the RAMO.12,15 Thus, in the bis-amine complexes, the irreversible redox properties are associated with increased sulfur character in the RAMO, leading to S-centered rather than Ni-centered

As suggested by the above observation, the tendency toward S-centered redox reactions in Ni(II) thiolate complexes with an N2S2 ligand donor atom set is modulated by the nature of the N donors. Various NiN2S2 complexes have been synthesized in which the N-donor ligands were varied from bis-amine to amine/amide to bis-amidate, and their redox and air stability were examined.11,15,42−44 In the bis-amine complexes (e.g., Ni(bmmp-dmed), the redox chemistry is irreversible and sulfur centered but the complexes are stable to air oxidation.12 Bis-amide N-donor ligand environments (e.g., Ni(emi))2−) result in reversible Ni-centered redox chemistry; however, the sulfur ligands are susceptible to air-oxidation.12 Model complexes with a mixed amine/amide N-donor ligand environment feature Ni(II) centers that have reversible redox K


Article

Inorganic Chemistry

Conversion of the Cys2 amidate to a secondary amine was explored in a semisynthetic variant of NiSOD, and the resulting Ni active site that was minimally structurally perturbed (XAS) could access the Ni(III) state (EPR) but had only ∼1% of WT-NiSOD catalytic activity.45 Here, we report that alteration of the His1 N-terminal amine by the addition of Ala to the N-terminus also results in a loss of SOD activity as well as the loss of a stable Ni(III) species in the resting enzyme, clearly establishing the role of mixed amine/ amidate ligation in producing enzyme activity. In principle, the addition of Ala to the N-terminus could result in the formation of a Ni center with either amide/ amidate or bis-amidate coordination. XAS was used to show that the Ni center in Ala0-NiSOD is in the correct locus to bind to both Cys residues in the Ni site and that the remaining ligands are N/O-donors consistent with ligation of two backbone N-donors. However, EXAFS data are not suitable to characterize the types of N-donor ligands that are present. To assess whether Ala0-NiSOD features amide/amidate or bisamidate coordination, DFT geometry optimizations were performed to generate models for each alternative, and the resulting models were evaluated on the basis of our spectroscopic data within the framework of DFT and TD-

Figure 13. Geometry-optimized active site models of S = 1/2 Ni(III)bound WT and Ala0 NiSOD: (a) WTox and (b) Ala0ox. Asterisks (*) denote the atoms that were frozen during geometry optimizations.

oxidation.12 The presence of mixed amine/amide N-donors in NiSOD was proposed to promote metal-centered oxidation without activating the S-donors toward oxidation.10

Figure 14. DFT-calculated MO diagram and isosurface plots of relevant filled (blue/white) and empty (red/white) β-MOs of Ala0ox. Inset: unpaired spin density plot, where white and yellow lobes show positive and negative spin density, respectively. Note that the HOMO/LUMO gap is not drawn to scale. L


Article

Inorganic Chemistry DFT computations. The prominent absorption feature at ∼22,300 cm−1 displayed by Ala0-NiSOD under anaerobic conditions is nicely reproduced in the TD-DFT-computed absorption spectrum of Ala0red, leading to the conclusion that this species most likely features an amide/amidate N-donor set rather than a bis-amidate ligand environment. In addition to the loss of SOD activity and altered spectroscopic features, Ala0-NiSOD proved to be unstable in air, leading to oxidation of the protein to a mixture of species containing disulfide and S-oxygenated complexes as revealed by ESI-MS. Similar oxidation products to all of these have been observed for a set of peptide models of the WT active site (albeit using KO2 as an oxidant),16 whereas the WT enzyme does not undergo oxidation in air, suggesting a role for the protein in protecting the thiolates from oxidation. There are also precedents in synthetic model chemistry involving oxidation by O2 of planar, low-spin Ni(II) dithiolates leading to disulfide complexes16,46 and to sulfoxygenates including sulfenate and sulfinate complexes.30,32 The observation that the sulfoxygenate NiSOD products feature only even numbers of O atoms is consistent with a dioxygenase-like mechanism incorporating both oxygen atoms from a single O2 unit,31 and the observation of a bis-sulfenate suggests addition of O2 across the cis-S donors followed by intramolecular O−O bond cleavage, as has been observed in a model system that features cis-thiolate coordination.47 Oxidation of the thiolate S-donor ligands to a disulfide ligand or to sulfoxygenates would be expected to decrease the σ-donor ability of the S-donor ligand and favor a high-spin Ni electronic configuration. Indeed, the MCD spectrum obtained from air-oxidized Ala0-NiSOD shows negatively signed, temperature-dependent features at ∼16,800 and ∼28,000 cm−1 that are remarkably similar to those displayed by the active site Cys → Ser NiSOD single and double variants.25 These Cys variants contain a high-spin (S = 1) Ni(II) ion in a six-coordinate ligand environment with likely coordination from the Ser oxygen(s) and solvent molecule(s) in addition to ligation from the terminal His1 amine, a nitrogen from the His1 imidazole ring, and the deprotonated backbone amide of Cys2.25 The EPR spectra of anaerobic Ala0-NiSOD and the mixture of oxidation products also lack the signature rhombic spectrum of low-spin d7 Ni(III). A comparison of the HOMO of Ala0red and of the β-LUMO of the hypothetical Ni(III)-bound Ala0ox model, which is the formally singly occupied orbital, to the corresponding MOs of WTred and WTox provides further insight into the origin of the unique properties of Ala0 NiSOD (Figure 15). Substitution of an amine with a weakly coordinating amide in the reduced state (Ala0red) causes a relatively insignificant decrease in the amount of Ni 3d-orbital character in the HOMO from ∼28 to 26%. Similarly, the total S 3p orbital contribution to this orbital remains essentially constant, though the individual contribution from SCys2 that coordinates trans to Neq,His1 increases by 5.9% and that from SCys6 decreases by 5.6% (Figure 15a). A comparison of the βLUMOs of the oxidized counterparts shows that the amount of Ni 3d-orbital character in this orbital decreases slightly from ∼73% in WTox to ∼70% in Ala0ox. Although the β-LUMO of WTox contains predominantly Ni dz2 orbital character, this MO has contributions from several Ni 3d orbitals in the case of Ala0ox (Figure 15b). Moreover, the sulfur 3p orbital contribution to the β-LUMO increases from ∼8% in WTox to ∼13% in Ala0ox and, in the latter model, this “hole

Figure 15. DFT-computed compositions and isosurface plots of the (a) HOMO of WTred and Ala0red and (b) β-LUMO of WTox and Ala0ox.

character” (and thus the unpaired spin density in the occupied spin-up counterpart) is localized entirely on SCys2. This increase in SCys2 3p orbital contribution to the β-LUMO likely explains why Ala0-NiSOD is unstable in air despite the similarities in the MO compositions of WTred and Ala0red. This conclusion is supported by a previous spectroscopic and computational study of a one-electron oxidized derivative of a NiN2S2 complex possessing a bis-amidate N-donor set.10 On the basis of DFT calculations for model systems of NiSOD, it has been proposed that protonation of the metalbound cysteine thiolates48 and the electronic preference for metal-centered redox chemistry10 protect the S-donors from oxidation by the O2•− substrate or the O2/H2O2 products. The results obtained in the present study provide compelling evidence that the mixed amine/amide ligation in NiSOD plays an important role by adjusting the redox properties of the active site to optimize Ni center reduction/oxidation while protecting the S-donors from oxidation by the O2•− substrate or the O2/H2O2 products. The role of the His imidazole ligand in Ala0-NiSOD is less clear. The DFT studies support an unligated, His-off imidazole, whereas the best EXAFS fits are obtained when it is included as a ligand, though the His-off model also provided a good fit. Five-coordinate Ni(II) complexes can adopt either square pyramidal or trigonal bipyramidal geometries, and the nickel ion can be either highspin (S = 1) or low-spin (S = 0).49 It has been demonstrated that low-spin complexes are typically formed utilizing low electronegative donor ligands such as P, As, and C, and highspin complexes are formed utilizing high electronegative donor ligands such as O and N with S-donor ligands found in both low- and high-spin complexes.50 A similar five-coordinate Hison complex was observed as an intermediate in an inner-sphere reaction pathway for NiSOD addressed using DFT to explore the potential energy surface.51 However, in this case, the complex was calculated to be high-spin, and the dismutation reaction was proposed to involve protonation of one Cys ligand.

■

CONCLUSIONS From an analysis of our spectroscopic data within the framework of DFT and TD-DFT computations, we propose that the Ala0 insertion to the N-terminus of NiSOD results in a M


Inorganic Chemistry

■

change from the amine/amidate coordination of the Ni(II) ion to an amide/amidate N-donor environment. Substitution of the N-terminal amine in the WT enzyme with an amide in the Ala0-NiSOD variant promotes Cys-thiolate oxidation, highlighting the necessity of the mixed amine/amidate ligand environment in the WT enzyme to promote Ni-based and reversible oxidation. The DFT-computed MO compositions for WT-NiSOD and Ala0red reveal small but significant differences in the composition of the HOMO and lead us to propose that the reason why Ala0-NiSOD is unstable in air is the altered electronic structure of its one-electron oxidized form, most importantly, the composition of β-LUMO. In the case of Ala0ox, this orbital contains significantly more sulfur 3p character than its counterpart of WTox, which renders the Ni− S bond more readily susceptible to S-based oxidation.

■

REFERENCES

(1) Sheng, Y.; Abreu, I. A.; Cabelli, D. E.; Maroney, M. J.; Miller, A.F.; Teixeira, M.; Valentine, J. S. Superoxide Dismutases and Superoxide Reductases. Chem. Rev. 2014, 114 (7), 3854−3918. (2) Barondeau, D. P.; Kassmann, C. J.; Bruns, C. K.; Tainer, J. A.; Getzoff, E. D. Nickel superoxide dismutase structure and mechanism. Biochemistry 2004, 43 (25), 8038−8047. (3) Wuerges, J.; Lee, J. W.; Yim, Y. I.; Yim, H. S.; Kang, S. O.; Carugo, K. D. Crystal structure of nickel-containing superoxide dismutase reveals another type of active site. Proc. Natl. Acad. Sci. U. S. A. 2004, 101 (23), 8569−8574. (4) Choudhury, S. B.; Lee, J.-W.; Davidson, G.; Yim, Y.-I.; Bose, K.; Sharma, M. L.; Kang, S.-O.; Cabelli, D. E.; Maroney, M. J. Examination of the Nickel Site Structure and Reaction Mechanism in Streptomyces seoulensis Superoxide Dismutase. Biochemistry 1999, 38 (12), 3744−3752. (5) Bar-Or, D.; Curtis, G.; Rao, N.; Bampos, N.; Lau, E. Characterization of the Co(2+) and Ni(2+) binding amino-acid residues of the N-terminus of human albumin. An insight into the mechanism of a new assay for myocardial ischemia. Eur. J. Biochem. 2001, 268 (1), 42−48. (6) Chung, K. C.; Cao, L.; Dias, A. V.; Pickering, I. J.; George, G. N.; Zamble, D. B. A high-affinity metal-binding peptide from Escherichia coli HypB. J. Am. Chem. Soc. 2008, 130 (43), 14056−14057. (7) Dias, A. V.; Mulvihill, C. M.; Leach, M. R.; Pickering, I. J.; George, G. N.; Zamble, D. B. Structural and biological analysis of the metal sites of Escherichia coli hydrogenase accessory protein HypB. Biochemistry 2008, 47 (46), 11981−91. (8) Higgins, K. A.; Chivers, P. T.; Maroney, M. J. Role of the Nterminus in determining metal-specific responses in the E. coli Ni- and Co-responsive metalloregulator. J. Am. Chem. Soc. 2012, 134 (16), 7081−93. (9) Iwig, J. S.; Leitch, S.; Herbst, R. W.; Maroney, M. J.; Chivers, P. T. Ni(II) and Co(II) sensing by Escherichia coli RcnR. J. Am. Chem. Soc. 2008, 130 (24), 7592−606. (10) Fiedler, A. T.; Bryngelson, P. A.; Maroney, M. J.; Brunold, T. C. Spectroscopic and Computational Studies of Ni Superoxide Dismutase: Electronic Structure Contributions to Enzymatic Function. J. Am. Chem. Soc. 2005, 127, 5449−5462. (11) Mathrubootham, V.; Thomas, J.; Staples, R.; McCraken, J.; Shearer, J.; Hegg, E. L. Bisamidate and mixed amine/amidate NiN2S2 complexes as models for nickel-containing acetyl coenzyme A synthase and superoxide dismutase: an experimental and computational study. Inorg. Chem. 2010, 49 (12), 5393−5406. (12) Shearer, J.; Dehestani, A.; Abanda, F. Probing variable amine/ amide ligation in Ni(II)N2S2 complexes using sulfur K-edge and nickel L-edge X-ray absorption spectroscopies: implications for the active site of nickel superoxide dismutase. Inorg. Chem. 2008, 47 (7), 2649−2660. (13) Shearer, J.; Zhao, N. [Me4N](Ni(II)(BEAAM)): a synthetic model for nickel superoxide dismutase that contains Ni in a mixed amine/amide coordination environment. Inorg. Chem. 2006, 45 (24), 9637−9639. (14) Fiedler, A. T.; Brunold, T. C. Spectroscopic and computational studies of Ni3+ complexes with mixed S/N ligation: implications for the active site of nickel superoxide dismutase. Inorg. Chem. 2007, 46 (21), 8511−8523. (15) Mullins, C. S.; Grapperhaus, C. A.; Kozlowski, P. M. Density functional theory investigations of NiN2S2 reactivity as a function of nitrogen donor type and N-H center dot center dot center dot S hydrogen bonding inspired by nickel-containing superoxide dismutase. JBIC, J. Biol. Inorg. Chem. 2006, 11 (5), 617−625. (16) Tietze, D.; Sartorius, J.; Seth, B. K.; Herr, K.; Heimer, P.; Imhof, D.; Mollenhauer, D.; Buntkowsky, G. New insights into the mechanism of nickel superoxide degradation from studies of model peptides. Sci. Rep. 2017, 7, 17194. (17) Bryngelson, P. A.; Arobo, S. E.; Pinkham, J. L.; Cabelli, D. E.; Maroney, M. J. Expression, reconstitution, and mutation of

ASSOCIATED CONTENT

S Supporting Information *

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.inorgchem.8b01499. Mutagenic primers (Table S1) and first coordination sphere EXAFS fits (Table S2), MW determination by SEC for WT- and Ala0-NiSOD (Figure S1), additional MS data including ESI-MS data for WT-NiSOD (Figure S2), and MS/MS analysis of peptides from the pepsin digest of Ala0-NiSODox (Figure S3) (PDF)

■

Article

AUTHOR INFORMATION

Corresponding Author

*Phone: 413-545-4876. Fax: 413-545-4490. E-mail: [email protected]. ORCID

Thomas C. Brunold: 0000-0001-6516-598X Michael J. Maroney: 0000-0002-5598-3038 Notes

The authors declare no competing financial interest.

■

ACKNOWLEDGMENTS This work was supported by grants from the National Science Foundation (CHE-1111462 to M.J.M) and the National Institutes of Health (GM 64631 to T.C.B.). S.D. acknowledges the NIH (Grant T32-GM008505) for support and the NSF (Grant CHE-0840494) for computational resources. The U.S. Department of Energy, Division of Materials Sciences and Division of Chemical Sciences, supported XAS data collection at the National Synchrotron Light Source (NSLS) at Brookhaven National Laboratory. The National Institutes of Health supports beamline X3B (formerly X9B) at NSLS. Pulse radiolysis studies were carried out at the Accelerator Center for Energy Research at BNL, which is supported by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences & Biosciences under contract DE-SC0012704. Mass spectrometry data were obtained at the University of Massachusetts−Amherst Mass Spectrometry Center. We thank Dr. Steve Eyles for his help with experimental design and data analysis. N


Article

Inorganic Chemistry recombinant Streptomycescoelicolor NiSOD. J. Am. Chem. Soc. 2004, 126 (2), 460−1. (18) Herbst, R. W.; Guce, A.; Bryngelson, P. A.; Higgins, K. A.; Ryan, K. C.; Cabelli, D. E.; Garman, S. C.; Maroney, M. J. Role of conserved tyrosine residues in NiSOD catalysis: a case of convergent evolution. Biochemistry 2009, 48 (15), 3354−69. (19) Ryan, K. C.; Johnson, O. E.; Cabelli, D. E.; Brunold, T. C.; Maroney, M. J. Nickel superoxide dismutase: structural and functional roles of Cys2 and Cys6. JBIC, J. Biol. Inorg. Chem. 2010, 15 (5), 795− 807. (20) Leitch, S.; Bradley, M. J.; Rowe, J. L.; Chivers, P. T.; Maroney, M. J. Nickel-specific response in the transcriptional regulator, Escherichia coli NikR. J. Am. Chem. Soc. 2007, 129 (16), 5085−5095. (21) Padden, K. M.; Krebs, J. F.; MacBeth, C. E.; Scarrow, R. C.; Borovik, A. S. Immobilized metal complexes in porous organic hosts: Development of a material for the selective and reversible binding of nitric oxide. J. Am. Chem. Soc. 2001, 123 (6), 1072−1079. (22) Carr, C. E.; Musiani, F.; Huang, H.-T.; Chivers, P. T.; Ciurli, S.; Maroney, M. J. Glutamate Ligation in the Ni(II)- and Co(II)Responsive Escherichia coli Transcriptional Regulator. Inorg. Chem. 2017, 56, 6459−6476. (23) Ankudinov, A. L.; Ravel, B.; Rehr, J. J.; Conradson, S. D. Realspace multiple-scattering calculation and interpretation of x-rayabsorption near-edge structure. Phys. Rev. B: Condens. Matter Mater. Phys. 1998, 58 (12), 7565−7576. (24) Neese, F. ORCA 2.9: An ab Initio, Density Functional, and Semiempirical program package; Mülheim an der Ruhr, Germany, 2010. (25) Johnson, O. E.; Ryan, K.; Maroney, M. J.; Brunold, T. C. Spectroscopic and Computational Investigation of Three Cys-to-Ser Mutants of Nickel Superoxide Dismutase: Insight into the Roles Played by Cys2 and Cy6 Active-Site Residues. JBIC, J. Biol. Inorg. Chem. 2010, 15, 777−793. (26) Schafer, A.; Horn, H.; Ahlrichs, R. Fully Optimized Contracted Gaussian Basis Sets for Atoms Lithiium to Krypton. J. Chem. Phys. 1992, 97, 2571−2577. (27) Becke, A. D. Density-Functional Thermochemistry. 3. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648−5652. (28) Lee, C.; Yang, W.; Parr, R. G. Development of the ColleSalvetti Correlation-Energy Formula into a Functional of the Electron-Density. Phys. Rev. B: Condens. Matter Mater. Phys. 1988, 37, 785−789. (29) Schrodinger, L. The PyMOL Molecular Graphics System, ver 1.3. (30) Kumar, M.; Colpas, G. J.; Day, R. O.; Maroney, M. J. Ligand oxidation in a nickel thiolate complex: a model for the deactivation of hydrogenase by oxygen. J. Am. Chem. Soc. 1989, 111 (21), 8323−5. (31) Mirza, S. A.; Pressler, M. A.; Kumar, M.; Day, R. O.; Maroney, M. J. Oxidation of nickel thiolate ligands by dioxygen. Inorg. Chem. 1993, 32 (6), 977−87. (32) Grapperhaus, C. A.; Darensbourg, M. Y. Oxygen Capture by Sulfur in Nickel Thiolates. Acc. Chem. Res. 1998, 31 (8), 451−459. (33) Colpas, G. J.; Maroney, M. J.; Bagyinka, C.; Kumar, M.; Willis, W. S.; Suib, S. L.; Baidya, N.; Mascharak, P. K. X-Ray Spectroscopic Studies of Nickel-Complexes, with Application to the Structure of Nickel Sites in Hydrogenases. Inorg. Chem. 1991, 30 (5), 920−928. (34) Xia, J.; Dong, J.; Wang, S.; Scott, R. A.; Lindahl, P. A. EXAFS, EPR, and Electronic Absorption Spectroscopic Study of the.alpha. Metallo Subunit of CO Dehydrogenase from Clostridium thermoaceticum. J. Am. Chem. Soc. 1995, 117 (27), 7065−70. (35) Andruniow, T.; Kozlowski, P. M.; Zgierski, M. Z. Theoretical analysis of electronic absorption spectra of vitamin B12 models. J. Chem. Phys. 2001, 115 (16), 7522−7533. (36) Dupont, C. L.; Neupane, K.; Shearer, J.; Palenik, B. Diversity, function and evolution of genes coding for putative Ni-containing superoxide dismutases. Environ. Microbiol. 2008, 10 (7), 1831−43. (37) Deutsch, E.; Root, M. J.; Nosco, D. L. In Adv. Inorg. Bioinorg. Mech.; Sykes, A. G., Ed.; Academic Press: New York and London, 1982; Chapter 6, Vol. 1, p 269.

(38) Adzamli, I. K.; Deutsch, E. Kinetics and Mechanism of the Reaction of Hydrogen Peroxide with Thiolato Complexes of Cobalt(III) and Chromium(III). Comments on the Nucleophilicity of Coordinated Sulfur and the Stability of Coordinated Sulfenic Acids. Inorg. Chem. 1980, 19, 1366−73. (39) Can, M.; Armstrong, F. A.; Ragsdale, S. W. Structure, Function, and Mechanism of the Nickel Metalloenzymes, CO Dehydrogenase, and Acetyl-CoA Synthase. Chem. Rev. 2014, 114 (8), 4149−4174. (40) Lubitz, W.; Ogata, H.; Rudiger, O.; Reijerse, E. Hydrogenases. Chem. Rev. 2014, 114 (8), 4081−4148. (41) Campecino, J. O.; Maroney, M. J. Reinventing the Wheel: The NiSOD Story. RSC Metallobio Ser. 2017, 10, 170−199. (42) Neupane, K. P.; Shearer, J. The influence of amine/amide versus bisamide coordination in nickel superoxide dismutase. Inorg. Chem. 2006, 45 (26), 10552−10566. (43) Mullins, C. S.; Grapperhaus, C. A.; Frye, B. C.; Wood, L. H.; Hay, A. J.; Buchanan, R. M.; Mashuta, M. S. Synthesis and Sulfur Oxygenation of a (N3S)Ni Complex Related to Nickel-Containing Superoxide Dismutase. Inorg. Chem. 2009, 48 (21), 9974−9976. (44) Hatlevik, O.; Blanksma, M. C.; Mathrubootham, V.; Arif, A. M.; Hegg, E. L. Modeling carbon monoxide dehydrogenase/acetyl-CoA synthase (CODH/ACS): a trinuclear nickel complex employing deprotonated amides and bridging thiolates. JBIC, J. Biol. Inorg. Chem. 2004, 9 (2), 238−246. (45) Campecino, J. O.; Dudycz, L. W.; Tumelty, D.; Berg, V.; Cabelli, D. E.; Maroney, M. J. A Semisynthetic Strategy Leads to Alteration of the Backbone Amidate Ligand in the NiSOD Active Site. J. Am. Chem. Soc. 2015, 137 (28), 9044−9052. (46) Kumar, M.; Day, R. O.; Colpas, G. J.; Maroney, M. J. Ligand Oxidation in a Nickel Thiolate Complex. J. Am. Chem. Soc. 1989, 111 (15), 5974−5976. (47) Grapperhaus, C. A.; Darensbourg, M. Y.; Sumner, L. W.; Russell, D. H. Template Effect for O2 Addition across cis-Sulfur Sites in Nickel Dithiolates. J. Am. Chem. Soc. 1996, 118, 1791−1792. (48) Szilagyi, R. K.; Bryngelson, P. A.; Maroney, M. J.; Hedman, B.; Hodgson, K. O.; Solomon, E. I. S. K-edge X-ray absorption spectroscopic investigation of the Ni-containing superoxide dismutase active site: New structural insight into the mechanism. J. Am. Chem. Soc. 2004, 126 (10), 3018−3019. (49) Busch, D. H.; Barefield, E. K. Low-spin five-coordinate nickel(II) complexes containing a macrocyclic tetradentate amine. Inorg. Chem. 1971, 10 (6), 1216−1220. (50) Sacconi, L.; Bertini, I. Low- and high-spin five-coordinate cobalt(II) and nickel(II) complexes with tris(2diphenylphosphinoethyl)amine. J. Am. Chem. Soc. 1968, 90 (20), 5443−5446. (51) Pelmenschikov, V.; Siegbahn, P. E. M. Nickel superoxide dismutase reaction mechanism studied by hybrid density functional methods. J. Am. Chem. Soc. 2006, 128 (23), 7466−7475.

O


Amide Ligation in Nickel Superoxide

Recommend Documents