Molecular docking analysis of selected Clinacanthus nutans constituents as xanthine oxidase, nitric oxide synthase, human neutrophil elastase, matrix metalloproteinase 2, matrix metalloproteinase 9 and squalene synthase inhibitors
Radhakrishnan Narayanaswamy1, Azizul Isha1, Lam Kok Wai2, Intan Safinar Ismail3
1 Laboratory of Natural Products, Institute of Bioscience, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia
2 Faculty of Pharmacy, Universiti Kebangsaan Malaysia, Jalan Raja Muda Abdul Aziz, 50300 Kuala Lumpur, Malaysia
3 Laboratory of Natural Products, Institute of Bioscience; Department of Chemistry, Faculty of Science, Universiti Putra Malaysia, 43400 Serdang, Selangor, Malaysia
|Date of Web Publication||10-Feb-2016|
Intan Safinar Ismail
Laboratory of Natural Products (LHS), Institute of Bioscience, Universiti Putra Malaysia, 43400 Serdang, Selangor
Source of Support: None, Conflict of Interest: None
| Abstract|| |
Background: Clinacanthus nutans (Burm. f.) Lindau has gained popularity among Malaysians as a traditional plant for anti-inflammatory activity. Objective: This prompted us to carry out the present study on a selected 11 constituents of C. nutans which are clinacoside A–C, cycloclinacoside A1, shaftoside, vitexin, orientin, isovitexin, isoorientin, lupeol and β-sitosterol. Materials and Methods: Selected 11 constituents of C. nutans were evaluated on the docking behavior of xanthine oxidase (XO), nitric oxide synthase (NOS), human neutrophil elastase (HNE), matrix metalloproteinase (MMP 2 and 9), and squalene synthase (SQS) using Discovery Studio Version 3.1. Also, molecular physicochemical, bioactivity, absorption, distribution, metabolism, excretion, and toxicity (ADMET), and toxicity prediction by computer assisted technology analyzes were also carried out. Results: The molecular physicochemical analysis revealed that four ligands, namely clinacoside A–C and cycloclinacoside A1 showed nil violations and complied with Lipinski's rule of five. As for the analysis of bioactivity, all the 11 selected constituents of C. nutans exhibited active score (>0) toward enzyme inhibitors descriptor. ADMET analysis showed that the ligands except orientin and isoorientin were predicted to have Cytochrome P4502D6 inhibition effect. Docking studies and binding free energy calculations revealed that clinacoside B exhibited the least binding energy for the target enzymes except for XO and SQS. Isovitexin and isoorientin showed the potentials in the docking and binding with all of the six targeted enzymes, whereas vitexin and orientin docked and bound with only NOS and HNE. Conclusion: This present study has paved a new insight in understanding these 11 C. nutans ligands as potential inhibitors against XO, NOS, HNE, MMP 2, MMP 9, and SQS.
Keywords: Clinacanthus nutans, clinacoside B, cycloclinacoside, isoorientin, isovitexin, shaftoside
|How to cite this article:|
Narayanaswamy R, Isha A, Wai LK, Ismail IS. Molecular docking analysis of selected Clinacanthus nutans constituents as xanthine oxidase, nitric oxide synthase, human neutrophil elastase, matrix metalloproteinase 2, matrix metalloproteinase 9 and squalene synthase inhibitors. Phcog Mag 2016;12, Suppl S1:21-6
|How to cite this URL:|
Narayanaswamy R, Isha A, Wai LK, Ismail IS. Molecular docking analysis of selected Clinacanthus nutans constituents as xanthine oxidase, nitric oxide synthase, human neutrophil elastase, matrix metalloproteinase 2, matrix metalloproteinase 9 and squalene synthase inhibitors. Phcog Mag [serial online] 2016 [cited 2021 May 12];12, Suppl S1:21-6. Available from: http://www.phcog.com/text.asp?2016/12/45/21/176111
- Isovitexin and isoorientin (Clinacanthus nutans constituent) showed potentials in the docking and binding with all of the six targeted enzymes (xanthine oxidase [XO], nitric oxide synthase [NOS], human neutrophil elastase [HNE], matrix metalloproteinase [MMP 2 and 9], and squalene synthase [SQS])
- Moreover, clinacoside B (C. nutans constituent) exhibited the least binding energy for the target enzymes except for XO and SQS
- Interestingly, all of the selected ligands from C. nutans showed the potential to dock and bind with HNE.
| Introduction|| |
Clinacanthus nutans (Burm. f.) Lindauis a small shrub, belongs to the family Acantaceae and native to tropical Asia. In Malaysia, it is commonly known as Sabah snake grass, presently it is cultivated to meet the huge market demand. It has been used traditionally to treat snake bites, dysentery, diabetes, fever, skin infections, and burns., Recent studies have shown that C. nutans have the significant biological activities such as antiviral, antioxidant, anti-inflammatory, and anti-bacterial activities., However, only two enzymes which are myeloperoxidase (MPO) and elastase inhibition activity of C. nutans have been reported so far.
C. nutans has been chemically studied  wherein the important chemical constituents isolated from C. nutans are clinacoside A–C, cycloclinacoside A1, shaftoside, vitexin, orientin, isovitexin, isoorientin, lupeol, and β-sitosterol. Therefore, these 11 chemical constituents were selected to be evaluated in this study on the docking behavior of xanthine oxidase (XO), nitric oxide synthase (NOS), human neutrophil elastase (HNE), matrix metalloproteinase (MMP 2 and 9), and squalene synthase (SQS) with investigation on the enzymes putative binding sites using Discovery Studio Version 3.1. (Accelrys, USA).
| Materials and Methods|| |
Chemical structures of the 11 selected ligands namely (i) shaftoside (CID442658); (ii) vitexin (C ID5280441); (iii) orientin (CID5281675); (iv) isovitexin (CID162350); (v) isoorientin (CID114776); (vi) lupeol (CID259846), and (vii) β-sitosterol (Chemspider ID192962) were retrieved from PubMed (www.pubmed.com) and Chemspider (www.chemspider.com) compound databases. Unavailable three dimensional structures of clinacoside A–C and cycloclinacoside A1 were generated using ACD (www.acdlabs.com/download/chemsk.html).
Target protein identification and preparation
The three-dimensional structures of the XO (Protein Data Bank (PDB) ID: 1FIQ with resolution of 2.50 Å), NOS (PDB ID: 4NOS with resolution of 2.30 Å), HNE (PDB ID: 1H1B with resolution of 2.00 Å), MMP 2 (PDB ID: 1QIB with resolution of 2.80 Å), MMP 9 (PDB ID: 4H1Q with resolution of 1.59 Å), and SQS (PDB ID: 3ASX with resolution of 2.00 Å) were obtained from the Research Collaborator for Structural Bioinformatics PDB (Anonymous, www.rcsb.org). A chain of all of the proteins (except for XO, where C chain) was preprocessed separately by deleting other chains (B, C, and D), ligand, as well as the crystallographically observed water molecules (water without hydrogen bonds).
Molecular descriptors calculation
Molinspiration online database was used for all the selected twelve ligands to calculate thirteen descriptors (www.molinspiration.com) which are log P, polar surface area, molecular weight, number of atoms, number of O or N, number of OH or NH, number of rotatable bonds, volume, drug likeness including G protein coupled receptors ligand, ion channel modulator, kinase inhibitor and nuclear receptor ligand, and the number of violations to Lipinski's rule.
Absorption, distribution, metabolism, excretion, and toxicity, and toxicity prediction by computer assisted technology test
Both of absorption, distribution, metabolism, excretion and toxicity (ADMET) and toxicity prediction by computer assisted technology (TOPKAT) tests were performed using Discovery Studio® 3.1 (Accelrys, San Diego, USA). ADMET analysis was performed using human intestinal absorption, aqueous solubility, blood brain barrier, cytochrome P4502D6 (CYP2D6), plasma protein binding, and hepatotoxicity descriptors. As for the TOPKAT analysis, aerobic biodegradability (AB), Ames mutagenicity, ocular irritancy, skin irritancy, skin sensitization, and oral toxicity in rat (LD50 in g/Kg of body weight) descriptors were used.
Docking studies were carried out on the crystal structures of XO, NOS, HNE, MMP 2, MMP 9, and SQS retrieved from PDB using the CDOCKER protocol under the protein-ligand interaction section in Discovery Studio® 3.1 (Accelrys, San Diego, USA). In general, CDOCKER is a grid-based molecular docking method that employs CHARMM force fields. A protein was firstly held rigid while the ligands were allowed to flex during the refinement. Two hundred random ligand conformations were then generated from the initial ligand structure through the high temperature molecular dynamics followed by random rotations, refinement by grid-based (GRID I) simulated annealing, and a final grid-based or full force field minimization. In this experiment, the ligand was heated to a temperature of 700 K in 2000 steps and the cooling steps were set in 5000 steps to 300 K with the grid extension set to 10 Å. Hydrogen atoms were added to the structures, and all ionizable residues were set at their default protonation state at a neutral pH. For each ligand, top 10 ligand binding poses were ranked according to their CDOCKER energies, and the predicted binding interactions were then analyzed from which the best among the 10 ligand binding poses was chosen and carried out in situ ligand minimization using standard protocol.
| Results and Discussion|| |
Lipinski's rule of five is a computational tool used to predict the solubility and permeability of chemical compounds based on their molecular properties. According to Lipinski's rule, poor absorption or permeation is more likely when log P value >5; molecular weight >500; more than 5 hydrogen bond donors; more than 10 hydrogen bond acceptors, and greater than 15 rotatable bonds. In the present study, clinacoside A–C and cycloclinacoside A1 showed no violation. On other hand, orientin and isoorientin showed two violations, while shaftoside showed three violations as given in [Table 1].
|Table 1: Molecular physicochemical descriptors analysis on 11 ligands using Molinspiration online software tool|
Click here to view
Bioactivity score is another computational approach used to determine whether or not a particular molecule is similar to the known drugs in their molecular properties and structure features. According to the bioactivity score, if >0 is active; if (−5.0 to − 0.0) is moderately active and if <−5.0 is inactive. In the present study, all of the 11 C. nutans constituents showed active score (>0) toward enzyme inhibitors descriptor. However, for other descriptors, these compounds exhibited active to moderate active scores with none showing inactive score (<−5.0), as shown in [Table 2].
|Table 2: Bioactivity score of 11 ligands using Molinspiration online software tool|
Click here to view
Optimizing desirable ADMET characteristics is now recognized as an important approach for drug discovery  as many drug candidates often fail to comply with the ADMET profiles. Many pharmaceutical companies opt for using ADMET profiling with in some cases used it as an alternative for the wet screening method. [Table 3] shows the ADMET profile of the 11 ligands wherein shaftoside, vitexin, orientin, isovitexin, and isoorientin are predicated to have hepatotoxic effect compared to all other ligands. Most of the ligands except orientin and isoorientin were predicted to have CYP2D6 inhibition effect. The toxicity profile of the 11 ligands as depicted in [Table 4] shows that cycloclinacoside A1 is nondegradable toward AB nature and all of the ligands were predicted to have ocular/eye irritancy effect in humans.
XO is the key regulatory enzyme in purine metabolism which catalyzes the oxidation of hypoxanthine to xanthine and then to uric acid. Hyperuricemia is a metabolic disorder characterized by elevated level of serum uric acid. Allopurinol and febuxostat are the two XO inhibitors widely used for the treatment of hyperuricemia. However, these drugs are associated with adverse effects such as gastrointestinal, hepatic, renal, and allergic reactions. In this study, 11 selected constituents of C. nutans were evaluated on the docking behavior of XO. The docking studies and binding free energy calculations as in [Table 5] show isoorientin to be having the highest interaction energy (−49.80 kcal/mol) with that of XO. In contrast, cycloclinacoside A1 showed the least interaction energy (−25.60 kcal/mol) compared to other ligands. However, clinacoside B, shaftoside, vitexin, orientin, lupeol and β-sitosterol were the six ligands which could not dock with XO due to the general poor binding phenomenon. Interestingly, isoorientin showed interaction with Molybdenum-Oxygen-Sulfur (MOS) complex which is the key component in XO. Lin et al. reported isovitexin as a competitive inhibitor of XO and suggested that the bulky sugar moiety of isovitexin has been responsible for the interaction with XO. Isovitexin and isoorientin isolated from Biophytum umbraculum, an African medicinal plant were recently found to be the XO inhibitors.
|Table 5: The interaction energy analysis of 11 ligands with that of XO using Discovery Studio® 3.1|
Click here to view
NOS are a family of enzymes that catalyze the production of NO from L-arginine. There are three isoforms of NOS in mammals, of which, two are constitutive and one is an inducible type. NO is an important cellular signaling molecule playing a vital role in various cellular processes. [Table 6] shows the docking studies and binding free energy calculations in which the highest interaction energy (−61.31 kcal/mol) was exhibited by shaftoside with that of NOS. Clinacoside B, on the contrary, showed the least interaction energy (−38.67 kcal/mol). Possibly due to the general poor binding phenomenon, lupeol and β-sitosterol could not dock with NOS. In this study, six ligands namely clinacoside B, clinacoside C, shaftoside, vitexin, orientin, and isovitexin exhibited interaction with the Glu377 amino acid residue of NOS as shown in [Table 6]. Lin et al. reported that isovitexin suppressed lipopolysaccharide mediated NOS in mouse macrophage cells.
|Table 6: The interaction energy analysis of 11 ligands with that of NOS using Discovery Studio® 3.1|
Click here to view
As for the docking studies and binding free energy calculations with HNE, shaftoside and clinacoside B exhibited the highest (−43.89 kcal/mol) and the least interaction energy (−30.15 kcal/mol), respectively, compared to other ligands [Table 7]. In this study, six ligands namely clinacoside A, clinacoside B, cycloclinacoside A1, orientin, isovitexin, and lupeol exhibited interaction with the Ser195 amino acid residue of HNE as shown in [Table 7]. Lupeol as HNE inhibitor was also reported by Mitaine-offer et al. Wanikiat et al. reported that acetone extract of C. nutans exhibited dose-dependent inhibition of N-formyl-methionyl-leucyl-phenylalanine induced HNE release.
|Table 7: The interaction energy analysis of 11 ligands with that of HNE using Discovery Studio® 3.1|
Click here to view
The maximum interaction energy in the docking studies and binding free energy calculations with that of MMP 2 was exhibited by isovitexin (−58.59 kcal/mol), whereas clinacoside B showed the least interaction energy of − 46.90 kcal/mol. In contrast, the other six ligands; clinacoside A, shaftoside, vitexin, orientin, lupeol and β-sitosterol, could not dock with MMP 2 which might be due to the general poor binding phenomenon. Four ligands namely clinacoside B, clinacoside C, isovitexin, and isoorientin exhibited interaction with the Glu202 amino acid residue of MMP 2 as shown in [Table 8].
|Table 8: The interaction energy analysis of 11 ligands with that of MMP 2 using Discovery Studio® 3.1|
Click here to view
MMP 9 is another target protein which its docking studies and binding free energy calculations showed isoorientin having the maximum interaction energy of − 60.45 kcal/mol. Clinacoside B, on the other hand, showed the very least interaction energy (−43.57 kcal/mol). Four of other ligands which are shaftoside, vitexin, orientin, and lupeol could not dock with MMP 9. This might be due to the general poor binding phenomenon as suggested by Akdogan et al. Among the ligands, clinacoside A, cycloclinacoside A1, isovitexin, and isoorientin exhibited interaction with Leu188 amino acid residue of MMP 9 as shown in [Table 9]. Although C. nutans has been known to have anti-inflammatory activity, until at the present, there is no available reported investigation for its MMP 2 and 9 inhibitory activity.
|Table 9: The interaction energy analysis of 11 ligands with that of MMP 9 using Discovery Studio® 3.1|
Click here to view
Finally the docking studies and binding free energy calculations with that of SQS in which isoorientin exhibited the largest interaction energy (−49.07 kcal/mol) while cycloclinacoside A1 showed smallest interaction energy (−34.16 kcal/mol) [Table 10]. Shaftoside, vitexin, orientin, lupeol, and β- sitosterol were the ligands which could not dock with SQS, which were suggested to be due to the general poor binding phenomenon. Other three ligands namely clinacoside A–C exhibited interaction with Asp80 amino acid, the residue of SQS. To our best of knowledge, there is no report yet available with regard to these ligands' SQS inhibitory activity.
|Table 10: The interaction energy analysis of 11 ligands with that of SQS using Discovery Studio® 3.1|
Click here to view
| Conclusion|| |
In the present study, it was found that isovitexin and isoorientin has the potential to dock and bind with all of the six targeted enzymes, whereas vitexin and orientin failed to dock and bind with four enzymes except NOS and HNE. As for β-sitosterol, it could dock and bind with only HNE and MMP 9. Interestingly, all of the selected ligands from C. nutans showed the potential to dock and bind with HNE. Hence, it is strongly suggested that the results of this present study has paved better understanding of these 11 ligands of C. nutans as potential XO, NOS, HNE, MMP 2, MMP 9, and SQS inhibitors in relation to the prevention of associated disorders of hyperuricemia, wound healing, and hyperlipidemia.
The author (R.N) would like to thank the Research Management Center of Universiti Putra Malaysia) for the Post-Doctoral financial support.
Financial support and sponsorship
Financial supported by Research Management Center, Universiti Putra Malaysia, Malaysia.
Conflicts of interest
There are no conflicts of interest.
| References|| |
Shim S, Ismail A, Khoo B. Perspective and insight on Clinacanthus nutans
Lindau in traditional medicine. Int J Integr Biol 2013;14:7-9.
Sakdarat S, Shuyprom A, Pientong C, Ekalaksananan T, Thongchai S. Bioactive constituents from the leaves of Clinacanthus nutans
Lindau. Bioorg Med Chem 2009;17:1857-60.
Cheeptham N, Towers GH. Light-mediated activities of some Thai medicinal plant teas. Fitoterapia 2002;73:651-62.
Roosita K, Kusharto CM, Sekiyama M, Fachrurozi Y, Ohtsuka R. Medicinal plants used by the villagers of a Sundanese community in West Java, Indonesia. J Ethnopharmacol 2008;115:72-81.
Ho SY, Titew WP, Priya M, Mohamed SA, Gabriel AA. Phytochemical analysis and antibacterial activity of methanolic extract of Clinacanthus nutans
leaf. Int J Drug Dev Res 2013;5:349-55.
Teshima KI, Kaneko T, Ohtani K, Kasai R, Lhieochaiphant S, Picheansoonthon C, et al
. Sulfur-containing glucosides from Clinacanthus nutans
. Phytochemistry 1998;48:831-5.
Dampawan P, Huntrakul C, Reutrakul V, Raston CL, White AH. Constituents of Clinacanthus nutans
and the crystal structure of Lup-20 (29)-ene-3-one. J Sci Soc Thailand 1977;3:14-26.
Wu G, Robertson DH, Brooks CL 3rd
, Vieth M. Detailed analysis of grid-based molecular docking: A case study of CDOCKER-A CHARMm-based MD docking algorithm. J Comput Chem 2003;24:1549-62.
Bhikshapathi DV, Shravan Kumar Y, Madhusudan Rao Y, Kishan V. Borrelidin: A prospective drug. Indian J Biotechnol 2010;9:18-23.
Radhakrishnan N, Lam KW, Abas F, Ismail IS. Molecular docking of curcumin analogues as human neutrophil elastase (HNE) inhibitors. Banglad J Pharmacol 2014;9:77-82.
Verma A. Lead finding from Phyllanthus debelis
with hepatoprotective potentials. Asian Pac J Trop Biomed 2012;2:S1735-7.
Segall MD, Beresford AP, Gola JM, Hawksley D, Tarbit MH. Focus on success: Using a probabilistic approach to achieve an optimal balance of compound properties in drug discovery. Expert Opin Drug Metab Toxicol 2006;2:325-37.
Darvas F, Keseru G, Papp A, Dormán G, Urge L, Krajcsi P. In silico
and ex silico
ADME approaches for drug discovery. Curr Top Med Chem 2002;2:1287-304.
Hodgson J. ADMET – turning chemicals into drugs. Nat Biotechnol 2001;19:722-6.
Borges F, Fernandes E, Roleira F. Progress towards the discovery of xanthine oxidase inhibitors. Curr Med Chem 2002;9:195-217.
Khanna S, Burudkar S, Bajaj K, Shah P, Keche A, Ghosh U, et al.
Isocytosine-based inhibitors of xanthine oxidase: Design, synthesis, SAR, PK and in vivo
efficacy in rat model of hyperuricemia. Bioorg Med Chem Lett 2012;22:7543-6.
Akdogan ED, Erman B, Yelekci K. In silico
design of novel and highly selective lysine-specific histone demethylase inhibitors. Turk J Chem 2011;35:523-42.
Sathisha KR, Khanum SA, Chandra JN, Ayisha F, Balaji S, Marathe GK, et al.
Synthesis and xanthine oxidase inhibitory activity of 7-methyl-2-(phenoxymethyl)-5H-[1,3,4]thiadiazolo[3,2-a] pyrimidin-5-one derivatives. Bioorg Med Chem 2011;19:211-20.
Lin CM, Chen CS, Chen CT, Liang YC, Lin JK. Molecular modeling of flavonoids that inhibits xanthine oxidase. Biochem Biophys Res Commun 2002;294:167-72.
Pham AT, Nguyen C, Malterud KE, Diallo D, Wangensteen H. Bioactive flavone-C-glycosides of the African medicinal plant Biophytum umbraculum
. Molecules 2013;18:10312-9.
Singh SP, Konwar BK. Molecular docking studies of quercetin and its analogues against human inducible nitric oxide synthase. Springerplus 2012;1:69.
Lin CM, Huang ST, Liang YC, Lin MS, Shih CM, Chang YC, et al.
Isovitexin suppresses lipopolysaccharide-mediated inducible nitric oxide synthase through inhibition of NF-kappa B in mouse macrophages. Planta Med 2005;71:748-53.
Mitaine-Offer AC, Hornebeck W, Sauvain M, Zèches-Hanrot M. Triterpenes and phytosterols as human leucocyte elastase inhibitors. Planta Med 2002;68:930-2.
Wanikiat P, Panthong A, Sujayanon P, Yoosook C, Rossi AG, Reutrakul V. The anti-inflammatory effects and the inhibition of neutrophil responsiveness by Barleria lupulina
and Clinacanthus nutans
extracts. J Ethnopharmacol 2008;116:234-44.
| Authors|| |
Intan Safinar Ismail, is an Associate Professor at the Department of Chemistry, Universiti Putra Malaysia (UPM), Malaysia since 2005. At the Institute of Biosciences (IBS), UPM, she is also heading the Laboratory of Natural Products, where she actively engaged in research activities. Her research interest is to hunt molecules from natural source having biologically active properties by using metabolomics/reductive approach.
[Table 1], [Table 2], [Table 3], [Table 4], [Table 5], [Table 6], [Table 7], [Table 8], [Table 9], [Table 10]