Digital Commons @ Trinity Digital Commons @ Trinity

: A disulﬁde-bridged peptide containing two Ni 2 + binding sites based on the nickel superoxide dismutase protein, {Ni 2 (SOD mds )} has been prepared. At physiological pH (7.4), it was found that the metal sites are mononuclear with a square planar NOS 2 coordination environment with the two sulfur-based ligands derived from cysteinate residues, the nitrogen ligand derived from the amide backbone, and a water ligand. Furthermore, S K-edge X-ray absorption spectroscopy indicated that the two cysteinate sulfur atoms ligated to nickel are each protonated. Elevation of the pH to 9.6 results in the deprotonation of the cysteinate sulfur atoms, and yields a binuclear, cysteinate bridged Ni 22 + center with each nickel contained in a distorted square planar geometry. At both pH = 7.4 and 9.6, the nickel sites are moderately air sensitive, yielding intractable oxidation products. However, at pH = 9.6, {Ni 2 (SOD mds )} reacts with O 2 at an ~3.5-fold faster rate than at pH = 7.4. Electronic structure calculations indicate that the reduced reactivity at pH = 7.4 is a result of a reduction in S(3p) character and deactivation of the nucleophilic frontier molecular orbitals upon cysteinate sulfur protonation.


Introduction
Nickel is an essential biological co-factor found at the active-sites of a number of microbial metalloenzymes and proteins (Chart 1) [1][2][3][4][5].Broadly divided into redox active and non-active nickel metalloproteins, it has been recognized that the majority of known redox active nickel containing metalloenzymes contain cysteinate sulfur ligation to nickel.Cysteinate ligation appears necessary to poise nickel-based one-electron redox couples so as to be accessible under physiological conditions.It has also been demonstrated that several redox inactive nickel transport and regulatory proteins also possess cysteinate ligands to Ni 2+ [6][7][8][9][10][11].
An interesting feature of the nickel-thiolate moiety is its ability to support ligand protonation without subsequent protonolysis [12][13][14].This is most often observed in (near) square planar Ni 2+ centers where the nucleophilic HOMO possesses significant S(3pπ) character, which effectively act as S-based lone-pairs.To date, two nickel containing metalloenzymes, nickel iron hydrogenase [NiFe]H 2 ase and nickel containing superoxide dismutase (NiSOD), have been demonstrated to possess at least one cysteinate ligand that becomes protonated, forming a Ni-S(H + )-Cys moiety under physiological conditions [15,16].Concerning the role of the Ni-S(H + )-Cys moiety in biochemical reactions, it has been proposed that these moieties can behave as proton donors/acceptors and sources of formal hydrogen atoms [17][18][19][20][21][22].However, their exact role(s) in biological reactions is currently unknown.
NiSOD is a nickel containing homohexameric metalloenzyme that disproportionates O 2 − into H 2 O 2 and O 2 by cycling between reduced Ni 2+ and oxidized Ni 3+ oxidation states [23][24][25].Each monomer contains a mononuclear nickel site that is coordinated by two cis-cysteinate sulfur atoms As stated above, NiSOD itself has been shown to possess at least one Ni-S(H + )-Cys moiety in its reduced form [15].Although not without controversy [33], we have provided strong evidence based on sulfur K-edge X-ray absorption studies that like NiSOD itself, NiSOD metallopeptide based mimics possess a Ni-S(H + )-Cys moiety at physiological pH as well [18,22].Studies have also suggested that the pK a of the Ni-S(H + )-Cys proton within these active-sites is ~8.5, and can become reversibly deprotonated at high pH (>9.0) [22].
In addition to mimicking the NiSOD active site, derivatives of the NiSOD inspired metallopeptides are capable of not only mimicking NiSOD, but also mimicking the active-site of cobalt containing NHase, and coordinating Cu 2+ [34,35].This inspired us to further derivatize a NiSOD metallopeptide mimic, SOD m1 (SOD m1 = (SOD mds = HCDLP-CGVYDA-PA), in order to generate different metal-site structures.The intent of such studies is not necessarily to generate NiSOD biomimetic metallopeptides, but instead to probe different metal coordination environments within a biologically derived scaffold.Herein, we present work on a nickel metallopeptide, {Ni 2 II (SOD mds )} (SOD mds = (T a CDLP-CGVYDA-PA) 2 , where T a is a 2-mercaptoacetate group).It will be demonstrated that at pH = 7.4 this metallopeptide possesses two mononuclear Ni 2+ sites that support the formation of the Ni-S(H + )-Cys moiety.Furthermore, we will show that the metallopeptide forms a dinuclear cysteinate bridged Ni 2 2+ center upon elevation of the pH to 9.6.Lastly, it will be demonstrated that Ni-S-Cys protonation protects the metallopeptide from oxidative damage by O 2 .The protection of the metallopeptide against oxidative damage will be rationalized in terms of an alteration of the electronic structure of the nickel-site upon protonation rendering the thiolate sulfurs relatively inert towards oxidation.

General Considerations
All manipulations were performed under an N 2 /H 2 (97:3) atmosphere in a COY anaerobic chamber.Fmoc/O t Bu protected amino acids and resins were obtained from Advanced Chemtech (Louisville, KY, USA).All other reagents obtained from commercial suppliers were of the highest purity available, and used as received.Analytical and semi-preparative reverse-phase HPLC were performed using Waters X-Bridge C-18 analytical (4.6 × 150 mm; 5 µm) and semi-preparative (30 × 150 mm; 5 µm) columns on a Waters Deltaprep 600 equipped with a photodiode array detector (detection wavelength set to 254 nm) (Waters Technology Corporation, Milford, MA, USA).Mass spectrometry was performed on either a Bruker Microflex MALDI-TOF mass spectrometer,(Bruker, Billerica, MA, USA) a ThermoFinnegan LCQ Deca XP ESI-MS,(Thermo Instrument Systems, Waltham, MA, USA) or a Waters Micromass 20 ESI-MS operating in positive ion mode (Waters Technology Corporation, Milford, MA, USA).NMR spectra were obtained on a 400 MHz Varian VNMRS NMR spectrometer (Agilent, Santa Clara, CA, USA).All chemical shifts (δ) are referenced to the residual protio solvent peak.Electronic absorption spectra were obtained on either a JASCO J-1500, CARY 50 or CARY 5000 UV-vis-NIR spectrometer (Agilent, Santa Clara, CA, USA).Circular dichroism spectra were obtained on a JASCO J-1500 spectropolarimeter (JASCO Inc, Easton, MD, USA).The simultaneous deconvolutions of the CD and electronic absorption spectra were performed using an in-house-written procedure for Igor Pro version 6 and 8 (Wavemetrics, Lake Oswego, OR, USA).Infrared spectra were collected using a Thermo Nicolet 6700 FTIR spectrometer with a diamond crystal ATR (Thermo Instrument Systems, Waltham, MA, USA).X-band EPR spectra were obtained on a Bruker EMXPlus EPR spectrometer equipped with a closed-cycle He cryostat (Bruker, Billerica, MA, USA).

Preparation of S-triphenylmethyl-thioglycolic Acid
S-triphenylmethyl-thioglycolic acid (T A -trityl) was prepared by a modification of the procedure on Martinage et al. [36].Briefly, 1.5 g (5.77 mmol) of triphenylmethanol was added to a 25 mL TFA solution of thioglycolic acid (400 µL, 5.77 mmol) and stirred for 5 h under argon at room temperature.The TFA was removed under vacuum, and the resulting orange solid was washed three times with toluene followed by three times with hexanes resulting in a white solid, which was analytically pure (1.46 g, 81% yield). 1 H NMR (CDCl 3 , 400 MHz): δ 7.30 (m, 15H), 3.04 (s, 2H).

Preparation of {Ni 2 (SOD mds )}
The peptide SOD mds (T a CDLP-CGVYD-PA) was prepared on an AAPTec Focus XC-2RV peptide synthesizer or by manual solid-phase peptide synthesis using HBTU/HOBt coupling strategies on a 0.12 mmole scale with alanine loaded Wang resin using a five-fold excess of activated protected amino acid.Following the coupling trityl protected Cys1 and removal of the Fmoc group, T A -trityl was coupled to the N-terminus using standard HBTU/HOBt coupling strategies.Global peptide deprotection and peptide cleavage from the resin was performed under N 2 using a cleavage cocktail comprised of 84.5:5:5:5:2.5 (TFA:phenol:water:thioanisole:EDT) over the course of 12 h.Following removal of the cleavage solution by vacuum on a Schlenk line, the resulting glassy product was washed four times with cold freshly distilled diethyl ether.HPLC and MALDI-TOF studies of the crude peptide mixture demonstrated the complete formation of the disulfide bridged dimer.The resulting crude peptide was subsequently purified by preparative HPLC (9:1 water:acetonitrile-4:6 water:acetonitrile over the course of 30 min; rt = 14.6 min) resulting in the pure disulfide bridged dimer (m1S 3 ) 2 (14% yield).MALDI-TOF MS: [SOD mda /Na] + exp: 2473.3 m/z; calcd: 2473.9 m/z.Solutions of SOD mds in 50 mM NEM buffer (pH 7.4 or 9.5) were prepared and 2.0 equiv of NiCl 2 (added from a pH 7.0 50 mM stock solution) per peptide were then added to solutions of the (m1S 3 ) 2 .The number of free thiol groups per peptide was verified using an Ellman's assay compared to the peptide concentration as determined by the combined absorbance of the Y residue and disulfide moiety (combined ε = 1,525 M −1 cm −1 at 278 nm) [37].ESI-MS data were obtained by injecting an air-free solutions of the pH 7.4 or 8.6 metallopeptide into the mass spectrometer using an air-tight syringe (ESI-MS: pH 9.

Determination of Ni-S(H + )-Cys pK a
Solutions of {Ni 2 (SOD mds )} were formed in 50 mM NEM buffer at a pH of 6.5.To these, aliquots of NaOH or HCl (0.5 M) were then added to the solution and the pH measured using an Orion ® micro-pH electrode (ThermoFisher Scientific, Waltham, MA, USA) and the electronic absorption measured following each addition.The resulting pH titration curve was constructed by monitoring the change in absorbance at λ = 320 nm, where the largest difference in absorbance between the two species is observed.

Kinetics of the Air Oxidation of {Ni 2
II (m1S 3 ) 2 } Air was bubbled for five minutes through the circular dichroism (CD) spectrum of anaerobically prepared solutions of 1.0 mM {Ni 2 (SOD mds )} at pH 7.4 or 9.5 (50 mM NEM buffer).A CD spectrum of these solutions was subsequently taken every 10 min for 12 h with the sample continuously exposed to air.The oxidation kinetics was modeled using pseudo-first-order reaction kinetics using KinTek Explorer v 5.2.Second order rate constants are reported per nickel site-one at pH = 9.6 and two at pH = 7.4.From the deconvoluted spectra, kinetic traces are reported in the decay of the {Ni 2 (SOD mds )} starting materials.

Nickel K-Edge X-ray Absorption Spectroscopy
Nickel K-edge X-ray absorption spectroscopic data were collected on the HXMA beamline (wiggler insertion device operating at 1.5 T) at the Canadian Light Source (Saskatoon, SA, Canada).Solutions of {Ni 2 (SOD mds )} (1.0 mM in 1:1 50 mM NEM buffer:glycerol at a pH of 7.4 or 9.5) were injected between Kapton tape windows in aluminum sample holders and quickly frozen in liquid nitrogen.Data were collected at 20 K with sample temperatures maintained using an Oxford liquid He flow cryostat.Light was monochromatized using a Si(220) double crystal monochromator, which was detuned 50% for harmonic rejection, and focused using a Rh mirror.Spectra were obtained in fluorescence mode using a 32-element solid-state Ge detector on both lines with a 3-micron cobalt filter placed between the sample and detector, and spectra were calibrated against the first inflection point of Ni-foil (8333 eV), which was simultaneously recorded with the metallopeptide data.Data were collected in 10 eV steps from 8133-8313 eV (1 s integration time per point), 0.3 eV steps from 8313-8363 eV (3 s integration time per point), 2 eV steps from 8363-8633 eV (5 s integration time per point), and 5 eV steps from 8633 eV-16 k (5 s integration time per point).Total fluorescence counts were maintained under 30 kHz, and a deadtime correction yielded no appreciable change to the data.The reported spectra represent the averaged spectra from five individual data sets.Prior to data averaging, each spectrum and detector channel was individually inspected for data quality.Data were subsequently processed and analyzed as previously reported using in-house written procedures for Igor Pro and FEFF 9.4 (University of Washington, Seattle, WA, USA) [18].

Sulfur K-Edge X-ray Absorption Spectroscopy
Solutions of {Ni 2 (SOD mds )} were prepared at a pH of either 7.4 or 9.6 (~1 mM in 50 mM NEM buffer) and injected into Lucite sample holders with polypropylene windows.Data were obtained at room temperature (~20 • C) on beamline X-19a at the NSLS (Upton, NY, USA) in a He purged sample chamber using a passivated implanted planar silicon (PIPS) detector.The photon energies were calibrated against the first inflection point of S 8 recorded before and after each sample; it was found that there was no detectable monochromator drift throughout the data collection.Data were obtained in 5 eV steps in the pre-edge region, 0.1 eV steps in the edge region, 1 eV steps in the near edge region, and 5 eV steps in the far edge region.The reported data represents the average of five individual scans.Following data averaging and a baseline was applied to each spectrum by fitting the pre-edge region to a polynomial function.This baseline was then subtracted from the whole spectrum.The region above the edge jump was then fit to a two-knot cubic spline, and the data normalized to the edge height.

Electronic Structure Calculation
Electronic structure calculations were performed using ORCA v 4.1.0(Max-Planck-Institut für Kohlenforschung, Mülheim a. d.Ruhr, Germany) [38].Unless otherwise stated, all calculations employed Ahlrichs' def2-tzvp basis set [39] on all atoms and the atom pairwise dispersion correction with Becke-Johnson damping to account for dispersive interactions [40,41].ORCA VeryTightSCF convergence criteria were used for the SCF cycles, with program defaults used for all other convergence criteria and settings.Geometry optimizations were performed at the BP86 level, and used the RI approximation and def2-tzvp/c auxiliary basis set [42].Single point calculations were performed at the PBE0 level and used the RIJCOSX approximation and def2-tzvp/j auxiliary basis set [43].
Nickel and sulfur K-edge X-ray absorption spectra were simulated using TD-DFT calculations (PBE0/def2-tzvp(-f) and the ZORA relativistic approximation) examining the first 25 transitions originating from each sulfur atom or the nickel atom.A Gaussian function to each transition (FWHM = 0.75 eV for S K-edge calculations and 1.2 eV for Ni K-edge calculations) followed by summing the individual transitions was used to produce the calculated spectra.A +36.7 eV (sulfur K-edge) or +171.3 eV (nickel K-edge) energy correction was applied to each transition.Atomic orbital population analyses were performed using a Löwdin population analysis.Isosurface plots were generated using Chimera v 1.13.1 [44].

Generation of {Ni 2 II (SOD mds )}
The apo-peptide SOD mds (SOD mds = (T a CDLP-CGVYDA-PA) 2 ) was prepared by standard solid phase peptide synthesis using Fmoc/ t Bu protection strategies.Cleavage of the peptide from the resin and subsequent global deprotection of the side-chain residues was effected by so-called reagent K. Reagent K is a peptide cleavage mixture that is used for peptides containing readily oxidizable residues such as Cys and Tyr, and contains a number of scavengers that dramatically reduce oxidative side reactions [45].Despite the use of reagent K, we found that the apo-peptide formed a disulfide bond upon cleavage from the resin; even workup under anaerobic conditions using rigorously purified diethyl ether yielded a crude peptide lacking monomeric disulfide free peptide as evidenced by mass spectrometry (MS) and an Ellman's assay (Figure 1). Figure 1 depicts an HPLC chromatogram of the crude product.LCMS examination of all components of the reaction mixture demonstrated that no monomeric peptide containing the N-terminal 2-mercaptoacetate group was produced.

Electronic Structure Calculation
Electronic structure calculations were performed using ORCA v 4.1.0(Max-Planck-Institut für Kohlenforschung, Mülheim a. d.Ruhr, Germany) [38].Unless otherwise stated, all calculations employed Ahlrichs' def2-tzvp basis set [39] on all atoms and the atom pairwise dispersion correction with Becke-Johnson damping to account for dispersive interactions [40,41].ORCA VeryTightSCF convergence criteria were used for the SCF cycles, with program defaults used for all other convergence criteria and settings.Geometry optimizations were performed at the BP86 level, and used the RI approximation and def2-tzvp/c auxiliary basis set [42].Single point calculations were performed at the PBE0 level and used the RIJCOSX approximation and def2-tzvp/j auxiliary basis set [43].
Nickel and sulfur K-edge X-ray absorption spectra were simulated using TD-DFT calculations (PBE0/def2-tzvp(-f) and the ZORA relativistic approximation) examining the first 25 transitions originating from each sulfur atom or the nickel atom.A Gaussian function to each transition (FWHM = 0.75 eV for S K-edge calculations and 1.2 eV for Ni K-edge calculations) followed by summing the individual transitions was used to produce the calculated spectra.A +36.7 eV (sulfur K-edge) or +171.3 eV (nickel K-edge) energy correction was applied to each transition.Atomic orbital population analyses were performed using a Löwdin population analysis.Isosurface plots were generated using Chimera v 1.13.1 [44].

Generation of {Ni2 II (SOD mds )}
The apo-peptide SOD mds (SOD mds = (T a CDLP-CGVYDA-PA)2) was prepared by standard solid phase peptide synthesis using Fmoc/ t Bu protection strategies.Cleavage of the peptide from the resin and subsequent global deprotection of the side-chain residues was effected by so-called reagent K. Reagent K is a peptide cleavage mixture that is used for peptides containing readily oxidizable residues such as Cys and Tyr, and contains a number of scavengers that dramatically reduce oxidative side reactions [45].Despite the use of reagent K, we found that the apo-peptide formed a disulfide bond upon cleavage from the resin; even workup under anaerobic conditions using rigorously purified diethyl ether yielded a crude peptide lacking monomeric disulfide free peptide as evidenced by mass spectrometry (MS) and an Ellman's assay (Figure 1). Figure 1 depicts an HPLC chromatogram of the crude product.LCMS examination of all components of the reaction mixture demonstrated that no monomeric peptide containing the N-terminal 2-mercaptoacetate group was produced.An Ellman's assay [38] demonstrated that SOD mds possessed four free thiol groups per dimeric apo-peptide.This is similar to what has been reported for all NiSOD metallopeptides generated to date.For these NiSOD inspired peptides lacking the N-terminal thiolate group, cleavage of the apo-peptide from the resin and subsequent aerobic work-up and purification yield purified monomeric peptides free of disulfide bonds between the two cysteinate sulfur atoms corresponding to Cys2 and Cys6 in the wild-type NiSOD sequence [26,27,30,34].Considering the only major modification between SOD mds and similar NiSOD inspired peptides generated to date is the presence of the 2-mercaptoacetate group, we suggest that the disulfide moiety within SOD mds results from the oxidative S-S bond formation between two N-terminal thiol groups from two different peptides.Furthermore, all identifiable monomeric peptides by LCMS resulting from incomplete/mis-coupling events lacked the N-terminal 2-mercaptoacetate group, lending further support of the N-terminal disulfide bond in SOD mds .This supposition will be shown to be further validated by computational modeling of the resulting nickel site (vide infra).
The addition of two equivalents of NiCl 2 to a pH 7.4 solution (50 mM N-ethylmorpholine, NEM) of SOD mds , forming {Ni 2 (SOD mds )}, yields a pinkish-tan solution.Raising the pH to 9.6 (50 mM NEM) causes a color change from the pinkish tan color to a more intensely colored reddish-brown solution.Although subtle changes in the far visible region of the electronic absorption spectra are noted at elevated solution pH, a much larger change is noted in the CD spectrum upon changing the solution pH from 7.4 to 9.6 (Figure 2).Monitoring the change in the electronic absorption spectrum vs. the change in pH shows a significant increase in the absorbance at 320 nm with two inflection points in the pH profile at a pH of 8.4 and 9.1.The change in the spectra in response to a change in pH is fully reversible.This suggests that there are two protonatable sites within {Ni 2 (SOD mds )}.We note that the Y9 phenolic group will become deprotonated at the higher pH value; however, the change in absorbance at 320 nm upon deprotonation of phenol will be minimal and not dramatically influence the intensity of the electronic absorption spectrum at 320 nm.Thus, we do not attribute the second deprotonation/protonation event to Y9 deprotonation.
Inorganics 2019, 7, x FOR PEER REVIEW 6 of 7 phase of a mixture of 0.1% TFA in water and 0.1% TFA in acetonitrile and a linear gradient of 9:1 water:acetonitrile-4:6 water:acetonitrile over the course of 30 min.Right: MALDI-TOF of the purified peptide SOD mds (* indicates [SOD mds ] + with the YDPA residues cleaved from the C-terminus of one of the monomers).
An Ellman's assay [38] demonstrated that SOD mds possessed four free thiol groups per dimeric apo-peptide.This is similar to what has been reported for all NiSOD metallopeptides generated to date.For these NiSOD inspired peptides lacking the N-terminal thiolate group, cleavage of the apopeptide from the resin and subsequent aerobic work-up and purification yield purified monomeric peptides free of disulfide bonds between the two cysteinate sulfur atoms corresponding to Cys2 and Cys6 in the wild-type NiSOD sequence [26,27,30,34].Considering the only major modification between SOD mds and similar NiSOD inspired peptides generated to date is the presence of the 2mercaptoacetate group, we suggest that the disulfide moiety within SOD mds results from the oxidative S-S bond formation between two N-terminal thiol groups from two different peptides.Furthermore, all identifiable monomeric peptides by LCMS resulting from incomplete/mis-coupling events lacked the N-terminal 2-mercaptoacetate group, lending further support of the N-terminal disulfide bond in SOD mds .This supposition will be shown to be further validated by computational modeling of the resulting nickel site (vide infra).
The addition of two equivalents of NiCl2 to a pH 7.4 solution (50 mM N-ethylmorpholine, NEM) of SOD mds , forming {Ni2(SOD mds )}, yields a pinkish-tan solution.Raising the pH to 9.6 (50 mM NEM) causes a color change from the pinkish tan color to a more intensely colored reddish-brown solution.Although subtle changes in the far visible region of the electronic absorption spectra are noted at elevated solution pH, a much larger change is noted in the CD spectrum upon changing the solution pH from 7.4 to 9.6 (Figure 2).Monitoring the change in the electronic absorption spectrum vs. the change in pH shows a significant increase in the absorbance at 320 nm with two inflection points in the pH profile at a pH of 8.4 and 9.1.The change in the spectra in response to a change in pH is fully reversible.This suggests that there are two protonatable sites within {Ni2(SOD mds )}.We note that the Y9 phenolic group will become deprotonated at the higher pH value; however, the change in absorbance at 320 nm upon deprotonation of phenol will be minimal and not dramatically influence the intensity of the electronic absorption spectrum at 320 nm.Thus, we do not attribute the second deprotonation/protonation event to Y9 deprotonation.Inorganics 2019, 7, 90 7 of 18

Nickel K-Edge X-ray Absorption Spectroscopy
The nickel K-edge X-ray absorption spectra of {Ni 2 II (SOD mds )} at pH 7.4 and 9.6 are depicted in

Nickel K-edge X-Ray Absorption Spectroscopy
The nickel K-edge X-ray absorption spectra of {Ni2 II (SOD mds )} at pH 7.4 and 9.6 are depicted in Figure 3.The XANES at both pH 7.4 and 9.6 are consistent with Ni II contained in a nomial square planar coordination environment; both display a weak Ni(1s → 3d) transition and a more intense low energy Ni(1s → 4pz) transition [46].However, the XANES spectra display notable differences at the two different solution pH values.At a pH of 7.4, the XANES region contains a poorly resolved, weak nominal Ni(1s → 3d) transition at 8333.4(3) eV and a higher energy more intense Ni(1s → 4pz) transition at 8337.3 (2) eV that is well separated from the edge.Raising the pH to 9.6 changes the overall shape of the edge indicating a change in coordination geometry about nickel.In addition, the Ni(1s → 3d) transition becomes more intense and red-shifts to 8331.0(2) eV while the Ni(1s → 4pz) transition becomes less intense, blue-shifts to 8338.2(4) eV, and is now poorly resolved from the edge.These observed changes indicate that a more rigorous square-planar coordination environment about nickel is generated at pH 7.4, which undergoes a distortion towards nominal D2d symmetry at pH 9.6 (vide infra).At pH 7.4, the EXAFS region of the Ni K-edge X-ray absorption spectrum is best modeled as a four-coordinate Ni-center with two Ni-S scatterers at 2.18 Å and two Ni-N/O scatterer at 1.91 Å.This model is consistent with a nickel center ligated by two cysteinate S atoms, one amidate nitrogen atom and a water ligand, with the N/O scatterers modeled in one shell.At pH 9.6, we see the loss of one N/O scatterer with the subsequent coordination of an addition sulfur ligand; the best fit to the data contained three Ni-S scatterers at 2.23 Å and one Ni-N scatterer at 1.89 Å.In addition, a vector for an outersphere Ni-Ni scatterer could be located at 3.25 Å.These data are consistent with the structural models depicted in Scheme 1; at pH 7.4, the data is consistent with each peptide containing At pH 7.4, the EXAFS region of the Ni K-edge X-ray absorption spectrum is best modeled as a four-coordinate Ni-center with two Ni-S scatterers at 2.18 Å and two Ni-N/O scatterer at 1.91 Å.This model is consistent with a nickel center ligated by two cysteinate S atoms, one amidate nitrogen atom and a water ligand, with the N/O scatterers modeled in one shell.At pH 9.6, we see the loss of one N/O scatterer with the subsequent coordination of an addition sulfur ligand; the best fit to the data contained three Ni-S scatterers at 2.23 Å and one Ni-N scatterer at 1.89 Å.In addition, a vector for an outersphere Ni-Ni scatterer could be located at 3.25 Å.These data are consistent with the structural models depicted in Scheme 1; at pH 7.4, the data is consistent with each peptide containing two mononuclear nickel sites that collapses into a binuclear cysteinate bridged Ni 2 center at high pH, the formation of which is well represented in small molecule Ni-thiolate chemistry [47,48].

Sulfur K-edge X-Ray Absorption Spectroscopy of {Ni2 II (SOD mds )}
The sulfur K-edge X-ray absorption spectra of {Ni2 II (SOD mds )} at pH 7.4 and 9.6 are depicted in Figure 4.At both pH values, there are higher energy edge features that are consistent with a disulfide bond, further demonstrating the presences of a bridging disulfide at both pH 7.4 and 9.6.At high pH, there is a pre-edge feature that possesses distinct asymmetry, and can be identified as terminal and bridging nickel-cysteinate S(1s → LUMO) transitions [15].At pH 7.4, this feature is lost while the edge broadens and gains intensity.This suggests that the S(1s → LUMO) transitions have been blueshifted into the edge, and is fully consistent with the EXAFS analysis of the metallopeptide at high and low pH; lowering the pH that induces the breaking of the bridging Ni-S Cys -Ni bonds, generating two mononuclear nickel centers.Furthermore, the sulfur K-edge spectrum at pH 7.4 is consistent with the formation of protonated coordinated Ni-S(H + )-Cys moieties at physiological pH; the reversible blue-shift of the S(1s → LUMO) transition upon lowering the pH is a hallmark of the formation of Ni-S(H + )-Cys moieties [15].This is further validated by time-dependent DFT calculations simulating the S K-edge X-ray absorption spectra (vide infra).

DFT Generated Structure for {Ni2 II (SOD mds )} and TD-DFT Calculations of Ni and S K-edge XAS Transitions
Based on the Ni and S K-edge XAS data, several computational models of the high and low pH forms of {Ni2 II (SOD mds )} were generated (Figure 5).One is a mononuclear nickel site that contains a formal Ni II center in an S2NO coordination environment with ligands derived from two thiolate sulfur atoms, an amidate nitrogen atom, and a water ligand.Models of this mononuclear site with variable thiolate sulfur atom protonation were also investigated.As computational models for the pH = 9.6 form of {Ni2 II (SOD mds )}, two dinuclear thiolate bridged Ni2 models were also generated.One possessed two bridging thiolates, one terminal thiolate and one amidate nitrogen per nickel center.The other dinuclear model possessed an identical Ni2S2 core, but with a disulfide bond derived from the 2-mercaptoacetate groups bridging the two ligand sets together.Geometry optimizations were 3.3.Sulfur K-Edge X-ray Absorption Spectroscopy of {Ni 2 II (SOD mds )} The sulfur K-edge X-ray absorption spectra of {Ni 2 II (SOD mds )} at pH 7.4 and 9.6 are depicted in Figure 4.At both pH values, there are higher energy edge features that are consistent with a disulfide bond, further demonstrating the presences of a bridging disulfide at both pH 7.4 and 9.6.At high pH, there is a pre-edge feature that possesses distinct asymmetry, and can be identified as terminal and bridging nickel-cysteinate S(1s → LUMO) transitions [15].At pH 7.4, this feature is lost while the edge broadens and gains intensity.This suggests that the S(1s → LUMO) transitions have been blue-shifted into the edge, and is fully consistent with the EXAFS analysis of the metallopeptide at high and low pH; lowering the pH that induces the breaking of the bridging Ni-S Cys -Ni bonds, generating two mononuclear nickel centers.Furthermore, the sulfur K-edge spectrum at pH 7.4 is consistent with the formation of protonated coordinated Ni-S(H + )-Cys moieties at physiological pH; the reversible blue-shift of the S(1s → LUMO) transition upon lowering the pH is a hallmark of the formation of Ni-S(H + )-Cys moieties [15].This is further validated by time-dependent DFT calculations simulating the S K-edge X-ray absorption spectra (vide infra).
Inorganics 2019, 7, x FOR PEER REVIEW 8 of 19 two mononuclear nickel sites that collapses into a binuclear cysteinate bridged Ni2 center at high pH, the formation of which is well represented in small molecule Ni-thiolate chemistry [47,48].
. Scheme 1. Interconversion of the structures of the nickel-sites of {Ni2 II (SOD mda )} at pH 7.4 and 9.6 and subsequent oxidative decomposition upon O2 exposure.

Sulfur K-edge X-Ray Absorption Spectroscopy of {Ni2 II (SOD mds )}
The sulfur K-edge X-ray absorption spectra of {Ni2 II (SOD mds )} at pH 7.4 and 9.6 are depicted in Figure 4.At both pH values, there are higher energy edge features that are consistent with a disulfide bond, further demonstrating the presences of a bridging disulfide at both pH 7.4 and 9.6.At high pH, there is a pre-edge feature that possesses distinct asymmetry, and can be identified as terminal and bridging nickel-cysteinate S(1s → LUMO) transitions [15].At pH 7.4, this feature is lost while the edge broadens and gains intensity.This suggests that the S(1s → LUMO) transitions have been blueshifted into the edge, and is fully consistent with the EXAFS analysis of the metallopeptide at high and low pH; lowering the pH that induces the breaking of the bridging Ni-S Cys -Ni bonds, generating two mononuclear nickel centers.Furthermore, the sulfur K-edge spectrum at pH 7.4 is consistent with the formation of protonated coordinated Ni-S(H + )-Cys moieties at physiological pH; the reversible blue-shift of the S(1s → LUMO) transition upon lowering the pH is a hallmark of the formation of Ni-S(H + )-Cys moieties [15].This is further validated by time-dependent DFT calculations simulating the S K-edge X-ray absorption spectra (vide infra).

DFT Generated Structure for {Ni2 II (SOD mds )} and TD-DFT Calculations of Ni and S K-edge XAS Transitions
Based on the Ni and S K-edge XAS data, several computational models of the high and low pH forms of {Ni2 II (SOD mds )} were generated (Figure 5).One is a mononuclear nickel site that contains a formal Ni II center in an S2NO coordination environment with ligands derived from two thiolate sulfur atoms, an amidate nitrogen atom, and a water ligand.Models of this mononuclear site with variable thiolate sulfur atom protonation were also investigated.As computational models for the pH = 9.6 form of {Ni2 II (SOD mds )}, two dinuclear thiolate bridged Ni2 models were also generated.One possessed two bridging thiolates, one terminal thiolate and one amidate nitrogen per nickel center.The other dinuclear model possessed an identical Ni2S2 core, but with a disulfide bond derived from the 2-mercaptoacetate groups bridging the two ligand sets together.Geometry optimizations were  Based on the Ni and S K-edge XAS data, several computational models of the high and low pH forms of {Ni 2 II (SOD mds )} were generated (Figure 5).One is a mononuclear nickel site that contains a formal Ni II center in an S 2 NO coordination environment with ligands derived from two thiolate sulfur atoms, an amidate nitrogen atom, and a water ligand (see Supplementary Materials).Models of this mononuclear site with variable thiolate sulfur atom protonation were also investigated.As computational models for the pH = 9.6 form of {Ni 2 II (SOD mds )}, two dinuclear thiolate bridged Ni 2 models were also generated.One possessed two bridging thiolates, one terminal thiolate and one amidate nitrogen per nickel center.The other dinuclear model possessed an identical Ni 2 S 2 core, but with a disulfide bond derived from the 2-mercaptoacetate groups bridging the two ligand sets together.Geometry optimizations were performed at the BP86/def2-tzvp level of theory with a dispersion correction; this functional generally yields computationally derived structures for transition metal complexes that are in excellent agreement with experimental data [49,50].
performed at the BP86/def2-tzvp level of theory with a dispersion correction; this functional generally yields computationally derived structures for transition metal complexes that are in excellent agreement with experimental data [49,50].An examination of the metric parameters derived for the mononuclear Ni-sites (pH 7.4 monomeric model) demonstrates that a model with each thiolate sulfur becoming protonated and a water ligand bound to nickel is most consistent with the experimental parameters derived from EXAFS.The two Ni-S bond lengths are calculated to be 2.169 Å (trans to the water ligand) and 2.196 (trans to the amidate nitrogen); the Ni-N bond length is calculated to be 1.925 Å and the Ni-O bond length is calculated to be 1.889 Å.These are fully consistent with the EXAFS derived average Ni-S and Ni-N/O bond lengths.The mono-and unprotonated models yield Ni-ligand bond lengths that are less consistent with the EXAFS data than the doubly protonated model.However, when one considers the inherent error in EXAFS-derived bond lengths, the mono and unprotonated models could still be considered valid structural models for {Ni2 II (SOD mds )} based on bond-lengths alone.Nevertheless, the experimental vs. theoretical Ni and S K-edge studies strongly suggest that the doubly-protonated model is most valid (vide infra).
It was determined that inclusion of the N-terminal disulfide bridge is required to reproduce the metric parameters observed in the high pH form of {Ni2 II (SOD mds )}.When the disulfide moiety is omitted from the computational model, the geometry about the nickel-sites is nearly planar: the S-Ni-S-Ni dihedral angles are all ~3°.This planar geometry yields a Ni-Ni distance that is 0.An examination of the metric parameters derived for the mononuclear Ni-sites (pH 7.4 monomeric model) demonstrates that a model with each thiolate sulfur becoming protonated and a water ligand bound to nickel is most consistent with the experimental parameters derived from EXAFS.The two Ni-S bond lengths are calculated to be 2.169 Å (trans to the water ligand) and 2.196 (trans to the amidate nitrogen); the Ni-N bond length is calculated to be 1.925 Å and the Ni-O bond length is calculated to be 1.889 Å.These are fully consistent with the EXAFS derived average Ni-S and Ni-N/O bond lengths.The mono-and unprotonated models yield Ni-ligand bond lengths that are less consistent with the EXAFS data than the doubly protonated model.However, when one considers the inherent error in EXAFS-derived bond lengths, the mono and unprotonated models could still be considered valid structural models for {Ni 2 II (SOD mds )} based on bond-lengths alone.Nevertheless, the experimental vs.
theoretical Ni and S K-edge studies strongly suggest that the doubly-protonated model is most valid (vide infra).It was determined that inclusion of the N-terminal disulfide bridge is required to reproduce the metric parameters observed in the high pH form of {Ni 2 II (SOD mds )}.When the disulfide moiety is omitted from the computational model, the geometry about the nickel-sites is nearly planar: the S-Ni-S-Ni dihedral angles are all ~3• .This planar geometry yields a Ni-Ni distance that is 0.1 Å longer than what is observed by EXAFS (3.374 Å).Inclusion of the disulfide bond forces the geometry around the nickel sites to become distorted away from planarity.Furthermore, the two nickel sites are no longer equivalent to one another with one center distorted more toward a tetrahedral-like coordination environment than the other (the above measured dihedral angles are now 33 • and 27 • ).This distortion brings the two Ni-Ni centers closer to one another by 0.1 Å, which reproduces the EXAFS derived Ni-Ni distance nicely (3.274 Å).Furthermore, the average Ni-S and Ni-N bond lengths are more consistent with the EXAFS-derived bond lengths.This lends further evidence for the presence of an N-terminal disulfide bond.Correlating the experimental Ni K-edge XANES at high and low pH with the time-dependent DFT (TD-DFT; PBE0/zora-def2-tzvp(-f)/ZORA) calculated Ni K-edge spectra of the computational models above supports the proposed structures (Figure 6).For the monomeric structures, we find that the nearly rigorous square planar coordination environment about Ni II yields Ni(1s → 3d) and Ni(1s → 4p z ) transitions that correlate well with the experimental data.As might be expected, the pH 9.6 model lacking the disulfide bridge resembles that of the pH 7.4 data, but is inconsistent with the pH 9.6 experimental data.Once the disulfide bridge is included in the computational model, the Ni(1s → 3d) transition red-shifts and increases in intensity while the Ni(1s → 4p z ) transition blue shifts and decreases in intensity.The reason for this is a consequence of 3d/4p mixing.In the D 2d distorted non-centrosymmetric coordination environment, the loss of the pseduo-inversion center allows for the Ni(1s → 3d) transition to acquire 4p character and the Ni(1s → 4p z ) transition to acquire 3d character.The result is an increase in the intensity of the Ni(1s → 3d) transition by gaining the dipole-allowed Ni(1s → 4p) character while the Ni(1s → 4p z ) transition losses intensity owing to an increase dipole-forbidden Ni(1s → 3d) character and a decrease in Ni(1s → 4p) character (Figure 6).lengths are more consistent with the EXAFS-derived bond lengths.This lends further evidence for the presence of an N-terminal disulfide bond.
Correlating the experimental Ni K-edge XANES at high and low pH with the time-dependent DFT (TD-DFT; PBE0/zora-def2-tzvp(-f)/ZORA) calculated Ni K-edge spectra of the computational models above supports the proposed structures (Figure 6).For the monomeric structures, we find that the nearly rigorous square planar coordination environment about Ni II yields Ni(1s → 3d) and Ni(1s → 4pz) transitions that correlate well with the experimental data.As might be expected, the pH 9.6 model lacking the disulfide bridge resembles that of the pH 7.4 data, but is inconsistent with the pH 9.6 experimental data.Once the disulfide bridge is included in the computational model, the Ni(1s → 3d) transition red-shifts and increases in intensity while the Ni(1s → 4pz) transition blue shifts and decreases in intensity.The reason for this is a consequence of 3d/4p mixing.In the D2d distorted noncentrosymmetric coordination environment, the loss of the pseduo-inversion center allows for the Ni(1s → 3d) transition to acquire 4p character and the Ni(1s → 4pz) transition to acquire 3d character.The result is an increase in the intensity of the Ni(1s → 3d) transition by gaining the dipole-allowed Ni(1s → 4p) character while the Ni(1s → 4pz) transition losses intensity owing to an increase dipoleforbidden Ni(1s → 3d) character and a decrease in Ni(1s → 4p) character (Figure 6).Calculation of the S K-edge X-ray absorption spectra using time-dependent DFT (TD-DFT; PBE0/zora-def2-tzvp(-f)/ZORA) also supports the structural assignments proposed above.The experimental S K-edge X-ray absorption spectrum of {Ni2 II (SOD mds )} obtained at pH 9.6 and the TD-DFT calculated S K-edge X-ray absorption spectrum using the disulfide bridged model match well, lending further credence to the validity of the proposed structural model (Figure 7).To determine the likely protonation state of the monomeric pH 7.4 form of {Ni2 II (SOD mds )}, the doubly-protonated, monoprotonated, and unprotonated monomeric models outlined above were examined.When one or both of the cysteinate sulfur atom(s) is unprotonated, a low energy pre-edge feature corresponding to the nominal S(1s → LUMO) transitions is produced in the calculated spectrum.As would be predicted, this pre-edge feature is weaker for the mono-protonated model relative to the unprotonated model owing to the fact that only one S(1s → LUMO) transition comprises this feature as opposed to two; the S(1s → LUMO) transition for a protonated thiolate sulfur atom will blue shift into the edge.It is therefore only the doubly-protonated model that reproduces the edge feature of the pH 7.4 experimental spectrum of {Ni2 II (SOD mds )}.The TD-DFT calculated S K-edge X-ray absorption spectrum for the doubly-protonated model has both S(1s → LUMO) transitions blue shifted by ~3.5 eV relative to the energy of the unprotonated sulfur atom S(1s → LUMO) transition, which reproduce the experimental spectrum well.Thus, we conclude that the pH 7.4 form of {Ni2 II (SOD mds )} possesses two protonated cysteinate sulfur atoms coordinated to nickel.Calculation of the S K-edge X-ray absorption spectra using time-dependent DFT (TD-DFT; PBE0/zora-def2-tzvp(-f)/ZORA) also supports the structural assignments proposed above.The experimental S K-edge X-ray absorption spectrum of {Ni 2 II (SOD mds )} obtained at pH 9.6 and the TD-DFT calculated S K-edge X-ray absorption spectrum using the disulfide bridged model match well, lending further credence to the validity of the proposed structural model (Figure 7).To determine the likely protonation state of the monomeric pH 7.4 form of {Ni 2 II (SOD mds )}, the doubly-protonated, monoprotonated, and unprotonated monomeric models outlined above were examined.When one or both of the cysteinate sulfur atom(s) is unprotonated, a low energy pre-edge feature corresponding to the nominal S(1s → LUMO) transitions is produced in the calculated spectrum.As would be predicted, this pre-edge feature is weaker for the mono-protonated model relative to the unprotonated model owing to the fact that only one S(1s → LUMO) transition comprises this feature as opposed to two; the S(1s → LUMO) transition for a protonated thiolate sulfur atom will blue shift into the edge.It is therefore only the doubly-protonated model that reproduces the edge feature of the pH 7.4 experimental spectrum of {Ni 2 II (SOD mds )}.The TD-DFT calculated S K-edge X-ray absorption spectrum for the doubly-protonated model has both S(1s → LUMO) transitions blue shifted by ~3.5 eV relative to the energy of the unprotonated sulfur atom S(1s → LUMO) transition, which reproduce the experimental spectrum well.Thus, we conclude that the pH 7.4 form of {Ni 2 II (SOD mds )} possesses two protonated cysteinate sulfur atoms coordinated to nickel.{Ni2 II (SOD mds )} reacts with O2 at both pH 7.4 and 9.6.The kinetics of O2 oxidation was followed by CD spectroscopy over the course of 12 h (Figure 8).It was found that the oxidation kinetics obey a pseudo-first order rate law under constant O2 concentration, and proceeds at a faster rate at high vs. low pH.Extraction of the second order rate constant for the oxidation reactions demonstrates that the reaction at pH 9.6 proceeds with a rate that is over 3.5-fold faster than at pH 7.4 (k = 1.8(3) × 10 −2 M −1 s −1 vs. 6.5(2) × 10 −2 M −1 s −1 ).Based on an analysis of the products formed during the oxidation reaction, we suspect that following the initial oxidation step there are multiple oxidation pathways leading to different final oxidation products.There are a number of items to note concerning the oxidation products formed at pH 7.4 and 9.6.First, under both pH conditions, at least one of the final soluble oxidation products contains nickel, and this product is identical by CD under both pH conditions.However, MS data of the solution and solid materials produced by O2 initiated oxidation indicates a complex mixture of unidentifiable   {Ni2 II (SOD mds )} reacts with O2 at both pH 7.4 and 9.6.The kinetics of O2 oxidation was followed by CD spectroscopy over the course of 12 h (Figure 8).It was found that the oxidation kinetics obey a pseudo-first order rate law under constant O2 concentration, and proceeds at a faster rate at high vs. low pH.Extraction of the second order rate constant for the oxidation reactions demonstrates that the reaction at pH 9.6 proceeds with a rate that is over 3.5-fold faster than at pH 7.4 (k = 1.8(3) × 10 −2 M −1 s −1 vs. 6.5(2) × 10 −2 M −1 s −1 ).Based on an analysis of the products formed during the oxidation reaction, we suspect that following the initial oxidation step there are multiple oxidation pathways leading to different final oxidation products.There are a number of items to note concerning the oxidation products formed at pH 7.4 and 9.6.First, under both pH conditions, at least one of the final soluble oxidation products contains nickel, and this product is identical by CD under both pH conditions.However, MS data of the solution and solid materials produced by O2 initiated oxidation indicates a complex mixture of unidentifiable There are a number of items to note concerning the oxidation products formed at pH 7.4 and 9.6.First, under both pH conditions, at least one of the final soluble oxidation products contains nickel, and this product is identical by CD under both pH conditions.However, MS data of the solution and solid materials produced by O 2 initiated oxidation indicates a complex mixture of unidentifiable products.Furthermore, the EPR spectra of the resulting products are silent down to 10 K, indicating that the Ni-site is in the formal Ni 2+ oxidation state.The resulting IR spectra of the produced solutions and solids showed no bands corresponding to S=O stretching frequencies, indicating that oxygen atom insertion reactions into the Ni-S moiety does not represent a major oxidation pathway.Instead, it is possible O 2 is initiating irreversible sulfur based ligand oxidation as has been observed in the work of Darensbourg, for example [51].Validating this possibility is the observation that the metallopeptide cannot be cleanly oxidized; attempts to chemically oxidize {Ni 2 II (SOD mds )} at pH 7.4 and 9.6 by a 3% hydrogen peroxide solution, ethanolic I 2 or MnO 4 − lead to the rapid bleaching of the solution and subsequent formation of unidentifiable tan insoluble aggregates, all of which yielded EPR and IR spectra consistent with the above formed from O 2 oxidation of the solutions.
A complex mixture of soluble and insoluble nickel containing products is also noted by Ni K-edge X-ray absorption spectroscopy.At both solution pH values, the Ni K-edge XANES no longer contains the prominent Ni(1s → 4p z ) transition, and is more reminiscent of a six coordinate Ni 2+ species (Figures S1 and S2).Because of the low signal to noise at high k, the EXAFS regions could only be simulated to k = 11 Å −1 The EXAFS region of the decomposition product obtained from air oxidation at pH 7.4 was best modeled as containing 1.4 Ni-S interactions (r = 2.22 Å) and 4.8 N/O interactions (r = 1.97 Å).The EXAFS region of the air oxidation product generated at pH 9.6 was best modeled with 0.6 Ni-S interactions (r = 2.24 Å), 4.3 N/O interaction (r = 1.94 Å), and 2.1 Ni-Ni interactions (r = 3.25 Å).In both cases, the resulting fitting statistics are poor with ε 2 values greater than 3.As this represents a mixture of soluble and insoluble compounds in a number of coordination environments, formulating likely structures about the nickel center(s) is not warranted based on the available data.
Electronic structure calculations suggest the reason for the increased stability of {Ni 2 (SOD mds )} at pH 7.4 vs. 9.6 results from the deactivation of the high-lying sulfur dominated Ni(3dπ)-S(3pπ) anti-bonding orbital upon protonation (Figure 9, Tables 1 and 2).For the monomeric species, the HOMO is identified as a nickel dominated Ni(3dπ)-N(2pπ) anti-bonding orbital.This is destabilized by 0.29 eV relative to the essentially non-bonding Ni(3d z2 ) orbital (HOMO-1).The HOMO-2 is a water O(2p) dominated O(2pπ)-Ni(3dπ) antibonding orbital followed by the Ni(3d xz ) dominated HOMO-3.Thus, none of the frontier MOs (FMOs) possess significant S(3p) character, rendering the sulfur atoms reasonably unreactive towards oxidative damage by O 2 .In contrast, deprotonation of the sulfur atoms of the mononuclear {Ni 2 II (SOD mds )} nickel-site dramatically alters the electronic structure of the complex.Electronic structure calculations reveal that the HOMO and HOMO-1 are significantly activated relative to the essentially non-bond Ni(3d z2 ) HOMO-2 by 0.73 and 0.42 eV, respectively.Furthermore, these two orbitals are S(3p) dominated S(3pπ)-Ni(3dπ) anti-bonding orbitals.Therefore, if the deprotonated monomeric form of {Ni 2 (SOD mds )} could be generated, we would predict it would be highly sensitive to O 2 damage owing to the activated S(3p) dominated HOMO.Given the slow rate of O 2 oxidation of {Ni 2 (SOD mds )} at pH = 7.4, it is likely that deprotonation of the Ni-S(H + )-Cys moieties are required for oxidative damage; as protonation is an equilibrium process, there will always be a small concentration of O 2 reactive unprotonated Ni-S-Cys bonds in solution.The disulfide bridged computational model was found to possess highly covalent Ni-ligand bonds.The HOMO, although containing less S(3p) than the HOMO of the unprotonated mononuclear computational model, still contains a significant amount of S(3p) character, and would therefore be predicted to be more susceptible to O 2 damage than the doubly-protonated mononuclear nickel-site.This is observed experimentally.

Discussion
In this study we have demonstrated the reversible formation of a dinuclear Ni

Discussion
In this study we have demonstrated the reversible formation of a dinuclear Ni2 II site within a peptide in response to pH.Dinuclear Ni2 II sites are not observed in other monomeric SOD metallopeptide-based mimics while the disulfide linked dimeric metallopeptide facilitates formation of a dinuclear Ni2 center.We speculate that dimer formation in {Ni2(SOD mds )} is facilitated by the close proximity of the two metal centers, which was enforced through a disulfide bond linkage near the individual nickel centers.
The driving force for conversion of the dinuclear Ni2 II site into two monomeric Ni II sites within {Ni2(SOD mds )} likely involves protonation of a terminal cysteinate sulfur atom.The nucleophilic HOMO of the dinuclear Ni2 II site possesses S(3p) character corresponding to the terminal thiolate sulfur coordinated to the more distorted nickel center, making it the likely protonation site.Protonation of that sulfur atom shortens the Ni-S(H + )-Cys bond length relative to the unprotonated model with a concomitant increase in the bridging Ni-S bond length trans to the protonation site (Figure 10).Association of water to the nickel center followed by breaking of the weakened Ni-S bond would lead to the eventual formation of two mononuclear Ni II centers.This study also gives insight into an additional mechanism of protection of Ni-thiolate bonds from oxidative damage against reactive oxygen species (ROSs).Nickel thiolates are susceptible to oxidative damage by O2 and H2O2 [52,53], yet the NiSOD active-site is robust against oxidative damage effected by such species.Possible explanations for protection of the NiSOD active-site against Ni-S oxidative damage by ROSs have been proposed, including electronic fine tuning of the Ni-S moiety via the mixed amine/amide coordination environment [54,55] and a fast rate of the O2 − disproportionation reaction relative to the rate of oxidation of the coordinate cysteinate sulfur atoms [45,48,49].In this study we have shown that oxidation of the Ni-S-Cys bond by O2 is slow for {Ni2(SOD mds )} at pH = 7.4, with a half-life of nearly 7 h under the reaction conditions investigated (ambient O2 concentration, 1.0 mM {Ni2(SOD mds )}).This is the result of both the deactivation and a significant reduction in S(3p) character to the nucleophilic FMOs upon protonation; the FMOs of the The disulfide bridge and methylene groups have been represented as small spheres and wires for clarity.
This study also gives insight into an additional mechanism of protection of Ni-thiolate bonds from oxidative damage against reactive oxygen species (ROSs).Nickel thiolates are susceptible to oxidative damage by O 2 and H 2 O 2 [52,53], yet the NiSOD active-site is robust against oxidative damage effected by such species.Possible explanations for protection of the NiSOD active-site against Ni-S oxidative damage by ROSs have been proposed, including electronic fine tuning of the Ni-S moiety via the mixed amine/amide coordination environment [54,55] and a fast rate of the O 2 − disproportionation reaction relative to the rate of oxidation of the coordinate cysteinate sulfur atoms [45,48,49].In this study we have shown that oxidation of the Ni-S-Cys bond by O 2 is slow for {Ni 2 (SOD mds )} at pH = 7.4, with a half-life of nearly 7 h under the reaction conditions investigated (ambient O 2 concentration, 1.0 mM {Ni 2 (SOD mds )}).This is the result of both the deactivation and a significant reduction in S(3p) character to the nucleophilic FMOs upon protonation; the FMOs of the doubly-protonated monomeric {Ni 2 (SOD mds )} computational model possesses little S(3p) character, while the HOMO and HOMO-1 of the unprotonated {Ni 2 (SOD mds )} computational model are both energetically activated and possess a large degree of S(3p) character.Thus, protonation will inherently protect the thiolate sulfur atoms against ROS induced oxidative damage.The dinuclear nickel site of {Ni 2 (SOD mds )} produced at pH = 9.6 was found to undergo oxidation at an increased rate relative to the mononuclear pH = 7.4 form.This is to be expected as the degree of terminal cysteinate S(3p) character to the more covalent anti-bonding FMOs has increased relative to the mononuclear pH = 7.4 doubly-protonated Ni form, rendering it more susceptible to oxidative damage.This study has also pointed to an additional role for the protonation of coordinated cysteinate sulfur atoms at metalloenzyme active-sites-poising the centers for reactivity via electronic fine tuning.As demonstrated above, the FMOs of the doubly-protonated mononuclear {Ni 2 (SOD mds )} are biased to the nickel-site, while those of the unprotonated mononuclear {Ni 2 (SOD mds )} computational model are biased to the thiolate sulfur atoms.Thus, one would expect that upon thiolate sulfur atom protonation reactivity would be shifted from the thiolate sulfur atoms to the nickel-site.For example, one could consider the active-site of [NiFe]H 2 ase, which has been shown to possess a terminal protonated cysteinate residue ligated to a D 2d distorted nickel center.It has been proposed that a key intermediate in the reactivity of [NiFe]H 2 ase is a formal Ni III -H species.Protonation of the cysteinate sulfur atom would bias the reactivity towards nickel, making the nickel-center the site susceptible to subsequent protonation events.Thus, one may envision a dual role for the Ni-S(H + )-Cys moiety in [NiFe]H 2 ase-gating of reactivity and proton donation to the hydride ligand.

Conclusions
A disulfide bridged metallopeptide has been prepared inspired by the metalloenzyme NiSOD.This metallopeptide contains two nickel centers in close proximity owing to a disulfide bridge between two peptide monomers.In response to pH, the mononuclear nickel-sites found at pH 7.4, which contain Ni II in an S 2 NO coordination motif reversibly form a dinuclear cysteinate Ni 2 II center at elevated pH (pH 9.6).The driving force for the interconversion of the dinuclear nickel center to two mononuclear nickel sites is proposed to be cysteinate S-atom protonation, which results in two coordinated protonated Cys S atoms at lower pH.It was shown that these Ni-S(H + )-Cys moieties reduce the O 2 initiated oxidative damage of the nickel-site, likely through the modulation of the electronic structure of the Ni-center rendering the S-atoms less nucleophilic upon protonation.This may have relevance in biological Ni systems, offering the Ni-S-Cys moiety protection against oxidative damage upon Ni-S(H + )-Cys formation.Furthermore, the modulation of the electronic structure of the Ni-site upon Ni-S(H + )-Cys formation suggests that reversible cysteinate sulfur atom protonation may be involved in the gating of biological reactivity at such metal-centers.

Figure 1 .
Figure 1.Left: Analytical HPLC chromatogram of the crude mixture resulting from the synthesis of SOD mds with a detector cutoff of 0 intensity units.Identifiable products are highlighted.A mobile

Figure 1 .
Figure 1.Left: Analytical HPLC chromatogram of the crude mixture resulting from the synthesis of SOD mds with a detector cutoff of 0 intensity units.Identifiable products are highlighted.A mobile phase of a mixture of 0.1% TFA in water and 0.1% TFA in acetonitrile and a linear gradient of 9:1 water:acetonitrile-4:6 water:acetonitrile over the course of 30 min.Right: MALDI-TOF of the purified peptide SOD mds (* indicates [SOD mds ] + with the YDPA residues cleaved from the C-terminus of one of the monomers).

Figure 2 .
Figure 2. Left: electronic absorption spectra (bottom) and CD spectra (top) of {Ni 2 II (SOD mds )} at pH 7.4 and 9.6.The solid black spectra represent the experimental data, the red dashed curves represent the individual transitions deconvolved from the spectra and the dashed black spectra represent the convolution of the individual transitions.Right: pH profile showing the change in absorbance at 320 nm vs. change in pH upon going from pH 7.4 to 9.6 (red circles) and pH 9.6 to 7.4 (blue squares).

Figure 3 .
Figure3.The XANES at both pH 7.4 and 9.6 are consistent with Ni II contained in a nomial square planar coordination environment; both display a weak Ni(1s → 3d) transition and a more intense low energy Ni(1s → 4p z ) transition[46].However, the XANES spectra display notable differences at the two different solution pH values.At a pH of 7.4, the XANES region contains a poorly resolved, weak nominal Ni(1s → 3d) transition at 8333.4(3) eV and a higher energy more intense Ni(1s → 4p z ) transition at 8337.3(2) eV that is well separated from the edge.Raising the pH to 9.6 changes the overall shape of the edge indicating a change in coordination geometry about nickel.In addition, the Ni(1s → 3d) transition becomes more intense and red-shifts to 8331.0(2) eV while the Ni(1s → 4p z ) transition becomes less intense, blue-shifts to 8338.2(4) eV, and is now poorly resolved from the edge.These observed changes indicate that a more rigorous square-planar coordination environment about nickel is generated at pH 7.4, which undergoes a distortion towards nominal D 2d symmetry at pH 9.6 (vide infra).

Figure 2 .
Figure 2. Left: electronic absorption spectra (bottom) and CD spectra (top) of {Ni2 II (SOD mds )} at pH 7.4 and 9.6.The solid black spectra represent the experimental data, the red dashed curves represent the individual transitions deconvolved from the spectra and the dashed black spectra represent the convolution of the individual transitions.Right: pH profile showing the change in absorbance at 320 nm vs. change in pH upon going from pH 7.4 to 9.6 (red circles) and pH 9.6 to 7.4 (blue squares).

Scheme 1 .
Scheme 1. Interconversion of the structures of the nickel-sites of {Ni 2 II (SOD mda )} at pH 7.4 and 9.6 and subsequent oxidative decomposition upon O 2 exposure.

3. 4 .
DFT Generated Structure for {Ni 2 II (SOD mds )} and TD-DFT Calculations of Ni and S K-Edge XAS Transitions

Figure 5 .
Figure 5. Computationally derived nickel site models of {Ni2 II (SOD mds )}.Metric parameters are provided next to the Ni-ligand bonds.
1 Å longer than what is observed by EXAFS (3.374 Å).Inclusion of the disulfide bond forces the geometry around the nickel sites to become distorted away from planarity.Furthermore, the two nickel sites are no longer equivalent to one another with one center distorted more toward a tetrahedral-like coordination environment than the other (the above measured dihedral angles are now 33° and 27°).This distortion brings the two Ni-Ni centers closer to one another by 0.1 Å, which reproduces the EXAFS derived Ni-Ni distance nicely (3.274 Å).Furthermore, the average Ni-S and Ni-N bond

Figure 5 .
Figure 5. Computationally derived nickel site models of {Ni 2 II (SOD mds )}.Metric parameters are provided next to the Ni-ligand bonds.

Figure 6 .
Figure 6.Comparison of the TD-DFT calculated Ni K-edge XANES with the experimental spectra for {Ni 2 (SOD mds )} at pH 7.4 (left) and 9.6 (right).

Figure 8 .
Figure 8. CD spectra following the air oxidation of {Ni2 II (SOD mds )} at pH 7.4 (left) and 9.6 (right), with the blue spectra representing the trace at t = 0 s, the red spectra representing the traces recorded every 600 s (10 min) over the course of 12 h, and the teal spectra represent the CD spectra of the solutions following 24 h of O2 exposure.The insets depict the kinetics traces highlighting the decay of {Ni2 II (SOD mds )} (blue trace) and best fit of the kinetic trace to a first order rate law.

Figure 8 .
Figure 8. CD spectra following the air oxidation of {Ni2 II (SOD mds )} at pH 7.4 (left) and 9.6 (right), with the blue spectra representing the trace at t = 0 s, the red spectra representing the traces recorded every 600 s (10 min) over the course of 12 h, and the teal spectra represent the CD spectra of the solutions following 24 h of O2 exposure.The insets depict the kinetics traces highlighting the decay of {Ni2 II (SOD mds )} (blue trace) and best fit of the kinetic trace to a first order rate law.

Figure 8 .
Figure 8. CD spectra following the air oxidation of {Ni 2 II (SOD mds )} at pH 7.4 (left) and 9.6 (right), with the blue spectra representing the trace at t = 0 s, the red spectra representing the traces recorded every 600 s (10 min) over the course of 12 h, and the teal spectra represent the CD spectra of the solutions following 24 h of O 2 exposure.The insets depict the kinetics traces highlighting the decay of {Ni 2 II (SOD mds )} (blue trace) and best fit of the kinetic trace to a first order rate law.

Table 1 .
Ni(3d), S(3p), N(2p), and water O(2p) Löwdin MO population analysis (%AO to MO) and energies (eV) relative to the Ni(3dz2) orbital for the LUMO through HOMO-3 of the computational models for doubly-protonated and unprotonated monomeric {Ni2(SOD mds )} computational models.Orbital compositions for the doubly-protonated model are given above the unprotonated model for each AO.

Table 1 .
Ni(3d), S(3p), N(2p), and water O(2p) Löwdin MO population analysis (%AO to MO) and energies (eV) relative to the Ni(3d z2 ) orbital for the LUMO through HOMO-3 of the computational models for doubly-protonated and unprotonated monomeric {Ni 2 (SOD mds )} computational models.Orbital compositions for the doubly-protonated model are given above the unprotonated model for each AO.LUMO HOMO HOMO-1 HOMO-2 HOMO- 2 II site within a peptide in response to pH.Dinuclear Ni 2 II sites are not observed in other monomeric SOD metallopeptide-based mimics while the disulfide linked dimeric metallopeptide facilitates formation of a dinuclear Ni 2 center.We speculate that dimer formation in {Ni 2 (SOD mds )} is facilitated by the close proximity of the two metal centers, which was enforced through a disulfide bond linkage near the individual nickel centers.The driving force for conversion of the dinuclear Ni 2 II site into two monomeric Ni II sites within {Ni 2 (SOD mds )} likely involves protonation of a terminal cysteinate sulfur atom.The nucleophilic HOMO of the dinuclear Ni 2 II site possesses S(3p) character corresponding to the terminal thiolate sulfur coordinated to the more distorted nickel center, making it the likely protonation site.Protonation of that sulfur atom shortens the Ni-S(H + )-Cys bond length relative to the unprotonated model with a concomitant increase in the bridging Ni-S bond length trans to the protonation site (Figure 10).Association of water to the nickel center followed by breaking of the weakened Ni-S bond would lead to the eventual formation of two mononuclear Ni II centers.

Figure 10 .
Figure 10.Structures of binuclear nickel-site models of {Ni2 II (SOD mds )} and {Ni2 II (SOD mds -S(H + )C1)}.The disulfide bridge and methylene groups have been represented as small spheres and wires for clarity.
[47,48]ics 2019, 7, x FOR PEER REVIEW 8 of 19 two mononuclear nickel sites that collapses into a binuclear cysteinate bridged Ni2 center at high pH, the formation of which is well represented in small molecule Ni-thiolate chemistry[47,48].. 3 1 trans to amidate nitrogen; 2 trans to water oxygen.

Table 2 .
Ni(3d), S(3p), and N(2p) Löwdin MO population analysis (%AO to MO) and energies (E, eV) relative to the HOMO for the LUMO+1 through HOMO-5 of the computational model for the disulfide bridged dinuclear {Ni 2 (SOD mds )} nickel site computational model.Orbital compositions for the doubly-protonated model are given above the unprotonated model for each AO.

Table 2 .
Cont. more distorted nickel site; 2 less distorted nickel site; 3 disulfide sulfur over Ni 1 ; 4 disulfide sulfur over Ni 2 ; 5 bridging thiolate sulfur; 6 bridging thiolate sulfur;7terminal thiolate sulfur ligated to Ni 1 ; 8 terminal thiolate ligated to Ni 2 ; 9 amidate nitrogen ligated to Ni 1 ; 10 amidate nitrogen ligated to Ni 2 . 1 Supplementary Materials:The following are available online at http://www.mdpi.com/2304-6740/7/7/90/s1,FiguresS1andS2:nickelK-edge XANES of the O 2 decomposition products of {Ni 2 II (SOD mds )}, Tables S1-S7: xyz files for computational models, TablesS8 and S9: alternative EXAFS models for the high and low pH forms of {Ni 2 II (SOD mds )}.Conceptualization, J.S.; methodology, B.K. and J.S..; formal analysis, B.K and J.S..; investigation, J.S., D.O. and B.K.; writing, B.K. and J.S.; supervision, J.S.; project administration, J.S.; funding acquisition, J.S.This research was funded by the National Science Foundation, grant number CHE-1565766, the National Institutes of Health, grant number R15-GM120641-01, and Trinity University.Molecular graphics and analyses performed with UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311.Nickel K-edge X-ray absorption spectroscopy was performed at the Canadian Light Source, which is supported by the Canada Foundation for Innovation, Natural Sciences and Engineering Research Council of Canada, the University of Saskatchewan, the Government of Saskatchewan, Western Economic Diversification Canada, the National Research Council Canada, and the Canadian Institutes of Health Research. Funding:Acknowledgments: