Potential Antioxidative and Anti-Hyperuricemic Components Targeting Superoxide Dismutase and Xanthine Oxidase Explored from Polygonatum Sibiricum Red.

Polygonatum sibiricum Red. (P. sibiricum) has been used as a traditional Chinese medicine with a wide range of pharmacology effects. However, the responsible bioactive compounds and their mechanisms of action concerning its antioxidative and anti-hyperuricemic activities remain unexplored. In this work, the antioxidant capacity of P. sibiricum was firstly evaluated with the 2,2-diphenyl-1-picrylhydrazyl (DPPH), 2,2’-azinobis-(3ethylbenzthiazoline)-6-sulfonic acid (ABTS) and ferric-reducing antioxidant power (FRAP) assays, from which the ethyl acetate (EA) fraction exhibited the highest DPPH, ABTS radical scavenging, and ferric-reducing capacities. Meanwhile, the EA fraction displayed the highest total phenolic and flavonoid contents among the four fractions. Next, the potential ligands from the EA fraction were screened out by bio-affinity ultrafiltration liquid chromatography-mass spectrometry (UF-LC-MS) with superoxide dismutase (SOD) and xanthine oxidase (XOD). As a result, N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, N-trans-feruloyltyramine were identified as potential SOD ligands, while N-cis-p-coumaroyltyramine was determined as potential XOD ligand. Additionally, these four ligands effectively interact with SOD and XOD in the molecular docking analysis, with binding energies (BEs) ranging from –6.83 to –6.51 kcal/mol, and the inhibition constants (Ki) from 9.83 to 16.83 μM, which were better than the positive controls. In conclusion, our results indicated that P. sibiricum has good antioxidative and anti-hyperuricemic activities, and its corresponding active ligands targeting SOD and XOD could be explored by the UF-LC-MS method.


Introduction
Polygonatum sibiricum Red. belonging to the Liliaceae family has been distributed in the temperate northern hemisphere countries, such as China, Japan, Korea, India, Russia, Europe, and North America [1,2]. In China, P. sibiricum is well-known as a traditional medicinal herb and functional food [3,4], as well as a health-improving substance [5], with a sweet fragrance and taste. P. sibiricum has been included in the "list of items that are both food and medicine" by the National Health Commission of the People's Republic of China. It was initially described as replenishing Qi and nourishing Yin, strengthening the spleen, and nourishing the lungs and kidneys in "Special Records of Famous Doctors". Thus far, it has been reported to have a variety of pharmacological applications and biological activities, including antioxidative and anti-hyperuricemic activities. Scientific reports have supported 17) was preserved in the herbarium of the Key Laboratory of Plant Germplasm Enhancement and Specialty Agriculture. The dried and pulverized rhizomes of P. sibiricum (1500 g) were extracted in ethanol (70%, 15 L) at room temperature for three times (2 days each time), and then concentrated to obtain the crude extract (CE, 150 g). The crude extract (150 g) was suspended in H 2 O (2.5 L) and successively partitioned with petroleum ether (2.5 × 3 L), ethyl acetate (2.5 × 3 L) and n-butanol (2.5 × 3 L) to yield PE (4.9 g), EA (3.4 g), n-Bu (14.7 g), and WA (86.8 g) fractions, respectively. The obtained samples were stored in airtight vials and kept at 4 • C for further tests.

Evaluation of the Antioxidant Capacity of P. sibiricum 2.3.1. DPPH Free Radical Scavenging Activity
The DPPH free radical scavenging activities of P. sibiricum samples were determined according to a previous study reported by Xu et al. [31] with slight modifications. Initially, 10 µL samples or the positive control solutions of vitamin C (46.875-3000 µM) were mixed with 190 µL of DPPH (100 µM) in a 96-well plate. The sample mixture was then incubated in the dark for 30 min at room temperature. Finally, the absorbance at 517 nm was measured with multifunctional microplate reader. Meanwhile, methanol was used as blank control in this assay, and all the samples and controls were tested in triplicate (n = 3). The DPPH free radical scavenging activity was determined using the formula: DPPH-free radical scavenging effect (%) = [(A C − A S )/A C ] × 100 (1) where A C and A S represent the absorbance value of the blank control and the tested sample or positive control, and the IC 50 value represents the 50% inhibition ratio of the DPPH free radicals.

ABTS Free Radical Scavenging Activity
The ABTS free radical scavenging activities of different P. sibiricum extracts were conducted following the method of Zhuang et al. [32] with a few minor modifications. Briefly, the working ABTS radical cation (ABTS + ) solution was prepared by mixing equal volumes of potassium persulfate (4.9 mM in H 2 O) and ABTS (7 mM in H 2 O), which are further incubated in the dark for 12 to 16 h. The ABTS + stock solution was then diluted with methanol to obtain an absorbance of 0.700 ± 0.03 at 734 nm. After that, 10 µL of appropriately diluted samples was combined with 190 µL ABTS + solution, and the absorbance of the mixture solution was recorded at 734 nm after incubation in the dark for 30 min. Vitamin C and methanol were utilized as positive and blank controls, respectively.
All samples and controls were tested in triplicate (n = 3). The results of ABTS scavenging activity were calculated as the DPPH scavenging method as described before (Equation (1)).

Ferric-Ion-Reducing Antioxidant Power (FRAP) Assay
The FRAP assay was performed using the method of Dordevic et al. [33] with modifications. The FRAP reagent (Fe 3+ -TPTZ solution) was composed of FeCl 3 ·6H 2 O (20 mM in H 2 O), TPTZ (10 mM in 40 mM HCl), and acetate buffer (300 mM, pH = 3.6) at a ratio of 1:1:10 (v/v/v), which was then stored at 37 • C before use. Next, 10 µL of properly diluted samples, 30 µL of distilled water, and 260 µL FRAP reagent were mixed in turn and incubated at 37 • C for 10 min. The absorbance was acquired at 593 nm by triplicate tests (n = 3). A calibration curve was established using FeSO 4 ·7H 2 O (62.5−2000 µM) with positive control (vitamin C). The activity of FRAP was expressed as mM Fe 2+ /g of the sample tested (mM Fe 2+ /g). The TPCs of PE, EA, n-Bu, WA, and CE extracts was determined (n = 3) by the Folin −Ciocalteu [34] method with modifications. In short, 20 µL of the diluted sample was mixed with 20 µL Folin−Ciocalteu reagent (in pure water, 25%, v/v) and the mixture was then incubated for 3 min. After that, 100 µL of sodium carbonate (Na 2 CO 3 , 200 µM) solution was added into the reaction system and incubated in the dark at room temperature for an hour. The absorbance was recorded at 760 nm. A calibration curve was established using gallic acid as the standard. The results of TPC were expressed in milligrams of GA equivalents (GAE) per gram of the sample (mg GAE/g sample).

Determination of Total Flavonoid Contents (TFC)
The TFCs of PE, EA, n-Bu, WA, and CE extracts were estimated (n = 3) using a previously described colorimetric method [35] with minor modifications. Briefly, 30 µL of appropriately diluted sample solution was mixed with 180 µL of distilled water and 20 µL of sodium nitrite solution (NaNO 2 , 5%, w/v). After incubation for 6 min, 40 µL of aluminum chloride solution (AlCl 3 , 10%, w/v) was added and incubated for a further 6 min. Then, 60 µL of sodium hydroxide solution (NaOH, 4%, w/v) was added and reacted for 15 min, and the absorbance of the mixture was determined at 510 nm. Rutin was used as the standard, and the results were represented as milligrams of rutin equivalent (RE) per gram of dry sample (mg RE/g sample).

Screening of the Potential Ligands of SOD and XOD with UF-LC-MS
Potential bioactive components with high relative binding affinity to SOD and XOD were screened by the UF-LC-MS procedures, which were carried out as previous studies [36,37]. Firstly, the EA fraction (5 mg) of P. sibiricum was dissolved with PBS (pH = 7.4, 995 µL with 5 µL DMSO) buffer solution and ultrasonicated for 30 min, which was subjected as the tested sample solution. 80 µL of sample solution (5 mg/mL) was incubated with SOD (0.2 U/µL) or XOD (0.1 U/µL) at 37 • C in the dark for 40 min. Secondly, the incubated solutions were transferred into 30 KDa cut-off ultrafiltration membranes, centrifuged at 10,000 rpm for 10 min at 25 • C, and immediately washed three times with PBS solution (pH = 7.4, 200 µL) through centrifugation to remove non-specific binding ligands. Thirdly, 200 µL acetonitrile (90%, v/v) was added and incubated for 10 min to release the compounds binding to SOD or XOD from the enzyme-ligand complexes, followed by centrifugation at 10,000 rpm for 10 min (n = 3). Finally, those ultrafiltrates were dried and reconstituted with 40 µL methanol for further analysis. In addition, the inactive enzyme group, denatured in boiling water (100 • C) for 10 min, was set up as negative control, and the treatment method is consistent with the active enzyme group. The flow rate was then 0.2 mL/min, and the injection volume was 20 µL. The following MS parameters were set: The capillary voltage (Vcap) was 3500 V and the fragmentor voltage was 175 V, respectively. The capillary temperature was 350 • C, and the drying gas flow rate was 8 L/min. The pressure in the nebulizer was set at 35 psi. The fixed collision energies were set as 10, 20, 40 and 60 V. The Mass Hunter workstation (Agilent) with a mass range of m/z 100-1500 was used to obtain profile data at a rate of one spectrum per second. The compounds were identified by comparing their retention times, parent ions, and mass fragments with references and databases.

Molecular Docking Study
The interaction mechanism between the potential ligands and corresponding target enzymes was further explored by molecular docking using AutoDock Tools 1.5.6 and the Discovery Studio 4.1 software based on the previous approach with minor modifications [37][38][39]. Firstly, the crystallized structures of the SOD (PDB 1CBJ) and XOD (PDB 1FIQ) were downloaded from RSCB Protein Data Bank (www.rcsb.org, accessed on 31 July 2022), and the 3D structures of the ligands with the lowest energy were established by ChemBio3D Ultra 12.0. Then, the 3D structures of ligands and receptors were processed by removing the water molecules, adding the hydrogen atoms, calculating the charge, and so on by AutoDock Tools. Subsequently, the docking active sites of SOD and XOD were optimized and obtained by Discovery Studio 4.1. Thereinto, the coordinates of the active sites of SOD and XOD were (X: 6.640; Y: 23.974; Z: 58.655) and (X: 28.671; Y: 29.977; Z: 101.417), respectively. Besides, the grid box was centered on the active sites of the receptors with a dimension size of 60 Å × 60 Å × 60 Å. Finally, molecular docking analysis between ligands and receptors was performed with 50 independent runs of the genetic algorithm by AutoDock Tools with other default parameters, and then the docking conformation was ranked according to the energy score.

Validation of Potential Ligands Activity by UF-LC-MS
To confirm the affinity of the potential ligands towards the target enzymes, the UF-LC-MS method was employed to estimate their relative IC 50 , and its procedure was the same as described before. Trolox, recognized for its strong antioxidant activity with an IC 50 = 2.82 mM when inhibiting SOD, was chosen as the positive control for SOD [40]. The potential ligands to SOD together with Trolox were subjected to UF-LC-MS for the determination of their corresponding BD values, therefore the relative IC 50 of potential ligands was calculated from these relative BD values ( Figure S1 and Table S1).

Statistical Analysis
All data in this work were expressed as mean ± standard deviation (SD) of triplicate measurements. The IC 50 values were calculated by plotting the percentages of scavenging activities or inhibition rates against the sample concentrations (six different concentration gradients in triplicate). Different software applications were used for statistical analysis including SPSS 25.0 (IBM Corp., New York, NY, USA), Origin 2021 (OriginLab Corporation, Northampton, MA, USA), and GraphPad Prism 5.0 (GraphPad Software Inc., San Diego, CA, USA).

Antioxidant Activities of P. sibiricum
Considering the complexity of chemical constituents and their diverse mechanisms of action, evaluating the antioxidant potential of a sample based on a single method is inappropriate. Thus, three different assays including DPPH, ABTS, and FRAP were employed in the current study to assess and compare the antioxidant activities of different P. sibiricum extracts [32]. As shown in Figure 1, the EA fraction displayed the strongest scavenging effect on DPPH and ABTS free radicals assays with the IC 50 values of 90.47 ± 3.17 µg/mL and 8.96 ± 0.21 µg/mL, respectively. This same fraction also showed the most prominent iron reducing ability among the five samples (PE, EA, n-Bu, WA, and CE), having a FRAP value of 2.59 ± 0.02 mM Fe 2+ /g. Moreover, in the three assays, the n-Bu fraction exhibited relatively higher antioxidant activities, followed by the PE fraction. However, the other extracts (WA, CE) showed modest antioxidant capacities. Based on the results from these three assays, the EA fraction exhibited the highest antioxidant activity as compared with the other four extracts. Hence, the EA fraction was selected for further research.
Considering the complexity of chemical constituents and their diverse mecha of action, evaluating the antioxidant potential of a sample based on a single meth inappropriate. Thus, three different assays including DPPH, ABTS, and FRAP wer ployed in the current study to assess and compare the antioxidant activities of differ sibiricum extracts [32]. As shown in Figure 1, the EA fraction displayed the strongest enging effect on DPPH and ABTS free radicals assays with the IC50 values of 90.47 µ g/mL and 8.96 ± 0.21 µ g/mL, respectively. This same fraction also showed the prominent iron reducing ability among the five samples (PE, EA, n-Bu, WA, and having a FRAP value of 2.59 ± 0.02 mM Fe 2+ /g. Moreover, in the three assays, the fraction exhibited relatively higher antioxidant activities, followed by the PE fra However, the other extracts (WA, CE) showed modest antioxidant capacities. Bas the results from these three assays, the EA fraction exhibited the highest antioxida tivity as compared with the other four extracts. Hence, the EA fraction was selecte further research.

Total Phenolic and Flavonoid Content
Numerous studies have supported the hypothesis that polyphenols and flavo in plants are natural antioxidants that can exhibit antioxidant effects by strongly capt free radicals like ROS. It was found that P. sibiricum possessed a good antioxidant ac and its EA fraction had the highest scavenging effect on DPPH and ABTS free rad The TPC and TFC of five P. sibiricum extracts were assessed using Folin-Ciocalteu an aluminum nitrate colorimetric methods. As shown in Table 1, TPC and TFC levels v between P. sibiricum CE and its four fractions. The TPC of the EA fraction (144.736 ± mg GAE/g) was the highest, followed by the n-Bu (19.143 ± 1.234 mg GAE/g) an (10.417 ± 0.899 mg GAE/g) fractions, and the CE has the weakest phenolic content of ± 0.014 mg GAE/g.
Similarly, the EA fraction presented the most abundant flavonoids (119.204 ± mg RT/g), followed by the n-Bu fraction (3.999 ± 0.381 mg RT/g), while the WA fra (0.078 ± 0.015 mg RT/g) had the lowest. Consequently, it is speculated that the gr antioxidant potentiality of the EA fraction could be attributed to its high levels of TP TFC.

Total Phenolic and Flavonoid Content
Numerous studies have supported the hypothesis that polyphenols and flavonoids in plants are natural antioxidants that can exhibit antioxidant effects by strongly capturing free radicals like ROS. It was found that P. sibiricum possessed a good antioxidant activity, and its EA fraction had the highest scavenging effect on DPPH and ABTS free radicals. The TPC and TFC of five P. sibiricum extracts were assessed using Folin-Ciocalteu and the aluminum nitrate colorimetric methods. As shown in Table 1, TPC and TFC levels varied between P. sibiricum CE and its four fractions. The TPC of the EA fraction (144.736 ± 6.419 mg GAE/g) was the highest, followed by the n-Bu (19.143 ± 1.234 mg GAE/g) and PE (10.417 ± 0.899 mg GAE/g) fractions, and the CE has the weakest phenolic content of 0.543 ± 0.014 mg GAE/g. Similarly, the EA fraction presented the most abundant flavonoids (119.204 ± 3.099 mg RT/g), followed by the n-Bu fraction (3.999 ± 0.381 mg RT/g), while the WA fraction (0.078 ± 0.015 mg RT/g) had the lowest. Consequently, it is speculated that the greatest antioxidant potentiality of the EA fraction could be attributed to its high levels of TPC and TFC.

Screening for SOD and XOD Ligands in P. sibiricum with UF-LC-MS
In previous studies, polysaccharides isolated from P. sibiricum were speculated as possible antioxidative and anti-hyperuricemic components on traditional animal experiments [28,29]. However, no substantial evidence for both its active components and their corresponding targets has been explored so far. To further explore the respective bioactive components in the EA fraction, the fast screening using bio-affinity ultrafiltration method with two targets (SOD and XOD) was applied to P. sibiricum [36,37]. The screened-out candidates might be regarded as potential SOD and XOD ligands contributing to the antioxidative and anti-hyperuricemic activities of the EA fraction. As shown in Figures 2 and 3, the components released from the binding complex and collected in the ultrafiltrates were analyzed by UPLC, in which 10 and 4 peaks with various binding abilities to SOD and XOD were observed. Herein, in order to further evaluate the affinity between enzymes and ligands, the binding degree (BD) was counted by the following equation: where Aa and Ab represent the peak area of the active and inactive enzyme groups, respectively. If peak areas in active group were greater than that of inactive group, the components were deduced as potential inhibitors [41].

Screening for SOD and XOD Ligands in P. sibiricum with UF-LC-MS
In previous studies, polysaccharides isolated from P. sibiricum were speculated as possible antioxidative and anti-hyperuricemic components on traditional animal experiments [28,29]. However, no substantial evidence for both its active components and their corresponding targets has been explored so far. To further explore the respective bioactive components in the EA fraction, the fast screening using bio-affinity ultrafiltration method with two targets (SOD and XOD) was applied to P. sibiricum [36,37]. The screened-out candidates might be regarded as potential SOD and XOD ligands contributing to the antioxidative and anti-hyperuricemic activities of the EA fraction. As shown in Figures 2 and  3, the components released from the binding complex and collected in the ultrafiltrates were analyzed by UPLC, in which 10 and 4 peaks with various binding abilities to SOD and XOD were observed. Herein, in order to further evaluate the affinity between enzymes and ligands, the binding degree (BD) was counted by the following equation: where Aa and Ab represent the peak area of the active and inactive enzyme groups, respectively. If peak areas in active group were greater than that of inactive group, the components were deduced as potential inhibitors [41]. The BD of potential ligands in the EA fraction targeting SOD and XOD are summarized in Table 2. For SOD, peak 10 possessed the highest binding degree (25.21%), followed by peak 6 (24.28%), peak 5 (22.93%), peak 9 (18.81%), peak 3 (17.00%), peak 7 (16.90%), peak 4 (14.86%), peak 12 (14.83%), and other lower binding ability peaks. For XOD, peak 7 exhibited a strong binding to XOD with the highest BD value of 39.72%, followed by peak 8 (9.83%), peak 3 (8.86%), and peak 4 (7.21%). In fact, the peak areas of some components in the inactive group exceeded that of the active group, with BDs < 0, were not considered as XOD ligands. When the inactive enzyme with denatured conformation was co-incubated with the sample solution, some ultrafiltration membrane would be blocked. Therefore, some non-specific small molecules were retained on the ultrafiltration membrane throughout the elution process. After dissociation, relatively larger components were obtained, resulting in greater final detected peak areas and negative BDs [41,42]. Based on this, peaks 5, 6, and 10 were considered to as potential SOD ligands, and peak 7 were presumed to be XOD potential ligands for further study.  The BD of potential ligands in the EA fraction targeting SOD and XOD are summarized in Table 2. For SOD, peak 10 possessed the highest binding degree (25.21%), followed by peak 6 (24.28%), peak 5 (22.93%), peak 9 (18.81%), peak 3 (17.00%), peak 7 (16.90%), peak 4 (14.86%), peak 12 (14.83%), and other lower binding ability peaks. For XOD, peak 7 exhibited a strong binding to XOD with the highest BD value of 39.72%, followed by peak 8 (9.83%), peak 3 (8.86%), and peak 4 (7.21%). In fact, the peak areas of some components in the inactive group exceeded that of the active group, with BDs < 0, were not considered as XOD ligands. When the inactive enzyme with denatured conformation was co-incubated with the sample solution, some ultrafiltration membrane would be blocked. Therefore, some non-specific small molecules were retained on the ultrafiltration membrane throughout the elution process. After dissociation, relatively larger components were obtained, resulting in greater final detected peak areas and negative BDs [41,42]. Based on this, peaks 5, 6, and 10 were considered to as potential SOD ligands, and peak 7 were presumed to be XOD potential ligands for further study.

Identification of SOD and XOD Ligands in P. sibiricum with UPLC-Q-TOF-MS/MS
The SOD and XOD ligands that were screened out above were identified using the high-resolution analytical instrument UPLC-Q-TOF-MS/MS in the negative-ion mode. The structures of these compounds were assigned by comparison with the MS/MS fragments and retention times reported in previous studies [43][44][45], as well as the corresponding standards compounds. The retention time (Rt), quasi-molecular ion ([M-H] − in negative ion mode), and characteristic fragments are listed in Table 2.
Firstly, the peak 10 was identified as N-trans-feruloyltyramine and taken as an example to explain fragmentation details of the homologous structures. Peak [44]. The structure of peak 10 was further determined by the retention time and secondary fragments of the corresponding standard. Peak 8 was identified as N-cis-feruloyltyramine possessing the same quasi-molecular ion and fragment ion as peak 10, which were further distinguished by the retention time and reference data about cis/trans isomers [44,45] ]was observed, which was deduced that an octopamine moiety in peaks 3 and 5 replaced the tyramine moiety in peaks 7 and 9 [44]. Hence, peaks 3 and 5 could be identified as N-cis-p-coumaroyloctopamine and N-trans-p-coumaroyloctopamine [45], which were already confirmed by the standard compound. Based on a similar principle, peaks 4 and 6, with the same [M-H]ion at m/z 328.1197, were detected as N-cis-feruloyloctopamine and N-trans-feruloyloctopamine, respectively, and the latter was confirmed by the corresponding standard. The C 2-3 bond saturation status and the type of substituted groups on the B-ring related to the fragmentation behaviors of homoisoflavonoids [46,47]. When the C 2-3 bond was saturated and the B-ring was replaced with a hydroxyl group, the predominant fragmentation was [M-H-CH 2 -B-ring+H] -, which was subsequently followed by the neu-  [48], while peak 13 was tentatively determined to be 5,7,2'-trihydroxy-8-methyl-4'-methoxyl homoisoflavanone [44]. Their structures are shown in Figure 4. The secondary mass spectrometry fragments of the three standard compounds with high binding degree to SOD enzyme as potential bioactive ligands and the retention time are shown in Figures S2-S5, which are completely consistent with our identification results. which was deduced that an octopamine moiety in peaks 3 and 5 replaced the tyramine moiety in peaks 7 and 9 [44]. Hence, peaks 3 and 5 could be identified as N-cis-p-coumaroyloctopamine and N-trans-p-coumaroyloctopamine [45], which were already confirmed by the standard compound. Based on a similar principle, peaks 4 and 6, with the same [M-H] -ion at m/z 328.1197, were detected as N-cis-feruloyloctopamine and N-trans-feruloyloctopamine, respectively, and the latter was confirmed by the corresponding standard. The C2-3 bond saturation status and the type of substituted groups on the B-ring related to the fragmentation behaviors of homoisoflavonoids [46,47]. When the C2-3 bond was saturated and the B-ring was replaced with a hydroxyl group, the predominant fragmentation was , which was unequivocally identified as 5,7,2',4'-tetrahydroxyl homoisoflavanone [48], while peak 13 was tentatively determined to be 5,7,2'-trihydroxy-8-methyl-4'-methoxyl homoisoflavanone [44]. Their structures are shown in Figure 4. The secondary mass spectrometry fragments of the three standard compounds with high binding degree to SOD enzyme as potential bioactive ligands and the retention time are shown in Figures S2-S4 and Figure S5, which are completely consistent with our identification results.

Molecular Docking
Molecular docking has been commonly performed to evaluate the interaction between the enzyme and the potential ligand and reveal the possible interaction modes by determining their docking energy, site of action, and the contributive key residues of the receptors [49,50].
In this work, according to the BD values of the potential ligands, the three compounds (peak 5, 6, 10) were docked with SOD and one compound (peak 7) to XOD. The binding energy (BE), inhibition constant (Ki), and hydrogen bonds are summarized in Table 3, and their best docking conformations within the binding sites are illustrated in Figure 5. Dithiocarbamate (DTC) and allopurinol (ALL) were set as positive controls towards SOD and XOD, correspondingly. For SOD, the three potential ligands, N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, and N-trans-feruloyltyramine, could smoothly enter the active pocket with low BE values. Thereinto, N-trans-feruloyltyramine (peak 10) exhibited the strongest affinity to SOD with the lowest binding energy (BE) of −6.83 kcal/mol and the inhibition constant (Ki) of 9.83 μM, followed by peak 6 (BE, −6.58 kcal/mol; Ki, 14.97 μM) and peak 5 (BE, −6.54 kcal/mol; Ki, 16.13 μM). Meanwhile, the BEs

Molecular Docking
Molecular docking has been commonly performed to evaluate the interaction between the enzyme and the potential ligand and reveal the possible interaction modes by determining their docking energy, site of action, and the contributive key residues of the receptors [49,50].
In this work, according to the BD values of the potential ligands, the three compounds (peak 5, 6, 10) were docked with SOD and one compound (peak 7) to XOD. The binding energy (BE), inhibition constant (Ki), and hydrogen bonds are summarized in Table 3, and their best docking conformations within the binding sites are illustrated in Figure 5. Dithiocarbamate (DTC) and allopurinol (ALL) were set as positive controls towards SOD and XOD, correspondingly. For SOD, the three potential ligands, N-trans-p-coumaroyloctopamine, N-transferuloyloctopamine, and N-trans-feruloyltyramine, could smoothly enter the active pocket with low BE values. Thereinto, N-trans-feruloyltyramine (peak 10) exhibited the strongest affinity to SOD with the lowest binding energy (BE) of −6.83 kcal/mol and the inhibition constant (Ki) of 9.83 µM, followed by peak 6 (BE, −6.58 kcal/mol; Ki, 14.97 µM) and peak 5 (BE, −6.54 kcal/mol; Ki, 16.13 µM). Meanwhile, the BEs and Kis of peak 5 (N-trans-p-coumaroyloctopamine), peak 6 (N-trans-feruloyloctopamine), and peak 10 (N-trans-feruloyltyramine) were lower than the positive control DTC (BE, −3.84 kcal/ mol; Ki, 1.52 mM), which showed that all the screened compounds have strong interactions with SOD. Particularly, the molecular docking results were in full compliance with the ultrafiltration screening results, with a relatively large BD value displaying a good affinity to SOD. Figure 5A revealed that N-trans-p-coumaroyloctopamine could form three hydrogen bonds with the active amino acid residue Val146, Val7 and Asp11. From Figure 5B, four hydrogen bonds were formed between N-trans-feruloyloctopamine and the active pocket of SOD (Val146, Val7, Asp11, Asn51). As Figure 5C showed that N-transferuloyltyramine could interact with SOD by forming four hydrogen bonds with the residues of Val146, Lys9, and Asp11. With regard to XOD, N-cis-p-coumaroyltyramine (peak 7) displayed a higher affinity to XOD with a lower BE of −6.51 kcal/mol, and the inhibition constant (Ki) of 16.83 µM, which was better than the positive drug ALL (BE, −5.08 kcal/mol; Ki, 189.72 µM). Three hydrogen bonds were formed between the amino acid residues of XOD (Gln144, Tyr1227, Ser1234) and N-cis-p-coumaroyltyramine in Figure 5D. Hydrogen bonds were the main contributor to stabilize the complex formed with these ligands-targets, moreover, there are some other forces including carbon hydrogen bond, pi-sigma, amide-pi stacked, pi-alkyl, and so on, which also contributed to the interactions between the two enzymes and those potential active ligands. In summary, the molecular docking study confirmed that N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, N-trans-feruloyltyramine, and N-cis-p-coumaroyltyramine interact well with SOD and XOD, respectively. Meanwhile, their BEs and Kis values were found better than the positive controls, which theoretically indicates that these compounds had a high inhibition effect on SOD and XOD. Hence, N-trans-p-coumaroyloctopamine, N-transferuloyloctopamine, and N-trans-feruloyltyramine were identified as potential ligands for SOD, and N-cis-p-coumaroyltyramine was determined as the potential ligand for XOD, whose activity was studied further. In summary, the molecular docking study confirmed that N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, N-trans-feruloyltyramine, and N-cis-p-coumaroyltyramine interact well with SOD and XOD, respectively. Meanwhile, their BEs and Kis values were found better than the positive controls, which theoretically indicates that these compounds had a high inhibition effect on SOD and XOD. Hence, N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, and N-trans-feruloyltyramine were identified as potential ligands for SOD, and N-cis-p-coumaroyltyramine was determined as the potential ligand for XOD, whose activity was studied further.

Antioxidant Capacity of Potential Ligands by UF-LC-MS
The relative IC50 of potential ligands with SOD was determined using UF-LC-MS, and the results are shown in Figure 6. The UF-LC-MS chromatogram is presented in Figure S1. As shown in Figure 6, compound 10 (N-trans-feruloyltyramine, relative IC50 = 1.52 mM) exhibits the highest SOD activity, followed by compound 6 (N-trans-feruloyloctopamine, relative IC50 = 1.81 mM) and compound 5 (N-trans-p-coumaroyloctopamine, relative IC50 = 2.06 mM). In addition, all three compounds showed better antioxidant activity than the positive control (Trolox). Their inhibition trend was consistent with the molecular docking and ultrafiltration screening results, with relatively large BD values displaying stronger inhibition to SOD.

Antioxidant Capacity of Potential Ligands by UF-LC-MS
The relative IC 50 of potential ligands with SOD was determined using UF-LC-MS, and the results are shown in Figure 6. The UF-LC-MS chromatogram is presented in Figure S1. As shown in Figure 6, compound 10 (N-trans-feruloyltyramine, relative IC 50 = 1.52 mM) exhibits the highest SOD activity, followed by compound 6 (N-trans-feruloyloctopamine, relative IC 50 = 1.81 mM) and compound 5 (N-trans-p-coumaroyloctopamine, relative IC 50 = 2.06 mM). In addition, all three compounds showed better antioxidant activity than the positive control (Trolox). Their inhibition trend was consistent with the molecular docking and ultrafiltration screening results, with relatively large BD values displaying stronger inhibition to SOD. In summary, the molecular docking study confirmed that N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, N-trans-feruloyltyramine, and N-cis-p-coumaroyltyramine interact well with SOD and XOD, respectively. Meanwhile, their BEs and Kis values were found better than the positive controls, which theoretically indicates that these compounds had a high inhibition effect on SOD and XOD. Hence, N-trans-p-coumaroyloctopamine, N-trans-feruloyloctopamine, and N-trans-feruloyltyramine were identified as potential ligands for SOD, and N-cis-p-coumaroyltyramine was determined as the potential ligand for XOD, whose activity was studied further.

Antioxidant Capacity of Potential Ligands by UF-LC-MS
The relative IC50 of potential ligands with SOD was determined using UF-LC-MS, and the results are shown in Figure 6. The UF-LC-MS chromatogram is presented in Figure S1. As shown in Figure 6, compound 10 (N-trans-feruloyltyramine, relative IC50 = 1.52 mM) exhibits the highest SOD activity, followed by compound 6 (N-trans-feruloyloctopamine, relative IC50 = 1.81 mM) and compound 5 (N-trans-p-coumaroyloctopamine, relative IC50 = 2.06 mM). In addition, all three compounds showed better antioxidant activity than the positive control (Trolox). Their inhibition trend was consistent with the molecular docking and ultrafiltration screening results, with relatively large BD values displaying stronger inhibition to SOD.

Conclusions
In the present study, the potential bioactive ligands with noteworthy antioxidative and anti-hyperuricemic activities together with their respective mechanisms of action were revealed for the first time from P. sibiricum. On the one hand, the antioxidant capacity of P. sibiricum was firstly evaluated with the DPPH, ABTS and FRAP assays, and the EA fraction exhibited the highest DPPH, ABTS radical scavenging and ferric reducing capacities with the IC 50 values of 90.47 ± 3.17 µg/mL, 8.96 ± 0.21 µg/mL and a FRAP value of 2.59 ± 0.02 mMFe 2+ /g. On the other hand, the EA fraction showed the highest total phenolic and flavonoid contents of 144.74 ± 6.42 mg GAE/g dw and 119.20 ± 3.10 mg RT/g dw among other extracts from P. sibiricum. Then, the bio-affinity ultrafiltration combining SOD and XOD with LC-MS/MS was used to screen out its potential bioactive antioxidative and anti-hyperuricemic compounds from P. sibiricum. As a result, 10 and 4 peaks with various binding abilities to SOD and XOD were explored based on their BD values. Peak 5 (N-trans-p-coumaroyloctopamine), peak 6 (N-trans-feruloyloctopamine), and peak 10 (N-trans-feruloyltyramine) were presumed to be potential ligands to SOD, and peak 7 (N-cis-p-coumaroyltyramine) was identified as the potential ligand to XOD. Additionally, molecular docking showed that these four compounds effectively interacted with SOD and XOD by hydrogen bonds and other varied interaction forces, and even were better than the corresponding positive controls (DTC and ALL). In summary, this study not only provides new evidence to support the antioxidative and anti-hyperuricemic pharmacological activities of P. sibiricum, but also offers the respective potential bioactive ligands from P. sibiricum for further study and application. More importantly, to the best of our knowledge, this is the first study to rapidly screen and identify the antioxidative and anti-hyperuricemic bioactive ligands from P. sibiricum targeting SOD and XOD through the UF-LC-MS method. These bioactive compounds have the potential to be applied in the development of P. sibiricum as a functional food or natural medicine against oxidization-related and hyperuricemia-related diseases in the near future.