Abstract: Cyclooxygenase (COX) plays a vital role in physiological process of inflammation. Therefore, we evaluated the inhibitory effect of phytocompounds derived from mangrove plants against COX protein using computational method. The 3D structure of avicenol, betulinic acid, heritonin, halprogin, lupeol, pyretrin, quercetin, rubrolide and triterpenoid obtained from PubChem database and docked against COX receptor by using Auto dock 4.0. The results indicated all the compounds have significant interaction and formation of hydrogen bonds. Lupeol acts as a potential inhibitor among the nine compounds, with binding energy -8.17 kcal mol-1 has formed two hydrogen bonds interaction at the residue THR212 and THR212. This study concludes number of effective bioactive compounds from mangrove will have great interest in pharmaceutical industry.
DUCTION
Inflammation is a protective response of body cells to various intrinsic or extrinsic stimuli such as irradiation, extreme temperature and toxic pathogens (Anderson and Borlak, 2008). Although, deregulated chronic inflammation can promote heart attack, paralysis and might contribute with atherosclerosis, hepatitis and rheumatoid arthritis (Chung et al., 2007). Non Steroidal Anti Inflammatory Drugs (NSAIDs) have therapeutic potential to rheumatoid arthritis and various types of inflammation. The target for NSAIDs is cyclooxygenase, a rate limiting receptor protein involved in the production and generation of inflammatory mediator such as prostaglandins from arachidonic acid and lipooxygenase. Particularly, COX-1 commonly known as housekeeping enzymes and constitutively expressed isoform in all tissues and associated with the production of prostaglandins (Selvam et al., 2005).
The mangrove ecosystem is one of the important coastal biodiversity. Generally, mangroves are woody plant or shrubs that grow in tropical, subtropical, estuaries, back waters, rivers. Commonly found in the Asian countries, islands of Indian Ocean, Arabian Sea, Bay of Bengal, South China and the Pacific (Qasim, l998). It contains array of primary and secondary metabolites such as alkaloids, flavonoids, terpenes, pheromones, tannins have toxicological, pharmacological and ecological importance (Bandaranayake, 2002). Specifically, phenolic compounds including flavonoid inhibit the expression of the inducible forms of COX, interleukins, tumor necrosis factor-α and adhesion molecules in inflammatory cells and tissues (Kim et al., 2004). Number of compounds was reported from mangroves which have therapeutic potential on diabetes, ulcer, inflammation, wound healing, antibiotics, cancer and cardiovascular disease by in vitro and in vivo methods (Silambarasan, 2012; Gurudeeban et al., 2013; Satyavani, 2013). Recently, Gurudeeban et al (2012a, b) evidently proved the DPPIV and alpha ketoglutarate inhibitory effect of mangrove based compounds using computational approaches, this might be minimize the usage of cost and experimental animals while entered in to in vivo studies. However, compounds with COX inhibitory effect are not validated by in silico approaches. Therefore, the present study aimed to evaluate the COX inhibitory effect of nine phyto compounds from mangrove plants using computational methods.
MATERIALS AND METHODS
Preparation of coordinate file: Three dimensional structure of cyclooxygenase was obtained from Protein Data Bank. The active sites of COX 1 protein was analysed by PDB sum. The 2D structure of compounds viz., avicenol A, betulinic acid, halprogin, heritonin, lupeol, pyretrin II, quercetin, rubrolide N and triterpenoid were retrieved from PubChem database (Table 1). The optimized ligand molecules were docked into distinguished model using Lig and Fit theory.
In silico analysis: Auto-dock Tools 4.0 was used to prepare, run and analyze the docking simulations. The pre calculated grid maps, one for each atom type present in the flexible molecules being docked and its stores the potential energy arising from the interaction with rigid macromolecules. The Lamarckian Genetic Algorithm (LGA) 23 was chosen search for the best conformers. During the docking process, a maximum of 10 conformers was considered. The population size was set to 150 and the individuals were initialized randomly. Maximum number of energy evaluation was set to 250000, maximum number of generations 50000, maximum number of top individual that automatically survived set to 1, mutation rate of 0.02, crossover rate of 0.8, step sizes were 0.2 Å for translations, 5.0 Å for quaternions and 5.0 Å for torsions. Cluster tolerance 0.5 Å, external grid energy 1000, max initial energy 0.0, max number of retries 10000 and 10 LGA runs were performed. Auto Dock was compiled and run under Windows XP operating system. Autodock results were analyzed to study the interactions and the binding energy of the docked structure. All the Auto Dock docking runs were performed in Intel CORETM i5, 64 bit operating system and 4GB RAM in Lenovo Win 7 PC.
RESULTS
The structure of protein viewed by PyMol and PDB sum predicted active sites are NAG 622, NAG 661 and ASN 68 (Fig. 1). The docking poses were ranked according to their docking scores and both the ranked list of docked ligands and their corresponding binding poses (Zhang et al., 2008). 10 docking runs were performed. Grid parameters were set as mentioned earlier and spacing between grid points was 0.375 Å. After the simulations were complete, the docked structures were analyzed and the interactions were observed.
Table 1: | List of ligand used for docking |
Fig. 1: | Active site of target receptor protein cyclooxygenase |
Table 2: | Molecular interactions of phyto compounds on cyclooxygenase receptor |
Hydrogen bond interactions and the binding distance between the donors and acceptors were measured for the best conformers. Distinct conformational clusters RMSD-tolerance and Van der Waals scaling factor was found to be 2.0 and 1.0 Å, respectively.
In silico studies: Molecular docking of mangrove derived bioactive molecules such as avicenol A, betulinic acid, heritonin, halprogin, lupeol, pyretrin, quercetin, rubrolide and triterpenoid into COX produced six clusters of conformers (Fig. 2). The binding energy, hydrogen bond interactions distance between hydrogen donor and acceptor atoms of each ligand represented in Table 2. .
Fig. 2(a-i): | Molecular interaction of cyclooxygenase receptor with mangrove derived bioactive molecules, (A) Avicenol A, (b) Betulinic acid, (c) Haloprogin, (d) Heritonin, (e) Lupeol, (f) Pyretrin II, (g) Quercetin, (h) Rubrolide N and (i) Triterpenoid |
Among the nine bioactive molecules lupeol had significant binding energy -8.17 kcal mol-1. The compounds interacted with receptor protein in the following amino acid residue LEU82, ARG388, ARG80, HIS388, TYR385, THR212, THR212, ASN 382, HIS 207, ASN 382, GLN 203, THR 206, GLU454, THR 212, ASN 382 and THR 212, respectively.
DISCUSSION
Mangrove species such as Rhizophora apiculata, Rhizophora mucronata, Excoecaria agallocha and Avicennia marina are also rich sources of phenolic, tannins, triterpenoids, alkaloids and flavonoid. In silico docking depends on the scoring functions and hydrogen bonds are used to predict their binding modes, affinities and orientation on the receptor protein (Duggan et al., 2010). In accordance to the present findings, quercetin the member of flavonoid exhibited inhibitory activity towards the COX with binding energy -5.57. Generally the flavonoids inhibited the phospholipase A2 and cyclooxygenase (COX) involved in inflammatory responses (Bitis et al., 2010). The docked ligands were selected based on docking energy and good interaction with the active sites residues and hydrogen bond donors (OH-NH group), not more than ten hydrogen bond acceptors, molecular weight under 500 g/mole, partition coefficient logP of less than 5 and rotatable bonds of less than 10 is taken as drug molecules and docking procedure is carried out (Lipinski et al., 2001). The mangrove derived triterpenoid and rubrolide-N showed good docking energy score on dihydrofolate reductase leading to the inhibition of growth of malarial parasites (Senthilraja et al., 2012). In the present finding, indicates the binding energy of -6.79 and -7.66 kcal/mol with COX 1 inhibitor by triterpenoid and rubrolide. Lupeol, Taraxerol and Betulinic acid were identified from the aerial roots of Avicennia marina (Mahera et al., 2011). We found that lupeol exhibited high energy and formed two hydrogen bonds interaction at the residue THR212 and THR212 with cluster RMSD 0.00 and reference RMSD 197.3 than compared to others. In conclusion, the binding mode of the tested nine compounds from halophytes inside the active site of COX-1 enzyme was predicted using in silico docking technique.
This study concludes that mangrove floral derived biomolecules will be an alternative source for inflammatory agent.
ACKNOWLEDGMENT
The authors are thankful to the authorities of Annamalai University for providing all support during the study period.