Abstract
Background
Inflammation-provoked disorders including cancer are arbitrated by cyclooxygenase-2 (COX-2). Celecoxib and niflumic acid are among the potent and selective inhibitors of this enzyme while aspirin (acetylsalicylic acid) and sodium salicylate are its non-selective and lesser potent inhibitors. Despite these proven studies, the comparative structural study of these selective and non-selective molecules at atomistic scale in complex state with COX-2 that may answer this differential inhibitory behavior has not been accomplished spotlighting the imperative need of additional research in this area. Thus, this study was framed to provide a strong explanation for the enigma of higher inhibitory activity of celecoxib-niflumic acid duo in comparison to aspirin and sodium salicylate towards COX-2.
Methods
A contemporary approach including advanced molecular docking against COX2, molecular dynamics of receptor-ligand complexes, simulation-trajectory-backed MMGBSA for different time points, radius of gyration (Rg) calculations, and e-pharmacophores approach was employed to attain a rational conclusion.
Results
Our findings demonstrated the higher binding affinity of celecoxib and niflumic acid over aspirin and sodium salicylate against COX-2. Although both selective and non-selective COX-2 inhibitors manifested nearly the same stability in the active site of this enzyme but the e-pharmocophoric features found in the case of selective inhibitors scored over non-selective ones. Thus, our findings excluded the differential stability to be the cause of stronger potency of selective inhibitors but attributed their potency to greater number of complementary features present in these inhibitors against the active site of inflammation engendering COX-2.
Similar content being viewed by others
Explore related subjects
Discover the latest articles and news from researchers in related subjects, suggested using machine learning.Data availability
All the data has been submitted as figures and tables with the manuscript.
References
Knights KM, Mangoni AA, Miners JO (2010) Defining the COX inhibitor selectivity of NSAIDs: implications for understanding toxicity. Expert Rev Clin Pharmacol 3:769–776
Ulrich CM, Whitton J, Yu JH, Sibert J, Sparks R, Potter JD et al (2005) PTGS2 (COX-2) -765G > C promoter variant reduces risk of colorectal adenoma among nonusers of nonsteroidal anti-inflammatory drugs. Cancer Epidemiol Biomark Prev Publ Am Assoc Cancer Res Cosponsor Am Soc Prev Oncol 14:616–619
Granström E (1984) The arachidonic acid cascade. The prostaglandins, thromboxanes and leukotrienes. Inflammation, 8 Suppl, S15–25
Williams CS, Mann M, DuBois RN (1999) The role of cyclooxygenases in inflammation, cancer, and development. Oncogene 18:7908–7916
Bhardwaj A, Kaur J, Wuest M, Wuest F (2017) In situ click chemistry generation of cyclooxygenase-2 inhibitors. Nat Commun 8:1
Zidar N, Odar K, Glavac D, Jerse M, Zupanc T, Stajer D (2009) Cyclooxygenase in normal human tissues–is COX-1 really a constitutive isoform, and COX-2 an inducible isoform? J Cell Mol Med 13:3753–3763
Minghetti L (2004) Cyclooxygenase-2 (COX-2) in inflammatory and degenerative brain diseases. J Neuropathol Exp Neurol 63:901–910
Sui H, Zhou S, Wang Y, Liu X, Zhou L, Yin P et al (2011) COX-2 contributes to P-glycoprotein-mediated multidrug resistance via phosphorylation of c-Jun at Ser63/73 in colorectal cancer. Carcinogenesis 32:667–675
Pang LY, Hurst EA, Argyle DJ (2016) Cyclooxygenase-2: a role in cancer stem cell survival and repopulation of cancer cells during therapy. Stem Cells Int 2016:2048731
Minghetti L (2007) Role of COX-2 in inflammatory and degenerative brain diseases. Subcell Biochem 42:127–141
Zhu X-T, Chen L, Lin J-H (2018) Selective COX-2 inhibitor versus non-selective COX-2 inhibitor for the prevention of heterotopic ossification after total hip arthroplasty: a meta-analysis. Med (Baltimore), 97
Hsu AL, Ching TT, Wang DS, Song X, Rangnekar VM, Chen CS (2000) The cyclooxygenase-2 inhibitor celecoxib induces apoptosis by blocking Akt activation in human prostate cancer cells independently of Bcl-2. J Biol Chem 275:11397–11403
Setiawati A, Setiawati A (2016) Celecoxib, a COX-2 Selective inhibitor, induces cell cycle arrest at the G2/M phase in hela cervical cancer cells. Asian Pac J Cancer Prev APJCP 17:1655–1660
Xie H, Gao L, Chai N, Song J, Wang J, Song Z et al (2009) Potent cell growth inhibitory effects in hepatitis B virus X protein positive hepatocellular carcinoma cells by the selective cyclooxygenase-2 inhibitor celecoxib. Mol Carcinog 48:56–65
Wang G, Li J, Zhang L, Huang S, Zhao X, Zhao X (2017) Celecoxib induced apoptosis against different breast cancer cell lines by down-regulated NF-κB pathway. Biochem Biophys Res Commun 490:969–976
Zakrocka I, Targowska-Duda KM, Wnorowski A, Kocki T, Jóźwiak K, Turski WA (2019) Influence of cyclooxygenase-2 inhibitors on kynurenic acid production in rat brain in vitro. Neurotox Res 35:244–254
Warner TD, Giuliano F, Vojnovic I, Bukasa A, Mitchell JA, Vane JR (1999) Nonsteroid drug selectivities for cyclo-oxygenase-1 rather than cyclo-oxygenase-2 are associated with human gastrointestinal toxicity: a full in vitro analysis. Proc Natl Acad Sci USA 96:7563–7568
Amann R, Peskar BA (2002) Anti-inflammatory effects of aspirin and sodium salicylate. Eur J Pharmacol 447:1–9
Bao W, Luo Y, Wang D, Li J, Wu X, Mei W (2018) Sodium salicylate modulates inflammatory responses through AMP-activated protein kinase activation in LPS-stimulated THP-1 cells. J Cell Biochem 119:850–860
Čeponytė U, Paškevičiūtė M, Petrikaitė V (2018) Comparison of NSAIDs activity in COX-2 expressing and non-expressing 2D and 3D pancreatic cancer cell cultures. Cancer Manag Res 10:1543–1551
Orlando BJ, Malkowski MG (2016) Substrate-selective inhibition of cyclooxygeanse-2 by fenamic acid derivatives is dependent on peroxide tone. J Biol Chem 291:15069–15081
Ganai SA (2021) Characterizing binding intensity and energetic features of histone deacetylase inhibitor pracinostat towards class I HDAC isozymes through futuristic drug designing strategy. Silico Pharmacol 9:18
Ganai SA, Srinivasan P, Rajamanikandan S, Shah BA, Mohan S, Gani M et al (2021) Delineating binding potential, stability of Sulforaphane-N-acetyl-cysteine in the active site of histone deacetylase 2 and testing its cytotoxicity against distinct cancer lines through stringent molecular dynamics, DFT and cell-based assays. Chem Biol Drug Des 98:363–376
Ganai SA, Mir MA, Shah BA, Qadri RA, Wani AH, Rajamanikandan S et al (2023) Evaluation of free radical quenching, anti-inflammatory activity together with anticancer potential of Lychnis coronaria and characterization of novel molecules from its extract through high resolution-liquid chromatography mass spectrometry coupled to structural biochemistry approach. J Biomol Structure Dyn 1–15
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S et al (2021) PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res 49:D1388–D1395
Kim S, Chen J, Cheng T, Gindulyte A, He J, He S et al (2020) PubChem in 2021: new data content and improved web interfaces. Nucleic Acids Res 49:D1388–D1395
LigPrep S (2021) LLC. Schrödinger Release 2021-1: LigPrep, Schrödinger, LLC, New York, NY, 2021. Schrödinger, LLC, New York
Madhavi Sastry G, Adzhigirey M, Day T, Annabhimoju R, Sherman W (2013) Protein and ligand preparation: parameters, protocols, and influence on virtual screening enrichments. J Comput Aided Mol Des 27:221–234
Farooq S, Ganai SA, Ganai BA, Mohan S, Uqab B, Nazir R (2021) Molecular characterization of lipase from a psychrotrophic bacterium Pseudomonas sp. CRBC14. Current genetics
Friesner RA, Murphy RB, Repasky MP, Frye LL, Greenwood JR, Halgren TA et al (2006) Extra precision glide: docking and scoring incorporating a model of hydrophobic enclosure for protein−ligand complexes. J Med Chem 49:6177–6196
Scott WRP, Hünenberger PH, Tironi IG, Mark AE, Billeter SR, Fennen J et al (1999) The GROMOS biomolecular simulation program package. J Phys Chem A 103:3596–3607
Schüttelkopf AW, van Aalten DM (2004) PRODRG: a tool for high-throughput crystallography of protein-ligand complexes. Acta crystallographica. Section D, Biological crystallography 60:1355–63
Hollingsworth SA, Dror RO (2018) Molecular dynamics simulation for all. Neuron 99:1129–1143
Berendsen HJC, Grigera JR, Straatsma TP (1987) The missing term in effective pair potentials. J Phys Chem 91:6269–6271
Hess B, Bekker H, Berendsen HJC, Fraaije JGEM (1997) LINCS: a linear constraint solver for molecular simulations. J Comput Chem 18:1463–1472
Durham E, Dorr B, Woetzel N, Staritzbichler R, Meiler J (2009) Solvent accessible surface area approximations for rapid and accurate protein structure prediction. J Mol Model 15:1093–1108
Rashid HU, Ahmad N, Abdalla M, Khan K, Martines MAU, Shabana S (2022) Molecular docking and dynamic simulations of cefixime, etoposide and nebrodenside A against the pathogenic proteins of SARS-CoV-2. J Mol Struct 1247:131296
Yoshino R, Yasuo N, Sekijima M (2019) Molecular dynamics simulation reveals the mechanism by which the influenza cap-dependent endonuclease acquires resistance against Baloxavir marboxil. Sci Rep 9:17464
Li J, Abel R, Zhu K, Cao Y, Zhao S, Friesner RA (2011) The VSGB 2.0 model: a next generation energy model for high resolution protein structure modeling. Proteins 79:2794–812
Salam NK, Nuti R, Sherman W (2009) Novel method for generating structure-based pharmacophores using energetic analysis. J Chem Inf Model 49:2356–2368
Ganai SA, Abdullah E, Rashid R, Altaf M (2017) Combinatorial in silico strategy towards identifying potential hotspots during inhibition of structurally identical HDAC1 and HDAC2 enzymes for effective chemotherapy against neurological disorders. Front Mol Neurosci 10:357
Hughes TS, Chalmers MJ, Novick S, Kuruvilla DS, Chang MR, Kamenecka TM et al (2012) Ligand and receptor dynamics contribute to the mechanism of graded PPARγ agonism. Structure (London, England : 1993). 20:139–50
Rai H, Barik A, Singh YP, Suresh A, Singh L, Singh G et al (2021) Molecular docking, binding mode analysis, molecular dynamics, and prediction of ADMET/toxicity properties of selective potential antiviral agents against SARS-CoV-2 main protease: an effort toward drug repurposing to combat COVID-19. Mol Divers
Baig MH, Sudhakar DR, Kalaiarasan P, Subbarao N, Wadhawa G, Lohani M et al (2014) Insight into the effect of inhibitor resistant S130G mutant on physico-chemical properties of SHV type beta-lactamase: a molecular dynamics study. PLoS One 9:e112456
Khan MT, Khan A, Rehman AU, Wang Y, Akhtar K, Malik SI et al (2019) Structural and free energy landscape of novel mutations in ribosomal protein S1 (rpsA) associated with pyrazinamide resistance. Sci Rep 9:7482
Ferreira de Freitas R, Schapira M (2017) A systematic analysis of atomic protein-ligand interactions in the PDB. MedChemComm 8:1970–81
Nicholls A, McGaughey GB, Sheridan RP, Good AC, Warren G, Mathieu M et al (2010) Molecular shape and medicinal chemistry: a perspective. J Med Chem 53:3862–3886
Borngräber S, Browner M, Gillmor S, Gerth C, Anton M, Fletterick R et al (1999) Shape and specificity in mammalian 15-lipoxygenase active site. The functional interplay of sequence determinants for the reaction specificity. J Biol Chem 274:37345–50
Kurumbail RG, Stevens AM, Gierse JK, McDonald JJ, Stegeman RA, Pak JY et al (1996) Structural basis for selective inhibition of cyclooxygenase-2 by anti-inflammatory agents. Nature 384:644–648
Carlson HA, Smith RD, Khazanov NA, Kirchhoff PD, Dunbar JB Jr, Benson ML (2008) Differences between high- and low-affinity complexes of enzymes and nonenzymes. J Med Chem 51:6432–6441
Acknowledgements
Rajamanikandan thanks Centre for Drug Discovery, Karpagam Academy of Higher Education, and Coimbatore for Molecular Dynamics. Malik FA sincerely thanks Rangaswamy R for providing the trial version of Schrödinger tool.
Author information
Authors and Affiliations
Contributions
Ganai SA designed the work. Ganai SA and Rajamanikandan did majority of the work. Shah BA, Lone A, Arwa F, and Malik FA also contributed significantly. Ganai SA and Rajamanikandan wrote the manuscript.
Corresponding author
Ethics declarations
Ethical approval
Not applicable.
Consent to participate
Not applicable to this study.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Springer Nature or its licensor (e.g. a society or other partner) holds exclusive rights to this article under a publishing agreement with the author(s) or other rightsholder(s); author self-archiving of the accepted manuscript version of this article is solely governed by the terms of such publishing agreement and applicable law.
About this article
Cite this article
Ganai, S.A., Rajamanikandan, S., Shah, B.A. et al. Comparative structural study of selective and non-selective NSAIDs against the enzyme cyclooxygenase-2 through real-time molecular dynamics linked to post-dynamics MM-GBSA and e-pharmacophores mapping. J Mol Model 29, 192 (2023). https://doi.org/10.1007/s00894-023-05603-7
Received:
Accepted:
Published:
DOI: https://doi.org/10.1007/s00894-023-05603-7