Virtual Screening Analysis and In-vitro Xanthine Oxidase Inhibitory Activity of Some Commercially Available Flavonoids

authors:

avatar Muthuswamy Umamaheswari , * , avatar Arumugam Madeswaran , avatar Kuppusamy Asokkumar

College of Pharmacy, Sri Ramakrishna Institute of Paramedical Sciences, Coimbatore-641044, Tamil Nadu, India.
Iranian Journal of Pharmaceutical Research : IJPR: Vol.12, issue 3; 317-323
article type: Original Article
received: November , 2011
accepted: February , 2012
DOI: https://doi.org/10.22037/ijpr.2013.1323

how to cite: Umamaheswari M, Madeswaran A, Asokkumar K. Virtual Screening Analysis and In-vitro Xanthine Oxidase Inhibitory Activity of Some Commercially Available Flavonoids. Iran J Pharm Res. 2013;12(3):e125675. https://doi.org/10.22037/ijpr.2013.1323.

Abstract

Allopurinol, the xanthine oxidase inhibitor, is the only drug available for the treatment of gout. We examined the xanthine oxidase inhibitory activity of some commercially available flavonoids such asepigallocatechin, acacatechin, myricetin, naringenin, daidzein and glycitein by virtual screening and in-vitro studies. The interacting residues within the complex model and their contact types were identified. The virtual screening analysis were carried out using AutoDock 4.2 and in-vitro xanthine oxidase inhibitory activity was carried out using xanthine as the substrate. In addition, enzyme kinetics was performed using LineweaverBurkplot analysis. Allopurinol, a known xanthine oxidase inhibitor was used as the standard. The docking energy ofglycitein was found to be -8.49 kcal/mol which was less than that of the standard (-4.47 kcal/ mol). All the selected flavonoids were found to exhibit lower binding energy (-8.08 to -6.03 kcal/ mol) than allopurinol. The docking results confirm that flavonoids showed greater inhibition of xanthine oxidase due to their active binding sites and lesser binding energies compared to allopurinol. This may be attributed to the presence of benzopyran ring in the flavonoids. In the xanthine oxidase assay, IC50 value of glycitein was found to be 12±0.86 μg/mL, whereas that of allopurinol was 24±0.28 μg/mL. All the remaining compounds exhibited IC50 values ranging between 22±0.64 to 62±1.18 μg/mL. In the enzyme kinetic studies, flavonoids showed competitive type of enzyme inhibition. It can be concluded that flavonoids could be a promising remedy for the treatment of gout and related inflammatory disorders. Further in-vivo studies are required to develop potential compounds with lesser side effects.

Introduction

Drug discovery and development is a complex, long term and interdisciplinary process. It is a multidimensional and sequential process that begins from target identification, lead discovery process, followed by lead optimization and pre-clinical in-vitro and in-vivo studies (1). Virtual screening of compound libraries has become a standard technology in modern drug discovery pipelines (2). Traditionally, drugs were synthesized from a variety of compounds and screened for its toxicity and biological activities and additionally examined for their pharmacokinetic profile. However, this process is generally time consuming (3). Structure based drug design is becoming a valuable and integral part of drug discovery process, which has been proven to be more effective than the ligand based drug design (4). Studies of interactions between protein domains and ligands are important in virtual screening analysis (5).

Virtual screening analysis can help in identifying drug targets via bioinformatics tools. They are used to analyze the target structures for possible binding sites, generation of candidate molecules, checking for their drug likeness, docking the molecules with the target, ranking them according to their binding affinities, and further optimization of the molecules to improve binding characteristics (6). Autodock 4.2 is a suite of automated docking tools. It usually starts with the definition of a binding site, in general at a restricted region of the protein. Autodock uses Monte Carlo and Simulated Annealing in combination with Genetic Algorithm which is used for global optimization (7).

Xanthine oxidase (XO) is a highly versatile enzyme that is widely distributed among different species from bacteria to man and within the various tissues of mammals. It is a member of group of enzymes known as molybdenum iron – sulphur flavin hydroxylases (8). It catalyses the oxidation of hypoxanthine to xanthine and then to uric acid, the final reactions in the metabolism of purine bases (9). The accumulation of uric acid in the body is responsible for several diseases and thus it plays a vital role in hyperuraecimia and gout (10). Inherited xanthine oxidase reductase (XOR) deficiency leads to xanthineuria and multiple organ failure syndrome caused by the accumulation of xanthine in different tissues (11).

Xanthine oxidase inhibitors (XOI) are much useful, since they possess lesser side effects compared to uricosuric and anti inflammatory agents (12). Allopurinol is the only clinically available XOI, which also suffers from many side effects such as hypersensitivity syndrome, Steven’s Johnson syndrome and renal toxicity. Thus, itis necessary to develop compounds with XOI activity with lesser side effects compared to allopurinol. Flavonoids and polyphenols have been reported to possess xanthine oxidase inhibitory activity (13).In addition, flavonoids also have anti inflammatory and antitumor properties (14). We thus began our work to look for virtual screening analysis and in-vitro xanthine oxidase inhibitory activity of some commercially available flavonoids.

Experimental

Softwares required

Python 2.7 - language was downloaded from www.python.com, Cygwin (a data storage) c:\program and Python 2.5 were simultaneously downloaded from www.cygwin.com, Molecular graphics laboratory (MGL) tools and AutoDock 4.2 was downloaded from studio visualizer 2.5.5 was downloaded from www.accelerys.com, Molecular orbital package (MOPAC), Chemsketch was downloaded from www.acdlabs.com. Online smiles translatory notation was carried out using cactus.nci.nih.gov/translate/.

Chemicals required

Allopurinol, xanthine, xanthine oxidase from bovine milk source and flavonoids such asepigallocatechin, acacatechin, myricetin, naringenin, daidzein, glyciteinwere purchased from Sigma Aldrich, USA. All other drugs and chemicals used in the study were obtained commercially and were of analytical grade.

Virtual screening analysis

The ligands such asepigallocatechin, acacatechin, myricetin, naringenin, daidzein, glycitein were built using Chemsketch and optimized using “Prepare Ligands” in the Autodock 4.2 for docking studies. Xanthine oxidase model from bovine milk source was downloaded from the RCSB protein data bank (Figure 1).

Xanthine oxidase from bovine milk source (3BDJ).

Xanthine oxidase from bovine milk source (3BDJ).

Figure 1

The optimized ligand molecules were docked into refined xanthine oxidase model using “Ligand Fit” in the Autodock 4.2 (15). These file preparations were carried out by plugin using scripts from the Autodock Tools package. The ligand score is an automated tool for protein-ligand docking that can define binding site, generate ligand conformations, dock each conformation, save the top docked structures (diverse poses) and apply scoring functions to each docked structure for the best binding mode. The binding sites for these molecules were selected based on the ligand-binding pocket of the templates. For each ligand, 10 poses were generated and scored using Autodock 4.2 scoring functions. Among these poses, the most suitable docking mode for flavone with a high score from consensus scoring functions was finally selected (16 - 18).

The binding position of ligands to the enzymexanthine oxidase (1 - Glycitein,2 - Naringenin, 3 - Daidzein,4 - Myricitein, 5 -Acacatechin,6 - Epigallocatechin and 7 - Allopurinol).

The binding position of ligands to the enzymexanthine oxidase (1 - Glycitein,2 - Naringenin, 3 - Daidzein,4 - Myricitein, 5 -Acacatechin,6 - Epigallocatechin and 7 - Allopurinol).

Figure 2

In-vitro xanthine oxidase inhibitory activity

The assay mixture consisted of 1mL of the test compound (5 – 100 μg/mL), 2.9 mL of phosphate buffer (pH 7.5) and 0.1 mL of xanthine oxidase enzyme solution (0.1 units/mL in phosphate buffer, pH 7.5), which was prepared immediately before use. After preincubation at 25 ºC for 15 min, the reaction was initiated by the addition of different concentration (5 – 100 μg/mL) of the substrate solution. The assay mixture was incubated at 25 ºC for 30 min. The reaction was stopped by adding 1 mL of 1 N HCl and the absorbance was measured at 290 nm using an UV spectrophotometer (19, 20). Allopurinol (5 – 100 μg/mL) was used as the standard. The percentage inhibition was calculated by,

Percentage inhibition ={ (A-B) – (C-D)}/(A-B)} X 100

where A is the activity of the enzyme without the compound, B is the control of A without the compound and enzyme, C and D are the activities of the compound with or without XO respectively. The assay was done in triplicate and IC50 values were calculated from the percentage inhibition (21).

Enzyme kinetics studies

Lineweaver – Burk plot analysis was performed to determine the mode of inhibition of flavonoids and compared with allopurinol.The assay was carried out in the presence or absence of flavonoids with varying concentrations of xanthine as the substrate, employing the xanthine oxidase assay as mentioned earlier. Lineweaver – Burk transformed values were plotted to determine the mode of enzyme inhibition (22, 23).

Results and Discussion

Virtual screening analysis

The virtual screening analysis was performed by the use of Autodock 4.2. The docking poses are ranked according to their docking scores and both the ranked list of docked ligands and their corresponding binding poses (24). In Figure 3, docked pose of xanthine oxidase enzyme (green) with glycitein ligand (blue) clearly demonstrates the binding positions of the ligand with the enzyme.

As shown in Table 1, flavonoids showed binding energy ranging from -8.49 kcal/mol to -6.03 kcal/mol. All the selected compounds had lesser binding energy when compared to the standard allopurinol(-4.47 kcal/mol). This proves that flavonoids consist of potential xanthine oxidase inhibitory binding sites when compared to the standard.

Docked structure of 3BDJ enzyme model (green) withglycitein ligand (blue).

Docked structure of 3BDJ enzyme model (green) withglycitein ligand (blue).

Figure 3

Table 1

Binding energies of the compounds based on their rank

COMPOUNDSBinding energies of the compounds based on their rank (kcal/mol)
12345678910
Epigallocatechin-6.03-5.75-5.57-5.40-5.38-5.23-5.19-4.81-4.97-4.93
Acacatechin-6.22-6.13-5.87-5.79-5.35-4.95-5.96-5.91-5.22-3.99
Myricetin-7.27-6.96-5.09-7.03-6.35-5.86-5.65-5.37-4.78-4.28
Naringenin-7.60-7.40-7.25-6.70-7.18-7.13-7.16-5.79-5.58-5.36
Daidzein-7.54-7.53-7.44-6.83-5.83-5.74-5.63-5.74-5.57-5.49
Glycitein-8.49-7.28-7.20-6.71-6.68-6.66-6.55-6.53-6.46-5.70
Allopurinol-4.47-4.47-4.46-4.46-4.45-4.20-4.09-4.09-3.99-3.87

In addition, two other parameters such asinhibition constant (Ki) and intermolecular energy were also determined. As shown in Table 2, flavonoids showed inhibition constantranging from 596.29 nM to 37.79 μM. All the selected compounds had lesser inhibition constantwhen compared to the standard (529.73 μM). Inhibition constant is directly proportional to the binding energy. We found a decrease in inhibition constant of all the selected flavonoids with a simultaneous decrease in the binding energy. When the binding energy of the compound decreases, there is an increase in activity. Thus, the xanthine oxidase inhibitory activity of the flavonoids were found to be higher compared to allopurinol.

Table 2

Inhibition Constant of the compounds based on their rank

COMPOUNDSInhibition Constant of the compounds based on their rank (μM, mM*, nM**)
12345678910
Epigallocatechin37.7960.8081.95110.66113.88146.80157.99300.11226.29242.72
Acacatechin27.7031.9849.9657.13119.08233.4442.5046.46149.991.19*
Myricetin4.677.94184.827.0222.1152.0372.20115.36315.69729.56
Naringenin2.663.774.8112.195.465.945.6657.2781.19117.43
Daidzein2.993.023.539.8953.2661.9574.3362.2882.3994.29
Glycitein596.29**4.635.3112.0712.7013.0215.9316.3518.5266.36
Allopurinol529.73534.14541.00541.30545.56830.851.01*1.01*1.18*1.45*

As shown in Table 3, flavonoids showed intermolecular energy ranging between -9.68 to -8.01 which was lesser when compared to the standard (-4.47). Intermolecular energy is also directly proportional to the binding energy. We found a decrease in intermolecular energy of all the selected compounds with a simultaneous decrease in the binding energy.This result further proved the xanthine oxidase inhibitory activity of all the selected flavonoids.

Table 3

Intermolecular energies of the compounds based on their rank

COMPOUNDSInter molecular energies of the compounds based on their rank
12345678910
Epigallocatechin-8.12-7.84-7.66-7.49-7.47-7.32-7.27-6.89-7.06-7.02
Acacatechin-8.01-7.92-7.66-7.58-7.14-6.74-7.75-7.70-7.01-5.78
Myricetin-9.36-9.05-7.18-9.12-8.44-7.94-7.74-7.46-6.86-6.37
Naringenin-8.80-8.59-8.45-7.90-8.37-8.32-8.35-6.98-6.77-6.55
Daidzein-8.43-8.43-8.33-7.72-6.73-6.64-6.53-6.63-6.47-6.39
Glycitein-9.68-8.47-8.39-7.90-7.87-7.86-7.74-7.72-7.65-6.89
Allopurinol-4.47-4.47-4.46-4.46-4.45-4.20-4.09-4.09-3.99-3.87

Docking ofallopurinol in themolybdopterin domain of xanthine oxidase showed the same position and interaction as that of salicylate binding except that Phe 1009 was far away from the bicyclic ring and a hydrogen bonding to Glu 802. We superimposed glycitein with salicylate revealing that bicyclic benzopyranone ring overlapped with salicylate ring, and the phenolic group of the selected flavonoids stretched to the space surrounding with several hydrophobic residues including Phe 1076, Phe 649, Leu 648, Leu 873, and Leu 1014. Several hydrogen bonds and electrostatic interactions were exposed including C7 hydroxyl binding to Glu 1261 via water, C5 hydroxyl close to the guanidium group of Arg 880 and C4 carbonyl interaction with hydroxyl group of Thr 1010. The probabilistic interactions between residues in the microenvironment of binding site to the benzopyran ring of flavonoids decisively prove the xanthine oxidase inhibitory activity by emphasizing its structural properties.

Based on the docking studies, the activity of the selected compounds was in order of

Lineweaver-Burk plot of inhibition of xanthine oxidase by allopurinol and catechin both showing competitive type of enzyme inhibition

Lineweaver-Burk plot of inhibition of xanthine oxidase by allopurinol and catechin both showing competitive type of enzyme inhibition

Figure 4

Table 4

In-vitro xanthine oxidase inhibitory activity of the selected compounds

COMPOUNDSPercentage Inhibition (Mean±SEM)
IC50μg/mL (Mean±SEM)
Concentration (μg/mL)
520406080
Epigallocatechin31.23±0.7344.53±0.9354.32±0.6474.53±1.2187.54±0.8736±0.64
Acacatechin32.00±0.7539.33±0.9252.33±0.6476.00±0.7885.33±1.2627±1.16
Myricetin28.43±0.4246.76±0.6558.45±0.5472.82±0.7686.72±0.8426±0.72
Naringenin33.26±0.6548.46±0.8264.73±0.4882.65±0.6592.76±0.4522±0.64
Daidzein32.87±0.5446.57±0.6262.37±0.6578.65±0.8387.54±0.7623±0.75
Glycitein29.64±0.7668.76±0.9375.46±0.9887.68±0.5495.63±0.7212±0.86
Allopurinol33.33±0.7038.00±0.8666.00±1.2076.67±1.4297.33±0.9229±0.28

Conclusion

The results of the present study clearly demonstrated the xanthine oxidase inhibitory activity of the selected flavonoids by virtual screening analysis and in-vitro assay. Virtual screening analysis is actually an added advantage to screen the xanthine oxidase inhibition. The nature of this inhibition, particularly the stronger effect of glycitein,naringenin, daidzeinandmyricetinthan the standard allopurinol, is interesting and merits its further characterization. Further investigations on the above compounds and in-vivo studies are necessary to develop potential chemical entities for the prevention and treatment of gout and related inflammatory disorders.

References

  • 1.

    Reichel A. The role of blood-brain barrier studies in the pharmaceutical industry. Curr. Drug Metab. 2006;7:183-203. [PubMed ID: 16472107].

  • 2.

    Seeliger D, Groot LB. Ligand docking and binding site analysis with PyMOL and Autodock/Vina. J. Comput. Aided Mol. Des. 2010;24:417-22. [PubMed ID: 20401516].

  • 3.

    Lin CN, Kuo SH, Chung MI, Ko FN, Teng CM. A new flavone C-glycoside and antiplatelet and vasorelaxing flavones from Gentiana arisanensis. J. Nat. Prod. 1997;60:851-53. [PubMed ID: 9287421].

  • 4.

    Hirashima A, Huang H. Homology modeling, agonist binding site identification, and docking in octopamine receptor of Periplaneta americana. Comput. Biol. Chem. 2008;32:185-90. [PubMed ID: 18430608].

  • 5.

    Powers R. Advances in nuclear magnetic resonance for drug discovery. Expert Opin. Drug Disc. 2009;4:1077-98.

  • 6.

    Breda A, Basso LA, Santos DS, Walter F. Virtual screening of drugs: score functions, docking, and drug design. Curr. Comp. Aided Drug Des. 2008;4:265-72.

  • 7.

    Morris MG, Goodsell SD, Halliday SR, Huey R, Hart EW, Belew KR, Olson JA. Automated docking using a lamarckian genetic algorithm and an empirical binding free energy function. J. Comput. Chem. 1998;19:1639-62.

  • 8.

    Symons CRM, Taiwo AF, Petersen LR. Electron addition to xanthine oxidase. An electron spins resonance study of the effects of ionizing radiation. J. Chem. Soc. 1989;85:4063-74.

  • 9.

    Zarepour M, Kaspari K, Stagge S, Rethmeier R, Mendel RR, Bittner F. Xanthine dehydrogenase AtXDH1 from Arabidopsis thaliana is a potent producer of superoxide anions via its NADH oxidase activity. Plant Mol. Biol. 2010;72:301-10. [PubMed ID: 19915948].

  • 10.

    Egwim EC, Vanchi MA, Egwin PO. Comparison of xanthine oxidase activities in cow and goat milks. Biokemistri. 2005;17:1-6.

  • 11.

    Borges F, Fernandes E, Roleira F. Progress towards the discovery of xanthine oxidase inhibitors. Curr. Med. Chem. 2002;9:195-217. [PubMed ID: 11860355].

  • 12.

    Umamaheswari M, Asokkumar K, Sivashanmugam AT, Remyaraju A, Subhadradevi V, Ravi TK. In-vitro xanthine oxidase inhibitory activity of the fractions of Erythrina stricta Roxb. J. Ethnopharmacol. 2009;124:646-8. [PubMed ID: 19467311].

  • 13.

    Lio M, Moriyama A, Matsumoto Y, Takaki N, Fukumoto M. Inhibition of xanthine oxidase by flavonoids. Agric. Biol. Chem. 1985;49:2173-6.

  • 14.

    Karrasch T, Kim JS, Jang IB, Jobin C. The flavonoid luteolin worsens chemical-induced colitis in NF-κBEGFP transgenic mice through blockade of NF-κB-dependent protective molecules. PLoS ONE. 2007;2:596-8.

  • 15.

    Venkatachalam CM, Jiang X, Oldfield T, Waldman M. LigandFit: A novel method for the shape-directed rapid docking of ligands to protein active sites. J. Mol. Graph. Model. 2003;21:289-307. [PubMed ID: 12479928].

  • 16.

    Sherman DH, Li S, Yermaliskaya VL, Kim Y, Smith AJ, Waterman RM, Podust ML. The structural basis for substrate anchoring, active site selectivity, and product formation by P450 PikC from Streptomyces venezuelae. J. Biol. Chem. 2006;281:26289-97. [PubMed ID: 16825192].

  • 17.

    Yasutake Y, Imoto N, Fujii T, Arisawa A, Tamura T. Crystal structure of cytochrome P450 MoxA from Nonomuraea recticatena (CYP105). Biochem. Biophys. Res. Commun. 2007;361:876-82. [PubMed ID: 17679139].

  • 18.

    Sugimoto H, Shinkyo R, Hayashi K, Yoneda S, Yamada M, Kamakura M, Ikushiro S, Shiro Y, Sakaki T. Crystal structure of CYP105A1 (P450SU-1) in complex with 1alpha,25-dihydroxyvitamin D3. Biochem. 2008;47:4017-27. [PubMed ID: 18314962].

  • 19.

    Sahgal G, Ramanathan S, Sasidharan S, Mordi MN, Ismail S, Mansor SM. In-vitro antioxidant and xanthine oxidase inhibitory activities of methanolic Swietenia mahagoni seed extracts. Molecules. 2009;14:4476-85. [PubMed ID: 19924080].

  • 20.

    Meda NTR, Lamien MA, Kiendrebeogo M, Lamien CE, Coulibaly AY, Mjllogo RJ, Nacoulma OG. In-vitro antioxidant, xanthine oxidase and acetylcholinesterase inhibitory activities of Balanites aegyptiaca (L.) Del. (Balanitaceae). Pak. J. Biol. Sci. 2010;13:362-8. [PubMed ID: 20836295].

  • 21.

    Werns SW, Grum CM, Ventura A, Hahn RA, Ho PP, Towner RD, Fantone JC, Schork MA, Lucchesi BR. Xanthine oxidase inhibition does not limit canine infarct size. Circulation. 1991;83:995-1005. [PubMed ID: 1999047].

  • 22.

    Chen Q, Groote R, Schonherr H, Vancso GJ. Probing single enzyme kinetics in real-time. J. Chem. Soc. 2009;38:2671-83.

  • 23.

    Henderson JF. Steady-state enzyme kinetics with high-affinity substrates or inhibitors. A statistical treatment of dose–response curves. J. Biochem. 1973;135:101-7.

  • 24.

    Pacher P, Nivorozhkin A, Szabo C. Therapeutic effect of xanthine oxidase inhibitors renaissance half a century after the discovery of allopurinol. Pharmacol. Rev. 2006;58:87-114. [PubMed ID: 16507884].

  • 25.

    Unno T, Sugimoto A, Kakuda T. Xanthine oxidase inhibitors from the leaves of Lagerstroemia speciosa (L). Pers. J. Ethnopharmacol. 2004;93:391-5.

  • 26.

    Priyatno LHA, Sukandar EY, Ibrahim S, Adnyna IK. Xanthine oxidase inhibitory activity of terpenoid and pyrrole compound isolated from Snake fruit (Salacca edulis Reinw.) cv. Bongkok. J. App. Sci. 2007;7:3127-30.

  • 27.

    Nguyen MTT, Awale S, Tezuka Y, Shi L, Tran LQ, Watanabe H, Kadota H. Xanthine oxidase inhibitors from the heart wood of Vietnamese Caesalpinia sappan. Chem. Pharm. Bull. 2005;53:984-8. [PubMed ID: 16079532].

  • 28.

    Hsu C. Antioxidant activity of extracts from Polygonum aviculare L. Biol. Res. 2006;39:281-8. [PubMed ID: 16874403].