
==== Front
Phys Chem Chem Phys
Phys Chem Chem Phys
CP
PPCPFQ
Physical Chemistry Chemical Physics
1463-9076
1463-9084
The Royal Society of Chemistry

39193659
d4cp01989j
10.1039/d4cp01989j
Chemistry
Proton-coupled electron transfer at a mis-metalated zinc site detected with protein charge ladders†‡
† This article is part of a collection celebrating the 85th birthday of Professor George Whitesides. A tribute to George Whitesides, written by the corresponding author of this study, is included at the end of this article.

‡ Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4cp01989j

Gonzalez Mayte a
Guberman-Pfeffer Matthew J. a
Koone Jordan C. a
Dashnaw Chad M. a
Lato Travis J. a
https://orcid.org/0000-0001-8265-5833
Shaw Bryan F. a
a Department of Chemistry and Biochemistry, Baylor University Waco TX USA bryan_shaw@baylor.edu

28 8 2024
11 9 2024
28 8 2024
26 35 2287022881
13 5 2024
12 8 2024
This journal is © the Owner Societies
2024
The Royal Society of Chemistry
https://creativecommons.org/licenses/by/3.0/ This article is licensed under a Creative Commons Attribution 3.0 Unported Licence. You can use material from this article in other publications without requesting further permissions from the RSC, provided that the correct acknowledgement is given.
Distinguishing proton-coupled electron transfer (PCET) from uncoupled electron transfer (ET) in proteins can be challenging. A recent investigation [J. C. Koone, M. Simmang, D. L. Saenger, M. L. Hunsicker-Wang and B. F. Shaw, J. Am. Chem. Soc., 145, 16488–16497] reported that protein charge ladders and capillary electrophoresis can distinguish between single ET, PCET, and two-proton coupled ET (2PCET) by directly measuring the change in protein net charge upon reduction/oxidation (ΔZET). The current study used similar methods to assess PCET in zinc-free, “double copper” superoxide dismutase-1 (4Cu-SOD1), where one copper is bound at the copper site of each monomer and one copper is bound at the bridging zinc site, resulting in a quasi-type III Cu center. At pH 7.4, the net charge (Z) of the 4Cu-SOD1 dimer was unaffected by reduction of all four Cu2+ ions, i.e., ΔZ4ET = −0.09 ± 0.05 per dimer (−0.02 ± 0.01 per copper atom). These values suggest that PCET is taking place at all four Cu atoms of the homodimer. Molecular dynamics and Poisson–Boltzmann calculations suggest that a metal-coordinating histidine at the zinc site (His71) is the proton acceptor. These data show how ligands of a naturally occurring zinc site can help facilitate PCET when the right redox metal is bound.

Distinguishing ET and PCET processes is challenging. This paper used “protein charge ladders” to show that PCET occurs at the non-redox “structural” zinc sight of SOD1, when a Cu ion displaces Zn.

National Science Foundation 10.13039/100000001 CHE: 2203441 Welch Foundation 10.13039/100000928 AA-1854 Baylor University 10.13039/100007492 Unassigned National Institutes of Health 10.13039/100000002 R01EY036371 pubstatusPaginated Article
==== Body
pmcIntroduction

A common method for determining whether an electron transfer (ET) reaction is coupled to proton transfer is to measure the redox potential (E) of the reaction as a function of pH.1–5 Stoichiometric proton-coupled electron transfer (PCET) reactions will yield a slope of approximately −60 mV pH−1 in Pourbaix-style plots of E vs. pH.6–8 Reactions involving electron transfer without proton coupling will yield horizontal lines. This method is not always suitable for metalloproteins because of pH-dependent conformational changes and the need for chemical mediators.9–11 Kinetic methods for detecting PCET also exist.12,13

A new method has been recently established for distinguishing ET and PCET in proteins.14 This new method uses “protein charge ladders” and capillary electrophoresis (CE) to measure how the protein's net charge (Z) fluctuates during redox cycling.14 The goal of the current study is to use this new method to determine if PCET can occur at a naturally occurring zinc site when a Cu2+ ion is bound at the site (instead of Zn2+) and to Cu1+.

A protein charge ladder is an electrostatic array of proteins – commonly prepared by acylation of surface lysine residues – that have systematically altered surface charge but similar shape.15 Originally developed by Whitesides and co-workers to quantify electrostatic effects in molecular recognition and protein net charge,16,17 a protein charge ladder of a metalloprotein can determine its degree of charge regulation upon electron transfer.14 When analyzed by CE, a charge ladder can quantify how the net electrostatic charge of the folded, solvated metalloprotein changes with metal oxidation state. Measuring charge regulation in redox reactions can discern whether a redox reaction involves PCET,18 ET,19 or even two-proton coupled electron transfer (2PCET), where single electron transfer is coupled to the transfer of two protons.14

Charge regulation during redox reactions involving metalloproteins can manifest as a net increase in the pKa of ionizable residues in response to a reduction in metal oxidation state (and vice versa).18–20 Charge regulation can explain why the change in net charge (ΔZET) of a metalloprotein upon single electron transfer falls between −1 ≤ ΔZ ≤ +1.19 Direct measurements of charge regulation during PCET and ET in proteins have been made for heme, Rieske, and type I and II copper metalloproteins using protein charge ladders and capillary electrophoresis.14,15,18,19 These studies show that protonation can occur: (i) stoichiometrically at one or two metal binding residues, with no detectable ΔZET of other residues,18 or (ii) sub-stoichiometrically (fractionally) across multiple non-coordinating residues that are distal to the active site.19

For example, in the [2Fe–2S] Rieske protein, protonation and ET appear to be perfectly stoichiometric; upon reduction of one Fe3+ to Fe2+, proton transfer occurs at zero, one, or two histidine residues as solvent pH increases from 5 to 10.14 Thus, at high pH, a Rieske protein can actually become more positively charged upon reduction of Fe3+, a consequence of 2PCET. With human Cu, Zn superoxide dismutase (hSOD1), PCET is also stoichiometric, involving proton transfer to a histidine that bridges Cu and Zn (His63), with no net change in protonation of other residues (Fig. 1).18 In the type I copper protein azurin, single electron transfer is coupled (to some degree) to the partial protonation of a handful of non-coordinating residues that are distal to the active site.19 It has been noted that the free energy associated with the protonation of these residues in azurin could account for a large fraction of the total reorganization energy measured by Gray and co-workers.21 This type of fractional, sub-stoichiometric protonation is also observed with myoglobin.19

Fig. 1 Electron transfer in the Cu, Cu derivative of Cu, Zn SOD1 (denoted CuCu, CuZn SOD1). Proton coupling with electron transfer at the canonical copper site is well established (His63 becomes protonated upon reduction of Cu2+ to Cu1+). Is a proton also coupled with electron transfer to copper bound at the canonical “structural” zinc site (CuZn)?

Materials and methods

Purification and characterization of SOD1

Wild-type (WT) superoxide dismutase 1 (SOD1) was recombinantly overexpressed in S. cerevisiae. SOD1 was purified as previously reported.22,23 Yeast stocks were made by transfecting Yep351-hSOD1 plasmids into EG118Δsod1 yeast. This was then followed by growing primary cultures to OD600nm ∼ 1.5 (approximately 36 hours) in YPD media and transferring into larger cultures for ∼7 days. The grown cells were spun down, lysed, and purified through ammonium sulfate precipitation followed by three chromatographic separations: hydrophobic interaction chromatography, ion exchange chromatography, and size exclusion chromatography. SOD1 purity was confirmed via SDS-PAGE following each chromatographic separation, and the final protein solution was characterized by mass spectrometry.

Demetalation of Cu, Zn SOD1

Purified SOD1 solutions were demetalated via sequential dialysis over 6 days in: (i) 100 mM ammonium acetate, 5 mM EDTA, pH 3.8; (ii) 100 mM ammonium acetate, 100 mM NaCl, pH 3.8; and (iii) 100 mM ammonium acetate, pH 5.5. The protein was in each buffer for ∼2 days with buffer replacement every ∼8 hours. To prevent metal contamination, all glassware was rinsed with 10 mM EDTA, followed by Milli-Q water. Inductively coupled plasma-mass spectrometry (7900 ICP-MS, Agilent Technologies) was used to verify full demetalation of SOD1 to be <0.08 equivalents of copper and zinc (Table S1, ESI‡).

Acetylation of apo-SOD1 to yield a “protein charge ladder”

The apo-SOD1 was acetylated with acetic anhydride as previously described to generate protein charge ladders.24 Briefly, the protein was washed into 100 mM HEPBS (pH 9) via centrifugal filtration to a final concentration of ∼150 μM. 100 mM acetic anhydride was made in 1,4-dioxane. Three molar equivalents of acetic anhydride were added to the protein solutions to generate a protein charge ladder. Extent of acetylation of SOD1 was confirmed through electrospray ionization-mass spectrometry (ESI-MS) using a ThermoFisher™ Discovery Orbitrap mass spectrometer. Protein samples were prepared by diluting with 0.2% formic acid and removing salt with a desalting column (Michrom BioResources, Inc., Auburn CA, USA) followed by eluting with 80% acetonitrile, 20% of 2% formic acid. Acetylated SOD1 was transferred into the appropriate running buffer for CE experiments through centrifugal filtration (4000 × g at 4 °C).

Remetalation of apo-SOD1 to 4Cu-SOD1 and reduction with sodium dithionite

Acetylated apo-SOD1 WT (i.e., the apo-SOD1 “protein charge ladder”) was titrated slowly with up to five molar equivalents of 50 mM of CuSO4 over 30 minutes. Solutions were incubated at 4 °C for a period of 2 days with metal binding and metal content verified with UV-vis and ICP-MS (Table S1, ESI‡). Protein charge ladders of reduced 4Cu-SOD1 were prepared by addition of 160 molar equivalents of sodium dithionite to the oxidized Cu–Cu SOD1 protein charge ladder.

Capillary electrophoresis

All CE experiments were performed using a Beckman P/ACE MDQ instrument with a bare fused-silica capillary. Electrophoresis was performed at 29 kV in corresponding buffers and the capillary was cooled to 22 °C to prevent Joule heating. The concentration of SOD1 protein charge ladders analyzed with CE was 30 μM SOD1 dimer. Dimethylformamide (DMF) was added as a neutral marker of electroosmotic flow (EOF). Capillary conditioning was performed in between experiments with washes of 0.1 M HCl (2 min), methanol (2 min), 0.1 M KOH (2 min), Milli-Q water (2 min), and running buffer (4 min). All protein concentrations were determined via UV-vis spectrometry ε280: 10 800 cm−1 M−1 for apo-SOD1 and ε265: 18 400 cm−1 M−1 holo-SOD1. For experiments involving protein charge ladders of reduced Cu, Cu SOD1, the running buffer used was deoxygenated via nitrogen bubbling for at least 3 hours prior to electrophoresis. All reduced protein solutions were analyzed with CE immediately following the addition of sodium dithionite.

Electrophoretic mobility was calculated by using eqn (1) to convert the migration times of unacetylated SOD1 and each rung of the charge ladder for each electropherogram.1

The equation represents the mobility of protein species defined by μ.

Computational methods overview

Nearly 3 μs of production-stage classical molecular dynamics, with and without variable protonation of titratable residues, and hundreds of Poisson–Boltzmann electrostatic energy calculations were performed for homodimeric CuCuCuZn-SOD1. The simulations comprised two stages: (1) evaluating the change in pKa for titratable residues in the second coordination sphere upon reduction of CuZn, (at this stage, the Cu was assumed to retain the H63, H71, H80, and D83 ligands in both redox states); and (2) assessing the energetic change upon protonating each of the four ligands to the reduced Cu in the Zn-binding site. All simulations were performed assuming the Cu in the Cu site was already reduced, and this Cu had lost coordination to both a water molecule and His63.

Metal site parameterization

A prerequisite for simulations of 4Cu-SOD1 was the derivation of molecular mechanics parameters for CuCu1+, CuZn2+, and CuZn1+. The parameters were derived with the assistance of the Metal Center Parameter Builder (MCPB)25 of the AmberTools22 Suite.26,27 Bonded parameters were obtained with the Seminario method28 from quantum mechanically (QM) optimized models including the Cu ion and coordinating ligands up-to-and-including the β-carbons. The QM optimization and frequency calculations were performed with the Becke, three-parameter, Lee–Yang–Parr (B3LYP) approximate density functional29,30 and a mixed basis set (LANL2DZ for Cu;31 6-31G(d) for H, C, N and O),32–34 as implemented in Gaussian 16 Rev. A.03.35 Merz–Singh–Kollman atomic partial charges for the Cu ion and coordinating ligands were obtained with the aforementioned model chemistry and the restrained electro-static potential (RESP) approach. Backbone (heavy and hydrogen) atoms of the coordinating ligands were held to the standard force field values (option 3c in MCPB). The metal center parameters were developed for use with the ff10 AMBER forcefield36 to be compatible with constant pH molecular dynamics simulations. The library and forcefield modification (frcmod) files needed for the simulations are provided as ESI‡ at the following GitHub repository: https://github.com/Mag14011/CuCuSOD1_SI_CompMethods.

Structure preparation

A dimer (chains A and F) of the highest (1.07 Å) resolution X-ray crystallographic structure for Cu–Zn human superoxide dismutase 1 (SOD1) available in the Protein Data Bank (accession code 2C9V)37 was selected as the starting point for the modeling in the present work. Preparatory steps included: (1) replacing Zn with Cu in the canonical Zn binding site, (2) replacing chain F with a superimposed copy of chain A (because there is a >3 Å break in chain F), (3) acetylating the N-terminus of each chain, (4) selecting the highest occupancy position when alternate coordinates were present, except for the Cu ion in the Cu site for which the low occupancy position corresponded to the reduced state, and (5) changing the residue names of metal-center ligands to match the above-described parametrization, and (6) changing the names of second coordination sphere ASP, GLU, and HIS residues to ASH, GLH, and HIP so these residues (as well as TYR and LYS for which a name change was not needed) can be titrated in constant pH molecular dynamics (CpHMD).

4Cu-SOD1 was simulated in the Cu1+–Cu2+ and Cu1+–Cu1+ redox states as a function of pH from −4 to 12 in increments of 2 pH units. 4Cu-SOD1 was also simulated in the Cu1+–Cu1+ state with each ligand (His63, His71, His80, and Asp83) to the Cu ion in the canonical Zn site either unprotonated or protonated. For the unprotonated cases, the bonded parameters between the ligand and Cu ion were deleted, but the electrostatic interaction was retained, essentially adopting a non-bonded or crystal-field-theory approximation. This procedure was necessary because it is not possible to simulate, at the classical mechanics level, the breaking of the coordination bond as the ligand becomes protonated.

tLEaP in the AmberTools22 package was used to add hydrogen atoms to the various SOD1 models. Each structure was placed at the center of a box of explicit water with at least a 15 Å buffer region to the boundary of the box. A sufficient number of counterions (6 and 8 Na+, respectively, for the Cu1+–Cu2+ and Cu1+–Cu1+ states) was added to achieve charge neutrality. The TIP3P water model38 and the monovalent ion parameters of Joung and Cheatham39 were used to model the solution state.

Molecular dynamics

Each solvated structure was subjected to 10 000 steps of steepest descent followed by 40 000 steps of conjugate-gradient minimization, all with a 10 kcal (mol Å2)−1 restraint on the heavy atoms of the protein backbone. Each system was subsequently heated in the NVT ensemble from 0 to 300 K at a rate of 0.3 K ps−1 and held at the final temperature for 3.0 ns. A 1.0 kcal mol−1 Å−2 restraint on the protein backbone was applied during the heating process and first 1.0 ns at the final temperature, but then reduced to 0.1 kcal mol−1 Å−2 for the remaining 2.0 ns. After thermalization, the density of each system was equilibrated under 1.0-bar pressure for 2.0 ns in the NPT ensemble. Production stage simulations were conducted in the NVT ensemble at 300 K. The lengths of the various production simulations are shown in (Table S3, ESI‡).

All NVT and NPT simulations (including the production-stage) employed periodic boundary conditions, the particle mesh Ewald40 treatment of electrostatic interactions with a direct sum cut-off of 10.0 Å, the SHAKE algorithm41,42 to rigidify bonds to hydrogen atoms, a Langevin thermostat with a collision frequency of 2 ps−1, and an integration timestep for the Langevin equation of motion of 2.0 fs. Pressure in NPT simulations was regulated with a Monte Carlo barostat having a relaxation time of 1.0 ps. PMEMD in its CPU and GPU43 implementations of the Amber22 package27 were used to perform the minimization and dynamical simulations, respectively.

Poisson Boltzmann electrostatic calculations

Configurations from the last 100 ns of each simulation in which a ligand to Cu1+ in the canonical Zn site was unprotonated or protonated were selected, stripped of the explicit water and ions, and submitted to the Poisson–Boltzmann surface area (PBSA) program of the AmberTools22 suite. The configurations were sampled every 1.0 ns, giving a total of 100 analyzed frames. The statistics of the total energy before and after protonating a given ligand were compared to assess which ligand is most likely to become protonated when the Cu in the Zn site is in the reduced (1+) state. The BASH script used for this post-simulation analysis is provided as ESI‡ in the following GitHub repository: https://github.com/Mag14011/CuCuSOD1_SI_CompMethods.

Results and discussion

In this paper, a pseudo-type III copper center was generated in homodimeric Cu, Zn-SOD1 by removing Zn2+ from the Zn site of each monomer and replacing it with Cu2+ (Fig. 1). This Cu–Cu derivative is denoted CuCuCuZn-SOD1 or 4Cu-SOD1.44–46 A lysine-acetyl protein charge ladder of this pseudo-type III copper derivative of SOD1 was created (Fig. 2) to determine whether it undergoes 2PCET, ET, or a mixture of the two processes upon reduction of all four Cu2+ ions in the 4Cu-SOD1 homodimer. Our previous studies of PCET in metalloproteins (with protein charge ladders) investigated the extent of PCET at naturally occurring redox sites that are properly metalated with redox active metals.14,18,19 In contrast, the current study uses protein charge ladders to determine whether PCET occurs at a “mis-metalated” non-redox active site (the zinc site of SOD1) when a redox active metal (copper) has displaced the natural metal (zinc).

Fig. 2 (A) Lysine-acetyl “protein charge ladders” of human wild type SOD1 synthesized by acetylation of surface lysine with acetic anhydride. (B) and (C) Raw electrospray ionization mass spectrum (B) and deconvoluted spectrum (C) for WT apo-SOD1 charge ladders. Note: the SOD1 homodimer dissociates under these ionization conditions (but remain dimeric during capillary electrophoresis).

Copper centers can typically be classified by stoichiometry, ligand type, and geometry of the metal site.47 Type III Cu centers are characterized by: (i) two copper ions (usually separated by ∼3.5 Å with a bridged substrate), (ii) a silent spectroscopic signature due to antiferromagnetic coupling,48,49 and (iii) coordination of both copper atoms by three histidine residues. In general, naturally occurring type III Cu proteins (e.g., multicopper oxidases) can perform twin PCET per se (i.e., PCET at both copper ions) by coupling proton transfer to an electron through bridging dioxygen substrates that become oxidized.50,51

The properly metalated Cu, Zn SOD1 protein is not a type III copper protein. Rather, it contains a unique type II copper center that shares a bridging histidine with a zinc atom.52 Type II copper centers consist of a single copper atom coordinated by four ligands (N or N/O) in a square planar geometry. The copper site in SOD1 is five coordinate in the Cu2+ oxidation state, coordinated by four histidine residues (one of which, His63, bridges both Cu and Zn sites) and one H2O (or O2− substrate). This copper site adopts a distorted square pyramid configuration (Fig. 1).49 The copper center transitions to a three coordinate distorted trigonal plane upon reduction – which involves PCET – wherein the bridging His63 is protonated upon reduction (Fig. 1).53 The zinc site contains three histidine and one aspartate that adopt a nearly tetrahedral coordination geometry when zinc is bound.54

The 4Cu-SOD1 derivative was prepared by removing all coordinated metal ions from the recombinant as-isolated human wild-type (WT) SOD1 dimer (using EDTA and pH 3.8 buffer). The removal of Cu and Zn was verified with inductively coupled plasma-mass spectrometry (ICP-MS; Table S1, ESI‡). Protein charge ladders of apo-SOD1 were then formed by acetylating its lysine residues with acetic anhydride, prior to remetalation with copper. Each acetylation increases the mass by 42 Da (Fig. 2A–C). Acetylation of ∼3 lysine residues on each SOD1 subunit (6 lysine per dimer) has a negligible effect on the structure of SOD1, according to previous analysis with hydrogen/deuterium exchange.55 It is important to remember that the SOD1 dimer dissociates during ESI-MS. The number of acetylated lysine inferred from the mass spectra correspond to each SOD1 monomer. Note that the recombinant SOD1 protein was expressed in S. cerevisiae, i.e., the N-terminus is also properly acetylated.

After acetylation of apo-SOD1 to form a lysine-acetyl protein charge ladder, four Cu2+ (per dimer) were titrated back into the lysine-acetyl charge ladder of apo-SOD1 dimer (Fig. 3), i.e., copper was added after lysine acetylation. Each metal binding site of the 4Cu-SOD1 derivative was occupied by copper, according to UV-vis spectroscopy (Fig. 3). Here, the sequential binding of Cu2+ to the SOD1 active site can be observed by the copper d–d transition at 500–800 nm (Fig. 3). Four molar equivalents of Cu2+ (per dimer) were added sequentially (1 eq. at a time) and monitored via UV-vis to ensure that Cu2+ bound to the active site (Fig. 3). Here, the expected increase in the d–d region was observed, indicating coordination Cu2+ (Fig. 3). ICP-MS verified copper stoichiometry (Table S1, ESI‡).

Fig. 3 UV-vis spectra of d–d band in a lysine-acetyl protein charge ladder of human WT SOD1 as four stoichiometric equivalents of Cu2+ (per dimer) are titrated into the metal-free homodimeric charge ladder over 48 h, pH 7.4 (blue traces). The red spectrum is of the copper-replete 4Cu-SOD1 protein after reduction with 160 stoichiometric equivalents of sodium dithionite (per SOD1 dimer).

CE and protein charge ladders were used to directly measure the net charge (Z) of 4Cu-SOD1 in both the oxidized Cu2+ and reduced Cu1+ state at multiple pH values.15,18,19,56 To ensure the maintenance of full copper metalation during CE, the solution of 4Cu-SOD1 protein charge ladder contained an extra molar equivalent of Cu2+ (i.e., five molar equivalents were added per dimer).

When separated by CE, a distribution of peaks (rungs) is observed for oxidized 4Cu-SOD1 (Fig. 4A). As expected, the electropherogram of the charge ladder shows a different distribution of peaks than the electrospray ionization mass spectrum. This difference is observed because electrospray ionization monomerizes SOD1 (complete dissociation of the SOD1 dimer), whereas the SOD1 dimer remains intact during CE. Thus, the third rung of the ladder (“Ac 2”, a dimer with two acetylations) can be comprised of either one monomer with zero acetylations bound to a monomer with two acetylation, or two monomers bound together, each with one acetylation.

Fig. 4 Use of capillary electrophoresis (CE) to measure change in net charge (ΔZET) of Cu–Cu SOD1 upon reduction of all four Cu2+ ions. (A) Electropherograms of protein charge ladders of oxidized (blue) and reduced (red) CuCuCuZn SOD1 at pH 7.4 (i.e. Cu2+, Cu2+ and Cu1+, Cu1+) are shown as overlays of 6 technical replicates. (B) Plot of the electrophoretic mobility (μ) of each rung versus the number of acetylated lysine(s) Ac(N). The x-intercept is equal to the quotient of the net charge (Z) of the unacetylated protein and the change in charge with each acetylation (i.e., ZAc(0)/ΔZAc). The charge of the reduced and oxidized proteins are indicated by Zred and Zox. (C) UV-vis spectra of a lysine-acetyl protein ladder of CuCuCuZn SOD1 used for CE analysis, before and after reduction of all four Cu2+ in the homodimer to Cu1+ using sodium dithionite and deoxygenated buffer.

The Z of the unmodified protein is determined by extrapolating the linear trend of electrophoretic mobility (μ) and acetylation number Ac(N) (Fig. 4B). In this paper, the first five to six rungs of each charge ladder were used to determine the Z of oxidized and reduced 4Cu-SOD1. CE only requires nanoliters of sample and has a run time of less than 15 minutes, which allows collection of multiple replicates (>6 technical replicates) to yield statistically significant values.18,19

The net charge of 4Cu-SOD1 does not change upon reduction at physiological pH

At pH 7.4, the net charge of the oxidized 4Cu-SOD1 dimer was measured to be Z = −6.35 ± 0.05 (Fig. 4A and B). This value differs slightly from a value measured previously for the properly metalated Cu, Zn SOD1 dimer (i.e., Z = −7.32 ± 0.17).18 This difference might be due to actual differences in net charge of the Cu–Cu and Cu–Zn SOD1 proteins (due to different residue pKa, or coordinated solvent/co-solvent ions) or it might be caused by the different buffer conditions. In the current study, 50 mM 3-morpholinopropane-1-sulfonic acid buffer (MOPS) was used, whereas potassium phosphate was used in prior studies.15

The reduction of all Cu2+ ions in the copper loaded SOD1 was accomplished by the addition of 160 equivalents of sodium dithionite (Fig. 4A). Upon complete reduction, Z changed from −6.35 ± 0.05 to −6.44 ± 0.05, i.e., ΔZ of −0.09 ± 0.05 per dimer or −0.02 per copper (Fig. 4B). This value suggests that 4Cu-SOD1 undergoes complete charge regulation at pH 7.4 upon reduction of all four copper ions. Previous studies have confirmed – as does this one – that both copper atoms in each subunit can be reduced with dithionite.14,18,19 UV-vis spectra were obtained for the oxidized SOD1 samples prior to CE (Fig. 4C).18,57

Charge regulation in 4Cu-SOD1 is minimally affected by pH

Charge was also measured at pH 5.6 and 8.3 (the physiologically relevant pH range of SOD1) to assess effects of pH on charge regulation in 4Cu-SOD1.

Electrophoresis was performed at pH 5.6 in 10 mM sodium acetate (Fig. 5A). Analysis of the electropherograms showed that the Z of 4Cu-SOD1 was −3.92 ± 0.15 per dimer, and −4.30 ± 0.17 after reduction with sodium dithionite (ΔZ4ET = −0.38 ± 0.23 or ΔZET = −0.09 ± 0.23 per electron) (Fig. 5B). The data suggests that charge regulation weakened (i.e., magnitude of ΔZ became larger) but persisted at lower pH. The percent of charge regulation per electron was 91%, 7% less than the charge regulation observed at pH 7.4. The lower resolution of the electropherograms at pH 5.6 is caused (in part) by the slower electroosmotic flow from protonation of silanol groups (SiOH) in the capillary, making measurements at pH ≤ 5 difficult with CE (Fig. 5A).

Fig. 5 (A) Capillary electropherograms of protein charge ladders of oxidized (blue) and reduced (red) CuCuCuZn SOD1 at pH 5.6 are shown as overlays of 6 technical replicates. (B) The plot of the electrophoretic mobility (μ) of each rung versus the number of acetylated lysine Ac(N) at pH 5.6. (C) Electropherograms of protein charge ladders of oxidized (blue) and reduced (red) of CuCuCuZn SOD1 at pH 8.3 are shown as overlays of 6 technical replicates. (D) The plot of mobility of each rung versus Ac(N) for each rung at pH 8.3. Values of the charge of the reduced and oxidized proteins are indicated by Zred and Zox.

For experiments at pH 8.3, CE experiments were performed in 20 mM Tris-Gly (Fig. 5C). Charge regulation was not as pronounced at pH 8.3, with ΔZ4ET = −0.51 ± 0.31 (ΔZET = −0.13 ± 0.14 per electron) (Fig. 5D). The percent of charge regulation per electron was 87% (i.e., 11% less than the charge regulation taking place at pH 7.4 and 4% less than at pH 5.6) (Table 1). These values suggest that PCET is occurring at pH 8.3.

The measured net charge (ZCE) of hSOD1, and the change in net charge upon reduction (ΔZET) of all four Cu2+ ions in each dimer

Metal redox state	pH	Z CE a	ΔZa (Zox − Zred)	ΔZ per electron	% Charge regulation per electron	
Cu2+	5.6	−3.92 ± 0.15	−0.38 ± 0.23	−0.09 ± 0.23	91 ± 6	
Cu1+	5.6	−4.30 ± 0.17	
Cu2+	7.4	−6.35 ± 0.05	−0.09 ± 0.05	−0.02 ± 0.05	98 ± 1	
Cu1+	7.4	−6.44 ± 0.05	
Cu2+	8.3	−7.19 ± 0.06	−0.50 ± 0.14	−0.13 ± 0.14	87 ± 4	
Cu1+	8.3	−7.69 ± 0.30	
a Values are listed per dimer.

Note that electropherograms for reduced 4Cu-SOD1 proteins (Fig. 4 and 5) appear to exhibit slightly lower resolution compared to electropherograms of oxidized proteins. In general, peak resolution in CE can be dependent on three major factors: (i) the length of the sample “plug”; (ii) diffusion between the time of injection and detection; and (iii) interactions between the protein and the surface of the fused silica capillary.58–60 Adhesion of protein to the capillary is the likely cause of the lower resolution of reduced SOD1.

It is also important to comment on the error associated with replicate electropherograms of some of the protein charge ladders. There is, for example, larger variation in mobility values for reduced charge ladders in Fig. 5A compared to Fig. 4A. The variable mobility values in Fig. 5A do not lead to large variations in the final, calculated net charge. The data set of electropherograms in Fig. 4A (copper-reduced) consists of nearly superimposed electropherograms. Here the mobility of rung 2 (i.e., Ac(2)) varied from μ = 5.87 to 5.98 cm2 kV−1 min−1 across replicate measurements. This variation results in a small variation in Z, i.e., Z = −6.52 to −6.43. In contrast, the copper-reduced electropherograms of Fig. 5A are more variable. Here, rung 2 varied from μ = 4.93 to 5.27 cm2 kV−1 min−1. This results in larger (but still small) variation in Z, i.e., from Z = −4.26 to −4.51.

PCET during CuZn reduction: where is the proton going?

Our previous report found that properly metalated SOD1 (Cu2Zn2-SOD1) undergoes complete charge regulation upon ET.18 This complete charge regulation was due to PCET via protonation of a bridging histidine upon reduction of Cu2+.61,62 We expected the value for ΔZET of 4Cu-SOD1 to be approximately −2 per dimer, accounting for single ET at the two CuZn and PCET at each CuCu. That is, the CuZn site was not expected to participate in PCET but rather single, uncoupled ET.

Classic electrochemical studies of SOD1 support the hypothesis that Cu2+ bound to the zinc site engages in PCET.3 For example, the Pourbaix diagram of a copper SOD1 derivative where the zinc site of each monomer is empty (Cu2E2 SOD1) suggests proton coupling is occurring at pH > 7,3i.e., a pH range where Cu is known to migrate from the copper site to the empty zinc site.63 The slope of mV pH−1 at pH > 7 for this CuZn derivative is, however, only −18 mV pH−1.3

The likely occurrence of PCET at the zinc site raises the question: which amino acid residues are being protonated during reduction of CuZn? Is protonation occurring on one residue, or multiple residues? Are these residues at, near, or distal to the Zn site?

The residue responsible for protonation during PCET of CuCu in SOD1, His63, is thought to exist as a doubly deprotonated imidazolate.64,65 There are several possible sites that could participate in PCET at CuZn: (i) the bridging H63 imidazolate, which would become doubly protonated, (ii) a CuZn coordinated residue becoming protonated (H71, H80, or D83 might become protonated), or (iii) residue(s) distal to the active site becoming protonated as electrostatic interactions span ∼10 Å in physiological conditions and are longer in hydrophobic environments. Distal ionizable residues surrounding the CuZn site might adjust pKa.15 Within 20 Å of the CuZn there are a total of 46 ionizable residues which may contribute to diffuse charge regulation following ET (Fig. 6). There are 9 residues located from 4 Å to 10 Å (H46, K70, R79, D124, K136, E133, H48, R69, S134). Here, there is 1 located within 4 Å (K136), 2 within 5 Å (H46, K136), and 5 within 6 Å (H46, K70, R79, D124, K136). These distal ionizable residues may be located near the active site to participate in diffuse charge regulation, where pKa values of multiple residues change.

Fig. 6 Protic residues within 20 Å of the zinc site in SOD1 that might be potential proton acceptors during electron transfer (PDB: 2C9V). The residues that coordinate metals in the zinc site are shown in yellow (H63, H71, H80, and D83). Ionizable residues are shown within 8 Å (green), 12 Å (light green), and 20 Å (grey). The red dashed circle denotes the classical Debye radius in physiological buffer (10 Å, but which can be much longer in hydrophobic interiors).

The most reasonable hypothesis is that protonation is taking place at the CuZn coordinating ligands (H71, H80, D83) (Fig. 7). Since the reduction of Cu2+ to Cu1+ induces a change in coordination geometry (i.e., tetrahedral to trigonal planar), theoretically, it is possible that H71 could become uncoordinated from Cu(i) where it could participate in a hydrogen bonding network that involves D124.66

Fig. 7 Five possible sites of protonation during electron transfer to Cu2+ bound in the zinc site of SOD1. Four sites include: H80, D83, H71, H63 and a fifth “site” includes one or more residues in the second coordination sphere (e.g., D76, D124, or D125) or residues that are more distal to the metal site.

Determining which residues are participating in PCET at the CuZn site can be experimentally challenging. Therefore, a two-stage computational approach based on classical molecular dynamics simulations to make an approximation of which residues are becoming protonated upon reduction of CuZn was conducted.

First, the change in pKa upon reduction of the Cu2+ ion in the Zn binding site was assessed for the 10 Asp, 10 Glu, 2 His, and 11 Lys residues in the second coordination sphere of each subunit of 4Cu-SOD1 (Table S2, ESI‡). From this analysis only four residues were predicted to exhibit pKa changes ≥1.0 unit: D101 (1.4 in one subunit, but 0.7 in the other), D124 (2.6–2.9 units in both subunits), and D125 (1.0 in one subunit, but −0.4 in the other). Furthermore, the pKa's of all these residues remained at least 2 units below the lowest pH examined experimentally, thereby suggesting that second-sphere titratable residues are not the site of charge regulation (protonation) upon reduction of CuZn.

Given this null result, the second stage of computations involved assessing the energetics of protonating each of the four primary coordination sphere ligands to the CuZn site: H63, H71, H80, and D83. Dynamical trajectories were performed in the fully reduced state for 4Cu-SOD1 with only three of these four residues explicitly bonded to the Cu1+ ion in the Zn binding site; the fourth possible ligand was simulated as either unprotonated or protonated. Since it is not possible to simulate the breaking of a ligand–metal bond in a classical mechanical simulation, we adopted the crystal field theory approximation that the ligand–metal interaction was entirely electrostatic in nature. This assumption allowed us to ask the question: is it electrostatically more favorable for one of these four ligands to interact with the Cu1+ center or to become protonated?

Energetic assessments of configurations from these trajectories with the Poisson–Boltzmann surface area (PBSA) method gave the following rank ordering of the four primary coordination sphere residues in terms of increasing stability upon protonation when CuZn is reduced: D83 (+31 kcal mol−1) < H80 (−4 kcal mol−1) ≈ H63 (−6 kcal mol−1) < H71 (−28 kcal mol−1) (Fig. 8). We therefore propose that H71 is the site of proton uptake upon reduction of the Cu in the Zn binding site of SOD1.

Fig. 8 Change in relative energies of copper ligands (of the zinc site) upon protonation in the Cu1+ oxidation state. In this simulation, both CuZn and CuCu are in the +1 oxidation state. Energies represent computed averages over thermal ensembles of configurations sampled during molecular dynamics trajectories and assessed with Poisson–Boltzmann calculations with the protein explicitly included.

Conclusion

This study has shown that a naturally occurring non-redox “structural” metal site – the zinc site of Cu, Zn SOD1 – has the necessary structural and chemical elements to undergo PCET when a redox active metal (copper) is bound. This near zero value of ΔZET for this site suggests that the entire protein framework of the 4Cu-SOD1 dimer is insensitive to the redox state of its four coordinated copper ions. It is possible that the four protons that are coupled to ET in the copper replete homodimer might act as electrostatic shields, to prevent distal non-coordinating ionizable residues within the Debye radius from having pKa shifts upon ET. This insensitivity of net charge to metal oxidation state might be important for the diffusion-limited reduction and oxidation of superoxide radicals in the SOD1 protein.

With redox reactions involving metalloproteins, it remains analytically challenging to determine the outer sphere reorganization of protons, solvent, co-solvent and buffer, and the adjustment of side chain charge. This analytical gap makes it hard to answer “how” and “why” metalloproteins do or do not regulate their Z in response to ET (or couple proton(s) to electron transfer). Protein charge ladders can – at least – inform us about the net change in charged groups in protein ET reactions. The results of this study continue to demonstrate the utility of this technique for probing ET and PCET processes in metalloproteins. The current results for 4Cu-SOD1 adds another data point to the small but growing trend of ΔZET and protein Z for metalloproteins (Fig. 9). The results suggest that: (i) proteins engaging in PCET can have wide ranges of net negative charge (from −5 to −10); (ii) proteins engaging in PCET and 2PCET are typically more negatively charged than proteins such as cytochrome c that engage in proton uncoupled, single ET. This trend points towards functional links between the net charge of metalloproteins and the type of redox reactions they perform in vivo.

Fig. 9 Plot of Zoxvs. ΔZET for all metalloproteins whose values of ΔZET and Z have been measured. * Denotes that the change in charge is expressed per 1e−.

George Whitesides tribute by Bryan F. Shaw

True story: George Whitesides paid for my wedding tux. In full disclosure, George bought the tuxedo for me to wear at “table one” when he was awarded the Priestley medal. I was a post-doc in the lab and we had drawn straws to see who would represent the lab at the ceremony. Buying an old tux was cheaper than renting a newer one. At the Priestley dinner, I sat between George's advisor, John D. Roberts and Robert's wife, Edith. After George left the table to receive the award, I couldn’t help myself. “Professor Roberts,” I asked, “what was George like as a grad student?” Professor Roberts – who won the Priestley two decades earlier – said in his enthusiastic tenor, “Well, you see, with George it was like I was the graduate student and he was the graduate advisor.” Feeling out of league, I changed the subject to marriage advice. Professor Roberts shared his unique wisdom. I wore that same tux a few months later when my fiancée Lizz and I were married. Soon after, our lives would take a tragic turn. Our first child, Noah, was born with aggressive tumors in both of his eyes (bilateral retinoblastoma). I remember leaving my desk in George's lab to meet Lizz at the ophthalmologist because “something is wrong with Noah's eyes.” Nothing has ever been the same. When I first told George the bad news, he consoled me and kindly asked “do you need any money?” We made it through treatment but Noah lost most of his vision. When I finally left Harvard to start my own lab at Baylor University, I modified my research focus. We started working on devices to detect eye disorders in children, and then on assistive technology to help students with blindness learn science. The tools that we make have helped children with eye disorders receive early diagnoses, and have helped them be included in science. Much of the courage to dare to do these different, un-chemical things, came from watching George and the amazing members of his lab. Seventeen years later, I am still married and Noah is now 16. Thank you George, for caring about everything and everyone, all the time, as much as you do. Happy 85th birthday!

Author contributions

MG: conceptualization, formal analysis, investigation, visualization, writing – original draft. MJG-P: data curation, investigation, writing – original draft. JCK: verification. CMD: visualization TJL: investigation. BFS: conceptualization, funding acquisition, methodology, resources, supervision, writing original draft. All authors contributed to writing, review and editing.

Data availability

The BASH script used for the computational analysis is provided as ESI‡ in the following GitHub repository: https://github.com/Mag14011/CuCuSOD1_SI_CompMethods.

Conflicts of interest

The authors declare no conflict of interest.

Supplementary Material

CP-026-D4CP01989J-s001

The authors acknowledge the National Science Foundation (CHE: 2203441), the National Institutes of Health (R01EY036371) and the Robert A. Welch foundation (AA-1854) for funding. MGP acknowledges the Office of Vice Provost research fellowship program at Baylor University for support.
==== Refs
Manner V. W. Mayer J. M. J. Am. Chem. Soc. 2009 131 9874 9875 10.1021/ja902942g 19569636
Azab H. A. Banci L. Borsari M. Luchinat C. Sola M. Viezzoli M. S. Inorg. Chem. 1992 31 4649 4655 10.1021/ic00048a037
Agarwal R. G. Wise C. F. Warren J. J. Mayer J. M. Chem. Rev. 2022 122 1482 10.1021/acs.chemrev.1c00791 34637292
Cui K. Soudackov A. V. Kessinger M. C. Xu J. Meyer G. J. Hammes-Schiffer S. J. Am. Chem. Soc. 2023 145 19321 19332 10.1021/jacs.3c05535 37611195
Armstrong F. A. Biochem. Soc. Trans. 1999 27 206 210 10.1042/bst0270206 10093735
Horvath S. Fernandez L. E. Appel A. M. Hammes-Schiffer S. Inorg. Chem. 2013 52 3643 3652 10.1021/ic302056j 23477912
Agarwal R. G. Coste S. C. Groff B. D. Heuer A. M. Noh H. Parada G. A. Wise C. F. Nichols E. M. Warren J. J. Mayer J. M. Chem. Rev. 2022 122 1 49 10.1021/acs.chemrev.1c00521 34928136
Armstrong F. A. Heering H. A. Hirst J. Chem. Soc. Rev. 1997 26 169 179 10.1039/CS9972600169
Frew J. E. Hill H. A. Eur. J. Biochem. 1988 172 261 269 10.1111/j.1432-1033.1988.tb13882.x 3280307
Hill H. A. Hunt N. I. Methods Enzymol. 1993 227 501 522 8255235
Armstrong F. A. Allen H. Hill O. Walton N. J. Q. Rev. Biophys. 1985 18 261 322 10.1017/S0033583500000366 3024201
Brenner S. Heyes D. J. Hay S. Hough M. A. Eady R. R. Hasnain S. S. Scrutton N. S. J. Biol. Chem. 2009 284 25973 25983 10.1074/jbc.M109.012245 19586913
Dempsey J. L. Winkler J. R. Gray H. B. Chem. Rev. 2010 110 7024 7039 10.1021/cr100182b 21082865
Koone J. C. Simmang M. Saenger D. L. Hunsicker-Wang L. M. Shaw B. F. J. Am. Chem. Soc. 2023 145 16488 16497 10.1021/jacs.3c03006 37486967
Gitlin I. Carbeck J. D. Whitesides G. M. Angew. Chem., Int. Ed. 2006 45 3022 3060 10.1002/anie.200502530 16619322
Gao J. Gomez F. A. Harter R. Whitesides G. M. Proc. Natl. Acad. Sci. U. S. A. 1994 91 12027 12030 10.1073/pnas.91.25.12027 7991578
Gao J. Mammen M. Whitesides G. M. Science 1996 272 535 537 10.1126/science.272.5261.535 8614800
Zhang A. Y. Koone J. C. Dashnaw C. M. Zahler C. T. Shaw B. F. Angew. Chem., Int. Ed. 2020 59 10989 10995 10.1002/anie.202001452 32212239
Zahler C. T. Zhou H. Abdolvahabi A. Holden R. L. Rasouli S. Tao P. Shaw B. F. Angew. Chem., Int. Ed. 2018 57 5364 5368 10.1002/anie.201712306 29451960
Marcus R. A. J. Chem. Phys. 1957 26 867 871 10.1063/1.1743423
Winkler J. R. Gray H. B. Chem. Rev. 1992 92 369 379 10.1021/cr00011a001
Shi Y. Rhodes N. R. Abdolvahabi A. Kohn T. Cook N. P. Marti A. A. Shaw B. F. J. Am. Chem. Soc. 2013 135 15897 15908 10.1021/ja407801x 24066782
Baumer K. M. Koone J. C. Shaw B. F. ACS Chem. Neurosci. 2020 11 304 313 10.1021/acschemneuro.9b00464 31895541
Shi Y. Mowery R. A. Shaw B. F. J. Mol. Biol. 2013 425 4388 4404 10.1016/j.jmb.2013.07.018 23871896
Li P. Merz, Jr. K. M. J. Chem. Inf. Model. 2016 56 599 604 10.1021/acs.jcim.5b00674 26913476
Case D. A. Cheatham III T. E. Darden T. Gohlke H. Luo R. Merz Jr K. M. Onufriev A. Simmerling C. Wang B. Woods R. J. J. Comput. Chem. 2005 26 1668 1688 10.1002/jcc.20290 16200636
Case H. M. A. D. A. , Belfon K. , Ben-Shalom I. Y. , Brozell S. R. , Cerutti D. S. , Cheatham, III T. E. , Cruzeiro V. W. D. , Darden T. A. , Duke R. E. , Giambasu G. , Gilson M. K. , Gohlke H. , Goetz A. W. , Harris R. , Izadi S. , Izmailov S. A. , Jin C. , Kasavajhala K. , Kaymak M. C. , King E. , Kovalenko A. , Kurtzman T. , Lee T. S. , LeGrand S. , Li P. , Lin C. , Liu J. , Luchko T. , Luo R. , Machado M. , Man V. , Manathunga M. , Merz K. M. , Miao Y. , Mikhailovskii O. , Monard G. , Nguyen H. , O’Hearn K. A. , Onufriev A. , Pan F. , Pantano S. , Qi R. , Rahnamoun A. , Roe D. R. , Roitberg A. , Sagui C. , Schott-Verdugo S. , Shen J. , Simmerling C. L. , Skrynnikov N. R. , Smith J. , Swails J. , Walker R. C. , Wang J. , Wei H. , Wolf R. M. , Wu X. , Xue Y. , York D. M. , Zhao S. and Kollman P. A. , Amber 2022, University of California, San Francisco, 2022
Seminario J. M. Int. J. Quantum Chem. 1996 60 1271 1277 10.1002/(SICI)1097-461X(1996)60:7<1271::AID-QUA8>3.0.CO;2-W
Becke A. J. Chem. Phys. 1993 98 5648 10.1063/1.464913
Stephens P. J. Devlin F. J. Chabalowski C. F. Frisch M. J. J. Phys. Chem. 1994 98 11623 11627 10.1021/j100096a001
Hay P. J. Wadt W. R. J. Chem. Phys. 1985 82 299 310 10.1063/1.448975
Ditchfield R. Hehre W. J. Pople J. A. J. Chem. Phys. 1971 54 724 728 10.1063/1.1674902
Hehre W. J. Ditchfield R. Pople J. A. J. Chem. Phys. 1972 56 2257 2261 10.1063/1.1677527
Hariharan P. C. Pople J. A. Theor. Chim. Acta 1973 28 213 222 10.1007/BF00533485
Frisch M. J. , Trucks G. W. , Schlegel H. B. , Scuseria G. E. , Robb M. A. , Cheeseman J. R. , Scalmani G. , Barone V. , Petersson G. A. , Nakatsuji H. , Li X. , Caricato M. , Marenich A. V. , Bloino J. , Janesko B. G. , Gomperts R. , Mennucci B. , Hratchian H. P. , Ortiz J. V. , Izmaylov A. F. , Sonnenberg J. L. , Williams-Young D. , Ding F. , Lipparini F. , Egidi F. , Goings J. , Peng B. , Petrone A. , Henderson T. , Ranasinghe D. , Zakrzewski V. G. , Gao J. , Rega N. , Zheng G. , Liang W. , Hada M. , Ehara M. , Toyota K. , Fukuda R. , Hasegawa J. , Ishida M. , Nakajima T. , Honda Y. , Kitao O. , Nakai H. , Vreven T. , Throssell K. , Montgomery, Jr. J. A. , Peralta J. E. , Ogliaro F. , Bearpark M. J. , Heyd J. J. , Brothers E. N. , Kudin K. N. , Staroverov V. N. , Keith T. A. , Kobayashi R. , Normand J. , Raghavachari K. , Rendell A. P. , Burant J. C. , Iyengar S. S. , Tomasi J. , Cossi M. , Millam J. M. , Klene M. , Adamo C. , Cammi R. , Ochterski J. W. , Martin R. L. , Morokuma K. , Farkas O. , Foresman J. B. and Fox D. J. , Gaussian 16, Revision A.03, Gaussian, Inc., Wallingford CT, 2016
Hornak V. Abel R. Okur A. Strockbine B. Roitberg A. Simmerling C. Proteins 2006 65 712 725 10.1002/prot.21123 16981200
Strange R. W. Antonyuk S. V. Hough M. A. Doucette P. A. Valentine J. S. Hasnain S. S. J. Mol. Biol. 2006 356 1152 1162 10.1016/j.jmb.2005.11.081 16406071
Jorgensen W. L. Chandrasekhar J. Madura J. D. Impey R. W. Klein M. L. J. Chem. Phys. 1983 79 926 935 10.1063/1.445869
Joung I. S. Cheatham III T. E. J. Phys. Chem. B 2008 112 9020 9041 10.1021/jp8001614 18593145
Darden T. York D. Pedersen L. J. Chem. Phys. 1993 98 10089 10092 10.1063/1.464397
Ryckaert J.-P. Ciccotti G. Berendsen H. J. J. Comput. Phys. 1977 23 327 341 10.1016/0021-9991(77)90098-5
Miyamoto S. Kollman P. A. J. Comput. Chem. 1992 13 952 962 10.1002/jcc.540130805
Gotz A. W. Williamson M. J. Xu D. Poole D. Le Grand S. Walker R. C. J. Chem. Theory Comput. 2012 8 1542 1555 10.1021/ct200909j 22582031
Valentine J. S. Pantoliano M. W. McDonnell P. J. Burger A. R. Lippard S. J. Proc. Natl. Acad. Sci. U. S. A. 1979 76 4245 4249 10.1073/pnas.76.9.4245 41239
Goto J. J. Zhu H. Sanchez R. J. Nersissian A. Gralla E. B. Valentine J. S. Cabelli D. E. J. Biol. Chem. 2000 275 1007 1014 10.1074/jbc.275.2.1007 10625639
Fee J. A. Briggs R. G. Biochim. Biophys. Acta, Protein Struct. 1975 400 439 450 10.1016/0005-2795(75)90200-7 169909
Andreini C. Banci L. Bertini I. Rosato A. J. Proteome Res. 2008 7 209 216 10.1021/pr070480u 17988086
Pirro F. La Gatta S. Arrigoni F. Famulari A. Maglio O. Del Vecchio P. Chiesa M. De Gioia L. Bertini L. Chino M. Nastri F. Lombardi A. Angew. Chem., Int. Ed. 2023 62 e202211552 10.1002/anie.202211552 36334012
Rich P. R. and Maréchal A. , in Comprehensive Biophysics, ed. E. H. Egelman, Elsevier, Amsterdam, 2012, pp. 72–93
Zhang W. Moore C. E. Zhang S. J. Am. Chem. Soc. 2022 144 1709 1717 10.1021/jacs.1c10948 35044761
Kipouros I. Solomon E. I. FEBS Lett. 2023 597 65 78 10.1002/1873-3468.14503 36178078
MacPherson I. S. Murphy M. E. Cell. Mol. Life Sci. 2007 64 2887 2899 10.1007/s00018-007-7310-9 17876515
Healy E. F. Flores R. Lynch V. M. Toledo S. J. Inorg. Biochem. 2020 210 111161 10.1016/j.jinorgbio.2020.111161 32619897
Tainer J. A. Getzoff E. D. Richardson J. S. Richardson D. C. Nature 1983 306 284 287 10.1038/306284a0 6316150
Shi Y. Abdolvahabi A. Shaw B. F. Protein Sci. 2014 23 1417 1433 10.1002/pro.2526 25052939
Colton I. J. Anderson J. R. Gao J. Chapman R. G. Isaacs L. Whitesides G. M. J. Am. Chem. Soc. 1997 119 12701 12709 10.1021/ja9723491
Ranguelova K. Bonini M. G. Mason R. P. Environ. Health Perspect. 2010 118 970 975 10.1289/ehp.0901533 20348042
Huang X. Coleman W. F. Zare R. N. J. Chromatogr. A 1989 480 95 110 10.1016/S0021-9673(01)84281-3
Ghosal S. Electrophoresis 2004 25 214 228 10.1002/elps.200305745 14743475
Grossman P. D. and Colburn J. C. , Capillary electrophoresis: Theory and practice, Academic Press, 2012
Tainer J. A. Getzoff E. D. Richardson J. S. Richardson D. C. Nature 1983 306 284 287 10.1038/306284a0 6316150
Fridovich I. Annu. Rev. Biochem. 1995 64 97 112 10.1146/annurev.bi.64.070195.000525 7574505
Valentine J. S. Pantoliano M. W. McDonnell P. J. Burger A. R. Lippard S. J. Proc. Natl. Acad. Sci. U. S. A. 1979 76 4245 4249 10.1073/pnas.76.9.4245 41239
Ellerby L. M. Cabelli D. E. Graden J. A. Valentine J. S. J. Am. Chem. Soc. 1996 118 6556 6561 10.1021/ja953845x
Valentine J. S. Doucette P. A. Potter S. Z. Annu. Rev. Biochem. 2005 74 563 593 10.1146/annurev.biochem.72.121801.161647 15952898
Seetharaman S. V. Winkler D. D. Taylor A. B. Cao X. Whitson L. J. Doucette P. A. Valentine J. S. Schirf V. Demeler B. Carroll M. C. Culotta V. C. Hart P. J. Biochemistry 2010 49 5714 5725 10.1021/bi100314n 20515040
