Microwave Saturation of Complex EPR Spectra and Free Radicals of Burnt Skin Treated with Apitherapeutic Agent

The effect of microwave power on the complex electron paramagnetic resonance spectra of the burn matrix after the therapy with propolis was examined. The spectra were measured with microwaves in the range of 2.2–79 mW. Three groups of free radicals were found in the damaged skin samples. Their spectral lines evolve differently with the microwave power. In order to detect these free radical groups, the lineshape of the spectra was numerically analysed. The spectra were a superposition of three component lines. The best fit was obtained for the deconvolution of the experimental spectra into one Gauss and two Lorentz lines. The microwave power changes also the lineshape of the spectra of thermally injured skin treated with the conventional agent—silver sulphadiazine. The spectral changes were different for propolis and for silver sulphadiazine. The number of individual groups of free radicals in the wound bed after implementation of these two substances is not equal. It may be explained by a higher activity of propolis than of silver sulphadiazine as therapeutic agents.


Introduction
During high temperature influence on skin, different types of free radicals are formed. The electron paramagnetic resonance (EPR) spectra may be used to solve the problem of the majority of free radical systems in the skin. Each group of free radicals is responsible for the EPR component line in the resultant resonance spectra [1][2][3]. The information about different free radicals in the thermally damaged skin is fragmentary. The knowledge about the types of free radicals in the skin is important for pharmacotherapy and alternative medicine [2][3][4][5].
The aim of this work was to check the usefulness of microwave saturation of EPR spectra to obtain information about the number of different groups of free radicals in the matrix of burnt skin during the therapy with propolis. The relative examination of the skin treated with the silver sulphadiazine was performed. The application of microwave effect on the lineshape of the EPR spectra in order to study free radicals in skin samples and free radicals in biological systems was discussed. The proposed innovatory analyses to skin samples are important to check the substances used in alternative medicine.

Experimental Animals. Ethics Committee of the Medical
University of Silesia in Katowice approved the study protocol, 2 Evidence-Based Complementary and Alternative Medicine based on the application of two 4-month-old, domestic pigs chosen as the experimental animals. The pig body weight at the beginning of the study was approximately 37.5 kg. Mentioned animal species has been proven to be useful for the assessment of wound healing phenomenon because of numerous analogies of human and pig skin, such as the thickness and composition of epidermis; structure and content of keratinous proteins; size orientation and distribution of the vessels; vascularization of the lower region of the follicle; the similar pattern of deep dermal burn healing [6,7]. The experimental animals were housed according to GLP standards of Polish Veterinary Law. Pigs general condition was estimated by regular control of behavior, body weight, and temperature. Thermal damages were inflicted according to Hoekstra et al. [6] standard model. The animals were divided into control and experimental group, each containing one pig. Control wound was treated with silver sulfadiazine to observe the healing process, occurring after application of the agent of choice in the local burn treatment [8], twice a day, from the 1st to the 5th day of the experiment. Experimental wounds were treated with propolis, twice a day, from the 1st to the 5th day of the experiment. Biopsies were taken five days after burn infliction.
After burns infliction, tissues were rinsed with an antiseptic solution. During the whole experiment, burns were covered with 55-75 mm layer of the propolis and silver sulfadiazine, twice a day. Subsequently, after the treatment with the mentioned therapeutic agents, the injuries were protected with a woven cotton material.

Sample Preparation to EPR Measurements.
Tissue material was mechanically purified, premicronized, and weighed. Then tissue samples were homogenized in acetone using a mechanical homogenizer (30000 rpm/min, 4 ∘ C, 30 min.) to obtain a homogeneous suspension. Subsequently received tissue homogenates were dried and reweighed. For EPR measurements experimental material was placed in the glass tubes. EPR signals were not detected for the empty glass tubes. The walls of the glass tubes were thin and their external diameter was of 3 mm. The mass of the samples located in the tubes was measured.

Detection of EPR
Spectra. An X-band (9.3 GHz) electron paramagnetic resonance spectrometer of Radiopan Firm (Poznan, Poland) was used for detection of free radicals in the skin samples. The EPR spectra were measured with a magnetic modulation of 100 kHz in the wide range of microwave power of 2.2-70 mW. The attenuations were changed from 0 dB to 15 dB. The microwave saturation effect in the EPR spectra, recorded at higher microwave powers, was obtained.

The Numerical
Analysis of the EPR Spectra. The lineshape of the first-derivative EPR spectra of wound bed treated with propolis and silver sulphadiazine was analysed. The experimental EPR spectra were deconvoluted into the component lines. The measured spectra were fitted by a superposition of theoretical lines with Gauss (G) and Lorentz (L) shape. The following sum of theoretical (G, L) lines were tested: GGG, LLL, GGL, and GLL.
The best fitting of the experimental spectra by the theoretical multicomponent lines was searched. The best result was the fitting with the lowest square error.
The parameters of the individual components in the EPR spectra such as amplitude ( ), linewidth (Δ pp ), and -factor were determined. The signal power for the component lines was calculated. The signal power is the percentage fraction of the individual line in the whole spectrum. It was calculated as the percentage fraction of the integral intensity of the individual line in the total integral intensity of the whole spectrum. The signal power gives information about the fraction of the individual type of free radicals in the tested sample.
The amplitude rises with increasing free radical contents in the sample linewidths depend on magnetic interactions and distances between free radicals. The -factor, in turn, is determined by the type of free radicals [9].
-Values were calculated from the EPR spectrum according to the formula [9,10] is Bohr magneton, and is resonance magnetic field. Microwave frequency (]) was measured by MCM101 recorder of EPRAD Firm (Poznan, Poland). The resonance values were obtained from the EPR lines.
Free radical concentration was measured by a double integration of the first-derivative EPR spectra and by a comparison to the integral intensity of the reference-ultramarine. The second reference-a ruby crystal-(Al 2 O 3 :Cr 3+ ) was also used.
The numerical analysis was performed for the signal which is given as a time series ( ), = 1 ⋅ ⋅ ⋅ represents the function values within the domain (−5, 5) and time step Δ = 10 −3 . The analysed signal is complex and cannot be described by a single basic function. For the analysis, the derivative of basic functions was used. The application is searching for the best match of function using function defined as a composition of basic functions. Before the approximation is performed, the signal is preprocessed by filtering and positioning of the reference point (0, 0).
For filtering purposes two methods were researched: the moving average filter [11] and the filter based on fast Fourier transformation [12]. The moving average filter removes the noise; however, the signal is deformed, and therefore finding maximal/minimal function values becomes impossible. The parameter defines a number of neighbors taken under consideration while averaging.
A reference point (0, 0) was found by using a and values which define the distance to move the filtered signal according to and axis. To find the reference values, two heuristic rules were defined: (a) the sum of values of samples separated by axis is close to 0, (b) the absolute value of the sum of functions from the left side of axis plus the absolute value of the sum of functions from the right side of should be maximized. Medicine   3   According to the heuristic rule, a simple linear optimization task was defined for  and  values: min

Evidence-Based Complementary and Alternative
All undefined values (outside the interval (−5, 5)) equal 0 and do not change the results. Equation (1a) is unaffected by value and therefore is estimated using linear approximation. In this case, the simplex method with bounds on the variables was used. Bounds are set empirically for to (−10, 10) and for to (−50, 50). Based on estimated value, the displacement is estimated according to (1b). The filtered and transformed function is processed by the estimation algorithm.
According to the prepared model of numerical analysis, the signal is approximated by the function. The function is a composition of derivatives of continuous elementary functions. The function was defined as follows: where 1 is weight of th function, is derivative of th function, and is quantity of elemental functions, ∈ .
The preliminary analysis of candidate elementary functions considered sinus, Lorentz, Gaussian, triangle, and quadratic functions. The accuracy was measured using root mean squared error (RMSE) [13]. The error was calculated as follows: where is the total number of the analyzed values. After preliminary research, the Gaussian and Lorentz function was selected as it offers the lowest error rate for the number of function equal or below 3. Therefore, the ( ) for = 1 ⋅ ⋅ ⋅ 3 was defined as The parameters define the Gaussian and Lorentz function representation. The parameter 2 defines mean/location parameter while 3 defines the variation/scale parameter.
The best fit was found by minimizing the difference between the real function and function modifying the matrix for given value and Gaussian/Lorentz functions: Boundary values were defined as for 1 , = 1 ⋅ ⋅ ⋅ parameters; the maximal signal value multiplied by a number of functions was defined for 2 , = 1 ⋅ ⋅ ⋅ parameters, which defines the position of function centers and should not exceed the analyzed interval (−5, 5) and 3 , = 1 ⋅ ⋅ ⋅ parameters to (0, 30) according to Gaussian and Lorentz functions requirements. To find the best match, the genetic algorithm [14] with Conjugate Gradient Method [15] was implemented. In the first phase, the algorithm finds a rough estimation of the parameters by modifying the values in matrix within given boundaries and then, for best 10 solutions, the Conjugate Gradient Method is executed to minimize the solution error. The population was set to 1000 and matrixes initial values were equally distributed within defined boundaries. Next, 100 epochs are run to find a rough estimation. The best solution from previous generation is moved to the next one while the remaining solutions undergo a mutation and a crossover. The mutation considers changing one of matrix elements by 1/100 of a parameter value range. The crossover is performed by a simple averaging of two samples. The probability of moving to the next epoch is based on the fitness function proposed by Chan and Sudhoff [14].

Results
The experimental EPR spectra of the skin of burned wounds were best fitted by the superposition of three lines. The shapes of the components were described by Gauss and Lorentz theoretical functions. One Gauss and two Lorentz lines were found in the spectral curves. Different free radical groups are responsible for the individual components in the spectra. The three groups of free radicals with Gauss and two Lorentz lines exist in the tested skin samples.
The components were obtained for the EPR spectra measured with microwave powers of 2.2, 7.0, 11.1, 22.2, 35.1, 55.3, and 70.0 mW. The number of the lines did not depend on the microwave power. The component lines of the EPR spectra of the thermally damaged skin cured with propolis and silver sulphadiazine were presented in Figures 1, 2, 3, 4, 5, 6, and 7. In Figures 1-7(a) the results of the numerical analysis for the skin after the therapy with propolis are shown, while in Figures 1-7 Table 1 for the injuries treated with propolis and, in Table 2, for the damaged skin treated with silver sulphadiazine. The amplitudes and linewidths change with the microwave power. These changes give way to easier detection of these lines in the resultant spectrum. All measured EPR lines are the broadest ones.
-Values were typical for free radicals paramagnetic centers; they were near 2.00. These values are characteristic for unpaired electrons in organic molecular units. A different effect of the microwave power on the component lines was observed for the wounded skin samples treated with propolis and silver sulphadiazine (Figures 1-7, Tables 1 and 2). It indicates a different composition of free radical system in these samples.
A lower free radical concentration (13.9 × 10 22 spin/g) characterizes the burn wound treated with propolis during 5 days rather than the injury treated with silver sulphadiazine (22.99 × 10 22 spin/g). The free radical concentrations of the individual types of free radicals in the wound therapeutically managed with propolis were also lower than for that treated with silver sulphadiazine.

Discussion
Different types of free radicals were found in the tested skin samples. Their EPR spectra reveal a multicomponent shape (Figures 1-7). The deconvolution of the EPR spectra into component lines of skin samples is difficult to perform because of high linewidths of the components (Tables 1 and  2). The individual lines overlap in the total spectrum and the resonance magnetic fields are similar to them. As a result, the individual lines are not visible in the EPR spectra and they should be numerically obtained. In this work the numerical procedures of deconvolution of the EPR spectra were additionally made easier by the use of the microwave saturation effect of the EPR lines. The microwave power increase led to the enhancement of the amplitude of the EPR lines because higher numbers of unpaired electrons come to higher excited energy levels [9]. However, this effect is observed only up to the maximal values of the microwave power and, after higher microwave powers, the saturation of the EPR lines occurs and its amplitude decreases with a former increase of microwave power. The microwave saturation effect-the decrease of amplitudes of the EPR line with increasing microwave power is the result of the inversion of distribution of unpaired electrons on the energy levels [16]. Higher amounts of unpaired electrons are located on higher energy levels rather than on ground levels. The microwave power increases, but the amount of electrons on the ground energy levels is too low. So, the EPR lines are quenched instead of being increased as it is expected. This effect is very useful to find the number of components in the complex EPR spectrum. The described above microwave saturation effect interacts on all types of free radicals in the sample, so it changes their EPR lines as was presented. However, the most important fact is that the microwave saturation effect appears at different microwave powers for the individual component lines. The majority of the microwave saturation effects in different groups of free radicals cause changes of the lineshape of their resultant EPR spectra. The numerical analysis of the EPR spectra, recorded with different microwave powers, and the comparison of their components make the determination of the individual lines easier. The described above information about the microwave saturation effect on EPR spectra was used in this work to analyse different groups of free radicals in burn wound after treatment with propolis. Our expectation is confirmed the three groups of free radicals and their three EPR lines were obtained by mathematical procedures. The number of the EPR lines is stable independently on the microwave power. (Figures 1-7). As it was expected, the parameters of the lines strongly changed with the microwave power (Tables  1 and 2 lines give information mainly about delocalized unpaired electrons. The microwave saturation effect was helpful in both kinds of tested samples, the injured skin treated with propolis and silver sulphadiazine. The electron paramagnetic resonance spectroscopic method, the microwave saturation effect, and the numerical analysis of the lineshape of the EPR spectra are proposed not only for the test of propolis influence on the skin but also they may be applied to examine free radicals in the drugs or tissues with several types of free radicals. The performed studies broaden our knowledge about the influence of propolis on free radicals in skin but the practical meaning of this study is equally important. The methods of the free radical examinations to alternative medicine purposes were presented.

Conclusions
The advanced spectroscopic method was applied to examine different groups of free radicals in thermally damaged skin treated with propolis. It was proved that the numerical analysis of the lineshape, continuously saturated by microwaves EPR spectra, is useful to determine the number of groups of free radicals in applied samples. These studies broaden our earlier numerical analysis of the EPR spectra measured only with one low microwave power. The microwave saturation of the complex EPR spectra is helpful to examine the multicomponent structure of free radicals in skin samples because of a different effect of the microwave power on the individual component spectral lines. The analysis for propolis was compared with the results of studies for the multicomponent spectra of the standard substance such as silver sulphadiazine. The EPR spectra were best fitted by a superposition of one Gauss and two Lorentz lines. The numerical analysis and microwave saturation of EPR spectra of burnt skin pointed out that the same number of groups of free radicals exists in the skin treated with propolis and silver sulphadiazine. However, the amount of free radicals of different groups in these samples is different. It is expected that it is the reason of a better therapeutic impact of propolis (containing compounds responsible for antioxidant, antimicrobial, antiinflammatory, antigenotoxic, antiangiogenic, and anticancer properties [17]) on the injured skin than of silver sulphadiazine, which was described in the other work [18]. These studies confirmed the usefulness of both the numerical analysis of the lineshape of EPR spectra and the microwave saturation of the spectra in the examination of complex system of free radicals in skin. The application of the tested methods may be broadened to other biological samples important in alternative medicine.