Anna
Barwinska-Sendra
a,
Arnaud
Baslé
a,
Kevin J.
Waldron
*a and
Sun
Un
*b
aInstitute for Cell and Molecular Biosciences, Newcastle University, Newcastle upon Tyne, NE2 4HH, UK. E-mail: kevin.waldron@newcastle.ac.uk
bDepartment of Biochemistry, Biophysics and Structural Biology, Institute for Integrative Biology of the Cell (I2BC), Université Paris-Saclay, CEA, CNRS UMR 9198, CEA-Saclay, Gif-sur-Yvette, F-91198, France. E-mail: sun.un@cea.fr
First published on 18th December 2017
The pathogenicity of Staphylococcus aureus is enhanced by having two superoxide dismutases (SODs): a Mn-specific SOD and another that can use either Mn or Fe. Using 94 GHz electron-nuclear double resonance (ENDOR) and electron double resonance detected (ELDOR)-NMR we show that, despite their different metal-specificities, their structural and electronic similarities extend down to their active-site 1H– and 14N–Mn(II) hyperfine interactions. However these interactions, and hence the positions of these nuclei, are different in the inactive Mn-reconstituted Escherichia coli Fe-specific SOD. Density functional theory modelling attributes this to a different angular position of the E. coli H171 ligand. This likely disrupts the Mn–H171–E170′ triad causing a shift in charge and in metal redox potential, leading to the loss of activity. This is supported by the correlated differences in the Mn(II) zero-field interactions of the three SOD types and suggests that the triad is important for determining metal specific activity.
Metalloproteins can precisely control the chemistry of their metal cofactors, and the SODs are a remarkable example of this phenomenon. Mn-, cam- and Fe-specific SODs (herein FeSOD) can all bind both metals. The active sites of these SODs, composed of a metal ion and its five ligands (three histidines, an aspartate, and a water molecule, (Fig. 1)) in a distorted trigonal bipyramidal arrangement, are virtually identical to within the average 2.0 Å resolution found for the published MnSOD and Mn-loaded camSOD crystallographic structures (see the ESI†). In spite of this, Mn- and FeSODs display highly metal-specific activity. This behavior is thought to be conferred by protein redox tuning of the metal,2,3 although there are few structural clues as to how this is accomplished.
Fig. 1 The structure of the interconnected metal sites of the two subunits of Mn/FeSODs, with carbon atoms in green, oxygen atoms in red, nitrogen atoms in blue, and the metal ion in purple. The residue numbering here, and in the text, is based on the familiar E. coli MnSOD (see the ESI†). |
Unlike that of the Fe(III)/Fe(II) couple (0.77 V), the standard electrode potential (E0) of the aqueous complexes of Mn(III)/Mn(II) (1.51 V) is outside of the range required for catalytic SOD activity (between −0.16 and +0.89 V).4 Hence, the protein environments of MnSOD and FeSOD modify their respective cofactor's E0 in a differential manner, and achieve this feat while having active site structures that are nearly identical, extending well beyond the primary ligand sphere. Since stabilization of the water ligand over an OH− will increase E0, favoring the M(II) ion with respect to M(III), it is likely that D167 and Q146, to which the water respectively donates and accepts a hydrogen-bond (Fig. 1), play a role in the metal redox tuning. The mutation of Q146 to glutamate in E. coli FeSOD reverses the hydrogen bonding direction to the water ligand and has been shown to increase E0 by more than 0.66 V.5 A similar range is observed for simple isostructural [M(4′-R-terpyridine)2]2+ complexes of Mn and Fe whose metal–ligand bond-lengths vary by less than 0.07 Å. Their E0 values are controlled by the electron-donating capacity of the R-group and, as a consequence, they can exhibit metal-specific and cambialistic ‘SOD-activity’ depending on R.6
Although these complexes are isostructural, their Mn forms are readily distinguished by their Mn(II) zero-field interactions.7 Importantly, this is also true for MnSOD, camSOD and FeSOD with Mn(II) in their active-sites (herein Mn(Mn)SOD, Mn(cam)SOD and Mn(Fe)SOD).8 The zero-field interaction is a magnetic property of paramagnetic species that have more than one unpaired electron and arises from magnetic spin–orbit- and spin–spin interactions of these electrons. It is characterized by its D and E values and its orientation with respect to the molecular frame. The sum |D| + E of the three SODs falls into distinct ranges. Mn(Mn)SODs have sums greater than 11.5 GHz, while they are less than 10.9 GHz for the catalytically inactive, Mn(II) reconstituted FeSODs, (Mn(Fe)SOD).8 Mn(cam)SODs and those carrying mutations that alter metal specificity, fall in-between.8 Based on the studies of inorganic model systems like those described above, it has been tempting to ascribe these observations to differences in charge polarization. If this is the case, the underlying cause is not discernible at 2.0 Å resolution of crystal structures. The positions of ligand protons relative to the metal centers, in particular those of the water–ligand, should provide valuable information regarding Mn(II) zero-field interaction of SODs and, more importantly, shed light on their metal specific activity.
To this end, using high-frequency electron paramagnetic resonance (EPR) techniques we have not only measured the Mn(II) zero-field interactions of the two active S. aureus proteins and the inactive Mn-substituted E. coli FeSOD (Mn(Fe)SOD), but also the hyperfine interactions between their Mn(II) centers and nearby proton and nitrogen nuclear spins. Each hyperfine interaction, An, is characterized by the size and orientation of its three principal tensor values (An,mm, m = x, y, z). These values are the sum of the isotropic (An,iso) and anisotropic dipolar (Tn,mm) contributions. An,iso has no explicit dependence on the structure, but is sensitive to the bonding between the metal and the nuclei. By contrast, the dipolar (Tn) contribution depends on the distance (rn) between the nucleus and the Mn(II) electron spin, and the angle (θn) which their inter-spin vector makes with B0. As illustrated in Fig. 2, rn and θn are explicitly dependent on the molecular structure. This makes hyperfine interactions an effective structural probe of the metal ligand sphere that can provide potentially finer details than crystallography.9 When the nuclei and Mn(II) are approximated as point-dipoles, as depicted in Fig. 2, the geometric dependence of An is easier to appreciate and is given by
(1) |
The point-dipole approximation is not completely appropriate because of the proximity of the 1H and 14N nuclei to the spatially large unpaired electron spin-density of Mn(II). Even so, eqn (1) does allow for simple estimations and trends. For protons surrounding the Mn(II) centers in SODs, their Mn(II) hyperfine couplings are dominated by the dipolar contribution and, hence, determined by their distance to the metal center. Using the SOD crystal structures and assuming standard hydrogen positions on ligands, the largest proton couplings will arise from the two protons on the water ligand with Tn > 3.5 MHz (with Aiso < 1 MHz), followed by the six imidazole protons (H2 and H5) that flank the three ligating histidine nitrogens with a Tn of 1.8–2.4 MHz (with Aiso < 0.1 MHz).9–11
The hyperfine interactions of the ligand protons and nitrogens are too small to be measured from the EPR spectra directly, but they can be measured using electron nuclear double resonances (ENDOR) and electron double resonance detected NMR (ELDOR-NMR). These techniques, described in detail elsewhere12,13 and in the ESI†, monitor the change in the EPR signal in response to the excitation of the spins of nuclei surrounding the metal ion, using either the allowed nuclear spin transitions (ENDOR, cyan transition in Fig. 2) or the forbidden simultaneous transition of both the electron and nuclear spins (ELDOR-NMR, magenta and purple transitions). Each nuclear spin will have two resonances at approximately:
(2) |
For S > 1/2 nuclei such as 14N, the nuclear quadrupolar interaction will further split the resonances (see the ESI†). In addition to single quantum (SQ, q = 1) transitions, multiple quantum transitions (q > 1) can be readily detected by ELDOR-NMR. These not only involve single nuclei, as is the case for nuclei like 14N that have a spin greater than 1/2, but also combinations of two or more nuclei, such as two spin 1/2 protons (q = 2), and even a DQ transition of a 14N and a SQ transition of a 1H (q = 3). Such multi-nuclei transitions can only arise from nuclei that are coupled to a common Mn(II), providing structural details not evident in the SQ and ENDOR spectra. One potential problem is that if the νn,NMR are small then the SQ and multiple-quantum resonances of different nuclei may overlap. Since νn,NMR are proportional to the applied magnetic field, we have exploited sufficiently high magnetic field of 3 T to partially alleviate this problem. In the following, we only consider the positive branch of the ELDOR-NMR spectra, and define νn′ = νobs + qνn,NMR, which can be subtracted from eqn (2) to obtain An,eff in a straightforward manner.
To relate the measured hyperfine couplings to the structure of the active-sites, we have used density functional theory (DFT) calculations on model supramolecular Mn(II) complexes based on the S. aureus crystal structures, as well as many other available SOD crystal structures. This combination of measurements and DFT calculations has allowed us to obtain a detailed picture of the electronic structure of the ligand sphere that suggests that the Mn–H171–E170′ triad may be as important as Q146 in determining the metal specificity of enzymatic activity of Mn- and Fe-dependent SODs.
Fig. 3 Graphic representation of Mn coordination in the active sites of the homodimers of S. aureus (A) MnSOD (ribbon in blue) and (B) camSOD (ribbon in teal). Mn is represented as purple spheres. The orange mesh represents the anomalous difference map rendered at 6.0σ and 15.0σ on the Mn ions at the active site of MnSOD and camSOD, respectively. Metal-coordinating ligand residues are represented as sticks with carbon atoms colored in grey, oxygens in red, nitrogens in blue and a water molecule represented as a red sphere (PDB accession codes: 5N56 for MnSOD and 5N57 for camSOD). |
Model | Owater–Mn–dihedral (°) | Mn-Interatomic distance (Å) | |||
---|---|---|---|---|---|
H171 | H81 | O–CO (D167) | OHOH | Oδ2,D167 | |
a Based on average of 57 Mn(Mn)SOD and 4 Mn(cam)SOD structures (see the ESI). b See text and Fig. 8 for details. | |||||
S. aureus | 30 ± 5 | −14 ± 3 | 19 ± 8 | 2.22 ± 0.01 | 2.09 ± 0.03 |
Mn(Mn)SOD | |||||
S. aureus | 25 ± 1 | −4 ± 2 | 20 ± 4 | 2.17 ± 0.04 | 1.94 ± 0.03 |
Mn(cam)SOD | |||||
Crystallographica | 29 ± 8 | −15 ± 10 | 16 ± 10 | 2.18 ± 0.12 | 2.01 ± 0.07 |
(Mn(Mn)SOD) | |||||
2.17 ± 0.12 | |||||
(Mn(cam)SOD) | |||||
GOb | −61 | −8 | 26.3 | 2.32 | 2.05 |
CDb | 33 | −8 | 15 | 2.30 | 2.04 |
CD(NH)b | 29 | −8 | 15 | 2.32 | 2.06 |
CD(OH)b | 29 | −8 | 15 | 2.33 | 2.07 |
Mn-Interatomic distance (Å) | |||||
---|---|---|---|---|---|
Nε2,H171 | Nε2,H81 | Nε2,H26 | H (water ligand) | H (water ligand) | |
Free | Hydrogen bound to D167 | ||||
S. aureus | 2.31 ± 0.11 | 2.15 ± 0.02 | 2.37 ± 0.01 | ||
Mn(Mn)SOD | |||||
S. aureus | 2.32 ± 0.02 | 2.33 ± 0.01 | 2.36 ± 0.01 | ||
Mn(cam)SOD | |||||
Crystallographica | 2.19 ± 0.07 | 2.17 ± 0.07 | 2.20 ± 0.07 | ||
GOb | 2.22 | 2.18 | 2.24 | 2.89 | 2.61 |
14N: (1.57, 1.65, 3.33) | 14N: (1.48, 1.55, 3.29) | 14N: (2.59, 2.61, 4.40) | (−1.99, −2.44, 7.78) | (−3.71, −4.53, 8.60) | |
1H2: (−2.10, −2.25, 4.39) | 1H2: (−2.07, −2.25, 4.05) | 1H2: (−2.11, −2.29, 4.44) | (−2.58, −2.58, 7.32) | (−3.95, −3.95, 8.35) | |
1H5: (−1.55, −1.65, 3.73) | 1H5: (−1.78, −1.87, 3.93) | 1H5: (−1.48, −1.57,3.57) | |||
(−1.63, −1.63, 3.71) | |||||
CDb | 2.20 | 2.19 | 2.25 | 2.97 | 2.59 |
CD(NH)b | 2.12 | 2.21 | 2.27 | 2.87 | 2.63 |
14N: (0.56, 0.61, 2.31) | 14N: (1.92, 1.99, 3.73) | 14N: (2.92, 2.94, 4.74) | (−2.42, −2.92, 7.40) | (−3.70, −4.51, 8.40) | |
1H2: (−1.90, −2.05, 3.87) | 1H2: (−1.98, −2.16, 4.02) | 1H2: (−1.98, −2.14, 4.15) | (−2.58, −2.58, 7.32) | (−3.95, −3.95, 8.35) | |
1H5: (−1.86, −1.99,4.39) | 1H5: (−1.73, −1.84, 3.89) | 1H5: (−1.57, −1.65, 3.66) | |||
(−1.63, −1.63, 3.71) | |||||
CD(OH)b | 2.10 | 2.22 | 2.27 | 2.82 | 2.64 |
14N: (0.33, 0.35, 2.03) | 14N: (1.95, 2.03, 3.75) | 14N: (2.96, 2.98, 4.77) | (−2.75, −3.304, 7.5) | (−3.66, −4.48, 8.30) | |
1H2: (−1.94, −2.10, 3.98) | 1H2: (−1.99, −2.16, 4.03) | 1H2: (−1.95, −2.11, 4.09) | |||
1H5: (−1.86, −2.00, 4.45) | 1H5: (−1.72, −1.82, 3.85) | 1H5: (−1.59, −1.66, 3.68) |
(3) |
Fig. 5 shows the 1H ELDOR-NMR spectra obtained at these two extremes and Fig. 6, the Davies121H ENDOR spectra taken at the Dyy,−3/2 magnetic-field position. The corresponding ν1H,NMR were 79, 185 and 154 MHz. The arrows in Fig. 4 indicate the specific magnetic-fields for each protein where the ELDOR-NMR and ENDOR spectra were obtained. ENDOR spectra at the three fields would have been ideal, since ENDOR resonances are typically sharper than the ELDOR-NMR ones.11,15 However, the intensity of the 1H ENDOR spectra taken at the Dzz,−5/2 and Dxx,−5/2 field positions was too low to be useful (Fig. 6). This was likely due to our inability to deliver sufficient radio-frequency power to excite the NMR transitions at ν1H,NMR of 79 and 185 MHz, which were far from 144 MHz for which our system is optimized. By comparison this was not the case for the Dyy,−5/2 field position, for which ν1H,NMR was 154 MHz.
Fig. 5 The 94 GHz 5 K SQ (bottom) and DQ (top) 1H ELDOR-NMR spectra of Mn(Mn)SOD (black), Mn(cam)SOD (green) and Mn(Fe)SOD (red) taken at the Dzz,−5/2 field positions. The SQ spectra taken at the Dxx,−5/2 field position are shown in the inset. The exact field positions are indicated by the arrows in Fig. 2. The Gaussian line shape simulations of the SQ spectra are also shown, with their colors corresponding to the measured spectra and their sums as dashed lines. The DFT hyperfine histograms of the GO (black) and CD(NH) (red) models are superimposed and their correspondence to the measured spectra is indicated by the dotted-lines. For the Dzz,−5/2 spectra, ν1H,NMR was 79 MHz and for the Dxx,−5/2, it was 185 MHz. See text and the ESI† for details. |
Fig. 6 The 94 GHz 5 K Davies 1H ENDOR spectra of Mn(Mn)SOD, Mn(Fe)SOD and Mn(cam)SOD taken at the Dyy,−3/2 magnet field positions indicated by the arrows in Fig. 2. The blue trace was obtained under the same conditions but at the D−xx,−5/2 magnetic-field position indicated by the black arrow in Fig. 2. The initial electron-spin inversion pulse was 200 ns, followed by a 16 μs radio-frequency pulse and standard spin-echo detection (12 and 24 ns pulses). The lower panel shows the calculated ENDOR spectra based on DFT hyperfine tensors obtained for the GO (red) and CD(NH) (black) models. The cyan traces show the calculated spectra obtained by symmetrizing the DFT hyperfine tensors and manually adjusting the Aiso values. The calculated spectra have been convolved using 150 kHz Gaussian and each arbitrarily scaled (see text and the ESI† for details). |
We were struck by the overall similarity of the 1H ELDOR-NMR and ENDOR spectra of the three proteins. The Dxx,−5/2 spectra (inset of Fig. 5) of Mn(Mn)SOD, Mn(cam)SOD and Mn(Fe)SOD were essentially identical and relatively featureless, extending asymmetrically about ν1H,NMR from −12.0 to 4.7 MHz, (or an An,eff range of −4.8 to 1.8 MHz) for the magnetic-field oriented along the Dxx direction (eqn (3)). In addition to the large matrix signal at ν − ν1H = 0, the Dzz,−5/2 spectra had two resolved features at 8 and 16 MHz (Fig. 5; see the ESI† for full spectra). Simulation based on eqn (2) and assuming Gaussian line-shapes showed that these features were unrelated and arose from two different protons with An,eff of 5.3 (2.2) and 6.4 (2.2) MHz (with the Gaussian widths in parentheses). The intensity pattern of the Dyy,−5/21H ENDOR spectra were also strikingly similar (Fig. 6). All three proteins had nearly identical resonances at −6.3 and −9.2 MHz. The tempting conclusion based on this close similarity of the 1H ELDOR-NMR and ENDOR spectra was that the positions of the five ligands at the level of the hydrogen atoms were the same for the three proteins.
This was not the case for Mn(Fe)SOD. As can be seen in Fig. 5, the positive half of all three DQ 1H ELDOR-NMR spectra is nearly identical extending to 30 MHz, twice that of the SQ spectra, demonstrating that in all three proteins the resolved high-frequency portion of the SQ spectra arose from two protons from a common Mn(II) center. The data lend themselves to two possible interpretations. The DQ high-frequency edge was defined by either: (1) two protons having the same Aeff of 6.4 MHz; or (2) one proton with Aeff = 5.3 MHz and the other 6.4 MHz. The resolution of the DQ spectra did not allow us to distinguish between these two possibilities. Mn(Fe)SOD also had at least one proton that was distinct. This could be seen in the positions of the low-frequency edges of DQ spectra. The low-frequency edge of Mn(Fe)SOD was 6 MHz lower, indicating that it had at least one pair of protons with a larger combined negative Aeff than the two active proteins. Both the 1H ENDOR and SQ ELDOR-NMR spectra of Mn(Fe)SOD were also distinct in the detail. The simulations of SQ spectra also required a third An,eff = 2.6 MHz component to explain the reproducibly lower amplitude of the Mn(Fe)SOD SQ 1H ELDOR-NMR spectrum in the 4–8 MHz region. The uniqueness of Mn(Fe)SOD was also evident in the structure of the ENDOR resonances at −3 MHz (orange region, Fig. 6). In the yellow regions, which arise from the smallest An,eff, all three proteins exhibited subtle but reproducible differences.
The close similarity between Mn(Mn)SOD and Mn(cam)SOD and the uniqueness of Mn(Fe)SOD extended to their 14N hyperfine interactions (Fig. 7). All three SODs had readily detectable SQ and DQ 14N ELDOR-NMR resonances. The large intense center peak (ν − ν14N = 0) obscured the low-frequency portions of SQ resonances, while the low-frequency regions of DQ resonances overlapped with the much larger SQ ones. Nonetheless, it was evident that the three proteins shared a common DQ feature, a partially resolved unequal doublet. The resonance at 20 MHz dictated another at 12 MHz since resonances occur pairwise (eqn (2)). Since, as in the case of radicals,16 the 1-spin DQ transition frequencies are purely hyperfine determined and are unaffected by 14N nuclear quadrupolar coupling, the doublet corresponded to an An,eff of 4 MHz, the largest 14N hyperfine coupling when the magnetic-field is along the Dzz direction. The unequal amplitudes of the DQ doublet likely arose from the overlap of the 12 MHz resonance with those arising from smaller Aeff. The three proteins also had a partially resolved SQ doublet with a separation ∼0.5 MHz smaller than 4 MHz. This may have been due to 14N nuclear quadrupolar contributions. What made Mn(Fe)SOD stand out was that its SQ and DQ intensity pattern was reproducibly different from those of the active proteins.
Fig. 7 The 94 GHz 5 K SQ (bottom) and DQ (top) 14N ELDOR-NMR spectra of Mn(Mn)SOD (black), Mn(cam)SOD (green) and Mn(Fe)SOD (red) taken at the Dzz,−5/2 field positions indicated by the arrows in Fig. 2. The DFT hyperfine histograms of the GO (black) and CD(NH) (red) models (ν14N,NMR = 5.7 MHz). |
Fig. 8 The DFT GO, CD(NH) and CD(OH) model structures. The numbers indicate the CM5 charges of each atom of the His171 ring with their hydrogen charges, if any, summed into them. The orientation of the ligands in Mn(Mn)SOD and Mn(cam)SOD is shown in Fig. 3 and 6. |
However, the measured ENDOR and ELDOR-NMR spectra were obtained relative to the Mn(II) zero-field interactions of the three proteins and not their molecular frame. Although we knew that the 1H ELDOR-NMR spectra in Fig. 5 corresponded to the case where the applied magnetic-field was parallel with the Dzz (θzf = 0), what was not known was how the zero-field interactions were orientated with respect to the molecular frame of the SODs. In simple ionic Mn(II) complexes, Dzz is often assumed to lie along the axis of greatest symmetry, so we initially placed it along the pseudosymmetry axis, defined by the Owater–Mn–NHis26 bonds, and carried out a brute-force search for the Euler angles, the angles of the three rotations that related the zero-field axes to the molecular frame, which yielded the best agreement between the DFT proton hyperfine couplings and ELDOR-NMR and ENDOR spectra. It is assumed that the DFT values were sufficiently accurate to reproduce the measured one. As will be seen, this appeared to the case. The search criteria were: (1) the A1H,eff values calculated from the DFT reproduce the three couplings found from the simulations of the ELDOR-NMR spectra; and (2) the DFT-based calculated Dyy,−5/2 ENDOR spectra match the −5 to −10 MHz region of the measured spectra. These requirements were fulfilled when Dzz was placed along the Nε,His26–Mn(II) bond and the Dxx along the Oδ2,Asp167–Mn bond as depicted in Fig. 2. The calculated ENDOR spectra based on DFT hyperfine tensors and orientation are shown in Fig. 6. The theoretical results, in the form of histograms that show frequency and transition probability for the magnetic-field oriented along the Dzz, are superimposed on the ELDOR-NMR spectra in Fig. 5 and 7 (see the ESI† for methodological details). In both cases, the agreement between the measured and calculated spectra was very good.
The line-shapes of the calculated ENDOR spectra of the water–ligand protons based on the DFT hyperfine tensors were more complex than the measured one. This appeared to be linked to the symmetry of the DFT tensors. Better agreement with the experimental data was obtained by averaging Txx and Tyy—that is by making the dipolar tensor axial (and similar to eqn (1)). This was consistent with single crystal 1H ENDOR measurements on the Mn(II) center in ConA by Carmieli and co-workers.9 They showed that the hyperfine tensors of water and histidine ligand protons were symmetric.9 For the SODs, small adjustments to the DFT Aiso were also needed. The biggest was for the Mn(Fe)SOD Hw proton (Fig. 6 and Table 1) of 0.4 MHz, while others needed less than 0.1 MHz. This resulted in good agreement between the calculated (cyan traces in Fig. 6) and measured ENDOR spectra. The DFT calculations apparently over-estimated the asymmetry of the hyperfine interactions.
As with the proton hyperfine interactions, the DFT derived 14N hyperfine tensors for GO and CD(NH) models accounted for all the salient features of the 14N ELDOR-NMR spectra of Mn(Fe)SOD, Mn(Mn)SOD and Mn(cam)SOD. To achieve this agreement required the inclusion of the nuclear quadrupolar coupling obtained from the DFT calculations. The largest hyperfine coupling of 4 MHz was consistent with that of Nε2,H26 (Table 1). Interestingly, 14N DQ histograms appear to account for the small difference in amplitude at 17 MHz (marked by an asterisk in Fig. 7) between Mn(Fe)SOD and the active proteins. Beyond this, the DFT calculations did not reveal the source of the lower amplitude of the Mn(Fe)SOD 14N ELDOR-NMR spectra in the 2–7 MHz region. A more detailed study using HYSCORE spectroscopy is planned to obtain more detailed information regarding Mn–N bonding.
The DFT calculations confirmed that the Aeff = 6.4 and 5.3 1H ELDOR-NMR components arose from the two water–ligand protons, the larger value from the solvent that forms hydrogen bonds with Oδ1,Asp167. The ENDOR resonances arising from these protons are respectively labelled as Hw′ and Hw in Fig. 6. The dipolar components of their hyperfine tensors were 4.1 MHz (22°, 2.62 Å) and 3.3 MHz (17°, 2.91 Å), the values in parentheses corresponding to the H–Mn(II)–Owater angles and H–Mn(II) distances from the DFT calculations. These were close to those measured for concanavalin-A, 3.5 MHz (19°, 2.77 Å) and 4.0 MHz (11°, 3.08 Å), (referred to as Hw1 and Hw2, respectively, in Carmieli et al.9 and the angles and distances coming from a neutron diffraction study17).
The calculated ELDOR-NMR transition probabilities were sensitive to geometry, which likely explained the small differences in amplitudes of SQ 1H spectra between the active proteins and Mn(Fe)SOD (Fig. 5).
As shown in Fig. 8, the ring plane of H171 in the GO DFT model is tilted −60° relative to the Owater–Mn–NHis26 axis compared to +30° in the crystal structures (and, by definition, in the CD(NH) DFT model). This larger dihedral angle causes the angular terms of the dipolar hyperfine couplings (eqn (1)) of H2,His171 and H5,His171 to become more negative leading to Aeff = −1.8 and −0.5 MHz, respectively, in the GO model, compared to −0.9 and +1.6 MHz in the CD(NH) model. It is for this reason that the DQ DFT histogram of the GO model extends to higher negative values than that of the CD models (Fig. 5). This semi-quantitatively matched what was seen in the experimental spectra. Such differences in the ring dihedral angle will manifest themselves most prominently in the orange region of the 1H ENDOR spectra in Fig. 6. The calculations showed that the extra feature in the Mn(Fe)SOD marked by the red arrow in Fig. 6 likely arose from the different positions of H2,His171.
In the same region of the 1H ENDOR spectra (Fig. 6), a green arrow marks an instance where none of the three proteins were the same. Based on the DFT calculations, it was tempting to assign this to the angular position of the ring of H26. In the S. aureus Mn(Mn)SOD, the ring plane of H26 relative to the Mn–Oδ1,Asp167 bond was 180°, while in Mn(cam)SOD it was 167°. A prominent resonance in the calculated spectra that coincides with this position arises from H5,His26. The DFT calculations show that a rotation of H26 about its Mn–Nε bond by 5 to 10° would be sufficient to explain the observed differences. The change in the 1H ELDOR-NMR and ENDOR would be on the order of 200 kHz. This would be difficult to discern in the former due to its broad resonances but readily apparent in the latter.
The calculations also showed that any differences in the Cβ proton hyperfine interactions would contribute to the yellow regions of the ENDOR spectra (Fig. 6). Since these regions were distinct for the three proteins, the positions of these protons, and by extension of the backbones, were unlikely to be the same for Mn(Mn)-, Mn(cam)SODs and Mn(Fe)SOD. This was consistent with the crystal structures of the two S. aureus proteins that showed that the positions of Cα of the ligands relative to the Mn(II) were indeed different.
The high level of agreement between the hyperfine measurements, DFT calculations, and the known crystallographic structures allowed us to derive details about the active sites of the three proteins.
The dihedral-constrained model (CD(NH) and CD(OH), Fig. 8) that we used had a complete H171, with its dihedral angles fixed to crystal structure values and, importantly, the carboxylate sidechain (modelled as an acetate) of E170′, a residue from the other homodimeric subunit. This carboxylate strongly forms hydrogen bonds with both the ring Nδ of H171 and its backbone amide nitrogen (Fig. 1). Energy minimization of these models yielded two energetically close structures: one with an elongated 1.1 Å Nδ,His165–H bond (CD(NH)) and the other, favored by only 0.5 kcal mol−1, with the proton transferred from the Nδ to the acetate. Both had a Nδ,His165–Oacetate distance of 2.63 Å indicative of a strong hydrogen bond (CD(OH)). Different attempts to globally optimize models that included E170′ generally led to E170′ disassociating itself from the complex. The hyperfine histograms of CD(NH) and CD(OH) were essentially the same.
The importance of E170 has already been established,18 but perhaps underestimated. The E. coli MnSOD protein in which E170 is mutated to an alanine contains only iron and is completely inactive even when reconstituted with manganese.18 Metal–histidine–carboxylate triads play a well-recognized important role in other metalloenzymes. In the cytochrome c peroxidase Fe(II)–His–Asp triad, a proton is shared by the histidine and aspartate, an arrangement that is readily perturbed by mutations and ionization of nearby residues.19 Our data and DFT calculations suggest that in Mn(Mn)- and Mn(cam)SOD, E170′ and H171 share a proton and the negative charge is spread over both, making the H171 imidazole-ring negative. The charge of the ring can range from −0.01 to as much as −0.42 (Fig. 8), based on the CM5 charge model which has been shown to provide useful charge distribution information for systems containing metals.20 This Mn–His–Glu triad is likely broken in Mn(Fe)SOD by the rotated position of the H171 ring, as evidenced by the shifted DQ 1H spectra. This leads to the negative charge being localized on the E170′ carboxylate sidechain, and gives the H171 ring a charge of +0.15. The altered position of the Mn(Fe)SOD H171 may be related to the subtle differences in the homodimer structure. Charge localization effects of the triad should closely parallel how the electron-donating capacity of the substituent in the [Mn(4′-R-terpy)2]2+ model affect zero-field interaction and Em. Since in Mn(Mn)SOD and Mn(cam)SOD, E170′ is predicted to have greater electron-donation into H171, these proteins should have a greater zero-field interaction and lower Em than Mn(Fe)SOD. This is indeed the case for the E. coli Mn(Mn)SOD and Mn(Fe)SOD proteins.2,3,8 An appealing physical interpretation of the correlation of |D| + E with activity also becomes apparent.8 Recalling that the zero-field interaction spans −D ± 3E (eqn (1)) in the equatorial plane of the active-sites, and that D < 0, |D| + E provides a measure of the zero-field interaction that only the direction of the electrostatic influence of the Nε,His171/E170′ pair is likely to dominate. The different zero-field interactions for Mn(Mn)SOD and Mn(cam)SOD suggest that there may be even more subtle differences in charge localization in these SODs to which even hyperfine spectroscopy is not sensitive, but that can be readily detected by high-field EPR.
An analysis of the structural data suggests a basis for redox control, and why catalysis by the Mn(Fe)SOD structure may be disrupted. Statistical analysis of the Mn and FeSOD structures (see the ESI†) shows that the metal–Oδ2,E170′ distances are: 6.97 ± 0.08 Å for MnSODs; 6.93 ± 0.07 Å for Mn and 6.92 ± 0.05 Å for Fe-containing cambialistic SODs; and 6.86 ± 0.07 Å for FeSOD. Although bordering on statistical significance, there is a trend towards longer metal–E170′ distances with Mn. Although the structure of Mn(Fe)SOD is yet to be published, if the tenth of an angstrom difference persists when Fe is replaced by Mn in FeSOD, this may be sufficient to disrupt the E170′–H171 hydrogen bond. Similarly, mutations that seemingly have little to do with redox tuning and that are far from the active site may cause volume changes that modify the E170′–H171 separation, resulting in a shift of catalysis towards ‘Mn’ or ‘Fe’-like behavior. This view provides an alternative or additional cause for the conversion of the P. gingivalis cambialistic SOD with high Mn activity to an FeSOD-like one when the G165T (E. coli numbering) mutation is introduced over 11 Å away from the active site. This change in specificity has been interpreted as a result of the displacement of W169 by the threonine, which in turn affected the Q146–water ligand interaction.21 Comparisons of the Q146–water–metal interactions are complicated by the different metals present in the 1.6 Å resolution structures, Fe in the wild type and Mn in the mutant (PDB 1UES and 1UER, respectively).21 Nonetheless in the crystal structures, although the distances from the metal to H171Nε are the same (2.23 ± 0.06 Å for Mn in the G165T mutant and 2.23 ± 0.03 Å for Fe in the wild-type), the metal to Oδ2,E177′ distances are different; 7.14 ± 0.09 Å for the cambialistic wild type and 7.02 ± 0.03 Å for the FeSOD-like mutant, following the same trend and exhibiting similar differences as those discussed above: the G165T mutation also shifts the Mn(II) zero-field interactions from cambialistic to Mn(Fe)SOD values.22 Hence, previous structural studies of the P. gingivalis camSOD enzyme appear to support our interpretation of our own data, which implicates the importance of the metal–histidine–E170′ triad in determining the metal-specificity of SOD catalysis.
Application of EPR techniques in combination with DFT calculations has provided us with a physical basis for understanding how the metal specific activities of SODs is determined at the molecular level. The metal specificity of the SODs are not only of fundamental interest, but play a critical role in determining how cells survive oxidative stress, including resisting the mammalian immune system's oxidative burst during infection by pathogens such as S. aureus.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7cp06829h |
This journal is © the Owner Societies 2018 |