General
Double beam UV-Visible Spectrometer (Elico SL 218), Bruker Advance-400 (400 MHz for 1H-NMR and 100 MHz for 13C-NMR), Agilent Cary 630 FTIR Spectrometer (Range: 4000 − 450 cm-1) and Agilent 6520 Q-TOF Mass Spectrometer with AGILENT 1200 HPLC System (ESI/APCI-HRMS & LC-HRMS) were used to record the spectra.
Plant material
The plant was collected from the Chauraas Campus, HNB Garhwal University, District Tehri (UK, India) and the identification of the plant was done by the expertise of the Department of Botany, FRI, Dehradun, Uttarakhand and herbarium is deposited (210/Dis./2021/Syst.Bot./Rev.Gen./4–5).
Isolation and characterization
The leaves of Nyctanthes arbor-tristis was dried and powdered (Dry weight = 2 kg) and then extracted with n-Hexane, Chloroform, Ethyl acetate and Methanol respectively with 1:5 (wt/vol.) ratio of plant material and the solvent by using Soxhlet extractor. The extracts were concentrated at reduced pressure (70kPa) and temperature by using rotatory evaporator (Roteva®, Medica instrument MFG.). The extracts; n-Hexane (72.30 gm), Chloroform (87.54 gm), Ethyl acetate (78.67 gm) and Methanol (76.43 gm) were stored for further experiment.The methanolic extract was adsorbed on 20 gm of silica gel and chromatographed on column packed with silica gel (100–120 mesh 1.4 Kg) and eluted with mixture of chloroform and methanol with increasing mixture polarity (chloroform-methanol, 98:2, 96:4 to 60:40). The eluent of chloroform-methanol mixture at 92:8 and 82:18 gives crude 1(light yellow color powdered crystal) and 2 (white color hydroscopic crystal). The crude 1 and 2 were concentrated, washed with 2N HCL and recrystallized with dissolving in hexane. The compound 1 (140 mg) with molecular formula C27H34O13 shows following spectral data: MP: 226–229; [α]D: -93˚ (MeOH); UV (MeOH): 206, 227, 300 and 308 nm; 1H NMR (400 MHz, CDCl3): δ 1.02 (3H, d, J = 7.0 Hz), 2.08 (1H, q, J = 7.0, 5.5 Hz), 2.33 (1H, td, J = 5.4, 2.8 Hz), 2.94 (1H, dd, J = 6.9, 5.4 Hz), 3.12–3.25 (2H, 3.18 (dd, J = 3.5, 2.7 Hz), 3.20 (dd, J = 3.5, 2.7 Hz)), 3.44 (1H, td, J = 4.5, 2.7 Hz), 3.58–3.87 (9H, 3.63 (t, J = 3.5 Hz), 3.72 (s), 3.77 (s), 3.81 (d, J = 4.5 Hz), 3.81 (d, J = 4.5 Hz)), 4.07 (1H, dd, J = 8.1, 7.0 Hz), 4.52 (1H, d, J = 2.7 Hz), 4.86 (1H, dd, J = 8.1, 6.9 Hz), 6.34–6.58 (2H, 6.39 (d, J = 2.8 Hz), 6.51 (d, J = 15.7 Hz)), 7.16 (2H, ddd, J = 8.8, 1.2, 0.5 Hz), 7.49 (2H, ddd, J = 8.8, 1.7, 0.5 Hz), 7.71 (1H, d, J = 15.7 Hz), 7.95 (1H, s); 13C NMR (100 MHz, CDCl3): δ 19.3 (1C, s), 40.0 (1C, s), 40.7–40.8 (2C, 40.7 (s), 40.7 (s)), 52.2 (1C, s), 56.0 (1C, s), 62.3 (1C, s), 71.1 (1C, s), 74.0 (1C, s), 76.7 (1C, s), 78.1 (1C, s), 80.2 (1C, s), 86.2 (1C, s), 96.8 (1C, s), 99.2 (1C, s), 109.7 (1C, s), 114.3 (2C, s), 115.1 (1C, s), 128.7 (2C, s), 130.3 (1C, s), 145.6 (1C, s), 154.1 (1C, s), 159.8 (1C, s), 166.8 (1C, s), 168.5 (1C, s); IR bands (KBr, MeOH): 3424, 2950, 1719, 1662, 1633, 1514, 1450, 1377, 1220, 1180, 1050, 877, 835, 750 and 591 cm− 1; Ms (m/z): 402, 226, 200, 178, 161, 156.
The compound 2 (163 mg) with molecular formula C26H32O14 shows following spectral data, MP: 226–229; [α]D: -69˚ (MeOH); 1H NMR (400 MHz, CDCl3): δ 1.02 (3H, d, J = 7.0 Hz), 2.08 (1H, quintd, J = 7.0, 5.5 Hz), 2.33 (1H, td, J = 5.4, 2.8 Hz), 2.94 (1H, dd, J = 6.9, 5.4 Hz), 3.12–3.25 (2H, 3.18 (dd, J = 3.5, 2.7 Hz), 3.20 (dd, J = 3.5, 2.7 Hz)), 3.44 (1H, td, J = 4.5, 2.7 Hz), 3.58–3.87 (6H, 3.63 (t, J = 3.5 Hz), 3.72 (s), 3.81 (d, J = 4.5 Hz), 3.81 (d, J = 4.5 Hz)), 4.07 (1H, dd, J = 8.1, 7.0 Hz), 4.52 (1H, d, J = 2.7 Hz), 4.88 (1H, dd, J = 8.1, 6.9 Hz), 6.39 (1H, d, J = 2.8 Hz), 6.57 (1H, d, J = 15.7 Hz), 6.77 (1H, dd, J = 8.4, 0.5 Hz), 7.28 (1H, dd, J = 1.9, 0.5 Hz), 7.63–7.80 (2H, 7.69 (dd, J = 8.4, 1.9 Hz), 7.74 (d, J = 15.7 Hz)), 7.95 (1H, s); 13C NMR (100 MHz, CDCl3): δ 19.3 (1C, s), 40.0 (1C, s), 40.7–40.8 (2C, 40.7 s), 52.2 (1C, s), 62.3 (1C, s), 71.1 (1C, s), 74.0 (1C, s), 76.7 (1C, s), 78.1 (1C, s), 80.2 (1C, s), 86.2 (1C, s), 96.8 (1C, s), 99.2 (1C, s), 109.7 (1C, s), 114.8 (1C, s), 115.1 (1C, s), 115.8 (1C, s), 128.2 (1C, s), 128.7 (1C, s), 145.6 (1C, s), 145.8 (1C, s), 146.3 (1C, s), 154.1 (1C, s), 166.8 (1C, s), 168.5 (1C, s); IR bands (KBr, MeOH): 3350, 2920, 1690, 1630, 1520, 1440, 1270, 1080, 950, 860, 810 and 770 cm− 1; Ms (m/z): 404, 226, 210, 178, 163, 156.
Cyclooxygenase inhibition assay
The repressive action of various extracts against COX-2 was calculated by the fair marketable quality COX-(human) inhibitor screening kit (Catalog No. 701070, 701080, Cayman Chemical, PLM Pvt. Ltd., New Delhi, India) based on the producer’s directions. The prepared reaction buffer contains 160 uL, 0.1 M Tris-HCl (pH 8.0), 5mM EDTA and 2 mM phenol. COX-1/2 (10 uL) was added in 10uL, 30uM solution of each extract in the presence of 10 uL Heme solution. All these test solutions were relaxed for 10 minutes at 37°C. Further, 10 uL arachidonic acid and 30 uL stannous chloride were added to test solution (sample) and blank solution (control) to start the reaction to produce the PGF2α by the reduction of PGH2 and relaxed exactly for 5 minutes at 37°C. The amount of produced PGF2α was calculated for sample and control by immunoassay method using the UV microplate reader at the wavelength 410 nm to calculate the percentage inhibition. The following equation was used to calculate the percentage inhibition:

Where [PGF2α]control is PGF2α produced in blank solution (without extracts) and [PGF2α]sample is PGF2α produced in the test sample solution (with extracts). The test compounds, dissolved in 5% DMSO, were evaluated at a concentration of 30 µM.
ADMET Calculation
In order to know drug likeliness of molecule/ligand PKCSM [39] and Swiss ADME [40] was used to predict the molecular weight, water partition coefficient (MLogP), H-bond donors, H-bond acceptor, total polar surface area (TPSA), Lipinski (drug-likeness), hepatotoxicity and Oral Rat Acute Toxicity (LD50) [41].
Molecular docking study
Geometry optimization of reported molecules
The designing of all compounds was done by using ACD Chemsketch and their optimization was done by choosing molecular mechanics (MM2) as a force field. These optimized compounds were used for docking. The following 21 reported phytochemicals of Nyctanthes arbor-tristis were used for the docking: abortitristoside A (1), abortitristoside B (2), abortitristoside C (3), abortitristosideD (4), abortitristoside E (5), nyctoside A (6), 6-beta-hydroxyloganin (7), nyctanthoside (8), rngyolne (9), astragalin (10), 7-o-trans-cinnamoyl-6b-hydroxyloganin (11), 6-o-trans-cinnamoyl-6b-hydroxyloganin (12), 6,7-di-o-benzoylnyctanthoside (13), desrhamnosylverbanscoside (14), arborside A (15), arborside A tetraacetate (16), arborside B (17), arborside B tetraacetate (18), arborside C (19), arborside C pentaacetate (20), arborside D hexaacetate (21) [10].
Protein Preparation
The groundwork of protein was done employing UCSF Chimera 1.11.2 in dock prep module, where the deletion of inherent ligand attached to crystal structure, hydrogen addition, water, solvents, replacing partial residues employing Dunbrack rotamer library, renovation of bromo-UMP to UMP, renovation of selenomethionine to methionine, alteration of methylselenyl-dCMP to CPM. AMBER.ff14SB force field were used to dispensed the charges on protein molecule [42]. Then the augmented protein was used for docking study.
Molecular Docking, Docking parameters & Post Docking modeling
The docking of all reference compounds and the selected phytochemicals of Nyctanthes arbor-tristis was performed using PyRx against COX-1 (PDB ID – 3KK6) and COX-2 (PDB ID – 3LN1). This software used generic algorithms for the docking. The defined parameters in PyRx for docking with Grid center (X = 26.0916, Y = 157.3507, Z = 10.1225 for COX-1 and X = 48.1391, Y = 33.8036, 25.3268 for COX-2), Dimension (X = 82.5649, Y = 89.4395, Z = 118.3433 Å for COX-1 and X = 133.5613, 89.4395, 118.3433 Å for COX-2), population size (n = 200), generations (g = 10) and number of solutions for each compound (s = 1). Based on the above set parameters, the compounds were screened. The post dock modeling of docked poses was executed by Discovery studio visualizer v3.5.The top four compounds were chosen by considering the lowest energy of docking/binding. The binding energy of a ligand to the protein is given by:

where EVDW means for Vander Waal energy, EH−BOND means for hydrogen bonding energy and EELECTROSTATIC means for electro statistic energy [43].