A Simulation Independent Analysis of Single- and Multi-Component cw ESR Spectra

1. Introduction

ESR spectroscopy is a useful and powerful tool for studying biological free radicals and transition metal cofactors in proteins [1,2,3,4,5,6,7]. Understanding the electronic structure and local environment in such systems provides important insights in catalytic mechanisms [8,9,10,11,12,13] and redox processes in biological systems [6,14,15,16]. However, extracting information from ESR spectra can be challenging, especially in the case of weak and poorly resolved coupling interactions between studied spins. Traditionally, ESR spectral analysis is carried out by the spectral simulations, followed by user interpretation, which not only introduce bias, but impairs the consistency of such analysis. The low-intensity signal components are difficult to analyze through standard spectral simulations in a variety of cases [17,18,19]. Accurate spectral analysis in such cases require quantum chemical computations or special techniques [20,21,22], but without sufficient expertise and knowledge in the field, the interpretations rely heavily on the researchers’ experience and intuition than the robustness of the method. An even more challenging problem, which occurs frequently, is the presence of more than one structurally similar molecules in a system. For closely resembling species, identifying multiple components by ESR alone could be a daunting task and the standard simulation tools are incapable of extracting such information from the spectra without user manipulation. Additionally, every standard ESR simulation methods requires user defined starting configuration for optimization [23,24], which could lead to overfitting and/or unintended manipulation. Hence, direct extraction of the relevant spin hamiltonian parameters, namely the g-factors and and hyperfine coupling constant (A) values from an experimental ESR spectrum would be ideal. The extracted parameter values can be used directly to interpret the electronic structure or if needed, can be further optimized by standard spectral fitting software.

The capability of the wavelet transform in picking periodic hyperfine pattern from poorly resolved ESR spectra was highlighted previously [25]. In a series of recent publications [26,27,28], we have presented an improved version of the wavelet transform based spectral analysis, which can decouple different frequency components in an ESR or nuclear magnetic resonance (NMR) spectrum, enhancing the resolution and providing an opportunity to extract spectral information or parameters in a selective and objective manner. We have modified the Noise Elimination and Reduction via Denoising (NERD) method [25,29], which is based on the discrete wavelet transform (DWT) method, for separation of the hyperfine lines in cw ESR spectra and extraction of spectral parameters initially [26]. Later, in analyzing ¹H NMR spectra of molecular mixtures, which features highly overlapped resonance lines due to the presence of scalar coupling between intramolecular protons, we recognized that spectral decomposition by the wavelet packet transform (WPT) is superior to DWT in separating the central frequencies from the multiplet structures encompassing them [27,28]. Using the same concept, a cw ESR spectrum is decomposed into its different frequency components by WPT and at an optimum level of decomposition, a detail component in the wavelet domain is used in extracting the hyperfine and/or superhyperfine structure, if any. The process has been explained with an illustrative example later in this work (cf. Method). In case of copper complexes, the analysis of the nitrogen-hyperfine structure is used in determination of both $g_{\perp }$ and hyperfine splitting by copper, $A_{\perp }\left(Cu\right)$.

The outcome of the method has been validated by comparing the extracted hyperfine coupling constant values from both partially resolved and unresolved ESR spectra of Tempol and Tempo, recorded in absence and in presence of oxygen, respectively. We demonstrate the consistency of the WPT based spectral analysis by using different wavelets in our analysis, namely Daubechies (Db6 and Db9) and Coiflet (coif3) wavelets. The method is applied across a wide range of copper complexes and both the copper- and nitrogen-hyperfine coupling constants are recovered from partially resolved and unresolved ESR spectra. The analysis didn’t use any prior knowledge about the structure of the complexes or user defined inputs in deriving the spectral parameters. The derived coordination geometry aligned with the structure predicted from analog studies, independent experiments and/or quantum computation, which further validates the results obtained by the WPT-based spectral analysis. The g- and A-values obtained were compared with the optimized parameters obtained from spectral fitting by EasySpin software for a set of selected cases. The comparison shows that our method remains unaffected by artifacts, such as saturation and the passage effect, which affects simulations or spectral fitting significantly. While the performance of both the methods were at par for well resolved spectra, unresolved spectral feature remained inaccessible to the spectral fitting strategy. In those later cases, to the best of our knowledge, the proposed WPT-based analysis is the only method that can separate the so-called hidden features from poorly resolved spectra and extract the relevant spectral parameters by a direct analysis of an ESR spectrum. In this work, along with describing the extraction of spectral parameters by WPT analysis of the ESR spectra of nitroxides and copper complexes, we demonstrated the efficiency of the method in identifying and analysing multi-component spectra.

2. Method

2.1. Overview of the Wavelet Packet Transform Theory

A continuous wavelet transform can be defined as [30],

$$F(\tau ,s)=\frac{1}{\sqrt{\left|s\right|}}{\int }_{-\infty }^{+\infty }f\left(\delta \right){\psi }^{*}\left(\frac{\delta -\tau }{s}\right)dt$$

(1)

where s is the inverse frequency (or frequency range) parameter, $\tau$ is the signal localization parameter, $\delta$ represents the chemical shift, $f\left(\delta \right)$ is the spectrum, $F(\tau ,s)$ is the wavelet-transformed signal at a given signal localization and frequency and ${\psi }^{*}\left(\frac{\delta -\tau }{s}\right)$ is the signal probing function called “wavelet”. Different wavelets are used to vary selectivity or sensitivity of adjacent frequencies with respect to signal localization. They are not dependent on a priori information of the signal or its characteristics.

Discrete wavelet transform (DWT) is expressed by two sets of wavelet components (Detail and Approximation) in the following way [30]:

$$\begin{array}{cc}\hfill D_j\left[n\right]& =\sum _{m=0}^{p-1}f\left[{\delta }_m\right]2^{\frac{j}{2}}\psi [2^j{\delta }_m-n]\hfill \end{array}$$

(2)

$$\begin{array}{cc}\hfill A_j\left[n\right]& =\sum _{m=0}^{p-1}f\left[{\delta }_m\right]2^{\frac{j}{2}}\varphi [2^j{\delta }_m-n]\hfill \end{array}$$

(3)

where $f\left[{\delta }_m\right]$ is the discrete input spectrum, p is the length of input signal $f\left[{\delta }_m\right]$, $D_j\left[n\right]$ and $A_j\left[n\right]$ are the Detail and Approximation components, respectively, at the $j^{th}$ decomposition level, and $\psi [2^j{\delta }_m-n]$ and $\varphi [2^j{\delta }_m-n]$ are wavelet and scaling functions, respectively. The maximum number of decomposition levels that can be obtained is N, where $N={log}_{2}p$, and $1\le j\le N$. The scaling and wavelet functions, at a decomposition level, are orthogonal to each other, as they represent non-overlapping frequency information. Similarly, wavelet functions at different decomposition levels are orthogonal to each other.

The Detail component $D_j\left[n\right]$ is the discrete form of Equation 1, where j and n are associated with s and $\tau$, respectively. The Approximation component $A_j\left[n\right]$ represent the remaining frequency bands not covered by the Detail components till $j^{\mathrm{th}}$ level. The signal $f\left[{\delta }_m\right]$ can be reconstructed using the inverse discrete wavelet transform as follows,

$$f\left[{\delta }_m\right]=\sum _{k=0}^{p-1}{A}_{j_0}\left[k\right]{\varphi }_{j_0,k}\left[{\delta }_m\right]+\sum _{j=1}^{j_0}\sum _{k=0}^{p-1}{D}_j\left[k\right]{\psi }_{j,k}\left[{\delta }_m\right]$$

(4)

where $j_0$ is the maximum decomposition level from which input signal needs to be reconstructed. Compared to that, both the Approximation and Detail components at each level are further decomposed into a set of Approximation and Detail components. A schematic diagram of DWT and WPT decomposition against increasing levels are shown for comparison in Figure 1 [27].

2.2. A Case Study: WPT Analysis of cw ESR Spectrum of Tempo

In this section we explain the details of WPT spectral analysis by using the cw ESR spectrum of Tempo as an example, shown in Figure 2. The ESR spectrum of Tempo is split into three major lines, corresponding to the hyperfine coupling of the ¹⁴N nucleus. Each of those lines are further split by the interaction of the ¹H nuclei in the molecule, giving rise to what we designate as the superhyperfine splitting. Our goal is to separate the superhyperfine or the fine structure of the Tempo spectrum from the rest. The first level of wavelet decomposition by the DB9 wavelet in the Figure 2 shows a pair of approximate (A1) and detail (D1) components. It can be seen that D1 comprises of noise and the entire spectrum is represented by A1. Hence, D1 and all the components that derive from D1 during successive decomposition are not used in the spectral analysis. Decomposition of A1 is continued till level-4, where the approximation component of AAA3 shows no superhyperfine splitting and correspondingly, the approximation component derived from the decomposition of DAA3 contains the superhyperfine splitting pattern. This is why, in this case, decomposition till level-4 is considered to be optimal in separating the hyperfine and superfine splittings of Tempo’s ESR spectrum. For cases, where the superhyperfine is not resolved or visible, the optimum level of decomposition is chosen to be 4 based on previous work [26].

3. Materials

A series of experimental X-band ESR spectra were used in our analysis and the corresponding molecular structures are given in Table 1. ESR spectra of Tempol and Tempo are very well studied and we presented their spectral analysis by WPT for validation purpose. Mutation of superoxide dismutase-1 (SOD1) is arguably correlated with the incidence of significant fractions of familial, and spontaneous cases of amyotrophic lateral schlerosis (ALS) disease [31,32]. For one of its mutants, SOD1:H48Q, cw ESR spectrum was collected at 9.26 GHz, at a temperature of 30 K and concentration in the range of 300 to 350 $\mu \mathrm{M}$ using 50 $\mu \mathrm{L}$ 0.4 mm glass capillaries. No cryoprotectant was added given the viscous nature of the sample. Due to the presence of a glutamine residue in place of a histidine in the mutant, the nitrogen-hyperfine structure along the $g_{\Vert }$ becomes partially resolved as a result of the reduction of the number of nitrogen in the copper coordination sphere. Cu-AHAHARA spectra were collected at 9.39 GHz for two different temperatures, 10 K and 100 K. ESR measurements of CuQu and CuQuA were done in DMSO solutions at 100 K with an ESR frequency of 9.32 GHz.

3.1. Experimental Section

3.1.1. Synthesis

• Copper quinoline (CuQuA): 160.0 mg (1.0 mM) 2-amino-8-quinolinol was taken in an 80.0 mL Schlenk flask with a magnetic stirrer bar and 15.0 mL methanol was added to make it a homogeneous yellow colour solution. 85.0 mg (0.49 mM) ${\mathrm{CuCl}}_2.2{\mathrm{H}}_2\mathrm{O}$ was taken separately with 10 mL MeOH in another Schlenk flask under the steam of ${\mathrm{N}}_2$. Then, the two methanolic solutions were charged under nitrogen flow. An immediate dark brown solution appeared after adding the metal salt to the ligand. Reducing the solvent with continuous ${\mathrm{N}}_2$/Ar flashing whitish deep brownish precipitate was observed. The reaction mixture was further stirred for two hours for completion with ${\mathrm{N}}_2$ purging. Then the precipitate was collected, washed with hexane and diethyl ether and dried under a high vacuum (yield= 137.0 mg, 76% w.r.t. 2-amino-8-quinolinol). XRD-suitable crystals (CCDC 2127156) were grown from slow diffusion from methanol/diethyl ether solution.
• Copper quinone (CuQu): Synthesized as reported above, where 8-hydroxyquinoline was utilized as the precursor ligand (yield= 120.0 mg).
• Copper AHAHARA: The AHAHARA peptide was synthesized by manual Fmoc solid-phase synthesis at elevated temperature using Amide Rink resin, Fmoc-protected amino acids and previously reported protocols [33].
• SOD1:H48Q: The recombinant SOD1:H48Q protein, replacing a histidine with glutamine, was expressed and purified as described in a previous work [34].

3.2. ESR Experiments

Oxygen-free samples of Tempol and Tempo were prepared by flame-sealing after a triple repeat of the freeze-thaw cycle under vacuum. Oxygen saturated samples were prepared by passing oxygen gas through the radical solution in water for 10 minutes and kept under oxygen atmosphere during measurements. Samples with an intermediate oxygen concentration were prepared and handled in air. The nitroxide concentration in all the samples was 100 $\mu \mathrm{M}$. The ESR spectra were recorded at 293K at a microwave (MW) frequency of 9.33 GHz, power of 0.1 mW, modulation frequency of 100 kHz and modulation amplitude of 0.1 G. The ESR spectrum of SOD1:H48Q was recorded at 30K at an MW frequency of 9.26 GHz, power of 0.06325 mW, modulation frequency of 100 kHz and modulation amplitude 4 G. The copper concentration was estimated to be 50 $\mu \mathrm{M}$. The ESR spectra of Cu-AHAHARA complex were acquired in Wilmad tubes using a Bruker Elexsys E500 EPR spectrometer equipped with a cryostat. Peptide stocks at pH 2 were prepared fresh by adding lyophilized solid (>90%) to 10 mM HCl until peptide concentration reached 1 mM. First, Cu(II) solution in water (1 mM, 75 $\mu$L) and peptide stock (1 mM in 10 mM HCl, 150 $\mu$L) were mixed and then buffer (91 mM Hepes, pH 8, 275 $\mu$L) was added to make a solution with 150 $\mu$M Cu(II) and 300 $\mu$M peptide (500 $\mu$L final volume). The mixture was incubated at room temperature for 2 hr, and after glycerol was added to a final concentration of 10%, transferred into an ESR tube and flash frozen in liquid nitrogen. The final pH of the sample was 7.6 as measured by a Spintrode electrode (Hamilton). ESR spectra were acquired at 100 K and 10 K using the following conditions: frequency 9.39 GHz, power 5 mW (D) or 2 mW (C), modulation frequency 100 kHz, modulation amplitude 8 G (10 K) and 4 G (100 K), time constant 163.8 ms. The ESR spectra of CuQu and CuQuA were collected using the Bruker cw ESR EMX spectrometer at the National Biomedical Resource for Advanced ESR Spectroscopy (ACERT). About 100 mg of the air-stable copper complexes were dissolved in 1 mL DMSO. The solutions were poured in 4 mm ESR tubes, frozen in liquid nitrogen and kept in liquid nitrogen Dewar until inserting them in the spectrometer. The spectra were recorded at 100 K, a microwave frequency of 9.316 GHz and attenuation of 30 dB for CuQu and 50 dB for CuQuA.

3.3. ESR Spectral Mix

For multi-component spectral analysis, four mixed spectra were calculated by mixing the ESR spectra of the compounds, CuQu and CuQuA, in proportions of (A) 2:1, (B) 4:1, (C) 10:1 and (D) 20:1 (Figure 9).

4. Results and Discussion

4.1. Validation of the Analysis

Interpretation of partially resolved ESR spectra by fitting has been the standard practice in the study of organic radicals and paramagnetic metal complexes. The approximations and the principles of such optimizations remain unknown to the users, whose interpretation of the results is largely dominated by the goodness of the fit and a priori knowledge about the structure or nature of the molecule under investigation. In contrast, a direct extraction of parameters, a process that lacks the visual aids, can raise queries about the accuracy of the outputs. That was our motivation to run the analysis on a series of X-band ESR spectra of Tempol and Tempo under varying oxygen concentration, shown in Figure 3. The superhyperfine lines are fully or partially resolved in absence of oxygen or when it is present in very low concentration. However, the features become invisible due to line broadening by increased oxygen concentration. Consequently, while the superhyperfine constant in the former cases could be obtained by simple visual inspection or spectral fitting, none of those are applicable in the latter. In all those cases, the WPT-based method recovered the superhyperfine structure from the corresponding ESR spectra and extracted the value of the coupling constant, illustrating both its advantage over the existing spectral analysis approaches and its consistency.

4.2. Robustness Against Wavelet Selection

An important and non-trivial variable in the WPT analysis of ESR spectra is the selection of wavelet or in other words, understanding the robustness of the analysis against the types of wavelets used in the analysis. We repeated the analysis shown in Figure 3 with three different types of wavelets; Coiflet-3 (Coif3), Daubechies-6 (Db6) and Daubechies-9 (Db9), and the results are summarized in Figure 4. The extracted superhyperfine constant values for Tempol and Tempo, $0.44\pm 0.01\mathrm{G}$ and $0.20\pm 0.01\mathrm{G}$, showed insignificant variation for all three wavelets, while producing the same number of splitting across all the analysis consistent with the molecular structure. This observation illustrated the robustness of the method, while further validating the accuracy of the extracted parameters. For the rest of the analysis in this work, we used Db9 wavelet.

4.3. Spectral Analysis of Partially Resolved ESR Spectra

We started our analysis by fitting the X-band cw ESR spectra of SOD1:H48Q recorded at 30K and Cu-AHAHARA, recorded at 10K and 100K by using the EasySpin software. The optimized simulations along with the experimental spectra are shown in Figure 5. For SOD1:H48Q, the EasySpin analysis yielded nitrogen-hyperfine constant values of 15.6 G (axial) and 9.0 G (parallel), while the simulations for Cu-AHAHARA spectra did not require any nitrogen-hyperfine parameter in fitting the spectra. It should be noted that the Cu-AHAHARA spectrum at 10K shows no nitrogen-hyperfine splitting; however, partial nitrogen-hyperfine splitting is evident for the spectrum at 100K. This apparent anomaly might have resulted from the partial saturation and passage effect [35,36] in case of the former, as evidenced by the behavior of the first integral of the spectrum, and consequently, the EasySpin fit for the 10K spectrum performed poorly compared to that of the 100K spectrum. In addition, the simulation presented for SOD1:H48Q in Figure 5[A] emphasized the potential presence of a second component, with a slightly shifted $g_{\Vert }$ and/or different hyperfine coupling constants.

In Figure 6, we presented the WPT analysis of the ESR spectra of (A) SOD1:H48Q, (B) Cu-AHAHARA at 10K and (C) Cu-AHAHARA at 100K. For all the analysis, the spectra were decomposed to level-4 using the DB9 wavelet. For SOD1:H48Q, our analysis yielded a nitrogen-hyperfine coupling constant of 14.1 G for the axial peaks in the range of 3200 to 3320 G. However, an analysis of the $g_{\Vert }$ splitting (Figure 6 [A.2]) revealed two different nitrogen-hyperfine splitting constants; 12.5 G along the dominant Cu parallel hyperfine splitting and 16.5 G along the minor component. Following that, we analyzed the splitting along the small axial peak between 3220 and 3400 G (not shown), which emphasized the presence of four nitrogens and a nitrogen-hyperfine constant of $A_N$ = 16.5 G. From this analysis, we could infer a spherical electron density distribution in the second component, suggesting a tetrahedral copper complex with a large hyperfine splitting of ∼205 G. Similar analysis for Cu-AHAHARA yielded nitrogen-hyperfine coupling constants of 13.9 G (10K) and 14.0 G (100K) along the $g_{\perp }$ signal component. The consistency of the results demonstrate the robustness of the WPT-based spectral analysis against artifacts, over-saturation effect in this case, which is a major advantage over the standard spectral fitting methods.

It should be noted that the copper coordination geometry of the Cu-AHAHARA complex could not be confirmed in a previous work from the ESR spectral analysis [33]. The WPT analysis presented in Figure 6 suggests strong overlap of nitrogen-hyperfine components in the $g_{\perp }$ region of the ESR spectra at 10K and 100K, while the former seems to contain spectral artifacts, making analysis by standard procedure error-prone. Hence, we analyzed the WPT components originated from the $g_{\Vert }$ region for the both the cases, which is summarized in Figure 7. In case of the spectrum collected at 10K, the WPT component originating from the ESR peak centered 2730 G showed nine evenly spaced lines at 12.5 G apart from each other, indicating an ${\mathrm{N}}_4$-coordination for the copper center in the Cu-AHAHARA complex. In addition, a smaller coupling constant in the $g_{\Vert }$ region in comparison to the value obtained in the $g_{\perp }$ region (13.9 G) suggests four equivalent nitrogens in an axial geometry. This interpretation aligns well with a series of independent experimental and theoretical studies conducted to elucidate the geometry of the complex previously [33]. For the spectrum recorded at 100K, a similar analysis in the $g_{\Vert }$ region did not reproduce the exact same results because of unresolved spectral overlapping in the WPT component. However, upon close inspection, we resolved half of the nitrogen-hyperfine splitting window for the $g_{\Vert }$ peaks centered at 2734 G and 2901 G. In this regard, it should be noted that only the even spacing between the peaks in a WPT component around a $g_{\Vert }$ peak was used as the criteria for recovering nitrogen-hyperfine splitting. The WPT analysis is highly accurate in recovering spectral information with respect to their location, but not necessarily the intensity. Factors like partial overlapping of resonance lines, residual noise and spectral artifacts affects the intensities of the peaks in a WPT component and hence, the recovered superfine splitting is unlikely to reproduce the relative intensity pattern expected for perfectly resolved spectra. The analysis of the Cu-AHAHARA ESR spectrum at 100K yielded the same nitrogen-hyperfine coupling constant of 12.5 G and suggested an ${\mathrm{N}}_4$ coordination for the copper center.

The WPT analysis for the X-band ESR spectra of two copper-nitrogen complexes, CuQu and CuQuA, are shown in Figure 8. The detail component in the optimal level of decomposition, which was 3 for both the spectra, revealed unresolved hyperfine structure due to nitrogens coordinated to the copper centers in the complexes. For example, Cu(II) splits the resonance line at $g_{\perp }$ into four lines, we will call them h-1, h-2, h-3 and h-4 for explanation purpose. Each of these lines split further due to the n interacting nitrogens into $m=(2\times n+1)$ lines. Given the small magnetic field window between h-1 and h-4, as well as the nitrogen-hyperfine coupling, most of the resonance lines cannot be resolved due to complete or partial overlapping. However, for most of the cases, it might be possible to identify the first $m/2$ lines originating from h-1 and the last $(m+1)/2$ lines originating from h-4, where the extent of overlapping of resonance lines is the least. Using this logic, from the detail component analysis for CuQu in Figure 8A, we identified 3 lines with a total separation of 53 G, which corresponded to an $m=5$ or 2 coupled nitrogens with h-1 and h-4 at 3203 G and 3314 G. Further, we calculated the nitrogen-hyperfine coupling constant to be $A\left(\mathrm{N}\right)=53/2$ G or 26.5 G, the Cu(II) hyperfine coupling $A\left(\mathrm{Cu}\right)=(3314-3203)/3$ G or 37 G and the $g_{\perp }=(3314+3203)/2$ G $\equiv 2.0427$. The calculated spin-Hamiltonian parameters are in good agreement with previously reported analogous ${\mathrm{N}}_2{\mathrm{O}}_2$-coordinated copper(II) complexes [37,38]. With a similar analysis for CuQuA, Figure 8B, we inferred that 3 nitrogens were coordinating with the Cu(II) center with $A\left(\mathrm{N}\right)=21.3$ G, while the $g_{\perp }$ and $A\left(\mathrm{Cu}\right)$ were calculated to be 2.0411 and 23.3 G. The original ${\mathrm{N}}_2{\mathrm{O}}_2$ coordination geometry for the copper center observed in CuQu complex can be expanded in CuQuA due to the presence of peripheral amine in the QuA ligand scaffold. The potential involvement of one amine group is supported by the EPR data, indicating an ${\mathrm{N}}_3{\mathrm{O}}_2$ coordination.

4.4. Analysis of Multi-Component ESR Spectra

It has already been shown in Figure 6 that the spectral analysis of the $g_{\Vert }$ components in the ESR spectrum of SOD1:H48Q revealed the presence of two components. In this section, we present further proof of how the WPT-based spectral analysis can be utilized efficiently in identifying multi-component spectra, even when the components are present in highly disproportionate amounts. A set of four mixed spectra was produced for the analysis by mixing the ESR spectra of CuQu and CuQuA in proportions of (A) 2:1, (B) 4:1, (C) 10:1 and (D) 20:1, shown in Figure 9. The WPT spectral analysis was conducted using Db9 wavelet at decomposition level of 4 and the WPT detail components for the spectral region between 2600 and 2775 G are shown in the insets of Figure 9. In all the cases, 5 lines with a splitting constant of 17 G were recovered between 2610 and 2670 G, which corresponds to CuQu, which is expected because of the abundance of the complex in the mixtures. However, a second component with a splitting constant of 16 G was identified in between 2700 and 2748 G. For the latter, only 4 or 3 lines were visible due to spectral overlap, but it can be separated from artifacts by the consistent positions of the peaks appearing the detail WPT component. It should be noted that this strategy clearly identifies a second minor component directly from the spectral analysis, in this case CuQuA, but it may not be possible to determine the structure of the component solely from the analysis due to spectral overlaps, leading to unrecoverable information loss. However, we like to emphasize that we used poorly resolved X-band ESR spectra of the individual complexes and given such constraints, we believe that it is a major achievement to detect a second component solely from ESR spectral analysis, even when the second component is present at as low a relative concentration as of 5%. In general, such findings from the WPT spectral analysis will help researchers decide on further structural analysis by employing high-frequency ESR and/or complementary techniques.

Figure 9. WPT analysis of ESR spectra of mixtures of CuQu and CuQuA, mixed in proportions of (A) 2:1, (B) 4:1, (C) 10:1 and (D) 20:1, respectively. The insets show the WPT component for the spectral region between 2600 G and 2775 G. Two components were identified in each of the spectral analysis from the difference in nitrogen-hyperfine splitting of 17 G (green double headed arrow) and 16 G (blue double headed arrow).

Figure 9. WPT analysis of ESR spectra of mixtures of CuQu and CuQuA, mixed in proportions of (A) 2:1, (B) 4:1, (C) 10:1 and (D) 20:1, respectively. The insets show the WPT component for the spectral region between 2600 G and 2775 G. Two components were identified in each of the spectral analysis from the difference in nitrogen-hyperfine splitting of 17 G (green double headed arrow) and 16 G (blue double headed arrow).

5. Conclusions

In this work, we have provided the complete recipe for a simulation-free analysis of cw ESR spectra using a wavelet packet transform-based algorithm and experimental X-band ESR spectra of organic free radicals (Tempol and Tempo) and metal complexes (copper (II) inorganic and biological complexes). We showcased three major accomplishments of the WPT-based spectral analysis over standard simulations: (1) direct extraction of nitrogen-hyperfine structure from poorly resolved/unresolved spectra, (2) insensitivity toward spectral artifacts and (3) identification of multi-component samples. Consistent extraction of superhyperfine coupling constant due to protons from both partially resolved and unresolved spectra for Tempol is presented, validating the accuracy of the new spectral analysis technique. The limitations of standard ESR spectral analysis are illustrated by running EasySpin simulations for partially resolved and poorly resolved spectra, with one of the cases containing some experimental artifacts. While the accuracy and efficiency of the simulations varied drastically in those cases, the WPT-based analysis extracted spectral parameters, revealing hyperfine and superhyperfine splittings not detectable by EasySpin, and the analysis was unperturbed by spectral artifacts. For the copper-nitrogen complexes, the resolved hyperfine structures confirmed the number of nitrogen atoms coordinating with copper centers as well as their coordination geometry. Finally, by close examination of the detail WPT components originating from the $g_{\Vert }$-regions of the spectra, the presence of two different components in case of SOD1:H48Q, a mutant of SOD1, has been confirmed along with their nitrogen-hyperfine coupling constants. The strength of the method in resolving multi-component spectra has been displayed further by analyzing four mixed spectra of two copper-nitrogen complexes, CuQu and CuQuA. The spectral analysis revealed the presence of the two components even when the proportion of the two components in the spectrum was 20:1. In addition, three different wavelets were used for selected cases, emphasizing the robustness of the method against the choice of wavelets. The comprehensive analysis presented in this work is expected to motivate broad adoption of the technique in analysing ESR spectra and the spectral parameters obtained from the analysis can be further optimized using simulations and/or other computational methods.