X-ray Crystallography and Electron Paramagnetic Resonance Spectroscopy Reveal Active Site Rearrangement of Cold-Adapted Inorganic Pyrophosphatase

Inorganic pyrophosphatase (PPase) catalyses the hydrolysis reaction of inorganic pyrophosphate to phosphates. Our previous studies showed that manganese (Mn) activated PPase from the psychrophilic bacterium Shewanella sp. AS-11 (Mn-Sh-PPase) has a characteristic temperature dependence of the activity with an optimum at 5 °C. Here we report the X-ray crystallography and electron paramagnetic resonance (EPR) spectroscopy structural analyses of Sh-PPase in the absence and presence of substrate analogues. We successfully determined the crystal structure of Mn-Sh-PPase without substrate and Mg-activated Sh-PPase (Mg-Sh-PPase) complexed with substrate analogue (imidodiphosphate; PNP). Crystallographic studies revealed a bridged water placed at a distance from the di-Mn centre in Mn-Sh-PPase without substrate. The water came closer to the metal centre when PNP bound. EPR analysis of Mn-Sh-PPase without substrate revealed considerably weak exchange coupling, whose magnitude was increased by binding of substrate analogues. The data indicate that the bridged molecule has weak bonds with the di-Mn centre, which suggests a ‘loose’ structure, whereas it comes closer to di-Mn centre by substrate binding, which suggests a ‘well-tuned’ structure for catalysis. Thus, we propose that Sh-PPase can rearrange the active site and that the ‘loose’ structure plays an important role in the cold adaptation mechanism.

In this paper, we initially examined the signal from first excited state (S = 1). The exchange coupling constant (J) was then determined from the temperature dependence of EPR signals from S = 1. Finally, the signal from S = 2 was derived from the Boltzmann curve of S = 2 and the Mn-Mn distance of Sh-PPase in solution was determined. X-ray crystallography and EPR spectroscopy were used for the analysis of the overall and active site structural change of Sh-PPase induced by the substrate. The findings reveal the cold adaptation mechanism in Sh-PPase.

Results and Discussion
X-ray crystal structure. The crystal structure of Mn-Sh-PPase in the absence of substrate was determined at 2.2 Å resolution (Fig. 2a). The M1 and M2 sites showed similar coordination spheres to other family II PPases (Fig. 2b,c). The M1 site in Mn-Sh-PPase was coordinated with one nitrogen of histidine, two oxygens of aspartates (D12 and D72), and a bridged water. Additionally, the M1 site was bound with an extra water molecule 8 . Coordination geometry of the M2 site in Mn-Sh-PPase was also similar with that of Bs-PPase and Sg-PPase with 5-coordinated bipyramidal geometry. However, the distances of a bridged water in the Sh-PPase from M1 and M2 (3.0 and 2.4 Å, respectively) were considerably longer compared with crystal structures of other di-Mn enzymes ( Table 1). In addition to the Mn-O distance, the Mn-Mn distance of Sh-PPase was also longer than arginase and catalase.
A previous crystallographic study of Bs-PPase showed that the addition of the fluoride ion in the PPase substitutes for the bridged water and dramatically reduces the reaction rate 9 . Therefore, we attempted to crystallise the PNP-bound form of Sh-PPase in the presence of sodium fluoride. However, the crystals obtained did not show any electron density for PNP, likely because Mn-Sh-PPase only weakly degrades PNP, even in the presence of sodium fluoride. In contrast, Sh-PPase activated by Mg 2+ ions (Mg-Sh-PPase) has much lower catalytic activity with PPi or PNP than Mn-Sh-PPase 13 , and we obtained well-diffracted crystals of Mg-Sh-PPase in the presence of PNP and sodium fluoride. X-ray diffraction data at 1.3 Å resolution were obtained. The electron density unambiguously showed that the PNP and four Mg 2+ atoms were bound to the di-metal centre of Sh-PPase. The overall structures of Sh-PPase with and without substrate were very similar with other family II PPases, and that the binding of substrate analogue induced the conformational change from the 'open' to 'closed' state ( Fig. 2a,d), as observed in Bs-PPase 7 . In the crystal structure of Mg-Sh-PPase with PNP, the bridged water between Mg 2+ ions at M1 and M2 is assumed to be replaced by a fluoride ion (Fig. 2e). An anomalous difference Fourier map indicated that M1 and M2 metal sites in the crystal structure of Mg-Sh-PPase contain a small fraction of metals other than Mg 2+ . As shown below, the EPR spectrum for apo Sh-PPase at 15 K showed an unexpected signal at g~4.3, which is characteristic of mono Fe 3+ ion 37 (see EPR results). This signal was considerably reduced by the addition of Mn to Sh-PPase, indicating that contaminated Fe 3+ was replaced by Mn 2+ . These observations suggest that Fe 3+ ion is contaminated in the active site of apo and Mg-Sh-PPase, even though all buffers were treated using ethylenediamine-N,N,N′,N′-tetraacetic acid (EDTA). The peak height of the anomalous difference map for M1 and M2 showed similar density level with those of sulphur atoms of Cys or Met residues. Based on the anomalous scattering factor of these elements (f ″ = 0.24 for S; f ″ = 1.6 for Fe) at an X-ray wavelength of 1.0 Å, M1 and M2 were occupied by Fe with an estimated content of 15% and 7%, respectively. Similarly, it was reported that M1 and M2 sites of the crystal structure of Bs-PPase with PNP and fluoride ion were also occupied by Fe/Mn in a 6/1 ratio 9 .
As mentioned above, a unique structural feature of Mn-Sh-PPase found in the present X-ray analysis is the longer distance between bridged water and two Mn atoms compared with that in other family II PPases (Table 1). In the crystal structure of Mg-Sh-PPase complexed with PNP, the site of bridged water was replaced by an oxygen atom of PNP. The active site of Mg-Sh-PPase complexed with PNP shows common coordination sphere, like other family II PPases 4,9 . The M1 site was converted to the 6-coordinated slight rhombic octahedral geometry by binding of PNP, and the M2 site was changed to a nearly axial octahedral coordination sphere (Fig. 2c,e). Based on the crystal structure of Bs-PPase complexed with PNP 9 , Fabrichniy and co-workers proposed that such a dynamic coordination change by binding of substrate must be why family II PPases prefer transition metals to Mg 2+ ions 8 . In addition, a bridged water in the crystal structure of Mg-Sh-PPase with PNP, as well as Bs-PPase, shows the formation of unique tri-metal coordination sphere (Fig. 2c,e), where it is bound with 3 metal ions (M1, M2 and M4). This unique tri-metal coordination should occur with all 3 lone pairs on the oxygen atom of a nucleophilic hydroxide anion (Fig. 2c). The tri-metal coordination needs to play an important role to hold a nucleophile at the suitable position and activate a bridged water for nucleophilic attack. We propose that the tri-metal coordination would be broken, most probably at a bond with M4, because of its longer distance, so that the oxygen atom of nucleophilic hydrogen oxide becomes able to attack the phosphorus of inorganic pyrophosphate. cW X-band epR measurements. Purified Sh-PPase has no transition metal in the active site (apo Sh-PPase). The apo Sh-PPase was incubated at 4 °C for 2 hours with the buffer containing 100 mM Tris-HCl (pH 7.0) and 15 mM MnCl 2 . Subsequently the excess Mn 2+ ions were removed by dialysis and buffer exchange using an ultrafiltration device. The continuous-wave (CW) X-band EPR spectrum showed a strong signal at g~2 with 6 hyperfine splitting (~9 mT) of a characteristic mono-nuclear Mn 2+ centre 38,39 (see Supplementary Fig. S1). The signals split by ~4.5 mT at ~290 mT (see Supplementary Fig. S1 inset) was also observed, which was characteristic of di-Mn 2+ centre (two equivalent I = 5/2) 18,19,[26][27][28] . The signal of mono-nuclear Mn 2+ was disappeared by additional buffer exchange of the sample by metal-free buffer using desalting spin column. Figure 3a shows CW X-band EPR spectra of apo-and Mn-Sh-PPase at 15 K. Addition of Mn ion to apo Sh-PPase diminished mono Fe 3+ signal at g~4.3 (Fig. 3a, asterisk). Alternatively, several hyperfine splitting with ~4.5 mT appeared in a wide range of the spectrum (Fig. 3b). This result confirmed that the active site of Sh-PPase has dominantly di-Mn centre and is consistent with our earlier ICP-mass results in which two metal binding sites of M1 and M2 are occupied by transition metals 12 . (2020) 10:4368 | https://doi.org/10.1038/s41598-020-61217-6 www.nature.com/scientificreports www.nature.com/scientificreports/ Determination of antiferromagnetic coupling. Figure 4a shows the temperature dependence of the EPR spectra for Mn-Sh-PPase. Since EPR spectra at g~2 were overlaid with signals contributed from several exited states of the 'ladder' di-Mn spin system, it was difficult to deconvolute the individual signals. EPR spectra at 2 K still showed a strong signal at ~350 mT with 6 hyperfine splitting of ~9 mT. The intensities of EPR signals www.nature.com/scientificreports www.nature.com/scientificreports/ for di-Mn centre at 2 K should be much smaller because the Boltzmann population was mostly in the ground state (S = 0), which is EPR-silent. This suggested that some of the mono-nuclear Mn 2+ ions remained in the sample, whose intensity are markedly increased at a low temperature. We propose that Mn 2+ ions that cannot be removed by non-specific binding to protein or mono-nuclear centre in the active site of Sh-PPase exist in Mn-Sh-PPase as a minor component.
We first analysed a well-isolated peak at the lowest field as indicated by arrow in Fig. 4a. The red circles in Fig. 4b are plots of the double integrated intensity of these signals with temperature. To obtain the accurate exchanged coupling constant (J), the absolute intensity was normalised appropriately. These signals reached the maximum intensity at ~6 K. If these signals would be the transition from S = 2 or higher spin system, the signals from S = 1 whose intensity reach the maximum at below 6 K must be appeared (Figs. 1 and 4b). No such signal, however, was observed below 6 K, except for the signals from mono-nuclear Mn 2+ as an impurity. These observations led us to suggest that this signal (marked by the arrow) was originated from S = 1 state. The Boltzmann curve for S = 1 was calculated by Eq. 3 and the best fitted curve was given an exchanged coupling constant, J = −0.85 ± 0.05 cm −1 (Fig. 4b. red line). This value was the smallest compared with other di-Mn enzymes ( Table 1) 21,28 . Most importantly, this is consistent with the X-ray results. In the crystal structure of Mn-Sh-PPase in the absence of substrate, the distances between bridged water and M1/M2 were considerably longer than other di-Mn enzymes, including other family II PPases. This longer metal-O (bridged) distances yield a considerably weak antiferromagnetic coupling in the di-Mn centre. Our observations reveal that the di-Mn centre of Mn-Sh-PPase in both single crystal and solution forms has longer distance with bridged water. This suggests that the active site of Mn-Sh-PPase has an intrinsically 'loose' structure. Cold-adapted enzymes, in general, is thought  www.nature.com/scientificreports www.nature.com/scientificreports/ to have flexible structure so that it can maintain the perturbation for the reactivity even at cold temperature. We conclude that the 'loose' structure of the active site in Mn-Sh-PPase allows binding of substrate and conformational change of active site required for the catalysis at low temperature, which is a key to cold adaptation mechanism of Sh-PPase.

Determination of axial ZfS parameters for S = 2 and di-Mn distance. EPR-derived parameters
were rarely directly combined with the structural information. Dismukes et al. developed a procedure to directly estimate the Mn-Mn distance from an empirically linear correlation between Mn-Mn distances and the axial ZFS parameters D (D 2 ) from S = 2 28 . As shown in Figs. 3 and 4a, however, it was difficult to assign signals from S = 2 due to the overlaying of several EPR transitions. To identify the signals from S = 2, the intensities of temperature dependence were plotted and compared with S = 2 Boltzmann curve with J = −0.85 cm −1 (Fig. 4b, blue line). The temperature dependence of the peak marked with the asterisk in Fig. 4a was the best fit with the curve of S = 2 ( Fig. 4b, blue squares). These signals were overlaid with some signals from other transitions of the di-Mn centre and mono-nuclear Mn impurity, so the intensities plotted in Fig. 4b were simply obtained from the peak height.
We simulated EPR spectra using spin Hamiltonian treated with single spin system, such as S = 1 and 2, as calculated using Eq. 2. The overall pattern of the EPR spectra could not reproduced by simulation from only S = 1 and 2 because experimental spectra were contributed from higher exited states of the di-Mn centre (S = 3, 4 and 5) and mono-Mn impurity. Since the estimation for the magnitude of D-value is sensitive to the peak position, we used EPR simulation to fit the peak position with experimentally assigned peaks of S = 1 and 2. The best fit simulation for S = 1 and S = 2 as shown in Fig. 5 and Supplementary Fig. S2 were given the axial ZFS parameters, |D 1 | = 0.248 ± 0.010 cm −1 and D 2 = −0.0575 ± 0.002 cm −1 , respectively. According to the linear correlation between Mn-Mn distances of X-ray crystal structures and the axial ZFS parameters for S = 2, the D 2 value was given to 3.57 ± 0.03 Å for the Mn-Mn distance of the Sh-PPase, which is comparable to most of other di-Mn enzymes and is within the range of linear correlation (Table 1, and see Supplementary Fig. S3) 21,28 . The EPR-derived Mn-Mn distance of Mn-Sh-PPase was comparable with the X-ray-derived distance of all family II PPases. In addition, we prepared two more samples of Mn-Sh-PPase and collected EPR spectra at 15 K (see Supplementary Fig. S4). The experimental errors were estimated for peak positions and D 2 for ~0.9 mT and 0.001 cm −1 , respectively. This error is corresponded to ~0.004 Å Mn-Mn distance, which is within a deviation of the empirically linear correlation between D 2 and Mn-Mn distance (±0.03 Å).
Substrate binding effect for Sh-PPase. Family II PPase only weakly converts the substrate analogue PNP to Pi 4,9 . Therefore, for the EPR measurement, enzyme concentration was decreased to ~100 µM to prevent consumption of PNP. All samples of Mn-Sh-PPase with substrate analogues were prepared by taking the sample from the same aliquot as that for Mn-Sh-PPase in the absence of substrate analogue and their concentration was decreased by addition of the buffer without any metals. Mn-Sh-PPase were then mixed with the solution with saturated PNP or methylendiphosphonic acid (PCP) in an EPR quartz tube and transferred to liquid nitrogen within 30 seconds. The CW X-band EPR spectrum at 15 K for Mn-Sh-PPase in the presence of PNP is denoted in red in Fig. 6a (red). The overall line-shape was slightly changed. In particular, the signal from S = 1 (Fig. 6a, arrow) was www.nature.com/scientificreports www.nature.com/scientificreports/ broadened and shifted to a higher magnetic field compared with Mn-Sh-PPase without substrate, indicating that the D 1 value became smaller. Importantly, the characteristic 45 mT hyperfine structure of the EPR spectrum at 40 K for the Mn-Sh-PPase with PNP was evident, even at a lower concentration of Mn-Sh-PPase (Fig. 6a, inset). This confirmed that the solution of Sh-PPase had a di-Mn centre in the active site in the presence of substrate, while the crystal structure was unresolved for Mn-Sh-PPase with substrate. The antiferromagnetic coupling constant of Mn-Sh-PPase with PNP was estimated by the signal intensity of S = 1, yielding J = −1.3 cm −1 , which was larger than the value for Mn-Sh-PPase without substrate (Fig. 6b, Supplementary S5a and Table 1). In addition, the signal from S = 2 of Mn-Sh-PPase complexed with PNP was shifted to ~9 mT lower magnetic field than that of substrate unbound enzyme (see Supplementary Fig. S6). This shift yielded a larger D 2 value, which corresponded to a slightly shorter Mn-Mn distance (3.52 Å) according to an empirically linear correlation (see Supplementary  Fig. S3). Since the Mn-Mn distance changed only marginally, the larger magnitude of the J value for Mn-Sh-PPase complexed with PNP suggested that a bridged water was closer to the di-Mn centre when substrate bound. In fact, the crystal structure of Mg-Sh-PPase with PNP displayed a shorter distance between M1/M2 and the bridged fluoride ion. This active site rearrangement is required for efficient reactivity. Our observations support the conclusion that the 'loose' structure of the active site is induced to a 'well-tuned' structure that is necessary for reactivity by binding of substrate. Hence, cold-adapted Sh-PPase is able to bind the substrate at a cold temperature due to the 'loose' active site. Once the substrate is bound, rearrangement of active site occurs simultaneously with the 'well-tuned' structure for nucleophilic attack.
The Southern Ocean is considered the largest high-nutrient low-chlorophyll (HNLC) region, with lower concentrations of minerals, such as transition metals. The 'loose' structure of Sh-PPase might be concerned with facilitating broad metal selectivity for Sh-PPase. This might confer a survival advantage on Shewanella sp. AS-11 in the metal-limited Southern Ocean.
In the sulphate-bound crystal structure of family II PPase from Streptococcus gordonii (Sg-PPase), sulphate ion bound in the substrate binding site and its conformation was in the 'closed' state like the substrate-bound form, www.nature.com/scientificreports www.nature.com/scientificreports/ indicating that sulphate ions might mimic the substrate or product in the active site 7 . To validate that sulphate ion can also bind to the substrate binding site of Sh-PPase, ammonium sulphate was mixed with Mn-Sh-PPase and EPR measurements were performed. Since ammonium sulphate was not the real substrate, higher concentration of enzyme could be used for EPR measurements. The EPR spectrum for Mn-Sh-PPase with sulphate ion was almost identical to that for PPase with PNP (Fig. 6a, grey). 45 mT hyperfine splitting at wide range of the spectrum was also observed, suggesting that the di-Mn centre was maintained. The antiferromagnetic coupling constant of Mn-Sh-PPase with sulphate was also estimated by the signal intensity of S = 1, yielding J = −1.2 cm −1 (Figs. 6b, S5b and Table 1) and the signal position of S = 2 (see Supplementary Fig. S6) was also identical with bound PNP. These observations suggested that sulphate ion bound correctly to the substrate binding site. We also performed EPR measurement for Mn-Sh-PPase complexed with another substrate analogue, PCP, which is used as the inhibitor for family II PPases 40 . The EPR spectrum for Mn-Sh-PPase with PCP clearly showed strong 6 hyperfine splitting characteristic of the mono-Mn 2+ ion, and no signal was observed from the di-Mn centre (Fig. 6a, orange). These observations suggested that binding of PCP induced the decomposition of the di-Mn centre. Thus, the inhibitory mechanism of PCP was different from that of PNP. In addition, the line-width of the EPR spectra of Mn-Sh-PPase with PNP and sulphate ion were somewhat broadened, suggesting that binding of substrate increased the perturbation of the di-Mn centre. Thus, we propose that the increased perturbation of the metal sites by substrate binding plays an important role for activation of bridged water for nucleophilic attack to the PPi substrate. conclusions X-ray crystallography and EPR spectroscopy were used to investigate a cold-adapted family II Sh-PPase from a psychrophilic bacterium that lives in the Southern Ocean. We successfully determined the crystal structure of Mn-Sh-PPase and Mg-Sh-PPase complexed with PNP. The overall structure consisted of N-and C-terminal domains, and the overall conformation was changed from the 'open' to 'closed' state by the binding of substrate, as occurs with other family II PPases. EPR analysis revealed that the active site in Sh-PPase with and without substrate had a di-Mn 2+ centre. The Mn-Mn distance was estimated to be 3.57 and 3.52 Å for Mn-Sh-PPase alone and complexed with substrate analogues, respectively. The most important observation was that the di-Mn centre of Mn-Sh-PPase in the absence of substrate experienced a considerably weak exchange coupling, which was increased by the binding of substrate analogues. In the X-ray crystal structure of Mn-Sh-PPase, the bridged water was further away from the di-Mn centre, and it was placed near the di-Mn centre in Mg-Sh-PPase with PNP. These observations suggest that the active site of Sh-PPase is a 'loose' structure and that the rearrangement of the active site upon binding of substrate generates a 'well-tuned' structure. Thus, we conclude that the 'loose' active site in Sh-PPase plays an important role for binding of the substrate at low temperature. This structural character is an advantage for the adaptation of Sh-PPase to a cold environment, such as the Southern Ocean, since the www.nature.com/scientificreports www.nature.com/scientificreports/ enzyme must maintain flexibility for catalysis at a cold temperature. Once the substrate is bound, the active site is changed to the 'well-tuned' structure that is necessary for the activation of nucleophilic water. Additionally, we propose that the 'loose' structure also facilitates the use of a wide range of metals, which confers a survival advantage on Shewanella sp. AS-11 to survive in the metal-limited condition of the NHLC region. In summary, X-ray crystallography and EPR spectroscopy allows viewing of the active site structure in Sh-PPase in both the single crystal and solution forms, and the rearrangement of the active site upon substrate binding. The present findings reveal that the 'loose' structure of Sh-PPase is important in the cold adaptation mechanism. However, free Mn 2+ ions still remained in the solution as indicated by a strong EPR signal with 6 hyperfine structure at g~2, which was characteristic of the mono-Mn 2+ ion (see Supplementary Fig. S1). To completely remove free or non-specifically bound Mn 2+ ions, the Mn-Sh-PPase protein was applied to the desalting spin column (APRO Science) that was equilibrated using a buffer comprised of 100 mM Tris-HCl (pH 7.5), 15 mM MgCl 2 , 50 mM KCl and 20 µM EDTA. After this treatment, the EPR spectrum was predominated by a complicated 11 hyperfine splitting characteristic of the di-Mn centre. The concentration of the enzyme was calculated from the absorbance at 280 nm (ε 280 nm = 23.3 mM/cm). PNP and PCP (both from Sigma-Aldrich) were used as an analogue of the substrate. As a note, Bs-PPase is also able to convert PNP, though its rate is much slower (~10 4 -10 5 times) than that for PPi 4,9 . Mn-Sh-PPase for substrate binding studies were used from the same aliquot as that for the unbound substrate, but the final protein concentrations were reduced to prevent the consumption of substrate analogues. Mn-Sh-PPase was transferred into a quartz EPR tube, mixed with analogues in the tube, and then frozen in liquid nitrogen within 30 seconds. X-ray data collection, phase determination and structure refinement. Crystals were soaked in the cryoprotectant containing reservoir solution with additional 30% glycerol. X-ray diffraction data were collected at SPring-8 beamline BL41XU using a wavelength of 1.0 Å. The diffraction images were indexed, integrated and scaled using X-ray detector software (XDS) 41 . The data collection and refinement statistics are shown in Table 2. The initial phase of the Mn-Sh-PPase crystal was obtained by the molecular replacement (MR) method. Crystal structure of family II PPase from Methanococcus jannaschii (PDB code 2EB0) was used as a search model in the Phaser program 42 . The initial phases of PNP and Mg-Sh-PPase were obtained by MR using the model of Mn-Sh-PPase and then automatically rebuilt using the ARP/wARP program 43 . The models were improved by multiple rounds of manual rebuilding using Coot 44 and a restrained refinement using Refmac5 45 in CCP4 46 package. epR spectroscopy. X-band EPR spectra were obtained using a Bruker ELEXSYS E580 spectrometer at the Analytical Research Centre for Experimental Sciences, Saga University. The spectra were acquired in the CW mode at ~9.597 GHz using an Oxford Instruments ESR 910 continuous helium flow cryostat. Typical experimental parameters were 2 mW microwave power, 100 kHz field modulation, 10 G modulation amplitude, ~24 G/s sweep rate and 164 ms time constant. epR analysis. The di-Mn 2+ centre has 'ladder' spin states from S = 0, 1, 2 … 5 with a weak antiferromagnetic exchange coupling (Fig. 1) 26,27 . This system can be described with spin Hamiltonian (Eq. 1) as follows:

Methods
where J is the isotropic exchange coupling constant between di-Mn ions; β is the electron Bohr magneton; g i , S i , D i , A i and I i are the g tensor, electron spin operator, axial ZFS parameter, hyperfine coupling constant and nuclear spin operator for each Mn 2+ ion, respectively. The first term is Heisenberg exchange coupling, the second term is dipolar exchange coupling, and the third term corresponds to Zeeman splitting, single ZFS, and hyperfine coupling for each Mn 2+ ion. The spin Hamiltonian diagonalization yields all energy levels. However, calculation of the energy levels of the di-Mn 2+ system (S 1 = S 2 = 5/2, I 1 = I 2 = 5/2) is complicated, because it requires the matrix diagonalization of 6 × 6 × 6 × 6 versus 6 × 6 × 6 × 6. Dismukes et al. reported an empirically linear correlation between the Mn-Mn distance from X-ray crystallography and EPR-derived ZFS parameter from the second excited state (D 2 ) ( Fig. 1) 28 . So, in this paper, the spin Hamiltonian was simplified as a single total spin system, such as S = 1, 2 … with two nuclear hyperfine coupling of two I = 5/2 (Eq. 2). In addition, the sign of the axial (D 2 ) and rhombic (E 2 ) ZFS parameters for S = 2 were assumed to be negative and zero, respectively. The axial ZFS parameters for S = 1 and 2 were obtained by the best fit of the peak position with EPR simulation generated by the EasySpin 5.2.20 program 47 operating in Matlab.

B g S D S S S E S S S A I S A I
The temperature dependence of EPR signal intensity was proportional to the weight of the Boltzmann population for total spin states. The Boltzmann curves for each total spin system were calculated by Eq. 3 as follows: where J is the isotropic exchange coupling, S is total spin and S i is each Mn spin. To obtain the J value, the intensities of the EPR signal with temperature were plotted and normalised appropriately, and the Boltzmann curve was calculated to fit the experimental plots.  Table 2. X-ray data collection and refinement statistics of Sh-PPase. a Values in parentheses are for the highestresolution shell. b R meas = Σ hkl {n/(n−1)} /2 Σ i | I i (hkl)-<I(hkl)>|/Σ hkl Σ i I i (hkl), where n is the multiplicity of reflection hkl, and <I(hkl)> is the average intensity of i observations. c R work = Σ hkl | F obs (hkl) -F calc (hkl)|/ Σ hkl F obs (hkl), where F obs and F calc are the observed and calculated structure factors, respectively. R free was calculated with 5% of the reflections. d RMSD, root mean square deviation.