Exploring the Antioxidant Properties of Caffeoylquinic and Feruloylquinic Acids: A Computational Study on Hydroperoxyl Radical Scavenging and Xanthine Oxidase Inhibition

Caffeoylquinic (5-CQA) and feruloylquinic (5-FQA) acids, found in coffee and other plant sources, are known to exhibit diverse biological activities, including potential antioxidant effects. However, the underlying mechanisms of these phenolic compounds remain elusive. This paper investigates the capacity and mode of action of 5-CQA and 5-FQA as natural antioxidants acting as hydroperoxyl radical scavengers and xanthine oxidase (XO) inhibitors. The hydroperoxyl radical scavenging potential was investigated using thermodynamic and kinetic calculations based on the DFT method, taking into account the influence of physiological conditions. Blind docking and molecular dynamics simulations were used to investigate the inhibition capacity toward the XO enzyme. The results showed that 5-CQA and 5-FQA exhibit potent hydroperoxyl radical scavenging capacity in both polar and lipidic physiological media, with rate constants higher than those of common antioxidants, such as Trolox and BHT. 5-CQA carrying catechol moiety was found to be more potent than 5-FQA in both physiological environments. Furthermore, both compounds show good affinity with the active site of the XO enzyme and form stable complexes. The hydrogen atom transfer (HAT) mechanism was found to be exclusive in lipid media, while both HAT and SET (single electron transfer) mechanisms are possible in water. 5-CQA and 5-FQA may, therefore, be considered potent natural antioxidants with potential health benefits.


Introduction
Coffee is one of the most popular beverages in the world, and it is known to provide several health benefits [1].It contains caffeine, which is a natural stimulant that can improve mental alertness, concentration, and mood [2].Moreover, coffee is rich in antioxidants, such as chlorogenic acid and caffeic acid, which can protect against oxidative stress and inflammation and may lower the risk of various chronic diseases, including type 2 diabetes, liver disease, and some types of cancer [3].Some studies also suggest that moderate coffee consumption may be associated with a lower risk of cognitive decline, depression, and mortality [4].
5-Caffeoylquinic (5-CQA) and 5-feruloylquinic (5-FQA) acids are two types of phenolic compounds found in abundance in coffee and other plants, which have been shown to possess various biological activities [5,6] (Figure 1).5-CQA is present in many dicotyledonous plants and is the most abundant phenolic acid in both Arabica and Robusta green coffee beans, accounting for 5-6% of dry beans [7].5-FQA, on the other hand, is relatively less abundant in green coffee beans (1%) and is found in many foods, such as black currant, carrot, highbush blueberry, and loquat [8].Both compounds are derivatives of quinic acid, which is a cyclic carboxylic acid.5-CQA, also known as chlorogenic acid, has a caffeic acid moiety attached to the quinic acid backbone.5-FQA has a similar structure but has a ferulic acid moiety attached to the quinic acid backbone.Both compounds contain multiple hydroxyl groups that make them potential antioxidants by donating electrons or hydrogen atoms to free radicals and reducing their reactivity.5-CQA has potent antioxidant properties and has been shown to protect against oxidative stress and inflammation [9,10].5-CQA may also help regulate blood sugar levels and improve insulin sensitivity [11].5-FQA has also antioxidant and anti-inflammatory properties [12].It has been shown to have neuroprotective effects and may help prevent cognitive decline [13,14].The beneficial effects of these compounds have been the subject of numerous biological studies.A study demonstrated the efficient absorption and rapid metabolism of 5-CQA, as well as some evidence of antioxidant activity in vivo [15].In addition, the antioxidant properties of six isomers of dicaffeoylquinic acid were analyzed, and the results showed that the isomers had varying antioxidant capacities, with some showing greater activity than others [16].The molecular mechanisms underlying the antioxidant activities of chlorogenic acid were examined in vitro, and it was found that the compound inhibited the activity of several key transcription factors involved in inflammation and oxidative stress and induced the expression of phase 2 detoxifying enzymes, which are important for protecting cells against oxidative damage [17].In addition, a study showed that a rosmarinic acid derivative has potent antioxidant activity that may have applications in the management of diabetes and hypertension [18].Despite the numerous investigations conducted, the precise mechanisms by which 5-CQA and 5-FQA exert their antioxidant activity, especially in their direct reactions with free radicals, remain unclear.Furthermore, there has been no comparison of the antioxidant potentials of these two compounds.Therefore, further studies are needed to fully understand the antioxidant capabilities of these biologically active molecules.
black currant, carrot, highbush blueberry, and loquat [8].Both compounds are derivatives of quinic acid, which is a cyclic carboxylic acid.5-CQA, also known as chlorogenic acid, has a caffeic acid moiety attached to the quinic acid backbone.5-FQA has a similar structure but has a ferulic acid moiety attached to the quinic acid backbone.Both compounds contain multiple hydroxyl groups that make them potential antioxidants by donating electrons or hydrogen atoms to free radicals and reducing their reactivity.5-CQA has potent antioxidant properties and has been shown to protect against oxidative stress and inflammation [9,10].5-CQA may also help regulate blood sugar levels and improve insulin sensitivity [11].5-FQA has also antioxidant and anti-inflammatory properties [12].It has been shown to have neuroprotective effects and may help prevent cognitive decline [13,14].The beneficial effects of these compounds have been the subject of numerous biological studies.A study demonstrated the efficient absorption and rapid metabolism of 5-CQA, as well as some evidence of antioxidant activity in vivo [15].In addition, the antioxidant properties of six isomers of dicaffeoylquinic acid were analyzed, and the results showed that the isomers had varying antioxidant capacities, with some showing greater activity than others [16].The molecular mechanisms underlying the antioxidant activities of chlorogenic acid were examined in vitro, and it was found that the compound inhibited the activity of several key transcription factors involved in inflammation and oxidative stress and induced the expression of phase 2 detoxifying enzymes, which are important for protecting cells against oxidative damage [17].In addition, a study showed that a rosmarinic acid derivative has potent antioxidant activity that may have applications in the management of diabetes and hypertension [18].Despite the numerous investigations conducted, the precise mechanisms by which 5-CQA and 5-FQA exert their antioxidant activity, especially in their direct reactions with free radicals, remain unclear.Furthermore, there has been no comparison of the antioxidant potentials of these two compounds.Therefore, further studies are needed to fully understand the antioxidant capabilities of these biologically active molecules.Xanthine oxidase (XO) is an enzyme found in the liver and other tissues that plays a crucial role in purine metabolism.It catalyzes the oxidation of hypoxanthine to xanthine and xanthine to uric acid, generating reactive oxygen species (ROS) as by-products.High levels of XO activity and ROS production have been associated with oxidative stress, inflammation, and various diseases, including cardiovascular disease and diabetes [19].Several components of coffee, including caffeine and chlorogenic acid, have been shown to inhibit XO activity and reduce ROS production [20,21].The XO inhibition potential of 5-CQA, in particular, has been the subject of several experimental studies; however, to our knowledge, 5-FQA has not yet been considered.For example, Arshad et al. investigated the XO inhibitory activity of eight phenolic compounds, including 5-CQA, using in vitro and computational methods, which revealed that these compounds possess varying levels of XO inhibitory activity, potentially lowering uric acid levels and combating hyperuricemia [21].On the other hand, Wan et al. explored the mechanisms by which eight Xanthine oxidase (XO) is an enzyme found in the liver and other tissues that plays a crucial role in purine metabolism.It catalyzes the oxidation of hypoxanthine to xanthine and xanthine to uric acid, generating reactive oxygen species (ROS) as by-products.High levels of XO activity and ROS production have been associated with oxidative stress, inflammation, and various diseases, including cardiovascular disease and diabetes [19].Several components of coffee, including caffeine and chlorogenic acid, have been shown to inhibit XO activity and reduce ROS production [20,21].The XO inhibition potential of 5-CQA, in particular, has been the subject of several experimental studies; however, to our knowledge, 5-FQA has not yet been considered.For example, Arshad et al. investigated the XO inhibitory activity of eight phenolic compounds, including 5-CQA, using in vitro and computational methods, which revealed that these compounds possess varying levels of XO inhibitory activity, potentially lowering uric acid levels and combating hyperuricemia [21].On the other hand, Wan et al. explored the mechanisms by which eight caffeoylquinic acids (CQAs) inhibit XO in vitro, showing that diCQAs have higher inhibitory activity than monoCQAs due to hydrophobic interaction and hydrogen bonding and that CQAs preferentially bind to the flavin adenine dinucleotide region of XO [22].In another work, Mohamed Isa et al. studied the XO inhibitory activity of the methanol extract of Plumeria rubra Linn flowers, which contain large amounts of phenol and flavonoids, showing that the extract exhibits strong XO inhibitory activity in vitro and can reduce serum uric acid levels in rats without toxicity.Despite these numerous studies on the potential inhibition of XO by 5-CQA, the exact mechanism of the interaction between this molecule and the enzyme is still unclear.Furthermore, given the structural similarity between 5-CQA and 5-FQA, it is reasonable to assume that 5-FQA also has potent XO inhibitory activity.Therefore, further research is needed to determine the XO inhibitory potential of 5-CQA and 5-FQA.
In this context and in order to explore in depth the antioxidant potential of 5-CQA and 5-FQA as well as to compare their capacities, this study focused on the evaluation of the reactivity toward the hydroperoxyl radical (HOO • ) in physiological environments and on the mechanism of XO inhibition.The HOO • radical is a moderately reactive ROS with a half-life of about a few seconds in biological cells, and it is the simplest of the peroxyl radicals that can cause significant damage to cells [23,24].Due to these characteristics, it is a suitable target for examining antiradical mechanisms [25].In this study, the main antiradical mechanisms were examined taking into account the influence of pH and physiological environments at all possible sites of 5-CQA and 5-FQA.All the work was carried out using DFT calculations at the M06-2X/6-311++G(d,p) level, which is one of the most accurate methods for describing radical reactions in solution.XO inhibition, on the other hand, was examined by blind docking and molecular dynamics simulations.These in silico approaches are useful methods for studying the interactions between a ligand and a target protein and can provide valuable information on the binding mechanism and stability of the complex.The findings of this study show that 5-CQA and 5-FQA are potent radical scavengers, with 5-CQA being more potent than 5-FQA.Both compounds were found to be more potent than the common antioxidants Trolox and BHT under physiological conditions.Furthermore, both compounds inhibit XO enzymes in the same way as quercetin.

Quantum Chemistry Calculations
Phenolic compounds exhibit their antiradical action through three distinct processes: HAT (hydrogen atom transfer), RAF (radical adduct formation), and SET (single electron transfer) [26][27][28][29][30]. HAT involves the transfer of a hydrogen atom from the antioxidant to the free radical, and it can occur in both polar and lipid environments.The Gibbs free energy of the reaction describes this one-step process.On the other hand, RAF is a one-step process where free radicals combine with antioxidants, forming the product [antioxidantfree radical] • [31].The Gibbs free energy of the reaction describes this process.Unlike HAT and RAF, SET involves an electron transfer and only occurs in polar environments [32].The phenolic OH bonds of the antioxidant can exist in both protonated and deprotonated forms in these environments, enabling electron transfer from the neutral or anionic state of the antioxidant.However, SET from an undissociated phenol is often not feasible, making it significant only for deprotonated phenols.The Gibbs free energy of electron transfer describes the dominant step in this pathway.These mechanisms are represented by the following equations: In this work, density functional theory (DFT) computations were carried out using the Gaussian 09 version E.01 software [33].The M06-2X functional, which is highly parametrized and includes empirical exchange-correlation, was used in combination with the 6-311++G(d,p) basis set, as it is known to be one of the most accurate approaches for calculating rate constants for radical reactions in solution [34][35][36].To account for solvation effects in physiological media, the SMD model was employed [37].Proton affinity (PA) values of OH groups were computed as described in the literature [38].The pKa values have been calculated according to the literature using the following equation [39]: ∆G sol RTln(10) where ∆G sol is the Gibbs free energy of the reaction in solution, R is the gas constant, and T is the temperature (298.15K).The quantum mechanics-based test for overall free radical scavenging activity (QM-ORSA) methodology was employed for the kinetic study, which was validated experimentally [25, 40,41].The rate constants (k) were calculated at 298.15 K using conventional transition state theory (TST) [42][43][44][45][46].
where σ, κ, kB, h, ∆G = are reaction symmetry numbers, tunneling corrections, Boltzmann constant, Planck constant, and Gibbs free energy of activation, respectively [47][48][49].The tunneling correction that refers to the ratio between quantum mechanics and classical mechanical rates of barrier crossing was calculated using the Eckart barrier [47].
The ratio of the probability of a specific pathway to the total probability of all possible pathways, branching ratio (Γ), was determined using the equations below: where k i and k overall are the rate constants of a specific reaction path and the sum of the rate constants of all reaction paths, respectively.The Gibbs free energy of activation (∆G = SET ) for the single electron transfer (SET) mechanism was predicted using the Marcus theory approach [50], with λ approximated as λ ≈ ∆E SET + ∆G 0 SET .
where ∆G 0 SET is the Gibbs free energy of the reaction.Finally, the apparent rate constants (k app ) were corrected using the Collins-Kimball theory to account for diffusion limitation [51].

Blind Docking Consensus Procedure
We applied the blind docking (BD) consensus approach to study the interactions between both ligands and XO.BD is a method that scans the whole protein surface finding potential hotspots for a specific ligand [52].First, we obtained the crystallized structure of the bovine xanthine oxidase described in the Protein Data Bank (PDB), whose code is 3nvy (https://www.rcsb.org/structure/3nvyaccessed on 15 September 2021).The PDB file obtained was processed using Maestro tools: Protein Preparation Wizard and System Builder.The ligand was processed with LigPrep, the tool of Maestro designed to prepare small molecules.The format of the structure obtained was mol2.All the mol2 files generated by Maestro used the force field OPLS3e [53].Protein and ligands were processed with ADT to obtain a pdbqt file.The force field of the pdbqt file was Gasteiger.BD calculations were performed using AutoDock Vina [54] and Lead Finder [55].The grid box size of all runs had measurements of 30 × 30 × 30 Angstroms as the default parameter.All the runs, including the BD consensus, were performed using an alpha version 1.0 of metascreener software (https://github.com/bio-hpc/metascreener,accessed on 9 July 2023) The scoring function of AD and LF considered the following parameters: Lennard-Jones interactions term (LJ), hydrogen bonds (H-bonds), electrostatic interactions, hydrophobic stabilization, entropic penalty due to the number of rotatable bonds, and internal energy of the ligand.Maestro software (Maestro-Desmond Interoperability Tools, Schrödinger, New York, NY, USA, 2020) was also used to calculate the interactions between ligand and protein residues.

Molecular Docking Simulations
The most relevant pose of each ligand obtained in the BD was selected to run molecular dynamic simulations (MD).Three MDs were carried out using Maestro-Desmond software version 2020-4 (Desmond Molecular Dynamics System, D. E. Shaw Research, New York, NY, 2020.Maestro-Desmond Interoperability Tools, Schrödinger, New York, NY, USA, 2020).Both complexes created were immersed in a box filled with water molecules using the simple point charge (SPC) scheme.The measurements of the box were 10 × 10 × 10 Å. Sodium ions were added to neutralize charges.Ions of chlorine and sodium were added to obtain a final NaCl concentration of 0.150 M. Periodic boundary conditions were used, and a cutoff of 9 Å was established for van der Waals interactions.The particle mesh Ewald (PWE) method with a tolerance of 10-9 was used in the electrostatic part.The energy minimization was carried out by 2000 steps using the steepest descent method with a threshold of 1.0 kcal/mol/Å.The NPT simulations were realized at 300 K with the Nose-Hoover algorithm [56], and the pressure was maintained at 1 bar with the Martyna-Tobias-Klein barostat [57].The OPLS3e force field was used in all the MDs and the duration of both MDs was 100 ns.

Acid-Base Equilibrium at Physiological pH
Understanding the behavior of phenolic groups at physiological pH is critical in studying free radical scavenging properties in water.While the pKa values of 5-CQA are available in the literature [58], the same cannot be said for 5-FQA.To address this, we followed a procedure outlined in a published paper to calculate the pKa of 5-FQA [39].We included the obtained pKa values and PA values of the OH groups of both compounds in Figure 2. Our findings indicate that for both 5-CQA and 5-FQA, the COOH group is the preferred deprotonation site, followed by the phenolic OH group at position 13, which has the lowest PA value.In water, the mono-anionic form dominates for both compounds, with a non-negligible amount of the dianionic form present.The neutral form is absent under these conditions, and therefore, we only considered the mono-anionic and dianionic forms in our study.However, in the lipid environment, only the neutral form exists, and we considered only this form in our thermodynamic and kinetic studies.

Thermodynamic Evaluation of the Antiradical Mechanisms
Figure 3 shows the ΔG values of the reaction between 5-CQA and 5-FQA with the HOO • radical in both aqueous and lipid media.In water, the HAT reaction of the nonaromatic OH groups was found to be energetically unfavorable, as indicated by positive ΔG values (31.9-45.1 kcal/mol).This suggests that these OH groups are not effective for antiradical action.Similarly, the RAF reaction at both C8 and C9 positions also showed positive ΔG values, making this mechanism less likely to occur.However, the RAF reaction at the C8 position was almost isergonic (2.1-2.7 kcal/mol), indicating that this process can take place.As a result, the kinetic study also considered the RAF mechanism at C8. On the other hand, the aromatic OH groups had negative ΔG values for both 5-CQA and 5-FQA (−5.8 to −13.0 kcal/mol), suggesting that they are the most likely sites for the HAT process.In terms of the SET mechanism, both forms showed positive ΔG values.However, the dianionic form exhibited the lowest ΔG values (2.0-4.8 kcal/mol vs. 29.3-29.8kcal/mol), indicating that this mechanism could also be possible for this form.
When considering pentyl ethanoate, where only the neutral form exists, the HAT reaction was found to be thermodynamically favorable for the phenolic OH groups (−1.4 to

Thermodynamic Evaluation of the Antiradical Mechanisms
Figure 3 shows the ∆G values of the reaction between 5-CQA and 5-FQA with the HOO • radical in both aqueous and lipid media.In water, the HAT reaction of the nonaromatic OH groups was found to be energetically unfavorable, as indicated by positive ∆G values (31.9-45.1 kcal/mol).This suggests that these OH groups are not effective for antiradical action.Similarly, the RAF reaction at both C8 and C9 positions also showed positive ∆G values, making this mechanism less likely to occur.However, the RAF reaction at the C8 position was almost isergonic (2.1-2.7 kcal/mol), indicating that this process can take place.As a result, the kinetic study also considered the RAF mechanism at C8.

Thermodynamic Evaluation of the Antiradical Mechanisms
Figure 3 shows the ΔG values of the reaction between 5-CQA and 5-FQA with the HOO • radical in both aqueous and lipid media.In water, the HAT reaction of the nonaromatic OH groups was found to be energetically unfavorable, as indicated by positive ΔG values (31.9-45.1 kcal/mol).This suggests that these OH groups are not effective for antiradical action.Similarly, the RAF reaction at both C8 and C9 positions also showed positive ΔG values, making this mechanism less likely to occur.However, the RAF reaction at the C8 position was almost isergonic (2.1-2.7 kcal/mol), indicating that this process can take place.As a result, the kinetic study also considered the RAF mechanism at C8. On the other hand, the aromatic OH groups had negative ΔG values for both 5-CQA and 5-FQA (−5.8 to −13.0 kcal/mol), suggesting that they are the most likely sites for the HAT process.In terms of the SET mechanism, both forms showed positive ΔG values.However, the dianionic form exhibited the lowest ΔG values (2.0-4.8 kcal/mol vs. 29.3-29.8kcal/mol), indicating that this mechanism could also be possible for this form.
When considering pentyl ethanoate, where only the neutral form exists, the HAT reaction was found to be thermodynamically favorable for the phenolic OH groups (−1.4 to On the other hand, the aromatic OH groups had negative ∆G values for both 5-CQA and 5-FQA (−5.8 to −13.0 kcal/mol), suggesting that they are the most likely sites for the HAT process.In terms of the SET mechanism, both forms showed positive ∆G values.However, the dianionic form exhibited the lowest ∆G values (2.0-4.8 kcal/mol vs. 29.3-29.8kcal/mol), indicating that this mechanism could also be possible for this form.
When considering pentyl ethanoate, where only the neutral form exists, the HAT reaction was found to be thermodynamically favorable for the phenolic OH groups (−1.4 to −6.9 kcal/mol).In contrast, the RAF reaction was an endogenous process for both compounds (4.4 and 4.8 kcal/mol).Despite this, the kinetic study also examined the RAF reaction at the C8 position.
Overall, the results indicate that the HAT mechanism is more likely to occur in the aromatic OH groups, while the RAF mechanism is less favorable.However, the RAF reaction at the C8 position is also possible.The SET mechanism was found to be less likely but could still be active for the dianionic form.

Reaction Kinetics under Physiological Conditions Effects of Polar Physiological Media
Using the QM-ORSA protocol [25], and based on the thermodynamic evaluation, we calculated the overall and individual rate constants for the reaction between 5-CQA and 5-FQA with the HOO • radical in water at physiological pH.The results are presented in Table 1, and the localized transition states can be seen in Figure 4. We found that both 5-CQA and 5-FQA reacted with high-rate constants of 2.68 × 10 8 M −1 s −1 and 2.28 × 10 7 M −1 s −1 , respectively.However, the reaction of 5-CQA is about 10 times faster than that of 5-FQA.Analysis of the data in Table 1 reveals that the SET mechanism is exclusive in the case of 5-FQA (Γ = 100%), while the HAT mechanism at 12OH plays a dominant role in the case of 5-CQA (Γ = 70%), despite the high-rate constant of the electron transfer process (7.95 × 10 7 M −1 s −1 ).This suggests that the deprotonated catechol moiety in water is highly reactive toward radical species.Similar findings were also observed for other catechol derivatives, such as 5-O-methylnorbergenin [59], quercetins [60], caftaric acid [61], and anthocyanidins [62].
Table 1.Calculated Gibbs free energy of activation (∆G = in kcal/mol), tunneling correction (κ), branching ratio (Γ in %), and rate constants (k in M −1 s −1 ) of the reaction of 5-CQA and 5-FQA with HOO • radical following the thermodynamically favorable antiradical mechanisms in water at physiological pH.

Effects of Lipid-Like Physiological Media
Since the SET process is not feasible in nonpolar media, it was not considered in the lipid environment, thus the study focused only on the RAF and HAT mechanisms.The results obtained from the reaction of 5-CQA and 5-FQA with the HOO • radical in pentyl ethanoate are presented in Table 2, and the localized TSs are shown in Figure 5.As we can see, RAF at the C8 position occurs with a low-rate constant that does not contribute to the overall rate constant for the two compounds studied (Γ = 0%).The HAT mechanism, on the other hand, is dominant under these conditions.The overall rate constants of 5-CQA and 5-FQA are 2.09 × 10 6 and 4.10 × 10 4 M −1 s −1 , respectively, reflecting the potent free radical scavenging activity of these compounds compared with typical antioxidants, such as Trolox (1.00 × 10 5 M −1 s −1 ) [63] and BHT (1.70 × 10 4 M −1 s −1 ) [64].These values are also higher than values for recognized antioxidants such as caftaric acid (k = 1.82 × 10 3 M −1 s −1 ) [61], carnosic acid (k = 5.70 × 10 3 M −1 s −1 ) [65], and cannabidiol (k = 2.60 × 10 3 M −1 s −1 ) [66], indicating that 5-CQA and 5-FQA are potent antioxidants in lipid media.As observed in water, the reactivity of 5-CQA is higher than that of 5-FQA, indicating that the former is a better antioxidant than the latter in both polar and lipid media.

Effects of Lipid-Like Physiological Media
Since the SET process is not feasible in nonpolar media, it was not considered in the lipid environment, thus the study focused only on the RAF and HAT mechanisms.The results obtained from the reaction of 5-CQA and 5-FQA with the HOO • radical in pentyl ethanoate are presented in Table 2, and the localized TSs are shown in Figure 5.As we can see, RAF at the C8 position occurs with a low-rate constant that does not contribute to the overall rate constant for the two compounds studied (Γ = 0%).The HAT mechanism, on the other hand, is dominant under these conditions.The overall rate constants of 5-CQA and 5-FQA are 2.09 × 10 6 and 4.10 × 10 4 M −1 s −1 , respectively, reflecting the potent free radical scavenging activity of these compounds compared with typical antioxidants, such as Trolox (1.00 × 10 5 M −1 s −1 ) [63] and BHT (1.70 × 10 4 M −1 s −1 ) [64].These values are also higher than values for recognized antioxidants such as caftaric acid (k = 1.82 × 10 3 M −1 s −1 ) [61], carnosic acid (k = 5.70 × 10 3 M −1 s −1 ) [65], and cannabidiol (k = 2.60 × 10 3 M −1 s −1 ) [66], indicating that 5-CQA and 5-FQA are potent antioxidants in lipid media.As observed in water, the reactivity of 5-CQA is higher than that of 5-FQA, indicating that the former is a better antioxidant than the latter in both polar and lipid media.In summary, the findings revealed that both 5-CQA and 5-FQA react with the H radical in physiological media at high-rate constants, with 5-CQA reacting faster tha FQA.Mechanistically, 5-CQA reacts primarily via the HAT mechanism in polar an pidic media, whereas 5-FQA appears to be more active in water via the SET mechan rather than the HAT mechanism.The overall rate constants of both compounds are c parable to or greater than those of reference antioxidants, suggesting that these c pounds are promising candidates for further study as potential antiradical agents in p iological media.

Blind Docking and Molecular Dynamics Simulations
The results of the BD calculation (Figure 6) show how the first three hotspots ha similar docking score (5-CQA: −11.21 kcal/mol to −9.94 kcal/mol, 5-FQA: −10.12 kcal to −9.47 kcal/mol).The first hotspot is the enzyme's active site, whose substrate is qu tin [67].Figure 7 shows the residues that interacted with both ligands.Regarding th teractions with residues, some of them are common with the quercetin complex: E interacted with both ligands and Glu802, Leu873, Arg880, and Phe914.Phe1009, Thr1 Val1011, and Leu1014 interacted only with 5-CQA.The average docking score of compound was as follows: −11.21 kcal/mol for the 5-CQA and -10.12 kcal/mol for th FQA.
The results of the MD show how the protein reached a stable conformation in complexes.The maximum RMSD value of the protein was 2.7 Å in both complexes garding the MD with the protein without ligand, the maximum RMSD value was 1 Moreover, the ligand stayed in a stable pose in the active site: 5-CQA had a maxim RMSD value of 4.5 Å, and 5-FQA had a value of 3.6 Å (Figure 8).
The main interactions of the 5-CQA along the trajectory were as follows: ionic i actions-Glu1261; hydrophobic interactions-Phe914, Phe1009, Val1011, and Ala1 and polar interactions-Gln767, Thr1010, Gln1040, Ala1078, Ser1080, and Ser1082 (Fi In summary, the findings revealed that both 5-CQA and 5-FQA react with the HOO • radical in physiological media at high-rate constants, with 5-CQA reacting faster than 5-FQA.Mechanistically, 5-CQA reacts primarily via the HAT mechanism in polar and lipidic media, whereas 5-FQA appears to be more active in water via the SET mechanism rather than the HAT mechanism.The overall rate constants of both compounds are comparable to or greater than those of reference antioxidants, suggesting that these compounds are promising candidates for further study as potential antiradical agents in physiological media.

Blind Docking and Molecular Dynamics Simulations
The results of the BD calculation (Figure 6) show how the first three hotspots have a similar docking score (5-CQA: −11.21 kcal/mol to −9.94 kcal/mol, 5-FQA: −10.12 kcal/mol to −9.47 kcal/mol).The first hotspot is the enzyme's active site, whose substrate is quercetin [67].Figure 7 shows the residues that interacted with both ligands.Regarding the interactions with residues, some of them are common with the quercetin complex: E1261 interacted with both ligands and Glu802, Leu873, Arg880, and Phe914.Phe1009, Thr1010, Val1011, and Leu1014 interacted only with 5-CQA.The average docking score of each compound was as follows: −11.21 kcal/mol for the 5-CQA and −10.12 kcal/mol for the 5-FQA.The results of the MD show how the protein reached a stable conformation in both complexes.The maximum RMSD value of the protein was 2.7 Å in both complexes.Regarding the MD with the protein without ligand, the maximum RMSD value was 1.8 Å.Moreover, the ligand stayed in a stable pose in the active site: 5-CQA had a maximum RMSD value of 4.5 Å, and 5-FQA had a value of 3.6 Å (Figure 8).The main interactions of the 5-CQA along the trajectory were as follows: ionic interactions-Glu1261; hydrophobic interactions-Phe914, Phe1009, Val1011, and Ala1078; and polar interactions-Gln767, Thr1010, Gln1040, Ala1078, Ser1080, and Ser1082 (Figure 9).The main interactions of the 5-FQA along the trajectory were as follows: ionic interactions-Arg912, Lys1045, Asp1191, and Glu1261; polar interactions-Thr1077, Ser1080, and Ser1082; and water interactions-Gly913 (Figure 10).These results show the high stability of 5-CQA and 5-FQA and Xanthine Oxidase complexes, and we could infer that both ligands could be inhibitors of the enzyme.

Figure 3 .
Figure 3. ΔG values calculated in kcal/mol for the reactions of 5-CQA and 5-FQA and their dissociated forms with the HOO • radical according to possible antiradical mechanisms under physiological conditions (water and pentyl ethanoate).HAT = hydrogen atom transfer, SET = single electron transfer, and RAF = radical adduct formation.

Figure 3 .
Figure 3. ΔG values calculated in kcal/mol for the reactions of 5-CQA and 5-FQA and their dissociated forms with the HOO • radical according to possible antiradical mechanisms under physiological conditions (water and pentyl ethanoate).HAT = hydrogen atom transfer, SET = single electron transfer, and RAF = radical adduct formation.

Figure 3 .
Figure 3. ∆G values calculated in kcal/mol for the reactions of 5-CQA and 5-FQA and their dissociated forms with the HOO • radical according to possible antiradical mechanisms under physiological conditions (water and pentyl ethanoate).HAT = hydrogen atom transfer, SET = single electron transfer, and RAF = radical adduct formation.

Figure 4 .
Figure 4. Optimized TS structures of the reaction of 5-CQA and 5-FQA with HOO • radical following HAT (hydrogen atom transfer) and RAF (radical adduct formation) mechanisms in water (W).

Figure 4 .
Figure 4. Optimized TS structures of the reaction of 5-CQA and 5-FQA with HOO • radical following HAT (hydrogen atom transfer) and RAF (radical adduct formation) mechanisms in water (W).

Figure 6 .
Figure 6.(A) Blinding the docking complex with all hotspots from the BD consensus.In (B), only the best three hotspots are represented.Illustrative images were prepared using the Pymol software version 2.5.0.

Figure 6 .
Figure 6.(A) Blinding the docking complex with all hotspots from the BD consensus.In (B), only the best three hotspots are represented.Illustrative images were prepared using the Pymol software version 2.5.0.

Figure 7 .
Figure 7. (A) Blinding the docking complex with hotspot number 1 resulted from the BD consensus.Illustrative images were prepared using the Pymol software version 2.5.0.(B) Two-dimensional protein-ligand interaction diagram generated between xanthine oxidase and 5-CQA (left) and 5-FQA (right), using the Ligand Interaction script in Maestro (Schrödinger, Inc., www.schrodinger.comaccessed on 20 June 2022).

Figure 7 .
Figure 7. (A) Blinding the docking complex with hotspot number 1 resulted from the BD consensus.Illustrative images were prepared using the Pymol software version 2.5.0.(B) Two-dimensional protein-ligand interaction diagram generated between xanthine oxidase and 5-CQA (left) and 5-FQA (right), using the Ligand Interaction script in Maestro (Schrödinger, Inc., www.schrodinger.comaccessed on 20 June 2022).

Figure 8 .
Figure 8.The root-mean-square deviation (RMSD) obtained for xanthine oxidase (blue (A) and orange (B)) bound to 5-FQA (red) or 5-CQA (purple) plotted with respect to the initial pose.The RMSD of the xanthine oxidase WT is represented in both plots (A,B) in green.

Figure 8 .
Figure 8.The root-mean-square deviation (RMSD) obtained for xanthine oxidase (blue (A) and orange (B)) bound to 5-FQA (red) or 5-CQA (purple) plotted with respect to the initial pose.The RMSD of the xanthine oxidase WT is represented in both plots (A,B) in green.