Discovery of andrographolide hit analog as a potent cyclooxygenase

Scientific Reports volume 13, Article number: 8147 (2023) Cite this article

433 Accesses

1 Altmetric

Metrics details

Cyclooxygenase-2 (COX-2) is the key enzyme responsible for the conversion of arachidonic acid to prostaglandins that display pro-inflammatory properties and thus, it is a potential target protein to develop anti-inflammatory drugs. In this study, chemical and bio-informatics approaches have been employed to find a novel potent andrographolide (AGP) analog as a COX-2 inhibitor having better pharmacological properties than aspirin and rofecoxib (controls). The full amino acid sequenced human Alpha fold (AF) COX-2 protein (604AA) was selected and validated for its accuracy against the reported COX-2 protein structures (PDB ID: 5F19, 5KIR, 5F1A, 5IKQ and 1V0X) followed by multiple sequence alignment analysis to establish the sequence conservation. The systematic virtual screening of 237 AGP analogs against AF-COX-2 protein yielded 22 lead compounds based on the binding energy score (< − 8.0 kcal/mol). These were further screened out to 7 analogs by molecular docking analysis and investigated further for ADMET prediction, ligand efficiency metrics calculations, quantum mechanical analysis, MD simulation, electrostatic potential energy (EPE) docking simulation, and MM/GBSA. In-depth analysis revealed that AGP analog A3 (3-[2-[(1R,4aR,5R,6R,8aR)-6-hydroxy-5,6,8a-trimethyl-2-methylidene-3,4,4a,5,7,8-hexahydro-1H-naphthalen-1-yl]ethylidene]-4-hydroxyoxolan-2-one) forms the most stable complex with the AF-COX-2 showing the least RMSD value (0.37 ± 0.03 nm), a good number of hydrogen bonds (protein–ligand H-bond = 11, and protein H-bond = 525), minimum EPE score (− 53.81 kcal/mol), and lowest MM-GBSA before and after simulation (− 55.37 and − 56.25 kcal/mol, respectively) value compared to other analogs and controls. Thus, we suggest that the identified A3 AGP analog could be developed as a promising plant-based anti-inflammatory drug by inhibiting COX-2.

Cyclooxygenase (COX or prostaglandin G/H synthase) plays a vital role in the generation of biological mediators such as prostaglandins (PGs), thromboxane, and prostacylins from arachidonic acid. The COX enzyme exists as two isoforms: (i) the COX-1 isozyme is constitutively active and primarily involved in homeostatic functions like the regulation of platelet aggregation and gastric acidity; (ii) in contrast, the COX-2 isomeric form induces during pathological conditions such as inflammation, pain, and fever1,2,3. To develop anti-inflammatory and analgesic drugs, COX-2 is a viable and specific target and various COX-2 inhibitory drugs have been developed in the past few decades4. The COX-2 comprises 604 amino acid sequence length with three major domains: (i) the epidermal growth factor (EGF) domain (34–72 amino acid), (ii) the membrane binding domain (73–116 amino acid), and (iii) the catalytic domain which contains the COX and peroxidase active sites [UniProt- P35354]. The COX and peroxidase active sites are spatially distinct but functionally linked5. The proximity to Tyr371 and Ser516 are the critical catalytic amino acids for the COX-2 active site6. The valine residue at the 509 position in the COX-2 forms a hydrophobic (HYD) pocket at its active site and plays a key role in the selective induction of the enzyme and thus inflammatory response7,8.

The COX-2 catalyzes the reaction in two steps: (i) the first step is the COX reaction in which arachidonate is converted to PG-G2 which occurs in the HYD channel of the core of the protein and (ii) the second step is the peroxidase reaction wherein PG-G2 is reduced to PG-H2 (an important precursor of prostacyclin and expressed during inflammation) that takes place at the heme-containing active site located near the protein surface9. Non-steroidal anti-inflammatory drugs (NSAIDs), especially coxibs, and aspirin, reduce inflammation by inhibiting the PG synthesis via covalent interaction with Ser-516 residue at the COX-2 active site10,11. Although these drugs are safe and effective, they can cause gastrointestinal toxicity, which has limited their use in highly sensitive patients. Designing of selective COX-2 inhibitors including celecoxib and rofecoxib (coxibs), inhibit COX-2 instead of COX-1 in order to provide pain relief and anti-inflammatory properties and prevent the gastric discomforts experienced with non-selective NSAIDs12,13,14. Rofecoxib is a highly selective NSAID due to the absence of carboxylic acid (thus less GI irritant) and the presence of two large aromatic cycles linked to a central heterocycle ring15. Beside the selectivity and efficacy, the high cost of the drug, potential adverse cardiovascular effects, and the increased risk of ischemic stroke, use of COX-2 inhibitors is controversial16,17,18. Hence, there is an urgent demand for a therapeutic drug against COX-2 with minimal or no side effects preferably from the natural origin19.

Secondary metabolites such as alkaloids, terpenoids, stilbenes, flavonoids, saponins, and fatty acids have been demonstrated to have COX-2 inhibitory activity20. Amongst these, andrographolide (AGP), which is obtained from the Andrographis paniculata plant, has been proven as an important natural anti-inflammatory agent along with various other pharmacological benefits21,22,23. Andrographolide is a natural electrophilic covalent inhibitor, derived from the terpenoid pathway and privileged with the structural features which are necessary for binding with multiple targets24,25,26. Structurally, AGP contains a complex 3D architecture with a high number of sp3-hybridized carbons, rich oxygen content, low nitrogen content, a polycyclic structure with aliphatic side chains, and no aromatic rings. Moreover, the chemical properties of andrographolide are strengthened by α-Alkylidene-β-hydroxy-γ-butyrolactone moiety which is a specific electrophilic fragment that enables it to irreversibly and covalently modify the target protein with thiol-containing amino acid residue25. This provides substantially to the binding affinity with the target protein and presents andrographolide as an efficient covalent inhibitor26,27.

A few research groups have recently attempted to understand the andrographolide signaling pathway targets which are directly or indirectly involved in the inflammatory reaction28. Tran et al. (2020) reported the anti-inflammatory activity of the andrographolide and its few derivatives in a lipopolysaccharide-induced acute lung injury model and emphasized that α–Alkylidene-β-hydroxy-γ-butyrolactone and ent-labdane moieties are vital for anti-inflammatory activity25. Dai et al. (2011) reported that AGP exerts its enhanced anti-inflammatory effect by decreasing serum inducible nitric oxide synthase (iNOS) activity, NO production, and PGE2 production29. Recently, Wang et al. (2019) reported that AGP and its derivative effectively inhibit the lipopolysaccharide-induced expression of COX-2 in murine macrophages30. Similarly, other research groups have reported that the andrographolide demonstrates anticancer activity by inhibiting the expression of COX-2 protein31,32. Although efforts have been paid to prove the anti-inflammatory activity of andrographolide, the pharmacokinetic parameters, i.e. poor solubility, and low bioavailability are the limiting factors for its successful drug development33.

In this study, systematic computational analysis is carried out to find out a novel natural plant-based andrographolide analog as an anti-inflammatory agent against COX-2 enzyme which has better potency, efficacy, and selectivity than aspirin and rofecoxib (controls). Molecular mechanics techniques of cheminformatics and quantum mechanical approaches are employed for detailed investigation. The full sequence length of the human AF-COX-2 protein structure is used after the validation against the available 3D crystal structure. Further, multiple sequence alignment was performed for the amino acid sequence length of the protein to find out conserved regions for biological activity. A total of 237 AGP analogs are used for screening followed by docking studies for the hit analogs, ADMET parameters prediction, ligand efficiency metrics calculations, in-depth quantum mechanical analysis, and MM/GBSA rescoring studies. Further, the stability of the hit AGP analogs complex was assessed by the molecular docking simulation and electrostatic potential energy (EPE) docking simulation analysis. Outcomes of the present study help in the development of plant-based COX-2 inhibitor as a future anti-inflammatory candidates.

The 3D structure of crystallized COX-2 (5F19, 5KIR, 5F1A, 5IKQ, and 1V0X) protein was retrieved from the PDB database while the predicted structure of human COX-2 protein, containing 604 amino acids, was obtained from the AlphaFold (AF) (ID: AF-P35354-F1)34,35,36,37. The protein structure was validated by SAVES v6.0 (https://saves.mbi.ucla.edu/), ProSA, and ProQ38,39,40,41,42. For the accuracy of the AF-COX-2 human protein structure, the validation results were compared with other available crystallized COX-2 PDB structures, i.e. 5F19, 5KIR, 5F1A and 5IKQ, and predicted COX-2 model (1V0X). Before any interaction analysis, the AF-COX-2 protein structure was refined by 3Drefine (https://3drefine.mu.hekademeia.org/)43. Subsequently, possible binding pockets for the AF-COX-2 protein were predicted by CASTp 3.0 and visualized by Chimera44,45.

To verify the presence of the conserved amino acids in the AF-COX-2 protein, multiple sequence alignment (MSA) was performed using the PRALINE program46. For the same, 10 different mammal species namely Homo sapiens, Rattus norvegicus, Mus musculus, Ovis aries, Bos Taurus, Oryctolagus cuniculus, Cavia porcellus, Equus caballus, Neovison vison, and Gallus gallus were selected and amino acid sequences of their respective COX-2 protein were retrieved from Uniprot (http://www.uniprot.org/). The resulting aligned amino acid sequences were colored based on the conservation index (0 to10), wherein the highly conserved sequences of amino acid amongst species are responsible for the particular biological function of the protein47.

A total of 237 AGP analogs were taken from the PubChem database that consists of the common α,β-unsaturated γ-lactone moiety which is responsible for the biological activity by forming adduct with the residues during the biological interaction25,26. Then, the analogs were converted into a autodock PDBQT format for virtual screening by Open Babel48,49. PyRx (AutoDock Vina) software was used for the site-directed virtual screening of the AGP analogs against the AF-COX-2 protein50. Further, based on the interacting amino acid residues present for the co-crystal aspirin and rofecoxib in the 5F19 and 5KIR PDB structure respectively, amino acids namely His75, Arg106, Val330, Tyr334, Val335, Tyr341, Tyr371, Trp373, Arg499, Phe504, Glu510, Ser516, and Leu517 were chosen as active site residues in the AF-COX-2 protein for setting up the grid for the virtual screening35,36. The grid box was set with search space center 4.76, 5.07, and − 0.12 at x, y, and z, respectively, with 8 exhaustiveness, while the dimensions (Å) were set for the grid at x, y, and z with a value of 29.54, 29.74, and 27.37 each. For shortlisting the hit AGP analogs, the cutoff value for binding energy (BE) was set to be lesser than − 8.00 kcal/mol.

The hit AGP analogs obtained from virtual screening were submitted for molecular docking analysis along with AGP, aspirin, and rofecoxib as controls. AutoDock4.2 was utilized for the site-directed docking of these compounds with AF-COX-251. Since the AF protein structure of COX-2 is a predicted protein structure, co-crystal ligands, i.e. aspirin and rofecoxib of the existing crystallized structure of COX-2 protein 5F19 and 5KIR, respectively, were taken as reference compounds for the docking analysis. Further, amino acids i.e. His75, Arg106, Val330, Tyr334, Val335, Tyr341, Tyr371, Trp373, Arg499, Phe504, Glu510, Ser516, and Leu517 were chosen as the binding site residues based on the aspirin and rofecoxib interacting residues in the 5F19 and 5KIR PDB structures. The gird box for docking was generated with 0.375 Å spacing and dimension along the x, y, and z coordinates with the values of 84, 88, and 94, respectively. The grid center was set to the centroid of the binding site residues with the dimension of 4.305, 4.394, and 1.611, respectively, alongside the x, y, and z–Axis. Further, GA (genetic algorithm) was employed to calculate the docking parameters with 25000000 number of energy evaluations, 27000 number of generations, and 2.0 Å RMSC (root mean square cluster) tolerance. Lamarckian genetic algorithm was used to perform docking simulation. The AutoGrid and AutoDock were utilized to perform molecular docking and to investigate the best binding pose of the ligands in protein–ligand complexes. Finally, the ligand's conformers with the lowest free binding energy were selected for further analysis. The output files obtained from the docking study were visualized and analyzed by ICM-Browser and Accelrys discovery studio visualizer52.

ADMETlab 2.0 was used to predict physicochemical, medicinal, absorption, distribution, metabolism, excretion, and toxicological parameters for AGP and shortlisted AGP analogs along with aspirin and rofecoxib53. It was also used to establish the drug-likeness of the AGP and shortlisted analogs.

Ligand binding efficiency metrics were computed in terms of inhibition constant (Ki), ligand efficiency (LE), ligand lipophilic efficiency (LLE), ligand efficiency scaling function (LE_Scale), fit quality (FQ), and lipophilicity corrected ligand efficiency (LELP) by following the Eq. (1–6)54,55,56,57,58.

Spartan 20 (wavefun.com) package was used to calculate quantum–mechanical (QM) descriptors in terms of the highest occupied molecular orbital (HOMO), lowest unoccupied molecular orbital (LUMO), HOMO–LUMO gap (HLG), and EPE parameters for assessing the molecular orbital properties, electrostatic potential, and elucidating various types of interactions. All the chemical compounds were optimized to compute equilibrium geometric variables in the gaseous state at the ground employing density functional theory (DFT) with Becke's three-parameter exchange potential59. Further, Lee Yang Parr (LYP) correlation functional, which is an exchange functional combining of Becke with a 6-31G* basis set, was used to calculate the ground state geometries60,61,62.

The native AF-COX-2 protein and its complexes with the hit AGP analogs (screened from the docking studies), AGP, aspirin, and rofecoxib were used to commence the MD simulation studies. GROMACS 2019.2 was used to carry out the MD simulation63. All the topological parameters for chemical compounds were obtained from the PRODRG server64. GROMOS96 54a7 force field was used to perform all the simulation studies65. Further, an explicit simple point charge (SPC) water model was used to solvate the system in a dodecahedron box type. The system was neutralized by the addition of counter-ions to sustain the protonated residual state at physiological pH 7.4. Subsequently, the initial steric hindrance was eliminated by the steepest descent energy method for 50,000 steps to minimize the system. Thereafter, NVT/NPT ensembles with a leap-frog integrator at 300 K temperature, and 1.0 bar pressure were used to equilibrate the system. The MD simulation was run for 100 ns for each equilibrated system with 5000 frames. The statistical MD simulation analysis was carried out through WebGRO (https://simlab.uams.edu/) by calculating the RMSD (root-means-square-deviation), RMSF (root-mean-square-fluctuation), the average area per residue, SAS (solvent–Accessible-surface) area, volume, and density, H-bonds (number of hydrogen-bonds in protein as well as between protein and ligand), and Rg (radius of gyration) value.

To assess the stability of adducts of AGP, hit analogs, aspirin, and rofecoxib with the selected peptide length of AF-COX-2, the EPE docking simulation study was carried out using Spartan 20. For the same, amino acid sequence length from 511 to 520 of AF-COX-2 (i.e. Val511, Gly512, Ala513, Pro514, Phe515, Ser516, Leu517, Lys518, Gly519, and Leu520) was selected as a peptide for forming the adducts with the compounds. For the same, first a model consisting of 10 amino acid residue peptide and one ligand was built. This was followed by energy minimization with the molecular mechanics option followed by the B3LYP/6-31G* single-point calculations. The energy of adduct was obtained directly from the output of the computations. Further, to calculate the electronic parameters, the data related to B3LYP/6-31G* basis set was taken as the number of unpaired electrons (multiplicity) = 0, total charge = neutral coupling, and coupling constant = empirical and stabilization energy was determined by deducing the sum of the energy of the individual constituents from the adduct energy. Then, EPE was calculated for each adduct and peptide which characterized the electron distribution on the surface of compounds66. Finally, the binding energy of complex compounds with peptide was calculated using the following Eq. (7).

The binding free energies in terms of molecular mechanics generalized born surface area (MM/GBSA) of AGP's hit analogs, AGP, aspirin, and rofecoxib complexes were calculated by Fast Amber Rescoring (FAR)67. The ff14SB force field was used for the small molecules, whereas the Generalized Amber Force Field2 (GAFF2) was employed for protein analysis68. Further, the AM1-BCC method was used to assign the partial charge in the analogs through an antechamber module of Amber69,70.

The full sequenced COX-2 protein structure AF-P35354-F1, comprising 1-604 amino acids, was used in this study, which is modeled by AlphaFold and uses a state-of-the–Art machine learning method for protein modeling. Before performing any computational analysis, this protein structure was validated by comparing the overall quality of the protein with the existing four crystallized COX-2 protein structures (PDB ID: 5F19, 5KIR, 5F1A, and 5IKQ) and a modeled structure (PDB ID: 1V0X). Supplementary Table S1 summarizes the results obtained from SAVES, ProSA, and ProQ for all the protein structures, which use ERRAT, VERIFY, PROVE, Whatcheck, Procheck, Z score, and LGscore to assess the quality of the protein by calculating the overall quality factor, compatibility of an atomic 3D model with its amino acid sequence, volumes of atoms, stereochemical parameters, and stereochemical quality. The AF-COX-2 protein structure successfully passed the VERIFY and PROVE with 97.74% overall quality factor with 90.9% residues in the most favored reason. In addition, z-score and LGscore values for the AF-COX-2 protein were found to be -8.97, and 7.639, respectively which indicates that AF-COX-2 model protein structure is of good quality40,42. Further, the overall quality of the AF-COX-2 protein structure was found to be the best among all other structures in terms of quality, compatibility, and stereochemistry parameters, and hence this 3D protein structure was used for further studies.

CASTp 3.0 predicted 96 possible binding pockets for the AF-COX-2 protein (Supplementary Table S2). Then, according to the shape, internal cavities, area, and volume values for the binding pockets, the amino acids from the second binding pocket (His75, Arg106, Phe191, Val214, Val330, Ile331, Tyr334, Val335, Leu338, Tyr341, Leu345, Lys346, Gln358, Asn361, Tyr371, Trp373, Lys454, Arg455, Met457, Arg499, Phe504, Glu510, Phe515, Ser516, Leu517, Leu520, Met521) were selected for further consideration. Next, binding site residues for the study were selected based on the interaction of rofecoxib and aspirin present in 5KIR and 5F19 PDB structure, respectively. Further, in the 5KIR PDB structure the binding site residues for rofecoxib were found to be Val344, Trp387, Phe518, Tyr385, 355, Ser530, Leu531, Arg120, 513, Glu524, and His90. Whereas, in the 5F19 PDB structure the aspirin binding site residues recognized as Tyr385, Ser530, Leu531, Glu524, Arg120, 513, and His 90. Therefore, to find out the more efficient compound than rofecoxib and aspirin, we have considered the binding site residues of both the compounds. Hence, we have selected the common binding residues (Tyr385, Ser530, Leu531, Arg120, and Glu524) along with the unique binding residues present on the rofecoxib and aspirn. Accordingly, we have considered His75, Arg106, Val330, Tyr334, Val335, Tyr341, Tyr371, Trp373, Arg499, Phe504, Glu510, Ser516, and Leu517 residues as binding site in this study for AF-COX-2 protein. Interestingly, CASTp 3.0 analysis revealed that all the selected binding site residues are also present in the predicted binding pocket for the AF-COX-2 and also may be important for ligand binding and protein inhibition71.

The AF-COX-2 protein has a modeled 3D structure with no co-crystallized ligand, therefore, binding site residue prediction in the refined protein structure was carried out by comparing it with the active site of the existing crystallized PDB structure of COX-2 present in the PDB database. For the same, 5F19 and 5KIR PDB structures of COX-2 were selected which have aspirin and rofecoxib as co-crystal ligands, respectively. Based on the binding sites of both the co-crystal ligands in the protein structure, His75, Arg106, Val330, Tyr334, Tyr341, Tyr371, Trp373, Arg499, Phe504, Glu510, Ser516, and Leu517 amino acids were selected as the active sites for binding interaction studies. To verify these residues as the prominent active site, multiple sequence alignment was performed using the COX-2 protein of 10 different mammalian species, which proved that except for Arg499, all the selected amino acids are fully conserved (Supplementary Fig. S1 and Table S3). This also signifies the importance of these amino acids in various biological functions and as a potential target for drug development47.

A set of 237 andrographolide analogs that consists of the common α,β-unsaturated γ-lactone moiety which is important for biological activity were screened against the full sequence validated AF-COX-2 protein to find out a potent plant based anti-inflammatory compound25,26. The site directed virtual screening was utilized to discover the hit AGP analogs by computational approach. It calculates the binding energy for all the subjected analogs by the results of protein–ligand interaction. Subsequently, the energetically more favorable analogs’ complexes were selected for further analysis. The results from the virtual screening analysis revealed that the lowest binding energy was found to be − 8.80 kcal/mol for PubChem ID: 132210370, while the highest binding energy was determined to be − 5.90 kcal/mol for PubChem ID:59876522 (Supplementary Table S4). Thereafter, based on the binding energy cutoff value (− 8.00 kcal/mol), 22 AGP analogs out of 237 were screened for further analysis. Therefore, these 22 AGP analogs were subjected to molecular docking studies to find out the potent COX-2 inhibitor along with AGP, and control drugs (aspirin and rofecoxib).

The binding energy values obtained from the molecular docking analysis were found to be in the range of − 9.35 kcal/mol (PubChem ID: 132210508) to − 3.2 kcal/mol (PubChem ID: 132217538) (Supplementary Table S5). In contrast, aspirin, rofecoxib, and AGP showed binding energy values of − 5.61 kcal/mol, − 8.17 kcal/mol, and − 7.95 kcal/mol, respectively. Out of 22 AGP analogs, 7 hit analogs, named A1–A7, were selected for further analysis based on the binding energy cutoff value of lesser than − 8.00 kcal/mol. In terms of binding energy, all seven hit analogs showed lesser value compared to aspirin and AGP, while five analogs (A1–A5) exhibited lesser value than rofecoxib (Table 1), indicating their potential towards COX-2 protein inhibition.

The binding pattern of AGP and its selected seven hit analogs with respect to the existing COX-2 protein inhibitors, i.e., aspirin and rofecoxib (controls), was examined to compare the molecular binding interaction and binding pose analysis (Table 2). This analysis reveals that Val335, Leu338, Ser339, Phe504, Val509, and Ala513 amino acids are common in the interactions of all the subjected compounds with AF-COX-2 protein. In addition, Trp373, Gly512, and Ser516 amino acids are shared by aspirin, AGP, and A1–A7 analogs, whereas Leu78, Val102, Arg106, Tyr341, Leu345, Leu517 residues are occupied by rofecoxib, AGP, and A1–A7. Further, among all the interacting amino acids, Arg106, Tyr371, Arg499, Met508, Val509, Glu510, Ala513, and Ser516 residues were involved in H-bond interaction, while Leu78, Val102, Arg106, Val335, Leu338, Tyr341, Leu345, Phe504, Met508, Val509, Ala513, Leu517 residues were engaged in the HYD interaction. AGP and hit analogs A1–A5 showed 3, 3, 4, 1, 3, and 4 H-bonds, respectively, while no H-bond was formed in the case of A6 and A7 analog. Intriguingly, it is evident from the binding pattern that AGP and A1–A7 analogs share the interacting amino acid residues of both aspirin and rofecoxib binding site of AF-COX-2 protein (Figs. 1 and 2), indicating their potential as COX-2 protein inhibitor with improved binding interaction. Further, the binding orientation of AGP and hit analogs display that oxalone moiety is faced towards the aspirin binding site, while the naphthalene/naphthol moiety is positioned towards the rofecoxib binding site (Supplementary Fig. S2).

Superimposition of aspirin, rofecoxib, AGP, and hit AGP analogs at the binding site of AF-COX-2. All the chemical compounds are represented as ball and stick model, protein as ribbon representation, and interacting amino acids as single letter code. Color code: light grey = AF-COX-2 protein, blue dotted surface = aspirin binding pocket, red dotted surface = rofecoxib binding pocket, blue = aspirin, red = rofecoxib, yellow = AGP, green = A1, orange = A2, magenta = A3, hot pink = A4, dark slate grey = A5, cyan = A6, and plum = A7. This figure is plotted using ICM-Browser52.

2D binding interactions of aspirin (a), rofecoxib (b), AGP (c), and AGP analogs (A1–A7) (d–j) with the binding site of AF-COX-2 protein. All the chemical compounds are depicted as ball and stick model; amino acids with their number are represented as a circle with three-letter code and a dashed line shows the interaction site of chemical compounds with amino acids. Color code: green = H-bonds, pink = HYD interactions.

The physicochemical, absorption, distribution, metabolism, excretion, and toxicological properties of AGP and A1–A7 analogs were predicted and compared with aspirin and rofecoxib to evaluate their drug-likeness (Table 3). Physicochemical properties such as molecular weight (MW), number of H-bond acceptors (nHA), number of H-bond donors (nHD), topological polar surface area (TPSA), and the logarithm of the n-octanol/water distribution coefficient (logP) were calculated for all the selected compounds. Further, medicinal chemistry properties, i.e. natural product-likeness score (NPscore), Lipinski rule, and Pfizer rule, were analyzed for all the compounds. Except for rofecoxib and A5 analog, all the compounds fulfilled the criteria for physicochemical and medicinal chemistry properties for drug-like candidates. Furthermore, the absorption properties such as Caco-2 (the human colon adenocarcinoma cell lines) permeability, Pgp-inhibitory activity (inhibitor of P-glycoprotein), human intestinal absorption (HIA range 0–0.3), and F30% (the human oral bioavailability 30%) were also computed for all the compounds. Except for rofecoxib and A6 analog, all the compounds fulfilled the criteria to be used as orally active drugs. In addition, distribution properties for example plasma protein binding (PPB), the volume of distribution (VD), blood–brain barrier (BBB) penetration, and fraction unbound in plasma (Fu) were predicted for all the studied molecules. The distribution properties of AGP and A1–A7 analogs were found to be in the acceptable range for a drug-like candidate. Moreover, the cumulative analysis of metabolism and excretion properties also favored the AGP and A1–A7 analog as potential drug candidates. Further, the toxicological properties assessment also predicted the non-toxic nature of the A1–A7 AGP analogs.

The ligand binding quality of the selected was calculated based on the molecular mass and lipophilicity with the context of the particular protein target. Typically, a low Ki value refers to high potency and it should be in the micromolar range for a molecule to be qualified as a hit or lead compound55. Ligand efficiency (threshold value = 0.3 kcal/mol/HA) gives an idea about how well a compound binds for its size with the protein56. Ligand lipophilic efficiency (threshold value = 3) is used to measure the affinity of a molecule towards lipophilicity which is the major factor for the promiscuity of the selective LLE and optimized compounds54,56,57. LE_Scale is a size-independent scaling function of ligand efficiency and is derived by fitting an exponential function to the maximal ligand efficiency values observed for a given heavy atom (HA). Fit quality (threshold value = 0.8), also called as scaled ligand efficiency, includes ligand efficiency and size into a single metric and results by dividing the observed LE with LE_scale55,56,58. Lastly, LELP (acceptable range = -10 to + 10) was calculated which presents a metric that indicates the price of ligand efficiency paid in lipophilicity54,55. Further, binding energy (BE) obtained from the docking studies (Table 1), HA and LogP acquired from physicochemical properties prediction (Table 3) were taken for the calculation of all the LE metrics. It can be observed from Table 4 that except for A5 and controls (aspirin and rofecoxib), all the screened analogs show the respective ligand efficiency metrics values which surpass the minimum thresholds and to be called as HIT analogs. Also, in terms of inhibition constant, analogs A1–A5 displayed a Ki value smaller than the controls (aspirin Ki = 75.604 µM and rofecoxib Ki = 0.995 µM) and AGP (Ki = 1.440 µM) indicating their greater binding affinity towards the inhibition of COX-2 protein than controls and AGP.

The frontier molecular orbitals (FMOs) such as HOMO, LUMO, and HLG were calculated to anticipate the electronic properties of the AGP, hit analogs, aspirin, and rofecoxib against AF-COX-2. In the case of AGP and A1–A7 analogs, HOMO and LUMO orbitals locate on the naphthalene/naphthol and oxalone moiety, respectively, whereas, in the case of aspirin, HOMO and LUMO orbitals localize on the benzoic acid part of the structure. In contrast, HOMO orbitals confine to the phenyl-furan moiety while the LUMO orbitals position on the methylsulfonylpheny-furan moiety of rofecoxib (Fig. 3). The HLG between two energetic states, i.e. HOMO and LUMO, explains the chemical reactivity of molecules where the lowest HLG value infers greater chemical reactivity and biological activity, and lesser stability of the compound by making the flow of electrons easier72. It also plays an important role in protein–ligand complex stabilization73. The HLG values were found to be in the range of 4.24 to 5.63 eV for all the investigated compounds with the order of A6 (5.63 eV) > A7 (5.51 eV) = A5 (5.51 eV) > aspirin (5.43 eV) > AGP (5.26 eV) > A4 (5.21 eV) > A2 (5.15 eV) > A3 (5.05 eV) > A1 (5.03 eV) > rofecoxib (4.24 eV). Based on the HLG value, both A1 and A3 show the most chemical reactivity, whereas A6 is the least reactive and most stable.

Illustration of HOMO and LUMO molecular orbitals of aspirin (a), rofecoxib (b), AGP (c), and AGP's hit analogs {A1–A7} (d–j) with their HLG (HOMO–LUMO-GAP). The upper structure shows LUMO while the lower structure depicted HOMO orbitals. Color code: red = most negative potential; blue = most positive potential.

The EPE map is used to predict the active sites for intermolecular interactions as a spatial electron density distribution over the chemical compounds74,75. Figure 4 shows the EPE map of aspirin, rofecoxib, AGP, and A1–A7 in which the blue color region represents the most positive potential that is responsible for nucleophilic interaction while the red color indicates the most negative potential that is the site for electrophilic interaction. It is evident from the EPE map analysis that acetyloxy and carboxyl groups in aspirin, methylsulfonyl and furan groups in rofecoxib, and oxalone and hydroxyl groups of naphthalene moiety in AGP and hit analogs are in the negative potential (red region) and thus these are the prominent site for the electrophilic interactions. In contrast, the positive potential (blue region), which is a potential site for nucleophilic interactions, lies around the hydrogen of the carboxyl group in aspirin, phenyl group of rofecoxib, and naphthalene or naphthol moiety in AGP and A1–A7 analogs. The positive EPE (kcal/mol) were found to be in the order of rofecoxib (36.24) < aspirin (59.54) < A7 (61.75) < A3 (62.15) < A4 (64.90) < AGP (65.08) < A1 (65.50) < A2 (65.79) < A6 (67.07) < A5 (68.07). The negative EPE (kcal/mol) trend was different and followed the order of A7 (− 48.98) > AGP (48.97) > A3 (48.31) > A6 (-47.84) > A4 (− 46.65) > A2 (− 46.53) > A1(− 45.97) > A5 (− 43.98) > rofecoxib (−43.51) > aspirin (− 41.70) (Fig. 4). As the negative value of EPE favors the stronger H-bond formation76,77, the AGP and its A1–A7 hit analogs showed better potential towards stronger H-bond interaction compared to the aspirin and rofecoxib.

Electrostatic potential energy maps of aspirin (a), rofecoxib (b), AGP (c), and hit AGP analogs {A1–A7} (d–j). The color code for EPE contours maps shows the most negative value in red to the most positive value in blue color; The order of potential is: blue (most + ve) > green > yellow > orange > red (most -ve).

Structurally, all the screened AGP analogs share the main chemical scaffold which consists of two cyclic groups (i.e. naphthalene/naphthol and oxalone) linked with ethyledine moiety. Bicyclic groups containing compounds manifest a myriad of biological activities78. These bicyclic groups containing drugs predominantly interact with proteins through hydrophobic, and aromatic residues. Further, it has also been demonstrated that the bicyclic compounds interact with the polar groups such as hydroxyl, and amides of the amino acid79. In this study, amongst the screened AGP analogs, there are two stereoisomers A1 (6S)/A3 (6R), and A2 (6S)/A4 (6R), whereas, analogs A5, A6, and A7 are distinct from each other. Compared to AGP, analogs A1 to A5 structurally differ at 5 and 6 positions of naphthalene moiety, while the rest of the chemical scaffold is the same. In place of naphthalene, analogs A6 and A7 possessed naphthol moiety and has a different group at the 3rd position of the naphthol ring. Further, molecular interaction analysis revealed that amongst screened analogs naphthalene moiety-containing compounds showed more affinity compared to naphthol towards the COX-2 protein in terms of binding energy, H-bonds and HYD interactions (Table 1 and 2). It is evident from Fig. 2 and Table 2 of the molecular interaction analysis that the oxolane moiety is mainly involved in the H-bond interaction while naphthalene/naphthol moiety mainly concerned with the HYD interactions. This was further confirmed by quantum mechanical descriptor analysis. Structurally, compared to AGP, analog A1–A4 has one less –OH group at the 5th position in the naphthalene moiety which increases the duration of action of these analogs which has been proved by ADMET prediction (Table 3). This further increases the HYD interaction of these types of analogs with the COX-2 protein and also increases the affinity with the receptor80. Therefore, based on the above relationship between the chemical structure features and the results obtained from the molecular docking and DFT analysis, we concluded that screened AGP analogs shared improved anti-inflammatory activity compared to AGP and the reference drugs (aspirin and rofecoxib).

MD simulation analysis was carried out for AGP and hit analogs’ complexes with AF-COX-2 protein along with appropriate controls {AF-COX-2 protein alone (neutral control) and its complex with aspirin and rofecoxib (positive control)}. The MD simulation trajectory analysis of AGP and hit analogs’ complex was compared with native AF-COX-2 protein, aspirin, and rofecoxib complex in terms of RMSD, RMSF, residual area, Rg, SAS area, SAS volume, SAS density, protein H-bonds, and protein–ligand complex H-bonds.

The RMSD measures the protein conformation by calculating the average distance between the backbone atoms of a protein during the simulation and assessing the structural stability. The average RMSD values for the AGP and hit analogs were found to be lower than the neutral control (native protein). In comparison with positive controls, the average RMSD value was observed to be lesser in analog A2, A3, and A4. The least value RMSD value was observed for the analog A3 indicating its ability to form a stable adduct with AF-COX-2 protein (Fig. 5a and Supplementary Table S6). The RMSF is calculated for the residual mobility of each protein–ligand complex as well as the native protein. No noticeable difference was found between the residual fluctuation profiles of all the complexes (Fig. 5b). This may be ascribed to the formation of H-bonds and van der Waal interactions between amino acid residues and ligands. As shown in Fig. 5b, the active binding site regions have lower fluctuations indicating their stability in the AF-COX-2 protein. Further, the average RMSF was observed to be lower in the A1, A3, A6, and A7 complexes than in the neutral control as well as the positive control (Supplementary Table S6). The individual RMSF for the interacting amino acid residues in each simulation complex system was computed and compared with the respective average RMSF value for the simulation system. Further, from the Supplementary Table S6, the average RMSF range for all the studied simulation systems was found to be 0.16 ± 0.10 to 0.21 ± 0.18 nm. Also, it is evident that except Arg106, Ser339, and Phe 504, rest of the interacting amino acid exhibited the RMSF value lower than 0.21 nm (which is highest average RMSF value amongst all the simulation system) (Supplementary Table S7). This means that the interacting binding residues are quite stable during the MD simulation. Moreover, the lower average RMSF value for the analogs A1, A2, A3, A6 and A7 than the protein demonstrate that the analogs are capable of forming stable interactions with the protein during MD simulation. This RMSF study suggested the ability of these four AGP analogs towards the formation of stable interaction with AF-COX-2 protein during the period of simulation.

Pictorial representation of MD simulation analysis results in terms of RMSD (a), RMSF (b), SAS area for each residue (c), and Rg (d) of the native AF-COX-2 protein and its complex with subjected ligands.

The SAS area for each residue during the simulation was also computed in order to calculate the area of each residue that is accessible for solvent (Fig. 5c)81. The SAS value per residue was found to be in the range of 0–2.49 nm2. It is evident that a few residues of the AF-COX-2 protein are not accessible for solvents such as Asn181 in rofecoxib and A2 complex, and Lys237 in A7 complex. Further, amongst all the studied analogs and control compounds, the average value of SAS area per residue was found to be highest for A3. Further, the maximum value of SAS area per residue (nm2) was found to be in the order of A4 (2.49) > A5 (2.44) > A3 (2.27) > A6 (2.22) > Rofecoxib (2.21) > A7 (2.19) > A2 (2.18) > A1 (2.11) > AGP (2.06) > aspirin (2.06) > protein (2.05) which signifies that analogs A3–A6 have the better accessibility for the solvent than the other complexes as well as native protein (Supplementary Table S8).

The Rg analysis was carried out to calculate the level of compactness in terms of protein folding and unfolding of the AF-COX-2 protein structure in the presence and absence of the ligands. The average Rg value for analog A3 and A6 was found to be equivalent to the Rg value of native protein (neutral control) but lesser than the positive control. This suggests a more stable and compact protein–ligand complex formation compared to other analogs and positive controls (Fig. 5d and Supplementary Table S8)82.

The Number of H-bonds in protein as well as between protein and ligand play a significant role in the stability of complexes, therefore, H-bond analysis was performed for native AF-COX-2 and its complexes with AGP, hit analogs, aspirin, and rofecoxib. It is evident from Fig. 6a that the analog A3 complex exhibited a maximum number of H bonds in protein, whereas, AGP and A3 complexes exhibited the maximum number of protein–ligand H-bonds as compared to other analogs and positive (Fig. 6b). This suggests the stronger and more efficient binding of the analog A3 with AF-COX-2 protein which is in good agreement with the results obtained from RMSD and RMSF analysis (Supplementary Table S6).

Graphical illustration of MD simulation studies of native COX-2 protein, and its complexes with aspirin, rofecoxib, AGP, and hit AGP analogs. (a) Number of H-bonds in protein, (b) number of H-bonds between protein and ligand, (c) SAS area, (d) SAS volume, and (e) SAS density.

SASA, the area of a complex surface that is accessible to a water solvent, is used to predict the level of conformational variations that occur through the interaction between proteins and ligands. Figure 6c shows the SASA plot for all the hit analogs, AGP, neutral, and positive control. The average SASA value for all the hit AGP analogs along with AGP was found to be more than the neutral control (native protein SASA = 273.57 nm2). Whereas, as compared to the positive control (rofecoxib and aspirin complex SASA = 273.13 and 279.12 nm2, respectively), AGP, A2, A3, A4, and A5 showed greater SASA values (i.e. > 279.12 nm2) (Supplementary Table S9). Amongst all the hit analogs and AGP, A3 exhibited the highest SASA value = 284.59 nm2. Conclusively, the A3 analog forms a relatively stable complex with AF-COX-2 protein after the interaction. Further, solvent–Accessible surface volume (Vsas), the volume enclosed by the center of a solvent probe rolling around the protein, was also calculated for all the studied ligands which signifies the stability of the complex system83. Typically, it measures the effect of forces on the protein surfaces exerted by solvents followed by the protein-solvent interactions. Also, Vsas is an alternative to refine the SASA term, in order to include the influence of the solvents’ effect on the protein's interior and also define the interaction between the protein and solvent84. It is an accurate and fast application to examine geometric volumes of the proteins. In addition, it can be used to calculate the volume that changes due to the interaction of protein–ligand complex with solvent85. Vsas along with SASA provides a better acquiescence of implicit and explicit nonpolar solvent forces on the protein and its effect on the protein folding86. As shown in Fig. 6d and Supplementary Table S9, AGP, A3, A4, and A5 analog complex showed a greater value of Vsas compared to all controls, while the highest value was observed for the A3 analog complex (119.24 nm\s3\n). These results indicate that the A3 analog forms the most stable complex with AF-COX-2 protein compared to all other analogs and controls. Furthermore, the SAS density, which is inversely proportional to the SASA, was determined to calculate the neighborhood density of burial HYD amino acids within the protein core that leads to protein folding87. The neighborhood density calculates the precise molecular and atomic quantities with coordinates for protein surface area accessible for solvent. It not only measures the hydrophobic effect of the amino acids density but also captures the electrostatic effect of the solvent on the protein folding and stability88. As demonstrated in Fig. 6e and Supplementary Table S9, the A3 analog possessed the least value of SAS density indicating better protein folding compared to other hit analogs and all controls after interaction with AF-COX-2 protein. All three solvent accessible surface analysis (area, volume, and density) together provide more emphasis on the MD simulation analysis related to the analogs’ complex stability in terms of protein folding and revealed more stability of analog A3 than the other studied compounds including the reference molecules. Overall, from MD simulation analysis, it can be inferred that A3 analog exhibit better protein stability than the other hit analogs, AGP, native protein, aspirin, and rofecoxib.

Moreover, the molecular interactions of the hit analogs were also analyzed after the MD simulation (Supplementary Fig. S3 and Table S10). Interestingly, except A7, the binding site for all the hit analogs were observed to be consistent before and after the simulation. Further, amino acid residues namely Val102, Arg106, Val335, Leu338, Ser339, Tyr341, Leu345, Phe504, Val509, Ala513, and Ser516 were found to be important residues before and after simulation as these are consistent before and after simulation. Further, a comparative analysis of before and after MD simulation in aspirin protein revealed that after simulation two hydrophobic interaction (Val509, Ala513) and two H-bonds (Arg106 and Arg499) are formed which were three and one respectively before the simulation. In case of rofecoxib, only three hydrophobic interactions (Leu338, Tyr371, and Trp373) were observed after the simulation, whereas before simulation one H-bond (residue) and six hydrophobic interactions (residue) were noted. Next, AGP protein complex exhibited six hydrophobic interactions without H-bond after simulation whereas three H-bond with six hydrophobic interaction were observed before simulation. Then, analog A1 protein complex showed two H-bonds with seven hydrophobic interactions after simulation while three H-bonds with five hydrophobic interactions were found before simulation. Analog A2 protein complex revealed nine hydrophobic interactions without any H-bonds after simulation whereas four H-bonds with six hydrophobic interactions were observed before simulation. Analog A3 protein complex displayed one H-bond with eight hydrophobic interactions after simulation while seven hydrophobic interactions along with one H-bond were observed before simulation. A4 protein complex exhibited one H-bond and seven hydrophobic interaction after simulation while three H-bonds with six hydrophobic interactions were observed before simulation. A5 protein complex revealed 2 H-bonds and seven hydrophobic interaction after simulation whereas 4 H-bonds with five hydrophobic interactions before simulation. A6 protein complex showed one H-bond with four hydrophobic interactions after simulation while only hydrophobic interactions were observed before simulation. Lastly, A7 protein complex revealed one H-bond with three hydrophobic interaction and only hydrophobic interactions were majorly involved before simulation (Supplementary Table S10). Overall, from the comparative analysis, it was observed that the hit analogs except the A7 remain close to the binding site identified by the molecular docking analysis which indicates the efficient binding stability of the complex of hit analogs with AF-COX-2. Also, analog A1, A3 and A5 complexes maintained their consistencies before and after simulation.

The EPE maps describe the distribution of electrons on the molecular surface, HOMO, and LUMO energy of the molecules. In this study, we have calculated the EPE, HOMO, and LUMO of the selected sequence of AF-COX-2 protein and its adduct with the AGP, hit AGP analogs, native protein, aspirin, and rofecoxib. Then, the binding energy of the adduct was computed, where the more negative value of binding energy indicates the more stability of the protein–ligand complex89. For this study, aspirin, rofecoxib, AGP, and A1–A7 AGP analogs were selected as ligands, while the protein sequences from 511 to 520 amino acids of AF-COX-2, i.e. Val511, Gly512, Ala513, Pro514, Phe515, Ser516, Leu517, Lys518, Gly519, Leu520 were selected as a peptide for docking simulation study by calculating EPE35. Later, the calculated binding energy values of the adducts were found to be in the range of − 53.81 kcal/mol (A3) to − 34.53 kcal/mol (aspirin) (Table 5). Further, all the AGP and all seven analogs exhibited lesser binding energy than the rofecoxib and aspirin, while the least binding energy value was observed for A3 suggesting the most stable protein–ligand complex.

The binding free energy of the interacting hit AGP analogs, AGP, aspirin, and rofecoxib to the AF-COX-2 protein was computed by MM/GBSA method for before and after simulation analysis. The results revealed that there is a subtle difference between the MM/GBSA score for before and after simulation analysis (Fig. 7). Also, It can be noted from Fig. 7 that amongst all the investigated compounds, the A3 analog showed the most negative binding free energy (before simulation = − 55.37 kcal/mol and after simulation = − 56.25 kcal/mol) in terms of MM/GBSA indicating its strong interaction with AF-COX-2 protein. The MM/GBSA score was found to be in the decreasing order of binding affinity A3 > A1 > A5 > Rofecoxib > A2 > A4 > A6 > A7 > AGP > aspirin. Further, it was observed that all the hit AGP analogs along with the AGP exhibited a more negative value of free energy compared to aspirin, whereas analogs A3, A1, and A5 showed the binding affinity better than Rofecoxib, aspirin, AGP, and other analogs. Moreover, analog A3, A1 and A5 complex with protein were observed to be more stable after the simulation analysis and also implies that the analog A3 showed the better binding towards the COX-2 protein as compared to the reference compounds and other studied analogs. The MM/GBSA results also corroborated the findings of docking, MD simulation, and DFT analysis, and establish that the A3 analog would be a promising compound to inhibit the AF-COX-2 activity and can be used as a promising natural anti-inflammatory drug.

The MM/GBSA energy score of rofecoxib (1), aspirin (2), AGP (3), and selected analogs A1–A7 (4–10, respectively) for before and after MD simulation.

Overall, our results reveals that andrographolide analog A3 is the most potent and effective plant based natural drug molecule against COX-2 protein. The efficacy of this analog was compared with other natural compounds as well as synthetic compounds for COX-2 inhibition. The binding energy for A3 from docking analysis was found to be − 8.56 kcal/mol which is better than virgin coconut oil derivatives (BE range for the derivatives: − 5.65 kcal/mol to − 7.58 kcal/mol), alliin (− 4.90 kcal/mol), pinoresinolm (− 8.38 kcal/mol), and syringaresinol (− 8.23 kcal/mol)90,91,92,93. Further, MD simulation as well as EPE docking simulation analysis were carried for the complex stability for A3. Moreover, the detailed comparative analysis results with the reference compounds supports that the identified analog A3 would be an efficient drug-like candidate against COX-2 protein. To the best of our knowledge, to date, the andrographolide analog A3 has not been explored for anti-inflammatory activity by inhibiting the COX-2. Hence, present work provides a foundation for the development of andrographolide analog A3 as a potential anti-inflammatory drug-like candidate acting via COX-2 inhibition. Though, we have employed a myriad of cheminformatic tools for identification and validation of analog A3 as a potential COX-2 inhibitor, further experimental authentication is required from the scientific community to expedite the drug development process.

In the present study, a promising plant-based anti-inflammatory compound against COX-2 was identified by employing a combined computational approach including molecular and quantum docking simulation studies. A validated full sequenced human AF-COX-2 protein structure was used for multiple sequence alignment to find out the active site residues for anti-inflammatory activity. The virtual screening of 237 AGP analogs against AF-COX-2 protein followed by interactive docking studies screened out 7 hit AGP analogs which showed their drug-like candidacy also from ADMET prediction analysis. Further, ligand efficiency metrics, DFT descriptors, i.e. HOMO, LUMO, HLG, and EPE established their good chemical reactivity. In addition, MD simulation, EPE docking simulation, and MM/GBSA study revealed that A3 analog (3-[2-[(1R,4aR,5R,6R,8aR)-6-hydroxy-5,6,8a-trimethyl-2-methylidene-3,4,4a,5,7,8-hexahydro-1H-naphthalen-1-yl]ethylidene]-4-hydroxyoxolan-2-one) forms the most stable complex with the AF-COX-2 and would be a promising natural anti-inflammatory agent. In near future, additional experimental validations are required to substantiate these findings.

All data generated or analyzed during this study are included in this article and its Supplementary Information files. Additional raw data or datasets generated and/or analyzed during the current study are available from the corresponding author on reasonable request.

Smith, W. L., DeWitt, D. L. & Garavito, R. M. Cyclooxygenases: Structural, cellular, and molecular biology. Annu. Rev. Biochem. 69, 145–182 (2000).

Article CAS PubMed Google Scholar

Vane, J. R. Inhibition of prostaglandin synthesis as a mechanism of action for aspirin-like drugs. Nat. New Biol. 231, 232–235 (1971).

Article CAS PubMed Google Scholar

Morita, I. Distinct functions of COX-1 and COX-2. Prostaglandins Other Lipid Mediat. 68–69, 165–175 (2002).

Article PubMed Google Scholar

Simon, L. S. Role and regulation of cyclooxygenase-2 during inflammation. Am. J. Med. 106, 37S-42S (1999).

Article CAS PubMed Google Scholar

Orlando, B. J. & Malkowski, M. G. Substrate-selective inhibition of cyclooxygeanse-2 by fenamic acid derivatives is dependent on peroxide tone. J. Biol. Chem. 291, 15069–15081 (2016).

Article CAS PubMed PubMed Central Google Scholar

Rouzer, C. A. & Marnett, L. J. Mechanism of free radical oxygenation of polyunsaturated fatty acids by cyclooxygenases. Chem. Rev. 103, 2239–2304 (2003).

Article CAS PubMed Google Scholar

Wong, E., Bayly, C., Waterman, H. L., Riendeau, D. & Mancini, J. A. Conversion of prostaglandin G/H synthase-1 into an enzyme sensitive to PGHS-2-selective inhibitors by a double His513→Arg and Ile523→Val mutation*. J. Biol. Chem. 272, 9280–9286 (1997).

Article CAS PubMed Google Scholar

Gierse, J. K. et al. A single amino acid difference between cyclooxygenase-1 (COX-1) and −2 (COX-2) reverses the selectivity of COX-2 specific inhibitors*. J. Biol. Chem. 271, 15810–15814 (1996).

Article CAS PubMed Google Scholar

Cipollone, F., Cicolini, G. & Bucci, M. Cyclooxygenase and prostaglandin synthases in atherosclerosis: Recent insights and future perspectives. Pharmacol. Ther. 118, 161–180 (2008).

Article CAS PubMed Google Scholar

Rouzer, C. A. & Marnett, L. J. Cyclooxygenases: Structural and functional insights. J. Lipid Res. 50, 29–34 (2009).

Article Google Scholar

Blobaum, A. L. & Marnett, L. J. Structural and functional basis of cyclooxygenase inhibition. J. Med. Chem. 50, 1425–1441 (2007).

Article CAS PubMed Google Scholar

Wallace, J. L. Selective cyclooxygenase-2 inhibitors: After the smoke has cleared. Dig. Liver Dis. 34, 89–94 (2002).

Article CAS PubMed Google Scholar

Marnett, L. J. The COXIB experience: A look in the rearview mirror. Annu. Rev. Pharmacol. Toxicol. 49, 265–290 (2009).

Article CAS PubMed Google Scholar

Hilário, M. O. E., Terreri, M. T. & Len, C. A. Nonsteroidal anti-inflammatory drugs: Cyclooxygenase 2 inhibitors. J. Pediatr. (Rio. J.) 82, 206–212 (2006).

Article Google Scholar

Ahuja, N., Singh, A. & Singh, B. Rofecoxib: An update on physicochemical, pharmaceutical, pharmacodynamic and pharmacokinetic aspects. J. Pharm. Pharmacol. 55, 859–894 (2010).

Article Google Scholar

Mukherjee, D., Nissen, S. E. & Topol, E. J. Risk of cardiovascular events associated with selective COX-2 inhibitors. JAMA 286, 954–959 (2001).

Article CAS PubMed Google Scholar

Funk, C. D. & FitzGerald, G. A. COX-2 inhibitors and cardiovascular risk. J. Cardiovasc. Pharmacol. 50, 470–479 (2007).

Article CAS PubMed Google Scholar

Chen, W. et al. New horizons in the roles and associations of COX-2 and novel natural inhibitors in cardiovascular diseases. Mol. Med. 27, 123 (2021).

Article CAS PubMed PubMed Central Google Scholar

Andersohn, F., Schade, R., Suissa, S. & Garbe, E. Cyclooxygenase-2 selective nonsteroidal anti-inflammatory drugs and the risk of ischemic stroke: A nested case-control study. Stroke 37, 1725–1730 (2006).

Article CAS PubMed Google Scholar

Attiq, A., Jalil, J., Husain, K. & Ahmad, W. Raging the war against inflammation with natural products. Front. Pharmacol. 9, 1–27 (2018).

Article Google Scholar

Lee, K. C., Chang, H. H., Chung, Y. H. & Lee, T. Y. Andrographolide acts as an anti-inflammatory agent in LPS-stimulated RAW264.7 macrophages by inhibiting STAT3-mediated suppression of the NF-κB pathway. J. Ethnopharmacol. 135, 678–684 (2011).

Article CAS PubMed Google Scholar

Yuan, L. et al. The semi-synthesis of novel andrographolide analogues and anti-influenza virus activity evaluation of their derivatives. Bioorganic Med. Chem. Lett. 26, 769–773 (2016).

Article CAS Google Scholar

Jiao, J. et al. Screening cyclooxygenase-2 inhibitors from Andrographis paniculata to treat inflammation based on bio-affinity ultrafiltration coupled with UPLC-Q-TOF-MS. Fitoterapia 137, 104259 (2019).

Article CAS PubMed Google Scholar

Singh, J., Petter, R. C., Baillie, T. A. & Whitty, A. The resurgence of covalent drugs. Nat. Rev. Drug Discov. 10, 307–317 (2011).

Article CAS PubMed Google Scholar

Tran, Q. T. N., Tan, D. W. S., Wong, W. S. F. & Chai, C. L. L. From irreversible to reversible covalent inhibitors: Harnessing the andrographolide scaffold for anti-inflammatory action. Eur. J. Med. Chem. 204, 112481 (2020).

Article CAS PubMed Google Scholar

Tran, Q. T. N., Tan, W. S. D., Wong, W. S. F. & Chai, C. L. L. Polypharmacology of andrographolide: Beyond one molecule one target. Nat. Prod. Rep. 38, 682–692 (2021).

Article CAS PubMed Google Scholar

Nguyen, V. S. et al. Specificity and inhibitory mechanism of andrographolide and its analogues as antiasthma agents on NF-κB p50. J. Nat. Prod. 78, 208–217 (2015).

Article CAS PubMed Google Scholar

Burgos, R. A., Alarcón, P., Quiroga, J., Manosalva, C. & Hancke, J. Andrographolide, an anti-inflammatory multitarget drug: All roads lead to cellular metabolism. Molecules 26, 5 (2021).

Article CAS Google Scholar

Dai, G. F. et al. Anti-inflammatory effect of novel andrographolide derivatives through inhibition of NO and PGE 2 production. Int. Immunopharmacol. 11, 2144–2149 (2011).

Article CAS PubMed Google Scholar

Wang, W. et al. Synthesis of new ent-labdane diterpene derivatives from andrographolide and evaluation of their anti-inflammatory activities. Eur. J. Med. Chem. 162, 70–79 (2019).

Article CAS PubMed Google Scholar

Peng, Y. et al. Andrographolide inhibits breast cancer through suppressing COX-2 expression and angiogenesis via inactivation of p300 signaling and VEGF pathway 11 medical and health sciences 1112 oncology and carcinogenesis. J. Exp. Clin. Cancer Res. 37, 1–14 (2018).

CAS Google Scholar

Liu, W. et al. Andrographolide potentiates PD-1 blockade immunotherapy by inhibiting COX2-mediated PGE2 release. Int. Immunopharmacol. 81, 106206 (2020).

Article CAS PubMed Google Scholar

Chen, M., Xie, C. & Liu, L. Solubility of andrographolide in various solvents from (288.2 to 323.2) K. J. Chem. Eng. Data 55, 5297–5298 (2010).

Article CAS Google Scholar

Tunyasuvunakool, K. et al. Highly accurate protein structure prediction for the human proteome. Nature 596, 590–596 (2021).

Article ADS CAS PubMed PubMed Central Google Scholar

Lucido, M. J., Orlando, B. J., Vecchio, A. J. & Malkowski, M. G. Crystal structure of aspirin–acetylated human cyclooxygenase-2: Insight into the formation of products with reversed stereochemistry. Biochemistry 55, 1226–1238 (2016).

Article CAS PubMed Google Scholar

Orlando, B. J. & Malkowski, M. G. Crystal structure of rofecoxib bound to human cyclooxygenase-2. Acta Crystallogr. Sect. Struct. Biol. Commun. 72, 772–776 (2016).

Article CAS Google Scholar

Orlando, B. J. & Malkowski, M. G. Substrate-selective Inhibition of cyclooxygeanse-2 by fenamic acid derivatives is dependent on peroxide tone*. J. Biol. Chem. 291, 15069–15081 (2016).

Article CAS PubMed PubMed Central Google Scholar

Colovos, C. & Yeates, T. O. Verification of protein structures: Patterns of nonbonded atomic interactions. Protein Sci. 2, 1511–1519 (1993).

Article CAS PubMed PubMed Central Google Scholar

Pontius, J., Richelle, J. & Wodak, S. J. Deviations from standard atomic volumes as a quality measure for protein crystal structures. J. Mol. Biol. 264, 121–136 (1996).

Article CAS PubMed Google Scholar

Wiederstein, M. & Sippl, M. J. ProSA-web: Interactive web service for the recognition of errors in three-dimensional structures of proteins. Nucleic Acids Res. 35, 407–410 (2007).

Article Google Scholar

Wallner, B. & Elofsson, A. Can correct protein models be identified?. Protein Sci. 12, 1073–1086 (2003).

Article CAS PubMed PubMed Central Google Scholar

Cristobal, S., Zemla, A., Fischer, D., Rychlewski, L. & Elofsson, A. A study of quality measures for protein threading models. BMC Bioinformatics 2, 5 (2001).

Article CAS PubMed PubMed Central Google Scholar

Bhattacharya, D., Nowotny, J., Cao, R. & Cheng, J. 3Drefine: An interactive web server for efficient protein structure refinement. Nucleic Acids Res. 44, W406–W409 (2016).

Article CAS PubMed PubMed Central Google Scholar

Tian, W., Chen, C., Lei, X., Zhao, J. & Liang, J. CASTp 3.0: computed atlas of surface topography of proteins. Nucleic Acids Res. 46, W363–W367 (2018).

Article CAS PubMed PubMed Central Google Scholar

Pettersen, E. F. et al. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612 (2004).

Article CAS PubMed Google Scholar

Simossis, V. A. & Heringa, J. PRALINE: A multiple sequence alignment toolbox that integrates homology-extended and secondary structure information. Nucleic Acids Res. 33, 289–294 (2005).

Article Google Scholar

Agarwal, S., Dixit, A. & Kashaw, S. K. Ligand and structure based virtual screening of chemical databases to explore potent small molecule inhibitors against breast invasive carcinoma using recent computational technologies. J. Mol. Graph. Model. 98, 107591 (2020).

Article CAS PubMed Google Scholar

Kim, S. et al. PubChem in 2021: New data content and improved web interfaces. Nucleic Acids Res. 49, D1388–D1395 (2021).

Article CAS PubMed Google Scholar

O’Boyle, N. M. et al. Open babel: An open chemical toolbox. J. Cheminform. 3, 33 (2011).

Article PubMed PubMed Central Google Scholar

Dallakyan, S. & Olson, A. J. Small-molecule library screening by docking with PyRx. In Chemical biology: Methods and protocols (eds Hempel, J. E. et al.) 243–250 (Springer, 2015). https://doi.org/10.1007/978-1-4939-2269-7_19.

Chapter Google Scholar

Morris, G. M. et al. AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J. Comput. Chem. 30, 2785–2791 (2009).

Article CAS PubMed PubMed Central Google Scholar

Abagyan, R., Totrov, M. & Kuznetsov, D. ICM—A new method for protein modeling and design: Applications to docking and structure prediction from the distorted native conformation. J. Comput. Chem. 15, 488–506 (1994).

Article CAS Google Scholar

Xiong, G. et al. ADMETlab 2.0: An integrated online platform for accurate and comprehensive predictions of ADMET properties. Nucleic Acids Res. 49, W5–W14 (2021).

Article CAS PubMed PubMed Central Google Scholar

Hopkins, A. L., Keserü, G. M., Leeson, P. D., Rees, D. C. & Reynolds, C. H. The role of ligand efficiency metrics in drug discovery. Nat. Rev. Drug Discov. 13, 105–121 (2014).

Article CAS PubMed Google Scholar

Onawole, A. T., Kolapo, T. U., Sulaiman, K. O. & Adegoke, R. O. Structure based virtual screening of the Ebola virus trimeric glycoprotein using consensus scoring. Comput. Biol. Chem. 72, 170–180 (2018).

Article CAS PubMed Google Scholar

Schultes, S. et al. Ligand efficiency as a guide in fragment hit selection and optimization. Drug Discov. Today Technol. 7, e157–e162 (2010).

Article CAS Google Scholar

Murray, C. W. et al. Validity of ligand efficiency metrics. ACS Med. Chem. Lett. 5, 616–618 (2014).

Article CAS PubMed PubMed Central Google Scholar

Reynolds, C. H., Bembenek, S. D. & Tounge, B. A. The role of molecular size in ligand efficiency. Bioorg. Med. Chem. Lett. 17, 4258–4261 (2007).

Article CAS PubMed Google Scholar

Becke, A. D. Density-functional thermochemistry. III: The role of exact exchange. J. Chem. Phys. 98, 5648–5652 (1993).

Article ADS CAS Google Scholar

Rassolov, V. A., Ratner, M. A., Pople, J. A., Redfern, P. C. & Curtiss, L. A. 6–31G* basis set for third-row atoms. J. Comput. Chem. 22, 976–984 (2001).

Article CAS Google Scholar

Gill, P. M. W., Johnson, B. G., Pople, J. A. & Frisch, M. J. The performance of the Becke–Lee–Yang–Parr (B–LYP) density functional theory with various basis sets. Chem. Phys. Lett. 197, 499–505 (1992).

Article ADS CAS Google Scholar

Stephens, P. J., Devlin, F. J., Chabalowski, C. F. & Frisch, M. J. Ab Initio calculation of vibrational absorption and circular dichroism spectra using density functional force fields. J. Phys. Chem. 98, 11623–11627 (1994).

Article CAS Google Scholar

Abraham, M. J. et al. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 1–2, 19–25 (2015).

Article ADS Google Scholar

Schüttelkopf, A. W. & Van Aalten, D. M. F. PRODRG: A tool for high-throughput crystallography of protein-ligand complexes. Acta Crystallogr. Sect. D Biol. Crystallogr. 60, 1355–1363 (2004).

Article Google Scholar

Schmid, N. et al. Definition and testing of the GROMOS force-field versions 54A7 and 54B7. Eur. Biophys. J. 40, 843–856 (2011).

Article CAS PubMed Google Scholar

Abdel-Mottaleb, M. S. A. & Abdel-Mottaleb, Y. Impact of magnesium, zinc, selenium, copper, and iodine food supplements on SARS-CoV, SARS-CoV-2 viruses and their adducts with human ACE2 enzyme: A Computational Based Investigation. Egypt. J. Chem. 64, 989–996 (2021).

Google Scholar

Wang, Z. et al. FarPPI: A webserver for accurate prediction of protein-ligand binding structures for small-molecule PPI inhibitors by MM/PB(GB)SA methods. Bioinformatics 35, 1777–1779 (2019).

Article CAS PubMed Google Scholar

Maier, J. A. et al. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 11, 3696–3713 (2015).

Article CAS PubMed PubMed Central Google Scholar

Jakalian, A., Jack, D. B. & Bayly, C. I. Fast, efficient generation of high-quality atomic charges. AM1-BCC model: II. Parameterization and validation. J. Comput. Chem. 23, 1623–1641 (2002).

Article CAS PubMed Google Scholar

Wang, E. et al. End-point binding free energy calculation with MM/PBSA and MM/GBSA: Strategies and applications in drug design. Chem. Rev. 119, 9478–9508 (2019).

Article CAS PubMed Google Scholar

Liang, J., Woodward, C. & Edelsbrunner, H. Anatomy of protein pockets and cavities: Measurement of binding site geometry and implications for ligand design. Protein Sci. 7, 1884–1897 (1998).

Article CAS PubMed PubMed Central Google Scholar

Banavath, H. N., Sharma, O. P., Kumar, M. S. & Baskaran, R. Identification of novel tyrosine kinase inhibitors for drug resistant T315I mutant BCR–ABL: A virtual screening and molecular dynamics simulations study. Sci. Rep. 4, 1–11 (2014).

Article Google Scholar

V, S. S., Palaka, B. K., Vel, G., Venkatesan, R., Ampasala, D. R. & Periyasamy, L. Identification of novel inhibitors of signal transducer and activator of transcription 3 over signal transducer and activator of transcription 1 for the treatment of breast cancer by in-silico and in-vitro approach. Process Biochem. 82, 153–166 (2019).

Article Google Scholar

Chaudhary, M. K. et al. Computational evaluation on molecular stability, reactivity, and drug potential of frovatriptan from DFT and molecular docking approach. Comput. Theor. Chem. 1191, 27–34 (2013).

Google Scholar

Noureddine, O., Issaoui, N. & Al-Dossary, O. DFT and molecular docking study of chloroquine derivatives as antiviral to coronavirus COVID-19. J. King Saud Univ. Sci. 33, 101248 (2021).

Article PubMed Google Scholar

Mottishaw, J. D., Erck, A. R., Kramer, J. H., Sun, H. & Koppang, M. Electrostatic potential maps and natural bond orbital analysis: Visualization and conceptualization of reactivity in Sangers reagent. J. Chem. Educ. 92, 1846–1852 (2015).

Article CAS Google Scholar

Uzzaman, M., Junaid, M. & Uddin, M. N. Evaluation of anti-tuberculosis activity of some oxotitanium(IV) Schiff base complexes; molecular docking, dynamics simulation and ADMET studies. SN Appl. Sci. 2, 1–11 (2020).

Article Google Scholar

Horton, D. A., Bourne, G. T. & Smythe, M. L. The combinatorial synthesis of bicyclic privileged structures or privileged substructures. Chem. Rev. 103, 893–930 (2003).

Article CAS PubMed Google Scholar

Hajduk, P. J., Bures, M., Praestgaard, J. & Fesik, S. W. Privileged molecules for protein binding identified from NMR-based screening. J. Med. Chem. 43, 3443–3447 (2000).

Article CAS PubMed Google Scholar

El-Haj, B. M. & Ahmed, S. B. M. Metabolic-hydroxy and carboxy functionalization of alkyl moieties in drug molecules: Prediction of structure influence and pharmacologic activity. Molecules 25, 1937 (2020).

Article CAS PubMed PubMed Central Google Scholar

Unal, M. A., Boyacioglu, B., Unver, H. & Elmali, A. Molecular simulation of PcCel45A protein expressed from Aspergillus nidulans to understand its structure, dynamics, and thermostability. J. Mol. Model. 25, 1. https://doi.org/10.1007/s00894-019-4175-4 (2019).

Article CAS Google Scholar

Lobanov, M. Y., Bogatyreva, N. S. & Galzitskaya, O. V. Radius of gyration as an indicator of protein structure compactness. Mol. Biol. 42, 623–628 (2008).

Article CAS Google Scholar

Chen, C. R. & Makhatadze, G. I. ProteinVolume: Calculating molecular van der Waals and void volumes in proteins. BMC Bioinformatics 16, 1–6 (2015).

Article Google Scholar

Kim, I. J. & Na, H. An efficient algorithm calculating common solvent accessible volume. PLoS ONE 17, e0265614 (2022).

Article CAS PubMed PubMed Central Google Scholar

Chen, C. R. & Makhatadze, G. I. ProteinVolume: Calculating molecular van der Waals and void volumes in proteins. BMC Bioinformatics 16, 101 (2015).

Article PubMed PubMed Central Google Scholar

Wagoner, J. A. & Baker, N. A. Assessing implicit models for nonpolar mean solvation forces: The importance of dispersion and volume terms. Proc. Natl. Acad. Sci. 103, 8331–8336 (2006).

Article ADS CAS PubMed PubMed Central Google Scholar

Durham, E., Dorr, B., Woetzel, N., Staritzbichler, R. & Meiler, J. Solvent accessible surface area approximations for rapid and accurate protein structure prediction. J. Mol. Model. 15, 1093–1108 (2009).

Article CAS PubMed PubMed Central Google Scholar

Weiser, J., Shenkin, P. S. & Still, W. C. Approximate solvent–Accessible surface areas from tetrahedrally directed neighbor densities. Biopolymers 50, 373–380 (1999).

3.0.CO;2-U" data-track-action="article reference" href="https://doi.org/10.1002%2F%28SICI%291097-0282%2819991005%2950%3A4%3C373%3A%3AAID-BIP3%3E3.0.CO%3B2-U" aria-label="Article reference 88" data-doi="10.1002/(SICI)1097-0282(19991005)50:43.0.CO;2-U">Article CAS PubMed Google Scholar

Abdel-Mottaleb, M. S. A. & Abdel-Mottaleb, Y. In search for effective and safe drugs against SARS-CoV-2: Part I] simulated interactions between selected nutraceuticals, ACE2 enzyme and S Protein simple peptide sequences. ChemRxiv https://doi.org/10.26434/chemrxiv.12155235.v1 (2020).

Article Google Scholar

Taidi, L., Maurady, A. & Britel, M. R. Molecular docking study and molecular dynamic simulation of human cyclooxygenase-2 (COX-2) with selected eutypoids. J. Biomol. Struct. Dyn. 40, 1189–1204 (2022).

Article CAS PubMed Google Scholar

Dhanjal, J. K. et al. Computational structure-based de novo design of hypothetical inhibitors against the anti-inflammatory target COX-2. PLoS ONE 10, e0134691 (2015).

Article PubMed PubMed Central Google Scholar

Jack, K. S., Asaruddin, M. R. B. & Bhawani, S. A. Pharmacophore study, molecular docking and molecular dynamic simulation of virgin coconut oil derivatives as anti-inflammatory agent against COX-2. Chem. Biol. Technol. Agric. 9, 73 (2022).

Article CAS Google Scholar

Sadeghi, M., Miroliaei, M., Fateminasab, F. & Moradi, M. Screening cyclooxygenase-2 inhibitors from Allium sativum L. compounds: In silico approach. J. Mol. Model. 28, 24 (2021).

Article PubMed Google Scholar

Download references

The authors thank the management of Vellore Institute of Technology, Vellore for providing the facilities to carry out this work.

School of Biosciences and Technology, Vellore Institute of Technology, Vellore, Tamil Nadu, 632014, India

Priyanka Jain & C. Sudandiradoss

Centre for Nanobiotechnology, Vellore Institute of Technology, Vellore, Tamil Nadu, 632014, India

Jitendra Satija

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

P.J.: Conceptualization and design of the study. Computational analysis, and writing the first draft. Dr. C.S.D.: Conceptualization, supervision, data-curation, review, and editing. J.S.: Conceptualization, review, and editing. All authors read and approved the final manuscript.

Correspondence to C. Sudandiradoss.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Jain, P., Satija, J. & Sudandiradoss, C. Discovery of andrographolide hit analog as a potent cyclooxygenase-2 inhibitor through consensus MD-simulation, electrostatic potential energy simulation and ligand efficiency metrics. Sci Rep 13, 8147 (2023). https://doi.org/10.1038/s41598-023-35192-7

Download citation

Received: 10 November 2022

Accepted: 14 May 2023

Published: 19 May 2023

DOI: https://doi.org/10.1038/s41598-023-35192-7

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.