Next Article in Journal
The Endothelial Centrosome: Specific Features and Functional Significance for Endothelial Cell Activity and Barrier Maintenance
Next Article in Special Issue
Polyphenols and Their Impact on Human Health
Previous Article in Journal
Recent Advances in Novel Compositions for Electrochemical Applications
Previous Article in Special Issue
Inhibition of Cancer Development by Natural Plant Polyphenols: Molecular Mechanisms
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 20
10.3390/ijms242015391
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inhibitory Effect of Quercetin on Oxidative Endogen Enzymes: A Focus on Putative Binding Modes

by
Stefania Olla
1,*,
Chiara Siguri
1,
Antonella Fais
2,
Benedetta Era
2,
Massimo Claudio Fantini
3 and
Amalia Di Petrillo
3,*
1
Istituto di Ricerca Genetica e Biomedica, Consiglio Nazionale delle Ricerche, Cittadella Universitaria di Monserrato, 09042 Monserrato, Italy
2
Department of Life and Environmental Sciences, University of Cagliari, 09042 Monserrato, Italy
3
Department of Medical Science and Public Health, University of Cagliari, 09042 Monserrato, Italy
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(20), 15391; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms242015391
Submission received: 26 September 2023 / Revised: 13 October 2023 / Accepted: 16 October 2023 / Published: 20 October 2023
(This article belongs to the Special Issue Polyphenols and Their Impact on Human Health: 2nd Edition)
Browse Figures
Review Reports Versions Notes

Abstract

:
Oxidative stress is defined as an imbalance between the production of free radicals and reactive oxygen species (ROS) and the ability of the body to neutralize them by anti-oxidant defense systems. Cells can produce ROS during physiological processes, but excessive ROS can lead to non-specific and irreversible damage to biological molecules, such as DNA, lipids, and proteins. Mitochondria mainly produce endogenous ROS during both physiological and pathological conditions. Enzymes like nicotinamide adenine dinucleotide phosphate oxidase (NOX), xanthine oxidase (XO), lipoxygenase (LOX), myeloperoxidase (MPO), and monoamine oxidase (MAO) contribute to this process. The body has enzymatic and non-enzymatic defense systems to neutralize ROS. The intake of bioactive phenols, like quercetin (Que), can protect against pro-oxidative damage by quenching ROS through a non-enzymatic system. In this study, we evaluate the ability of Que to target endogenous oxidant enzymes involved in ROS production and explore the mechanisms of action underlying its anti-oxidant properties. Que can act as a free radical scavenger by donating electrons through the negative charges in its phenolic and ketone groups. Additionally, it can effectively inhibit the activity of several endogenous oxidative enzymes by binding them with high affinity and specificity. Que had the best molecular docking results with XO, followed by MAO-A, 5-LOX, NOX, and MPO. Que’s binding to these enzymes was confirmed by subsequent molecular dynamics, revealing different stability phases depending on the enzyme bound. The 500 ns simulation showed a net evolution of binding for NOX and MPO. These findings suggest that Que has potential as a natural therapy for diseases related to oxidative stress.

1. Introduction

Quercetin (Que) is a 3,3′,4′,5,7-pentahydroxyflavone with the chemical formula C15H10O7. The presence of five hydroxyl groups in the structure and the pyrocatechol-kind of benzene ring make Que a strong anti-oxidant and a good scavenger of free radicals [1]. The structure of Que contains a keto carbonyl group, and the oxygen atom present on the first carbon, being basic, can form salts with acids. Furthermore, a dihydroxy group between the A ring, o-dihydroxy group B, C ring C2, C3 double bond, and 4-carbonyl are the active groups in Que (Figure S1A). Que’s biological activities have been largely attributed to these active phenolic hydroxyl groups and double bonds [2].
Que is highly present in daily human foods, such as onions, apples, red wine, and tea. It has various biological effects, including anti-oxidant, anti-cancer, anti-inflammatory, anti-viral, and anti-aging effects. Additionally, it has anti-aggregant and vasodilator effects, which help protect against cardiovascular diseases [3,4]. In humans consuming flavonoid-rich foods, the concentration range of Que plasma levels is 0.3–7.6 μM, mainly in the form of glucuronidated and sulfated metabolites [5].
Free radicals and reactive oxygen species (ROS) are produced in cells due to common metabolic processes or external sources [6]. Excessive oxidative stress disrupts the balance between the oxidation and anti-oxidation systems. This causes non-specific and irreversible damage to biological molecules, such as lipids, proteins, and DNA, leading to a loss of function [7].
In many chronic diseases, such as cancer, neurodegenerative diseases, and metabolic diseases, oxidative stress is often the primary trigger [8].
Possible sources of endogenous ROS include the nicotinamide adenine dinucleotide phosphate oxidase (NOX), xanthine oxidase (XO), lipoxygenase (LOX), myeloperoxidase (MPO), and monoamine-oxidase (MAO) enzymes [9]. Increased activation of these enzymes has been linked to oxidative stress as well as the onset and advancement of inflammatory diseases.
Notably, Que has exhibited promising protective effects by attenuating the expression and enzymatic activity of these endogenous oxidative enzymes [10,11,12,13,14,15].
NOX (EC 1.6.3.1) is a transmembrane glycoprotein that catalyzes the production of superoxide (O2−) from oxygen and NADPH. O2− reacts quickly to produce a burst of additional oxidants, including hydrogen peroxide (H2O2), which is typically further converted by MPO into hypochlorous acid (HOCl) and hydroxyl radical (•OH). For this reason, NOX is usually considered a major source of ROS and oxidative stress in eukaryotic cells [16].
XO (EC 1.17.3.2) is an iron-molybdenum flavoprotein, predominantly found in the cytoplasm of mammalian tissues, that plays a crucial role in purine catabolism [17]. It catalyzes the oxidation of hypoxanthine to xanthine and xanthine to uric acid (UA) with the production of a superoxide anion. The accumulation of UA has been shown to initiate the inflammatory process through the NLRP3 inflammasome and ROS production, contributing to inflammation-related tissue damage [18].
MPO (EC 1.11.2.2) is an enzyme found in the azurophilic/primary granules of neutrophils and, to a lesser extent, in monocytes. MPO catalyzes the oxidation of chloride and other halide ions in H2O2 to generate HOCl and other highly reactive products that mediate efficient anti-microbial action [19].
LOXs (EC 1.13.11.12) are a large monomeric protein family with non-heme, non-sulfur, iron cofactor-containing dioxygenases that catalyze the oxidation of polyunsaturated fatty acids (PUFA) in lipids containing a cis,cis-1,4-pentadiene into cell signaling agents [20]. The main products of LOXs are leukotrienes and lipoxins, which play important roles in several inflammation-related diseases, such as arthritis, asthma, cancer, and allergies.
MAO (EC 1.4.3.4) is a riboflavin protein distributed on the outer mitochondrial membrane that catalyzes the oxidative deamination of biogenic and xenobiotic amines, producing the corresponding aldehydes, hydrogen peroxide, and ammonia [21]. Two isoforms of MAO, MAO-A and MAO-B, differ in their substrate and tissue distributions. MAO-A preferentially deaminates serotonin and norepinephrine, while MAO-B metabolizes dopamine [22]. The generation of H2O2 via MAOs is reported to be a cytotoxic factor involved in oxidative stress and neurodegenerative disorders, such as Parkinson’s disease [23].
Due to their hydrophobic nature, Que and other phenolic compounds like curcumin exhibit a strong affinity for mitochondria and can be traced through autofluorescence. This unique property of Que may contribute to its specific effect on mitochondrial dysfunction [24].
It is also important to address the limitations of Que’s bioavailability, especially in its aglycone form, which has low solubility in water and gastrointestinal fluids, leading to precipitation and reduced absorption [25]. Various que derivatives have differing bioavailabilities, with que glycosides found in onions having the highest absorption rates. These challenges can be overcome through strategies like encapsulation in colloidal particles and crystal engineering, as demonstrated in recent articles [26,27].
Our comprehensive investigation delves into Que’s potential to target oxidant enzymes implicated in ROS production. To explore these aspects thoroughly, we conducted docking analyses and molecular dynamics (MD) simulations to validate the obtained data and calculate the MMGBSA binding free energy before and during the entire dynamics. The use of docking allowed us to explore binding interactions between Que and targets, providing insight into initial structural arrangements and potential binding modes. On the other hand, the MD simulations performed provided a dynamic view of the ligand-protein complex, allowing us to assess the stability of the binding and observe any conformational changes over time. By calculating the free energy of binding with MMGBSA, we quantified the thermodynamic aspects of the ligand-protein interaction.

2. Results

2.1. Quercetin Anti-Oxidant Activity: In Silico Approach

Density functional theoretical (DFT) approaches [28] at the B3LYP/6-311++G(d, p) level were applied for the optimization of Que, and revealed the presence of C2′H/C3OH, C3OH C4O, and C4O/C5O H bonds (Figure S1A). Moreover, optimized geometry (Table S1), electrostatic potential (ESP), average local ionization energy (ALIE), electron density (ED), highest occupied molecular orbital (HOMO), lowest un-occupied molecular orbital (LUMO), and natural bond orbitals (NBO) were calculated.
ESP analysis is a robust tool for revealing the chemical and reactivity attributes of molecular systems [29]. In our computations, the ESP values for Que range from −2.65 eV to 3.18 eV (Table 1).
These values provide insights into the distribution of electron-rich and electron-poor regions within the molecule, indicating potential sites of high and low electron density (Figure S1D).
ALIE, a pivotal property measuring the energy required to extract an electron from a specific location within the molecular system, reveals a range for Que from 8.84 eV to 16.81 eV (Table 1). Lower ALIE values signify regions where electrons are less tightly bound, suggesting sites with higher reactivity (Figure S1B).
The ED analysis involves determining the probability of locating an electron in the vicinity of a molecule [30]. This property contributes to understanding the electronic structure and behavior of Que. Our computations (Figure S1C) provide valuable insights into the spatial distribution of electron density around the molecule.
The HOMO and LUMO orbitals were visually presented through color-coded surfaces (Figure 1). The HOMO and LUMO energy levels were determined to be −5.68 eV and −1.94 eV, respectively, and their energy difference, known as the energy gap (ΔE), was found to be 3.74 eV (Table 1). This gap indicates Que’s capacity for engaging in free radical optimization mechanisms, with smaller ΔE values suggesting a greater propensity for effective participation in radical scavenging reactions [31]. In accordance with Koopman’s theorem (1934), the LUMO energy level corresponds to the electron affinity (EA = −ELUMO), indicating its ability to accept electrons during molecular interactions, and the HOMO energy level is associated with the ionization potential (IP = −EHOMO), signifying its role as an electron donor. For Que, the IP was calculated as 5.68 eV and the EA as 1.94 eV (Table 1).
Natural bond orbital (NBO) analysis is a robust method for studying intra-molecular and inter-molecular bonding, charge transfer, and conjugative interactions within molecular systems. Higher stabilization energy values indicate stronger electron donor-acceptor interactions, especially enhancing conjugative interactions [32]. The results of NBO studies, including natural population analysis, can be found in Table S2, while the second-order perturbation analysis of the Fock matrix is presented in Table S3.

2.2. Molecular Docking

Molecular docking was performed for Que and all selected endogenous oxidative enzymes. Table 2 shows the information regarding the five X-rays used. For all docking results obtained, the binding free energy of the receptor-ligand complex was predicted using MMGBSA calculations [33].

2.2.1. NAD(P)H Oxidase

The results of docking showed that Que establishes H-bonds with Asp 179 and Cys 242 and interacts with the different residues of NOX shown in Figure 2 with a Glide Gscore of −5.35 kcal/mol (Table 3). Subsequently, in this docking pose, the ligand-receptor binding energy was calculated, obtaining a MMGBSA dG bind of −4.82 kcal/mol (Table 3).

2.2.2. Xanthine Oxidase

The results of docking on XO showed a Glide Gscore of −8.87 kcal/mol and that Que establishes 4 four H-bonds with Glu 802, Ser 876, Arg 880, and Thr 1010, respectively, and two π-π interactions with Phe 914 and Phe 1009, as well as interactions with the residues shown in Figure 3. The best pose resulting from the docking simulation was then superimposed on the native Que, showing an RMSD of 0.111 Å (Figure S2) and a similar orientation of the binding mode. As before, the docking results were subjected to MMGBSA calculations, obtaining a free energy of −22.18 kcal/mol (Table 3).

2.2.3. Myeloperoxidase

The results of docking on MPO revealed a Glide Gscore of −4.91 kcal/mol and that Que is oriented within the active site of MPO with the A-ring and C-ring positioned over the iron heme, a crucial component of MPO catalytic activity. This orientation suggests a potential mechanism of inhibition of MPO activity by Que. Que interacts with Phe 147 and Phe 407 via π-π interactions and establishes two H-bonds with Glu 116 (Figure 4). Subsequently, the binding between Que and MPO obtained by docking was subjected to the MMGBSA calculation, obtaining a dG binding of −28.81 kcal/mol (Table 3).

2.2.4. Lipoxygenase

The results indicated that Que establishes several crucial interactions with amino acid residues within the 5-LOX binding site and has a Glide Gscore of −5.94 kcal/mol. Specifically, Que formed five H bonds with His 130, Glu 134, Thr 137, Arg 138, and Asp 166, and engaged in π-cation interactions with Arg 101 (Figure 5). Moreover, the binding pose was subjected to the MMGBSA calculation, obtaining a dG binding value of −37.42 kcal/mol (Table 3).

2.2.5. Monoamine Oxidase

The results indicated a Glide Gscore of −8.14 kcal/mol and that Que exhibited interactions with specific amino acid residues within the MAO-A active site. Notably, Que formed a H-bond with Ala 111 and engaged in a π-π interaction with Phe 208, as illustrated in Figure 6. Then, the obtained docking pose was subjected to the MMGBSA calculation with a positive free energy of 5.72 kcal/mol (Table 3).

2.3. Molecular Dynamics

MD studies were performed to simulate the dynamic trend and stability of the molecular Que-enzyme complexes obtained from the docking studies. In detail, the Desmond software (Schrödinger Release 2022-1) [39] was used to run all simulations for 500 ns by setting the Que-enzyme complex with the best Glide Gscore as the initial structure of each dynamic.

2.3.1. NADP(H) Oxidase

To study the stability of MD between Que and NOX, we calculated the trajectory RMSD, which is plotted in Figure 7A. All the frame trajectories are aligned on the reference frame (frame 0), and the average RMSD for the protein structure is 1.515 ± 0.114 Å, while Que averages 8.854 ± 1.427 Å. Throughout the 500 ns simulation, the protein remains relatively stable without significant variations, while Que stabilizes after approximately 82 ns.
To obtain detailed information on the binding energy and its variations during the simulation, MMGBSA calculations were performed on every 10th frame (see Section 4), with an average of −56.421 ± 8.635 kcal/mol. The frame with the lowest value is 4190 (419 ns), with −73.563 kcal/mol, while the frame 80 (80 ns) has the highest value at −22.723 kcal/mol.
During the MD simulation, Que forms different interactions with the binding site of the enzyme, some of which are stable, while most are intermittent. The most stable ones include the H-bonds with Ser 115 (observed 74% of the simulation time), Cys 133 (54%), and Leu 241 (43%) (Figure 7B and Table S4). The only interaction consistently observed in both docking and MD is an H-bond with Cys 242 (10%).
To gain deeper insights into the conformational changes of Que during the MD simulation, we performed trajectory clustering. The results of the five most populated clusters were submitted to MMGBSA calculations, and the results are shown in Table 4, while Figure S3A illustrates the positions of the clusters within the binding site, comparing them with the best Glide Gscore docking pose. Notably, the clusters exhibit different positions compared to the poses extracted during docking.

2.3.2. Xanthine Oxidase

The MD of XO and Que indicates that the protein increases steadily, reaching a peak after 150 ns of simulation before stabilizing (Figure 8A). In contrast, Que follows a stable trend until 330 ns, after which its RMSD begins to increase, forming a series of peaks and plateaus. The average RMSD for XO is 2.271 ± 0.261 Å, while Que exhibits an average RMSD of 7.627 ± 5.620 Å.
The average MMGBSA dG binding is −31.369 ± 6.173 kcal/mol. Frame 1 (0.1 ns) has the lowest value at −55.603 kcal/mol, while frame 3930 (393 ns) records the highest value at −10.938 kcal/mol.
The analysis of the interactions disclosed that Que forms stable interactions with Ala 1079 (H-bond, 40%), Glu 879 (water bridge, 38%), and Pro 1076 (water bridge, 33%) of NOX (Figure 8A) and maintains a H-bond with Glu 802 (21%), and two π-π stackings with Phe 1009 (16% and 10%), also present in the docking (Figure 8B).
We also clustered the trajectory script, and the results are provided in Table 5. The first cluster has an MMGBSA energy of −32.656 kcal/mol, while Figure S3B depicts the alignment of the five clusters with the docking pose. In this case, the docking position and the clusters are positioned in the same portion, but the five clusters are like each other in the binding while the docking position differs (Figure S3B).

2.3.3. Myeloperoxidase

The average RMSDs for the protein structures of MPO and Que are 1.352 ± 0.098 Å and 9.462 ± 3.813 Å, respectively. While the RMSD of MPO is generally stable throughout the simulation, that of Que tends to increase until it stabilizes after 214 ns (Figure 9A).
Subsequently, MMGBSA was calculated in frames with an average of −30.663 ± 7.097 kcal/mol, and the lowest value is observed at frame 2710 (271 ns) with −62.350 kcal/mol, while the highest is recorded at frame 160 (16 ns) with −14.895 kcal/mol. In this case, the results align with the MMGBSA dG binding values obtained from the docking pose.
Analysis of the binding interactions between Que and MPO showed that the only interaction between Que and the active site residues for at least 10% of the simulation time is an H-bond with Glu 116 (19%), which can also be observed in the docking simulation (Figure 9B).
Finally, all frames of the simulation were clustered, and the top five clusters were reported in Table 6 with the MMGBSA energy values that were calculated. An analysis of the binding modes of the 5 clusters shows that Que occupies MPO binding sites differently (Figure S3C). Clusters 2, 4, and 5 are located within the binding site, effectively replacing the space previously occupied by the docking pose, while clusters 1 and 3 gradually move out of the binding site.

2.3.4. Lipoxygenase

The RMSD analysis was performed, obtaining an average RMSD of 1.557 ± 0.103 Å for 5-LOX, while that of Que was 2.881 ± 0.441 Å. The protein’s RMSD trend is stable, while that of Que is quite stable throughout the simulation but shows a series of peaks at intervals of 29–35 ns, 175–210 ns, and 377–380 ns (Figure 10A).
Then, the MMGBSA calculations performed on the trajectory yielded an average MMGBSA dG binding of −30.663 ± 7.097 kcal/mol, ranging from −51.253 kcal/mol at frame 20 (2 ns) to −21.253 kcal/mol at frame 1910 (191 ns). These results do not deviate significantly from those obtained from the MMGBSA calculations performed on the docking poses.
The analysis of interactions between Que and the 5-LOX showed that the most stable interactions are two π-cation interactions with Arg 101 (69% and 35%), also found in docking studies, and a water bridge with Tyr 142 (41%). Conversely, while the interactions with residues Glu 134, Arg 138, and Asp 166 establish a series of H-bonds during the docking simulation, in the MD, all these interactions transform into water bridges (32%, 18%, and 16%, respectively) (Figure 10B).
In addition, a cluster analysis on all frames of the simulation pointed out that the five most populated clusters are all located in a region of the binding site also occupied by the reference docking pose, but the binding mode differs slightly (Figure S3D). Table 7 shows the results of the MMGBSA energy calculated for the five most populated clusters.

2.3.5. Monoamine Oxidase

The RMSD analysis shows high RMSD values for the protein compared to previous simulations. The RMSD of MAO-A, with an average of 5.430 ± 1.001 Å, displays instability with a high peak of approximately 180 ns, while in the final part, approximately 400 ns, it finds stability. The RMSD of Que is unstable in the early part of the simulation and finds its stability after approximately 179 ns, and its average is 3.039 ± 0.562 Å (Figure 11A).
The MMGBSA calculations executed on the trajectory displayed an average MMGBSA of −38.836 ± 4.158 kcal/mol, with a minimum of (−53.135 kcal/mol) and a maximum of (−26.992 kcal/mol). Such results exhibit a significant discrepancy compared to the MMGBSA dG binding value of the docking pose, in which it even assumes a positive value (see Table 3).
The analysis of Que binding with MAO-A showed that the most frequent bonds are two H-bonds with Thr 336 (96%) and Cys 323 (79%), a water bridge with Asn 181 (74%), and a π-π stacking interaction with Phe 208 (68%) also present in the docking (Figure 11B and Table S4).
In addition, cluster analysis showed that the five most populated clusters are all located in a region of the binding site also occupied by the reference docking pose (Figure S3E), and the MMGBSA energies are shown in Table 8.

3. Discussion

Que is a potent anti-oxidant in the flavonoid family due to the presence of a phenolic hydroxyl group and double bonds [40].
In recent years, several in silico studies have been conducted to elucidate the molecular basis of its anti-oxidant activity.
Although the different DFT studies cannot be directly compared due to the use of different software and levels of theory, the HOMO and LUMO orbitals carried out by us and Zheng et al. [34] presented similar results, with the HOMO orbitals predominantly located in the B and C ring and the LUMO orbitals distributed throughout the molecule. Reported values are good indicators of chemical reactivity, especially for aromatic systems [35,36], and the results suggest that Que is a highly effective free radical scavenger.
Que anti-oxidant effects may, in part, arise from its ability to inhibit these enzymes. In fact, it has been reported that Que’s ability to attenuate the expression and enzymatic activity of various endogenous oxidative enzymes, including NOX, XO, LOX, MPO, and MAO, both in vivo and in vitro [9,10,11]. To validate these observations, we conducted docking analyses and molecular dynamics simulations.
Regarding the NOX enzyme, the results of our molecular docking study raise interesting questions regarding the inhibitory potential of Que on NOX activity, especially for the MMGBSA result (−4.82 kcal/mol). These results warrant further investigation into the dynamics and binding kinetics of Que within the NOX binding pocket. The MD confirms that the binding between Que and NOX finds excellent stability after approximately 100 ns and remains stable throughout the remainder of the simulation. The new MMGBSA calculations also improved by obtaining an average of −56.421 kcal/mol. The difference is probably due to the change in binding compared to docking and the better fit of Que on NOX observed during the simulation, which is evident from the greater formation of strong interactions for most of the simulation (Table S4).
The previous studies, including X-ray structural analysis [37] and various docking simulations [38,41], consistently highlighted Que’s potential as an XO inhibitor. Building upon these promising insights, our own docking studies reaffirmed the previous findings supporting Que’s potential as an XO inhibitor. Moreover, MMGBSA energy in the docking pose yielded a substantial binding energy of −22.18 kcal/mol, indicating a discrete binding affinity confirmed by MD. The five most populated clusters also maintained binding in the same docking zone, and Que maintained excellent stability and some interactions for most of the simulation (Table S4), losing it slightly but then regaining it.
Regarding MPO, our results reaffirm the potential inhibitory role of Que in modulating MPO activity, as previously studied by Pereira and colleagues [42]. Que effectively interacts with MPO, primarily through H bonds with Glu 116. The Glide Gscore energy values were not optimal, while the MMGBSA of the docking pose yielded a moderate energy result. The MD study revealed changes in Que binding within the binding site, as evidenced by the ligand’s RMSD exhibiting a wide range, especially at the beginning of the dynamics, where Que undergoes modifications in its binding. This is further confirmed by the range of MMGBSA values obtained, ranging from −62.350 to −14.895 kcal/mol, and by the different binding poses found in the five most populated clusters.
Regarding lipoxygenase (5-LOX) enzymes, our results are like the docking results presented by Vyshnevska et al. [43]. The MD studies confirm the binding stability between Que and this enzyme with low RMSDs for enzyme and ligand, and the MMGBSA studies performed throughout the dynamics are comparable to docking. Furthermore, the five most populated clusters have bindings comparable to docking. These results contribute to the growing body of evidence supporting Que’s role as a potential natural agent for controlling inflammation and related disorders through 5-LOX inhibition.
Comparing our results to those reported by Larit et al. [44], we observed some differences in Que’s interactions with MAO-A. Larit et al. reported that Que interacted with multiple residues, including Ala 111, Ile 180, Asn 181, Phe 208, Gln 215, Thr 336, and Tyr 444. In contrast, our docking study on the MAO-A active site showed that Que primarily formed an H-bond with Ala 111 and a π-π interaction with Phe 208. An anomalous result was found for the energy of MMGBSA for docking. MD studies overturned this positive value and obtained an average MMGBSA energy of −38.836 kcal/mol. Dynamics analysis showed that Que stabilizes after approximately 180 ns of simulation and then maintains its stability. The interactions with Ala 111 and Phe 208 are also preserved during MD, additionally forming those with Thr 336, Asn181, Cys 323, and Tyr 444. Overall, these results suggest that Que has the potential to inhibit MAO-A, albeit with variations in the specific residues involved.

4. Materials and Methods

4.1. DFT Approach

Jaguar (Schrödinger Release 2022-1: Jaguar, Schrödinger, LLC, New York, NY, USA, 2022) [45] was utilized to conduct the DFT calculations [28] using the B3LYP (Becke’s three-parameter exchange potential and Lee-Yang-Parr correlation function) method [46] in conjunction with the 6-311++G** (d, p) basis set, which successfully recommended itself for similar tasks and objects [47,48]. The calculations included the determination of the following properties: HOMO, LUMO, their energy gap (ΔE), ESP, ALIE, NBO, and ED.

4.2. Protein Preparation

On the Protein Data Bank (RCSB PDB) [49], we selected 5 X-ray structures, which were prepared using the Protein Preparation Wizard tool (Schrödinger Release 2022-1: Protein Preparation Wizard) [50]. The preparation included the removal of crystallized molecules, filling in missing side chains and loops, determining bond orders using the Chemical Component Dictionary (CCD) database [51], and adding hydrogen atoms. Epik (Epik, Schrödinger, LLC, New York, NY, USA, 2021) was employed to generate the states of heteroatoms at pH 7.4 ± 2.0. Protein optimization utilized PROPKA [52], and heavy atoms were converged to achieve an RMSD of 0.30 Å.

4.3. Ligand Preparation

Conformational analysis of Que was conducted using the Quantomechanic (QM) Conformer and Tautomer Prediction tool within Jaguar. The most favorable conformation was subsequently prepared using LigPrep (Schrödinger Release 2022-1: LigPrep, Schrödinger, LLC, New York, NY, USA, 2022) with the OPLS4 force field. and pH range of 7.0 ± 2.0 while preserving the specified chirality.

4.4. Molecular Docking

Among various docking programs [53,54], Glide (Schrödinger Release 2022-1: Glide, Schrödinger, LLC, New York, NY, USA, 2022) [55,56] was selected for the docking simulations, and three poses per ligand were retained using a threshold of 0.50 kcal/mol, and the OPLS_2005 [57] force field was employed.
To create grids, the Receptor Grid Generator was utilized. The center of the grid, with dimensions of 46 × 46 × 46 Å, was determined based on the centroid of the residues within the binding site. 5 grids were generated, one for each enzyme.
The docking results were subjected to binding energy calculation with MMGBSA (Molecular Mechanics with Generalized Born and Surface Area solvation) [58] of Prime (Schrödinger Release 2022-1: Prime, Schrödinger, LLC, New York, NY, USA, 2022) [59,60] using VSGB [61] as the solvation model, and OPLS4 [62] as the force field.

4.5. Molecular Dynamics

MD simulations were conducted using Desmond (Schrödinger Release 2022-1: Desmond Molecular Dynamics System, New York, NY, USA, 2022 [39]. For each system, the pose with the best docking score was chosen, and a system setup was performed. The TIP3P (Transferable Intermolecular Potential with 3 Points) [63] solvent model was used, and the ligand-receptor complex was placed in an orthorhombic water box that extended 10.0 Å, and the box volumes were minimized and neutralized by adding ions (Na+ or Cl−). The OPLS4 force field was selected.
MD simulations were carried out for 500 ns, with trajectory recorded at 100 ps, in the NPT ensemble, maintaining a constant temperature (300.0 K) and pressure (1.01325 bar) using the Nosé-Hoover thermostat [64] and the Martyna-Tobias-Klein barostat [65] methods, respectively. Other parameters were kept at their default settings. The behavior and interactions between Que and the proteins were analyzed by utilizing the Simulation Interaction Diagram tool integrated into Desmond [39].
The trajectories obtained were submitted to a post-MMGBSA analysis using the thermal_mmgbsa.py script implemented in Desmond [39]. Out of 5001 frames, one frame was analyzed every 10 ns, resulting in a total of 500 frames analyzed.
Additionally, trajectory frames were clustered based on Que’s RMSD using Schrödinger’s trj_cluster.py script [66], a fully automated clustering procedure, and output files were generated for the 5 most populated clusters in each MD experiment. Each cluster was then submitted to MMGBSA calculations with Prime [59,60].

5. Conclusions

Que, a flavonoid widely distributed in fruits and vegetables, is known for its anti-microbial, anti-viral, anti-oxidant, and anti-inflammatory properties. Its anti-oxidant activity is mainly attributed to its intrinsic ability to neutralize free radicals due to the hydroxyl substitutions and the catechol-type B-ring. Additionally, it can inhibit the expression of pro-inflammatory and oxidant genes and endogenous oxidizing enzymes. The results of docking studies show that Que can interact with oxidative enzymes, making it a promising candidate for developing new anti-inflammatory agents. The study found that Que demonstrated the greatest effectiveness in molecular docking simulations with XO, followed by MAO-A, 5-LOX, NOX, and MPO, which are implicated in several inflammatory pathologies. To confirm these findings and explore the mechanisms of action further, MD simulation studies were conducted over a 500 ns period. These simulations confirmed binding and contributed to a more detailed understanding of the interactions between Que and the five oxidative enzymes. For the NOX and MPO simulations, clear variations in binding with Que were observed, while excellent binding energy values calculated with MMGBSA were observed for all simulations. In addition, all enzymes except MPO showed strong interactions preserved during the simulations. In fact, the only stable and weak interaction with MPO was with Glu116 (19%). These findings collectively emphasize the remarkable potential of Que as a versatile compound with significant implications for its development as a novel class of anti-inflammatory agents, provided strategies to enhance its bioavailability are explored and implemented.
However, it is crucial to acknowledge that Que has its limitations, particularly in terms of bioavailability. Its aglycone form exhibits low systemic exposure. This is an important consideration for the potential therapeutic use of Que. Novel delivery methods, including inclusion complexes, prodrugs, nanocrystals, microemulsions, liposomes, and phospholipid formulations, as well as polymer nanoparticles and micelles, hold significant promise and warrant further exploration to improve que’s bioavailability, medicinal value, and applicability.

Supplementary Materials

The supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ijms242015391/s1.

Author Contributions

A.D.P., S.O. and C.S. wrote the article and created tables and figures. M.C.F. and B.E. reviewed the article. A.D.P. and A.F. critically reviewed the article and supervised the project. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.M.; Zhang, Z.Y.; Wang, R.X. Protective Mechanisms of Quercetin against Myocardial Ischemia Reperfusion Injury. Front. Physiol. 2020, 11, 506938. [Google Scholar] [CrossRef]
  2. Wang, G.; Wang, Y.; Yao, L.; Gu, W.; Zhao, S.; Shen, Z.; Lin, Z.; Liu, W.; Yan, T. Pharmacological Activity of Quercetin: An Updated Review. Evid. Based Complement. Altern. Med. 2022, 2022, 3997190. [Google Scholar] [CrossRef] [PubMed]
  3. Di Petrillo, A.; Orrù, G.; Fais, A.; Fantini, M.C. Quercetin and Its Derivates as Antiviral Potentials: A Comprehensive Review. Phytother. Res. 2022, 36, 266–278. [Google Scholar] [CrossRef] [PubMed]
  4. Kim, J.K.; Park, S.U. Quercetin and Its Role in Biological Functions: An Updated Review. EXCLI J. 2018, 17, 856. [Google Scholar] [CrossRef]
  5. Lotito, S.B.; Frei, B. Consumption of Flavonoid-Rich Foods and Increased Plasma Antioxidant Capacity in Humans: Cause, Consequence, or Epiphenomenon? Free Radic. Biol. Med. 2006, 41, 1727–1746. [Google Scholar] [CrossRef] [PubMed]
  6. Phaniendra, A.; Jestadi, D.B.; Periyasamy, L. Free Radicals: Properties, Sources, Targets, and Their Implication in Various Diseases. Indian J. Clin. Biochem. 2015, 30, 11. [Google Scholar] [CrossRef] [PubMed]
  7. Liguori, I.; Russo, G.; Curcio, F.; Bulli, G.; Aran, L.; Della-Morte, D.; Gargiulo, G.; Testa, G.; Cacciatore, F.; Bonaduce, D.; et al. Oxidative Stress, Aging, and Diseases. Clin. Interv. Aging 2018, 13, 757–772. [Google Scholar] [CrossRef]
  8. Sharifi-Rad, M.; Anil Kumar, N.V.; Zucca, P.; Varoni, E.M.; Dini, L.; Panzarini, E.; Rajkovic, J.; Tsouh Fokou, P.V.; Azzini, E.; Peluso, I.; et al. Lifestyle, Oxidative Stress, and Antioxidants: Back and Forth in the Pathophysiology of Chronic Diseases. Front. Physiol. 2020, 11, 694. [Google Scholar] [CrossRef] [PubMed]
  9. Sarniak, A.; Lipińska, J.; Tytman, K.; Lipińska, S. Endogenous Mechanisms of Reactive Oxygen Species (ROS) Generation. Postep. Hig. I Med. Dosw. 2016, 70, 1150–1165. [Google Scholar] [CrossRef]
  10. Lu, N.; Sui, Y.; Tian, R.; Peng, Y.Y. Inhibitive Effects of Quercetin on Myeloperoxidase-Dependent Hypochlorous Acid Formation and Vascular Endothelial Injury. J. Agric. Food Chem. 2018, 66, 4933–4940. [Google Scholar] [CrossRef]
  11. Luo, M.; Tian, R.; Yang, Z.; Peng, Y.Y.; Lu, N. Quercetin Suppressed NADPH Oxidase-Derived Oxidative Stress via Heme Oxygenase-1 Induction in Macrophages. Arch. Biochem. Biophys. 2019, 671, 69–76. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, C.; Wang, R.; Zhang, G.; Gong, D. Mechanistic Insights into the Inhibition of Quercetin on Xanthine Oxidase. Int. J. Biol. Macromol. 2018, 112, 405–412. [Google Scholar] [CrossRef]
  13. Abdalla, F.H.; Cardoso, A.M.; Schmatz, R.; Gonçalves, J.F.; Baldissarelli, J.; Martins, C.C.; Zanini, D.; De Oliveira, L.S.; Da Costa, P.; Pimentel, V.C.; et al. Protective Effect of Quercetin in Ecto-Enzymes, Cholinesterases, and Myeloperoxidase Activities in the Lymphocytes of Rats Exposed to Cadmium. Mol. Cell. Biochem. 2014, 396, 201–211. [Google Scholar] [CrossRef] [PubMed]
  14. Da Silva, E.L.; Abdalla, D.S.P.; Terao, J. Inhibitory Effect of Flavonoids on Low-Density Lipoprotein Peroxidation Catalyzed by Mammalian 15-Lipoxygenase. IUBMB Life 2000, 49, 289–295. [Google Scholar] [CrossRef] [PubMed]
  15. Yoshino, S.; Hara, A.; Sakakibara, H.; Kawabata, K.; Tokumura, A.; Ishisaka, A.; Kawai, Y.; Terao, J. Effect of Quercetin and Glucuronide Metabolites on the Monoamine Oxidase-A Reaction in Mouse Brain Mitochondria. Nutrition 2011, 27, 847–852. [Google Scholar] [CrossRef]
  16. Panday, A.; Sahoo, M.K.; Osorio, D.; Batra, S. NADPH Oxidases: An Overview from Structure to Innate Immunity-Associated Pathologies. Cell. Mol. Immunol. 2014, 12, 5–23. [Google Scholar] [CrossRef] [PubMed]
  17. Bortolotti, M.; Polito, L.; Battelli, M.G.; Bolognesi, A. Xanthine Oxidoreductase: One Enzyme for Multiple Physiological Tasks. Redox Biol. 2021, 41, 101882. [Google Scholar] [CrossRef] [PubMed]
  18. Braga, T.T.; Forni, M.F.; Correa-Costa, M.; Ramos, R.N.; Barbuto, J.A.; Branco, P.; Castoldi, A.; Hiyane, M.I.; Davanso, M.R.; Latz, E.; et al. Soluble Uric Acid Activates the NLRP3 Inflammasome. Sci. Rep. 2017, 7, 39884. [Google Scholar] [CrossRef] [PubMed]
  19. Rehring, J.F.; Bui, T.M.; Galán-Enríquez, C.S.; Urbanczyk, J.M.; Ren, X.; Wiesolek, H.L.; Sullivan, D.P.; Sumagin, R. Released Myeloperoxidase Attenuates Neutrophil Migration and Accumulation in Inflamed Tissue. Front. Immunol. 2021, 12, 654259. [Google Scholar] [CrossRef]
  20. Singh, P.; Arif, Y.; Miszczuk, E.; Bajguz, A.; Hayat, S. Specific Roles of Lipoxygenases in Development and Responses to Stress in Plants. Plants 2022, 11, 979. [Google Scholar] [CrossRef]
  21. Dhiman, P.; Malik, N.; Sobarzo-Sánchez, E.; Uriarte, E.; Khatkar, A. Quercetin and Related Chromenone Derivatives as Monoamine Oxidase Inhibitors: Targeting Neurological and Mental Disorders. Molecules 2019, 24, 418. [Google Scholar] [CrossRef] [PubMed]
  22. Finberg, J.P.M.; Rabey, J.M. Inhibitors of MAO-A and MAO-B in Psychiatry and Neurology. Front. Pharmacol. 2016, 7, 340. [Google Scholar] [CrossRef]
  23. Uttara, B.; Singh, A.; Zamboni, P.; Mahajan, R. Oxidative Stress and Neurodegenerative Diseases: A Review of Upstream and Downstream Antioxidant Therapeutic Options. Curr. Neuropharmacol. 2009, 7, 65–74. [Google Scholar] [CrossRef] [PubMed]
  24. Schrödinger Software Release 2022-1|Macs in Chemistry. Available online: https://www.macinchem.org/blog/files/279f5ca3fa90d6bf7db20adc4b0afed0-2841.php (accessed on 17 May 2023).
  25. Guo, M.; Qin, S.; Wang, S.; Sun, M.; Yang, H.; Wang, X.; Fan, P.; Jin, Z. Herbal Medicine Nanocrystals: A Potential Novel Therapeutic Strategy. Molecules 2023, 31, 6370. [Google Scholar] [CrossRef] [PubMed]
  26. Kandemir, K.; Tomas, M.; McClements, D.J.; Capanoglu, E. Recent Advances on the Improvement of Quercetin Bioavailability. Trends Food Sci. Technol. 2022, 119, 192–200. [Google Scholar] [CrossRef]
  27. Cai, X.; Fang, Z.; Dou, J.; Yu, A.; Zhai, G. Bioavailability of Quercetin: Problems and Promises. Curr. Med. Chem. 2013, 20, 2572–2582. [Google Scholar] [CrossRef]
  28. Van Mourik, T.; Bühl, M.; Gaigeot, M.P. Density Functional Theory across Chemistry, Physics and Biology. Philos. Trans. R. Soc. A Math. Phys. Eng. Sci. 2014, 372, 20120488. [Google Scholar] [CrossRef]
  29. Gadre, S.R.; Suresh, C.H.; Mohan, N. Electrostatic Potential Topology for Probing Molecular Structure, Bonding and Reactivity. Molecules 2021, 26, 3289. [Google Scholar] [CrossRef]
  30. Cukrowski, I.; Dhimba, G.; Riley, D.L. A Molecular-Wide and Electron Density-Based Approach in Exploring Chemical Reactivity and Explicit Dimethyl Sulfoxide (DMSO) Solvent Molecule Effects in the Proline Catalyzed Aldol Reaction. Molecules 2022, 27, 962. [Google Scholar] [CrossRef]
  31. Islam, N. Investigation of Comparative Shielding of Morin against Oxidative Damage by Radicals: A DFT Study. Cogent Chem. 2015, 1, 1078272. [Google Scholar] [CrossRef]
  32. Weinhold, F.; Landis, C.R. Natural Bond Orbitals and Extensions of Localized Bonding Concepts. Chem. Educ. Res. Pract. 2001, 2, 91–104. [Google Scholar] [CrossRef]
  33. Gasteiger, J.; Marsili, M. A New Model for Calculating Atomic Charges in Molecules. Tetrahedron Lett. 1978, 19, 3181–3184. [Google Scholar] [CrossRef]
  34. Zheng, Y.Z.; Deng, G.; Liang, Q.; Chen, D.F.; Guo, R.; Lai, R.C. Antioxidant Activity of Quercetin and Its Glucosides from Propolis: A Theoretical Study. Sci. Rep. 2017, 7, 7543. [Google Scholar] [CrossRef] [PubMed]
  35. Politzer, P.; Murray, J.S.; Bulat, F.A. Average Local Ionization Energy: A Review. J. Mol. Model. 2010, 16, 1731–1742. [Google Scholar] [CrossRef] [PubMed]
  36. Sjoberg, P.; Murray, J.S.; Brinck, T.; Politzer’, P.; Politzer, P. Average Local Ionization Energies on the Molecular Surfaces of Aromatic Systems as Guides to Chemical Reactivity. Can. J. Chem. 2011, 68, 1440–1443. [Google Scholar] [CrossRef]
  37. Cao, H.; Pauff, J.M.; Hille, R. X-Ray Crystal Structure of a Xanthine Oxidase Complex with the Flavonoid Inhibitor Quercetin. J. Nat. Prod. 2014, 77, 1693–1699. [Google Scholar] [CrossRef]
  38. Enroth, C.; Eger, B.T.; Okamoto, K.; Nishino, T.; Nishino, T.; Pai, E.F. Crystal Structures of Bovinemilk Xanthine Dehydrogenase and Xanthine Oxidase: Structure-Based Mechanism of Conversion. Proc. Natl. Acad. Sci. USA 2000, 97, 10723–10728. [Google Scholar] [CrossRef]
  39. Bowers, K.J.; Chow, D.E.; Xu, H.; Dror, R.O.; Eastwood, M.P.; Gregersen, B.A.; Klepeis, J.L.; Kolossvary, I.; Moraes, M.A.; Sacerdoti, F.D.; et al. Scalable Algorithms for Molecular Dynamics Simulations on Commodity Clusters. In Proceedings of the 2006 ACM/IEEE Conference on Supercomputing, Tampa, FL, USA, 11–17 November 2006; p. 43. [Google Scholar] [CrossRef]
  40. Vo, Q.V.; Nam, P.C.; Thong, N.M.; Trung, N.T.; Phan, C.T.D.; Mechler, A. Antioxidant Motifs in Flavonoids: O−H versus C−H Bond Dissociation. ACS Omega 2019, 4, 8935–8942. [Google Scholar] [CrossRef]
  41. Bader, R.F.W.; Nguyen-Dang, T.T. Quantum Theory of Atoms in Molecules–Dalton Revisited. Adv. Quantum Chem. 1981, 14, 63–124. [Google Scholar] [CrossRef]
  42. Pereira, G.S.; Percebom, I.; Mendes, S.; Souza, P.S.S.; Diniz, L.F.A.; Costa, M.F.; Lopes, B.R.P.; Toledo, K.A. Quercetin Inhibits Neutrophil Extracellular Traps Release and Their Cytotoxic Effects on A549 Cells, as Well the Release and Enzymatic Activity of Elastase and Myeloperoxidase. Braz. J. Biol. 2022, 84, e252936. [Google Scholar] [CrossRef]
  43. Vyshnevska, L.; Severina, H.I.; Prokopenko, Y.; Shmalko, A. Molecular Docking Investigation of Anti-Inflammatory Herbal Compounds as Potential LOX-5 and COX-2 Inhibitors. Pharmacia 2022, 69, 733–744. [Google Scholar] [CrossRef]
  44. Larit, F.; Elokely, K.M.; Chaurasiya, N.D.; Benyahia, S.; Nael, M.A.; León, F.; Abu-Darwish, M.S.; Efferth, T.; Wang, Y.H.; Belouahem-Abed, D.; et al. Inhibition of Human Monoamine Oxidase A and B by Flavonoids Isolated from Two Algerian Medicinal Plants. Phytomedicine 2018, 40, 27–36. [Google Scholar] [CrossRef]
  45. Bochevarov, A.D.; Harder, E.; Hughes, T.F.; Greenwood, J.R.; Braden, D.A.; Philipp, D.M.; Rinaldo, D.; Halls, M.D.; Zhang, J.; Friesner, R.A. Jaguar: A High-Performance Quantum Chemistry Software Program with Strengths in Life and Materials Sciences. Int. J. Quantum Chem. 2013, 113, 2110–2142. [Google Scholar] [CrossRef]
  46. Becke, A.D. Density-functional Thermochemistry. III. The Role of Exact Exchange. J. Chem. Phys. 1993, 98, 5648–5652. [Google Scholar] [CrossRef]
  47. Leopoldini, M.; Russo, N.; Toscano, M. Gas and Liquid Phase Acidity of Natural Antioxidants. J. Agric. Food Chem. 2006, 54, 3078–3085. [Google Scholar] [CrossRef]
  48. Brovarets’, O.O.; Hovorun, D.M. A Never-Ending Conformational Story of the Quercetin Molecule: Quantum-Mechanical Investigation of the O3′H and O4′H Hydroxyl Groups Rotations. Appl. Sci. 2020, 10, 1147. [Google Scholar] [CrossRef]
  49. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef]
  50. Madhavi Sastry, G.; Adzhigirey, M.; Day, T.; Annabhimoju, R.; Sherman, W. Protein and Ligand Preparation: Parameters, Protocols, and Influence on Virtual Screening Enrichments. J. Comput. Aided Mol. Des. 2013, 27, 221–234. [Google Scholar] [CrossRef]
  51. Westbrook, J.D.; Shao, C.; Feng, Z.; Zhuravleva, M.; Velankar, S.; Young, J. The Chemical Component Dictionary: Complete Descriptions of Constituent Molecules in Experimentally Determined 3D Macromolecules in the Protein Data Bank. Bioinformatics 2015, 31, 1274–1278. [Google Scholar] [CrossRef]
  52. Søndergaard, C.R.; Olsson, M.H.M.; Rostkowski, M.; Jensen, J.H. Improved Treatment of Ligands and Coupling Effects in Empirical Calculation and Rationalization of p K a Values. J. Chem. Theory Comput. 2011, 7, 2284–2295. [Google Scholar] [CrossRef]
  53. Eberhardt, J.; Santos-Martins, D.; Tillack, A.F.; Forli, S. AutoDock Vina 1.2.0: New Docking Methods, Expanded Force Field, and Python Bindings. J. Chem. Inf. Model. 2021, 61, 3891–3898. [Google Scholar] [CrossRef]
  54. Das, S.; Shimshi, M.; Raz, K.; Nitoker Eliaz, N.; Mhashal, A.R.; Ansbacher, T.; Major, D.T. EnzyDock: Protein-Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes. J. Chem. Theory Comput. 2019, 15, 5116–5134. [Google Scholar] [CrossRef] [PubMed]
  55. Friesner, R.A.; Murphy, R.B.; Repasky, M.P.; Frye, L.L.; Greenwood, J.R.; Halgren, T.A.; Sanschagrin, P.C.; Mainz, D.T. Extra Precision Glide: Docking and Scoring Incorporating a Model of Hydrophobic Enclosure for Protein-Ligand Complexes. J. Med. Chem. 2006, 49, 6177–6196. [Google Scholar] [CrossRef] [PubMed]
  56. Friesner, R.A.; Banks, J.L.; Murphy, R.B.; Halgren, T.A.; Klicic, J.J.; Mainz, D.T.; Repasky, M.P.; Knoll, E.H.; Shelley, M.; Perry, J.K.; et al. Glide: A New Approach for Rapid, Accurate Docking and Scoring. 1. Method and Assessment of Docking Accuracy. J. Med. Chem. 2004, 47, 1739–1749. [Google Scholar] [CrossRef]
  57. Banks, J.L.; Beard, H.S.; Cao, Y.; Cho, A.E.; Damm, W.; Farid, R.; Felts, A.K.; Halgren, T.A.; Mainz, D.T.; Maple, J.R.; et al. Integrated Modeling Program, Applied Chemical Theory (IMPACT). J. Comput. Chem. 2005, 26, 1752–1780. [Google Scholar] [CrossRef]
  58. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA Methods to Estimate Ligand-Binding Affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
  59. Jacobson, M.P.; Friesner, R.A.; Xiang, Z.; Honig, B. On the Role of the Crystal Environment in Determining Protein Side-Chain Conformations. J. Mol. Biol. 2002, 320, 597–608. [Google Scholar] [CrossRef] [PubMed]
  60. Jacobson, M.P.; Pincus, D.L.; Rapp, C.S.; Day, T.J.F.; Honig, B.; Shaw, D.E.; Friesner, R.A. A Hierarchical Approach to All-Atom Protein Loop Prediction. Proteins Struct. Funct. Bioinform. 2004, 55, 351–367. [Google Scholar] [CrossRef]
  61. Li, J.; Abel, R.; Zhu, K.; Cao, Y.; Zhao, S.; Friesner, R.A. The VSGB 2.0 Model: A next Generation Energy Model for High Resolution Protein Structure Modeling. Proteins Struct. Funct. Bioinform. 2011, 79, 2794–2812. [Google Scholar] [CrossRef] [PubMed]
  62. Lu, C.; Wu, C.; Ghoreishi, D.; Chen, W.; Wang, L.; Damm, W.; Ross, G.A.; Dahlgren, M.K.; Russell, E.; Von Bargen, C.D.; et al. OPLS4: Improving Force Field Accuracy on Challenging Regimes of Chemical Space. J. Chem. Theory Comput. 2021, 17, 4291–4300. [Google Scholar] [CrossRef]
  63. Jorgensen, W.L.; Chandrasekhar, J.; Madura, J.D.; Impey, R.W.; Klein, M.L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  64. Hoover, W.G. Canonical Dynamics: Equilibrium Phase-Space Distributions. Phys. Rev. A 1985, 31, 1695. [Google Scholar] [CrossRef] [PubMed]
  65. Martyna, G.J.; Tobias, D.J.; Klein, M.L. Constant Pressure Molecular Dynamics Algorithms. J. Chem. Phys. 1994, 101, 4177–4189. [Google Scholar] [CrossRef]
  66. Frey, B.J.; Dueck, D. Clustering by Passing Messages between Data Points. Science 2007, 315, 972–976. [Google Scholar] [CrossRef] [PubMed]
Ijms 24 15391 g001
Figure 1. Molecular orbitals of Que and energy gap (ΔE). Blue indicates a positive phase, while red corresponds to a negative phase.
Figure 1. Molecular orbitals of Que and energy gap (ΔE). Blue indicates a positive phase, while red corresponds to a negative phase.
Ijms 24 15391 g001
Ijms 24 15391 g002
Figure 2. Ligand interaction diagram of Que docked in NOX. Purple arrows represent H-bonds.
Figure 2. Ligand interaction diagram of Que docked in NOX. Purple arrows represent H-bonds.
Ijms 24 15391 g002
Ijms 24 15391 g003
Figure 3. Ligand interaction diagram of Que docked in XO. Purple arrows represent H-bonds, whereas green lines represent π-π stackings. The MOS 1328 is the dioxothiomolybdenum (VI) ion.
Figure 3. Ligand interaction diagram of Que docked in XO. Purple arrows represent H-bonds, whereas green lines represent π-π stackings. The MOS 1328 is the dioxothiomolybdenum (VI) ion.
Ijms 24 15391 g003
Ijms 24 15391 g004
Figure 4. Ligand interaction diagram of Que docked in MPO. Purple arrows represent H-bonds, whereas green lines represent π-π stackings.
Figure 4. Ligand interaction diagram of Que docked in MPO. Purple arrows represent H-bonds, whereas green lines represent π-π stackings.
Ijms 24 15391 g004
Ijms 24 15391 g005
Figure 5. Ligand interaction diagram of Que docked in 5-LOX. Purple arrows represent H-bonds, and red lines represent π-cation interactions.
Figure 5. Ligand interaction diagram of Que docked in 5-LOX. Purple arrows represent H-bonds, and red lines represent π-cation interactions.
Ijms 24 15391 g005
Ijms 24 15391 g006
Figure 6. Ligand interaction diagram of Que docked in MAO-A. Purple arrows represent H-bonds, whereas green lines represent π-π stackings.
Figure 6. Ligand interaction diagram of Que docked in MAO-A. Purple arrows represent H-bonds, whereas green lines represent π-π stackings.
Ijms 24 15391 g006
Ijms 24 15391 g007
Figure 7. (A) RMSD analysis of Que (Lig fit on Prot, red) and NOX (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and NOX. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Figure 7. (A) RMSD analysis of Que (Lig fit on Prot, red) and NOX (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and NOX. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Ijms 24 15391 g007
Ijms 24 15391 g008
Figure 8. (A) RMSD analysis of Que (Lig fit on Prot, red) and XO (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and XO. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Figure 8. (A) RMSD analysis of Que (Lig fit on Prot, red) and XO (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and XO. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Ijms 24 15391 g008
Ijms 24 15391 g009
Figure 9. (A) RMSD analysis of Que (Lig fit on Prot, red) and MPO (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and MPO. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds.
Figure 9. (A) RMSD analysis of Que (Lig fit on Prot, red) and MPO (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and MPO. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds.
Ijms 24 15391 g009
Ijms 24 15391 g010
Figure 10. (A) RMSD analysis of Que (Lig fit on Prot, red) and 5-LOX (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and 5-LOX. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while red lines represent π-cation interactions.
Figure 10. (A) RMSD analysis of Que (Lig fit on Prot, red) and 5-LOX (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and 5-LOX. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while red lines represent π-cation interactions.
Ijms 24 15391 g010
Ijms 24 15391 g011
Figure 11. (A) RMSD analysis of Que (Lig fit on Prot, red) and MAO-A (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and MAO-A. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Figure 11. (A) RMSD analysis of Que (Lig fit on Prot, red) and MAO-A (Cα, blue). On the x-axis, the simulation time is represented in ns, while on the y-axis, the protein RMSD (left) and ligand RMSD (right) are plotted. (B) Detailed interactions between Que and MAO-A. Interactions happening for more than 10% of the simulation time are shown. Purple arrows represent H-bonds, while green lines represent π-π stackings.
Ijms 24 15391 g011
Table 1. Theoretical results of Que using B3LYP/6-311++G** method with Jaguar software (Schrödinger Release 2022-1).
Table 1. Theoretical results of Que using B3LYP/6-311++G** method with Jaguar software (Schrödinger Release 2022-1).
HOMO−5.68 eV
LUMO−1.94 eV
ΔE HOMO-LUMO3.74 eV
IP5.68 eV
EA1.94 eV
ESP−2.65 eV (min value) 3.18 eV (max value)
ALIE8.84 eV (min value) 16.81 eV (max value)
Table 2. Selected X-ray structures from the Protein Data Bank.
Table 2. Selected X-ray structures from the Protein Data Bank.
PDB IDResolution (Å)EnzymeOrganismCo-Crystallized Ligand(s)Ref.
2CDU1.80NOXFructilactobacillus sanfranciscensisFlavin-adenine dinucleotide (FAD), adenosine-5′-diphosphate (ADP)[34]
2Z5Y2.17MAO-AHomo sapiensFlavin-adenine dinucleotide (FAD), decyl(dimethyl)phosphine oxide (DCX), 7-methoxy-1-methyl-9h-beta-carboline (HRM)[35]
3NVY2.00XOBos taurusMolybdopterin (MTE), quercetin (QUE)[36]
5FIW1.70MPOHomo sapiens2-acetamido-2-deoxy-beta-D-glucopyranose (NAG), beta-D-mannopyranose (BMA)[37]
6NCF2.875-LOXHomo sapiens(3alpha,8alpha,17alpha,18alpha)-3-(acetyloxy)-11-oxours-12-en-23-oic acid (AFX)[38]
Table 3. Docking and Prime MMGBSA results for the interactions of Que with the analyzed enzymes. The results are all expressed in kcal/mol.
Table 3. Docking and Prime MMGBSA results for the interactions of Que with the analyzed enzymes. The results are all expressed in kcal/mol.
EnzymeGlide GscoreMMGBSA dG Bind
NOX−5.35−4.82
XO−8.87−22.18
MPO−4.91−28.81
5-LOX−5.94−37.42
MAO-A−8.145.72
Table 4. Calculated clusters derived from the trajectory of the MD simulation on the Que-NOX complex.
Table 4. Calculated clusters derived from the trajectory of the MD simulation on the Que-NOX complex.
ClusterN. of MembersMMGBSA dG Bind
162−47.260 kcal/mol
259−48.970 kcal/mol
357−55.649 kcal/mol
456−51.466 kcal/mol
556−54.908 kcal/mol
Table 5. Calculated clusters derived from the trajectory of MD simulation on the Que-XO complex.
Table 5. Calculated clusters derived from the trajectory of MD simulation on the Que-XO complex.
ClusterN. of MembersMMGBSA dG Bind
1151−32.656 kcal/mol
297−39.741 kcal/mol
396−36.597 kcal/mol
491−37.437 kcal/mol
586−28.251 kcal/mol
Table 6. Calculated clusters derived from the trajectory of MD simulation on the Que-MPO complex.
Table 6. Calculated clusters derived from the trajectory of MD simulation on the Que-MPO complex.
ClusterN. of MembersMMGBSA dG Bind
1108−30.793 kcal/mol
296−22.584 kcal/mol
394−23.547 kcal/mol
489−33.276 kcal/mol
582−29.198 kcal/mol
Table 7. Calculated clusters derived from the trajectory of MD simulation on the Que-5-LOX complex.
Table 7. Calculated clusters derived from the trajectory of MD simulation on the Que-5-LOX complex.
ClusterN. of MembersMMGBSA dG Bind
158−40.127 kcal/mol
255−41.856 kcal/mol
350−43.550 kcal/mol
450−35.036 kcal/mol
548−36.262 kcal/mol
Table 8. Calculated clusters derived from the trajectory of MD simulation on the Que-MAO-A complex.
Table 8. Calculated clusters derived from the trajectory of MD simulation on the Que-MAO-A complex.
ClusterN. of MembersMMGBSA dG Bind
1156−43.672 kcal/mol
2154−40.106 kcal/mol
3146−29.574 kcal/mol
4137−37.437 kcal/mol
5131−34.486 kcal/mol
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Olla, S.; Siguri, C.; Fais, A.; Era, B.; Fantini, M.C.; Di Petrillo, A. Inhibitory Effect of Quercetin on Oxidative Endogen Enzymes: A Focus on Putative Binding Modes. Int. J. Mol. Sci. 2023, 24, 15391. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms242015391

AMA Style

Olla S, Siguri C, Fais A, Era B, Fantini MC, Di Petrillo A. Inhibitory Effect of Quercetin on Oxidative Endogen Enzymes: A Focus on Putative Binding Modes. International Journal of Molecular Sciences. 2023; 24(20):15391. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms242015391

Chicago/Turabian Style

Olla, Stefania, Chiara Siguri, Antonella Fais, Benedetta Era, Massimo Claudio Fantini, and Amalia Di Petrillo. 2023. "Inhibitory Effect of Quercetin on Oxidative Endogen Enzymes: A Focus on Putative Binding Modes" International Journal of Molecular Sciences 24, no. 20: 15391. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms242015391

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop